微通道入口角度對微液滴生成過程的影響

李藝凡，夏國棟，王 軍

（北京工業(yè)大學(xué)環(huán)境與能源工程學(xué)院，北京 100124）

微通道入口角度對微液滴生成過程的影響

李藝凡，夏國棟，王 軍

（北京工業(yè)大學(xué)環(huán)境與能源工程學(xué)院，北京 100124）

為研究不同流體物性參數(shù)條件下微通道入口角度對微液滴生成過程的影響，采用流體體積（volume of fluid，VOF）法對錯流微通道內(nèi)液滴生成過程進行了三維數(shù)值模擬，對液滴生成過程中兩相壓差及液滴生成周期和尺寸進行了研究.結(jié)果表明:拉伸擠壓流型下，微液滴尺寸隨著通道入口角度增大呈現(xiàn)先減小后增大趨勢.表面張力較大或連續(xù)相黏度較小時，入口角度為90°的微通道有利于生成較短液滴；表面張力較小或連續(xù)相黏度較大時，入口角度為120°的微通道有利于生成較短液滴.對于固定的入口角度，微液滴尺寸隨著表面張力系數(shù)增大而增大，隨連續(xù)相黏度增大而減小.當(dāng)流體物性參數(shù)不同時，生成最短液滴的通道結(jié)構(gòu)不同.綜合考慮通道入口角度、表面張力系數(shù)和連續(xù)相黏度，提出了預(yù)測拉伸擠壓流型下量綱一液滴長度的經(jīng)驗關(guān)聯(lián)式，為微液滴系統(tǒng)的設(shè)計和應(yīng)用提供了指導(dǎo).

兩相流；微液滴；入口角度；壓差；流體物性參數(shù)

近年來，微加工技術(shù)的飛速發(fā)展使微流控系統(tǒng)在生物分析、化學(xué)合成及相關(guān)領(lǐng)域發(fā)揮越來越重要的作用.由于消耗樣品更少，混合速度更快，不易造成交叉污染等潛在優(yōu)勢，離散化微氣泡或微液滴系統(tǒng)在化學(xué)合成、生物鑒定、過程強化等方面得到了廣泛應(yīng)用［1-3］.

T型微通道是最簡單也是最早用于研究微液滴的通道結(jié)構(gòu)［4］.實際應(yīng)用中，往往需要精確控制微氣泡或者微液滴的尺寸，進而對氣泡或液滴進行操控.一些學(xué)者就微通道結(jié)構(gòu)尺寸、流體物性參數(shù)等對氣泡或液滴尺寸的影響進行了研究［5-12］，結(jié)果表明，微通道入口角度對微氣泡的尺寸有很大影響，通過改變?nèi)肟诮嵌瓤梢缘玫礁〉臍馀?，從而增強兩相熱質(zhì)傳遞［13-16］.但通道入口角度對微液滴或微氣泡的影響尚無定論，Tan等［13］在研究中發(fā)現(xiàn)通道入口角度為90°，即常規(guī)T型微通道時，能夠產(chǎn)生最短的氣泡，而侯璟鑫等［14］和Fries等［15］的研究結(jié)論表明，入口角度為銳角時產(chǎn)生的氣泡長度最短.

兩相壓力的變化對微液滴的生成有重要影響，在拉伸擠壓流型下，兩相壓差和表面張力在液滴生成過程中起主導(dǎo)作用［17-19］.微通道入口角度能夠改變流體的慣性碰撞方向，導(dǎo)致兩相壓差發(fā)生變化，從而對微液滴的生成周期和液滴尺寸產(chǎn)生顯著影響.當(dāng)流體物性參數(shù)不同時，微通道入口角度對微液滴尺寸的影響不同，但文獻中往往在特定流體物性參數(shù)下研究通道幾何結(jié)構(gòu)對液滴尺寸的影響，導(dǎo)致得到的結(jié)論不同甚至相反，對于微通道幾何結(jié)構(gòu)及流體物性參數(shù)的綜合作用尚缺乏系統(tǒng)深入的研究.

因此，本文首先對不同入口角度的微通道內(nèi)微液滴的形成過程進行三維數(shù)值模擬，深入分析微通道入口角度對微液滴生成過程中兩相壓差和表面張力的影響.進而對不同表面張力系數(shù)和連續(xù)相黏度條件下，微通道入口角度對微液滴尺寸的影響進行比較，為微液滴系統(tǒng)的設(shè)計和實際應(yīng)用提供指導(dǎo).

1 物理模型

圖1（a）為T型微通道三維結(jié)構(gòu)示意圖.其中水平通道（主通道）入口為入口1，支通道入口為入口2.α為兩相入口夾角，本文分別取30°、60°、90°、120°和150°，α=90°時為常規(guī)T型微通道.連續(xù)相由入口1進入主通道，離散相由入口2進入支通道.支通道長Ld=300 μm，主通道長Lc=3 400 μm.支通道寬wd=100 μm，主通道寬w=100 μm.支通道與主通道高度h相同，均為80 μm.圖1（b）為微液滴生成過程示意圖.圖中A、B兩點在z軸方向上均位于通道中心平面，點A在x軸方向距兩側(cè)通道壁面均為50 μm，與入口2的距離為250 μm.點B在y軸方向上距離兩側(cè)通道壁面均為50 μm且與入口1的距離為300 μm.點A和B的位置分別在離散相流體和連續(xù)相流體內(nèi)，通過研究這2點的壓力變化，分析液滴生成過程中兩相壓差.ra，up和ra，down分別表示兩相上游界面的軸向曲率半徑和下游界面的軸向曲率半徑（x-y平面內(nèi)）.

2 數(shù)值方法

本文采用FLUENT軟件中的VOF模型對微液滴的形成過程進行三維數(shù)值模擬.流動為層流，流體為不可壓縮牛頓流體，忽略重力作用，控制方程包括連續(xù)性方程、動量方程，分別為

每個網(wǎng)格單元中流體所占體積fvol=fVcell.式中:Vcell為計算單元的體積；f（0≤f≤1）為流體體積函數(shù).以油水兩相流動為例，f=1代表該網(wǎng)格單元中全部為水；f=0表示該網(wǎng)格單元中全部為油；0＜f＜1表示該網(wǎng)格單元中存在油水相界面.式（1）（2）中密度ρ和動力黏度系數(shù)μ等物性參數(shù)均為兩相物性參數(shù)的體積加權(quán)平均值，即

式中下標(biāo)c和d分別代表油相（連續(xù)相）和水相（離散相）.水的體積分?jǐn)?shù)為

采用連續(xù)表面力（continuum surface force，CSF）模型對微尺度條件下表面張力進行計算，表面張力源項可表示為

式中:Fs是表面張力源項，表示流場中表面張力作用的影響，存在于包含界面的控制單元內(nèi)；σ為兩相表面張力系數(shù)；κ為界面曲率，κ通過對相界面單位法向量求散度得到

在遠(yuǎn)離壁面的網(wǎng)格單元內(nèi)，相界面單位法向量為

在壁面附近的網(wǎng)格單元內(nèi)，對于給定的接觸角θ，計算單元內(nèi)相界面單位法向量為

式中:^nw表示壁面單位法向量；^тw表示壁面單位切向量.

微通道計算區(qū)域均采用六面體結(jié)構(gòu)化網(wǎng)格進行劃分.幾何界面重構(gòu)方法采用PLIC算法，壓力速度耦合采用PISO格式，采用PRESTO！算法對壓力項進行離散，采用二階迎風(fēng)格式對動量進行離散.2個入口均為速度入口邊界條件，出口相對壓力為零.時間步長應(yīng)根據(jù)計算結(jié)果的穩(wěn)定性及收斂性選取，本文為0.1 ms.本文以甘油的水溶液作為離散相，以硅油作為連續(xù)相.離散相和連續(xù)相的密度分別為1 162、984 kg/m3，離散相動力黏度為10.58 mPa·s，連續(xù)相動力黏度為28.6～228.6 mPa·s，表面張力系數(shù)為1.5～32.5 mN/m.選取網(wǎng)格數(shù)為2.96萬、7.75萬、9.99萬和16.9萬的常規(guī)T型微通道計算單元進行網(wǎng)格獨立性檢驗（α=90°，vc=0.83 mm/s，vd=0.69 mm/s，μc=68.6 mPa·s，σ=12.5 mN/m）.采用液滴在微通道中心平面內(nèi)的軸向最大長度L與主通道寬度w的比值定義液滴的量綱一尺寸.網(wǎng)格數(shù)為2.96萬、7.75萬、9.99萬的計算單元與網(wǎng)格數(shù)為16.9萬的計算單元得到的量綱一液滴尺寸比較，最大誤差分別為4.5％、2.1％和0.01％.因此，選取網(wǎng)格數(shù)為9.99萬（即每100 μm設(shè)置16個節(jié)點）進行數(shù)值計算，既保證了結(jié)果的精確性，又節(jié)省計算時間.

為了驗證數(shù)值方法的有效性，本文與Li等［18］采用相同條件（α=90°，vd=0.69 mm/s，vc=0.83、1.25、2.00、2.99、3.46 mm/s，ρd=1 162 kg/m3，ρc= 984 kg/m3，μd=10.58 mPa·s，μc=68.6 mPa·s，σ= 12.5 mN/m），對不同毛細(xì)數(shù)（Ca）下微液滴的量綱一尺寸進行計算，并與Li等［18］的實驗數(shù)據(jù)進行比較，如圖2所示.從圖2可以看出，在所研究的Ca范圍內(nèi)，數(shù)值模擬結(jié)果與文獻得到的結(jié)果吻合度較好，驗證了該數(shù)值模擬方法的有效性.

式中μc、vc、σ分別為連續(xù)相動力黏度、連續(xù)相流速和兩相表面張力系數(shù).

3 結(jié)果與討論

3.1微通道入口角度對微液滴生成過程的影響

當(dāng)連續(xù)相流速vc=0.83 mm/s，離散相流速vd= 0.69 mm/s時，毛細(xì)數(shù)較?。–a=0.004 6），T型微通道內(nèi)液液兩相流型為彈狀流，又稱拉伸擠壓狀流型，如圖3所示.彈狀流中，液滴占據(jù)整個主通道，液滴軸向長度大于主通道寬度，液滴呈彈狀或塞狀.微液滴的生成過程分為3個階段:第1階段為t=3.237 6～3.275 5 s，液滴生成后離散相向支通道內(nèi)收縮，而后逐漸恢復(fù)到初始狀態(tài)并準(zhǔn)備進入主通道，稱該階段為“準(zhǔn)備”階段；第2階段為t=3.275 5～3.390 0 s，離散相在主通道占據(jù)越來越大的空間，t=3.390 0 s時幾乎完全堵塞主通道，稱該階段為“發(fā)展”階段；第3階段為t=3.390 0～3.570 5 s，由于離散相堵塞主通道，連續(xù)相壓力增大，離散相形成“頸部”.“頸部”在兩相壓差作用下逐漸變細(xì)，同時離散相被向下游拉伸，t=3.570 5 s時離散相“頸部”幾乎與T型拐角處壁面接觸，最終離散相在“頸部”斷裂，稱這個階段為“斷裂”階段.微液滴生成周期為0.332 9 s，其中3個階段所需時間分別為TⅠ=0.037 9 s，TⅡ=0.114 5 s，TⅢ=0.180 5 s.

彈狀流時，連續(xù)相流速較低，連續(xù)相對離散相的黏性剪切作用較弱，由于微液滴生成過程中離散相堵塞主通道，因此兩相壓差在微液滴生成過程中起主導(dǎo)作用.兩相入口夾角會影響流體慣性碰撞的方向，進而影響微液滴生成過程.圖4、5分別為微通道入口角度α=30°和150°時微液滴的生成過程.如圖4所示，α=30°時，離散相“頸部”斷裂后，生成的液滴向下游運動.支通道與主通道交叉處的兩相界面靠近支通道右側(cè)壁面的部分向支通道收縮，同時交叉處兩相界面靠近支通道左側(cè)壁面的部分向主通道運動，逐漸進入主通道，這2個過程同時發(fā)生（如t=3.266 8～3.310 0 s），直到t=3.390 0 s時，離散相幾乎完全堵塞主通道.如圖5所示，當(dāng)α= 150°時，離散相“頸部”斷裂后，通道交錯處兩相界面靠近支通道左側(cè)壁面的部分向支通道收縮，同時兩相界面靠近支通道右側(cè)壁面的部分向主通道運動（如t=2.722 5～2.850 0 s）.當(dāng)t=2.946 2 s時，兩相界面幾乎與主通道側(cè)壁面接觸.α=30°和150°時，由于通道入口角度的變化，兩相相互作用方向改變，液滴生成過程的第1、2階段同時進行.

當(dāng)α=30°時，由于兩相入口夾角較小，兩相流體接近于并行流動，而α=150°時，兩相流體流動方向接近相反.α=30°時，液滴生成過程的前2個階段所需時間為123.2 ms，比T型微通道內(nèi)液滴生成的前2個階段所需時間短（TⅠ+TⅡ=152.4 ms）.α =150°時，液滴生成的前2個階段所需時間為223.7 ms，較T型微通道所需時間長.這是由于入口角度越大，連續(xù)相對離散相的阻礙作用及兩相流體的慣性碰撞作用越強，因此離散相進入主通道所需時間變長.α=30°、90°和150°時，液滴生成過程的“斷裂”階段所需時間分別為302.5、180.5、171.3 ms.可見，α越大，“斷裂”階段所需時間越短.有圖4、5可知，α=30°和150°時，液滴生成周期分別為425.7、395.0 ms，液滴生成周期越長，液滴尺寸越大.

圖6比較了不同入口角度α下離散相與連續(xù)相間壓差隨時間的變化.兩相壓差Δp=pd-pc，pd表示離散相壓力，pc表示連續(xù)相壓力，圖1中A、B兩點處壓力分別代表離散相壓力和連續(xù)相壓力.由圖6可知，液滴生成過程中，兩相間壓差呈周期性變化，壓差波動規(guī)律與微液滴生成過程相吻合.α越大，兩相壓差的平均值越大.α=30°時液滴生成周期最長，這與圖3～5結(jié)果一致.

已有研究表明，彈狀流時表面張力對微液滴生成有重要作用，表面張力與拉普拉斯壓差成正比，對液滴生成有阻礙作用［20］.微液滴生成過程中，作用在離散相下游界面的壓力為pdown=-σ（1/ra，down+ 1/rr，down）.式中:ra，down和rr，down分別為離散相下游界面的軸向曲率半徑（x-y平面內(nèi)）和徑向曲率半徑（y-z平面內(nèi)）；負(fù)號表示壓力作用方向沿x軸負(fù)方向.同樣地，作用在兩相上游界面的壓力為pup=σ （1/ra，up+1/rr，up）.式中:ra，up和rr，up分別為兩相上游界面的軸向曲率半徑和徑向曲率半徑，壓力的作用方向為x軸正方向.由于通道的高度h小于通道寬度w，界面徑向曲率半徑的變化被通道高度限制，因此rr，down≈rr，up≈h/2.從圖3～5可以看出，對于不同的α，微液滴生成過程中離散相下游界面的軸向曲率半徑基本相同，ra，down=w/2，即慣性碰撞方向基本不會影響下游界面的軸向曲率半徑.但隨著α增大，ra，up明顯減小.離散相兩端界面上的拉普拉斯壓差為ΔpL=pdown+pup=-σ（2/w-1/ra，up），從圖3～5可以看出，不同α?xí)r均有ra，up＞w/2，因此ΔpL的方向為x軸負(fù)方向，且α越小，ra，up越大，ΔpL越大.因此α較小時，表面張力對液滴生成的阻礙作用較大.同時，由圖6可知，α=30°時兩相壓差最小，即促使離散相頸縮的動力最小，因此，α=30°時“斷裂”階段所需時間比α=90°和α=150°長.此外，通道入口角度對離散相進入主通道的速度和離散相“頸部”收縮速度均產(chǎn)生影響，在上述3種結(jié)構(gòu)中，T型微通道內(nèi)液滴生成周期最短，液滴尺寸最小.

3.2表面張力變化時微通道入口角度對液滴尺寸的影響

微液滴尺寸對微尺度兩相流傳熱傳質(zhì)特性有重要影響.文獻中往往在特定的流體物性參數(shù)條件下研究通道幾何結(jié)構(gòu)對微液滴尺寸的影響，忽略了通道幾何尺寸和流體物性參數(shù)的綜合作用.圖7比較了兩相表面張力系數(shù)不同時，微液滴量綱一長度隨通道入口角度的變化（vc=0.83 mm/s，vd=0.69 mm/s，μc=68.6 mPa·s）.由圖可知，入口角度對微液滴尺寸的影響隨著表面張力系數(shù)的變化而不同.表面張力系數(shù)較大時（σ≥7.5 mN/m），隨著入口角度增大，微液滴尺寸逐漸減小，α=90°時液滴長度最小，之后液滴尺寸逐漸增大.當(dāng)σ=4.5、2.5 mN/m，液滴長度隨入口角度的變化也有類似規(guī)律，但α= 120°時液滴長度達到最小值.表面張力系數(shù)σ減小到2.0 mN/m時，微通道內(nèi)兩相流動處于彈狀流與射流的過渡期，α=30°和60°的微通道內(nèi)液滴長度大幅減小，而α≥90°的微通道內(nèi)液滴尺寸變化不大.σ=1.75 mN/m時，α=90°的微通道內(nèi)微液滴生成機制向射流轉(zhuǎn)變，液滴尺寸明顯減小.σ=1.5 mN/m時，不同入口角度的微通道內(nèi)液滴生成機制均完全轉(zhuǎn)變?yōu)樯淞?，離散相沿著下游方向被拉長并與連續(xù)相平行流動，離散相尾部流體周期性地脫落形成液滴.射流機制下，入口角度對液滴尺寸影響很小.隨著σ從32.5 mN/m減小到1.5 mN/m，α= 30°、60°、90°、120°、150°時液滴尺寸分別減小了58.9％、50.8％、43.3％、50.4％、51％，α=30°的微通道內(nèi)液滴尺寸變化最大，常規(guī)T型微通道內(nèi)（α= 90°）液滴尺寸受兩相表面張力系數(shù)影響最小.

連續(xù)相流速vc=0.83 mm/s，離散相流速vd= 0.69 mm/s，毛細(xì)數(shù)Ca=0.032 5，表面張力系數(shù)σ= 1.75 mN/m條件下，不同入口角度的微通道內(nèi)兩相流型如圖8所示.由圖8可知，液滴與離散相脫離的位置開始向支通道與主通道交叉點的下游移動，液滴與離散相斷裂處產(chǎn)生拖尾.α越小拖尾越明顯，液滴生成機制越接近射流機制，即通過減小通道入口角度，可使微通道內(nèi)流型提前轉(zhuǎn)變?yōu)樯淞? α=30°、60°、90°的微通道內(nèi)液滴尺寸較小且基本相等.α=120°、150°的微通道內(nèi)離散相拖尾不明顯，流型更接近于彈狀流，液滴尺寸較大.

3.3連續(xù)相黏度變化時微通道入口角度對液滴尺寸的影響

圖9為不同連續(xù)相黏度條件下微液滴尺寸隨通道入口角度的變化（μd=10.58 mPa·s，σ=12.5 mN/ m，vc=0.83 mm/s，vd=0.69 mm/s）.如圖9所示，由于α=30°的微通道內(nèi)兩相流動方向接近平行，兩相壓差較小，連續(xù)相黏度較?。é蘡=28.6 mPa·s）時，連續(xù)相的黏性剪切作用很弱，因此微通道中兩相平行流動，無微液滴產(chǎn)生.在α≥60°的微通道內(nèi)，隨著α增大，液滴尺寸先減小再增大，α=90°時取得最小值.μc=68.6、108.6、148.6 mPa·s時，5種微通道中均能產(chǎn)生液滴，液滴尺寸隨α的變化規(guī)律與前述相同.連續(xù)相黏度增大到188.6 mPa·s之后，液滴尺寸的最小值點向α增大的方向移動，α=120°時生成的液滴最小.雖然彈狀流時連續(xù)相黏性剪切作用較小，但從圖9可看出連續(xù)相黏度對液滴尺寸具有明顯影響.5種微通道中微液滴尺寸均隨連續(xù)相流體黏度的增大而減小，且連續(xù)相黏度越大，液滴尺寸減小的速率越慢.隨μc從28.6 mPa·s增大到228.6 mPa·s，α=30°的微通道內(nèi)兩相流型由平行流變?yōu)閺棤盍?；?60°、90°、120°、150°時液滴尺寸分別減小了28.9％、17.4％、24.5％、23.9％，α=60°的微通道內(nèi)液滴尺寸變化最大，常規(guī)T型微通道內(nèi)（α= 90°）液滴尺寸受連續(xù)相黏度影響最小.

當(dāng)μc=28.6 mPa·s，vc=0.83 mm/s，vd=0.69 mm/s，σ=12.5 mN/m，Ca=0.001 9時，比較不同入口角度的微通道內(nèi)兩相流型，結(jié)果如圖10所示.該條件下α=30°的微通道內(nèi)兩相流流型為平行流.其他結(jié)構(gòu)微通道內(nèi)均為彈狀流，液滴以拉伸擠壓的方式分散于連續(xù)相流體中，表面張力和兩相壓差在液滴生成過程中起主導(dǎo)作用.從圖10可明顯看出，隨著α增大，液滴尺寸先減小后增大.

綜上所述，微通道入口角度，兩相表面張力系數(shù)及連續(xù)相黏度均對微液滴尺寸有重要影響.綜合考慮以上3個參數(shù)，通過對彈狀流條件下的模擬數(shù)據(jù)進行擬合，得到了預(yù)測彈狀流下量綱一液滴長度的經(jīng)驗關(guān)聯(lián)式

由于σ=1.5、1.75和2.0 mN/m時，不同入口角度的微通道內(nèi)兩相流動逐漸由彈狀流向射流機制過渡.射流機制下液滴尺寸基本不受入口角度影響，因此液滴尺寸不符合式（11）.α=30°且μc= 28.6 mPa·s時，通道中無液滴生成，因此式（11）不適用于此工況.圖11比較了彈狀流條件下數(shù)值模擬得到的量綱一液滴長度與式（11）計算得到的量綱一液滴長度.從圖11可以看出，關(guān)聯(lián)式預(yù)測值與數(shù)值模擬結(jié)果的誤差基本在±10％以內(nèi).因此，在本文研究的參數(shù)范圍內(nèi)，式（11）能夠較準(zhǔn)確地預(yù)測量綱一液滴長度.

4 結(jié)論

1）彈狀流時，兩相壓差和表面張力在微液滴生成過程中起主導(dǎo)作用.改變微通道入口角度能夠使兩相壓差發(fā)生變化，影響離散相進入主通道和離散相頸縮所需時間，從而影響微液滴生成周期和液滴尺寸.

2）表面張力系數(shù)不同時，入口角度對微液滴尺寸的影響不同.拉伸擠壓機制下，微液滴長度隨入口角度增大呈現(xiàn)先減小后增大趨勢.表面張力系數(shù)在7.5～32.5 mN/m范圍內(nèi)時，入口角度90°的微通道內(nèi)液滴長度最短；表面張力系數(shù)為2.5、4.5 mN/ m，入口角度為120°時能得到最短液滴.隨著表面張力減小，入口角度較小的微通道內(nèi)更易形成射流，且射流機制下入口角度對微液滴尺寸影響很小.

3）連續(xù)相黏度越大，入口角度對微液滴尺寸的影響越小，且隨著連續(xù)相黏度增加，能夠生成最短液滴的通道結(jié)構(gòu)由常規(guī)T型微通道變?yōu)槿肟趭A角為120°的微通道.入口角度不變時，微液滴尺寸隨著連續(xù)相黏度增大而減小，且連續(xù)相黏度越大液滴尺寸減小的速率越小.

4）基于本文的參數(shù)變化范圍，考慮通道入口角度、表面張力系數(shù)及連續(xù)相黏度的變化，提出了預(yù)測彈狀流下量綱一液滴長度的關(guān)聯(lián)式.結(jié)果表明，該關(guān)聯(lián)式能夠較好地預(yù)測微液滴尺寸.

［1］TEH SY，LINR，HUNGLH，etal.Droplet microfluidics［J］.Lab on a Chip，2008，8（2）:198-220. ［2］SONG H，TICE J D，ISMAGILOV R F.A microfluidic system for controlling reaction networks in time［J］. Angewandte Chemie，2003，115（7）:792-796.

［3］CHE Z Z，WONG T N，NGUYEN N T.Heat transfer enhancement by recirculating flow within liquid plugs in microchannels［J］.International Journal of Heat and Mass Transfer，2012，55（7）:1947-1956.

［4］THORSEN T，ROBERTS R W，ARNOLD F H，et al. Dynamicpatternformationinavesicle-generating microfluidic device［J］.Applied Physics Letters，2001，86（18）:4163-4166.

［5］LIU H H，ZHANG Y H.Droplet formation in microfluidic cross-junctions［J］.Physics of Fluids，2011，23（8）: 082101.

［6］SHUI L L，MUGELE F，BERG A V D，et al.Geometrycontrolleddropletgenerationinhead-onmicrofluidic devices［J］.Applied Physics Letters，2008，93（15）: 153113.

［7］DESSIMOZ A L，CAVIN L，RENKEN A，et al.Liquidliquidtwo-phaseflowpatternsandmasstransfer characteristics in rectangular glass microreactors［J］. Chemical Engineering Science，2008，63（16）:4035-4044.

［8］LIU H H，ZHANG Y H.Droplet formation in a T-shaped microfluidic junction［J］.Journal of Applied Physics，2009，106（3）:034906.

［9］GUILLOT P，COLIN A.Stability of parallel flows in a microchannel after a T junction［J］.Physical Review E，2005，72（6）:066301.

［10］王維萌，馬一萍，陳斌.十字交叉微通道內(nèi)微液滴生成過程的數(shù)值模擬［J］.化工學(xué)報，2015，66（5）: 1633-1641. WANG W M，MA Y P，CHEN B.Numerical simulation of droplet generation in crossing micro-channel［J］. Journal of Chemical Industry and Engineering，2015，66 （5）:1633-1641.（in Chinese）

［11］DANG M H，YUE J，CHEN G W.Numerical simulation of Taylor bubble formation in a microchannel with a convergingshapemixingjunction［J］.Chemical Engineering Journal，2015，262:616-627.

［12］QIAN D Y，LAWAL A.Numerical study on gas and liquid slugs for Taylor flow in a T-junction microchannel ［J］.Chemical Engineering Science，2006，61（23）: 7609-7625.

［13］TAN J，DUL，XUJH，etal.Surfactant-free microdispersion process of gas in organic solvents in microfluidic devices［J］.AIChE Journal，2011，57 （10）:2647-2656.

［14］侯璟鑫，錢剛，周興貴.氣體入口角度和截面寬高比對微通道內(nèi)泰勒氣泡行為的影響［J］.化工學(xué)報，2013，64（6）:1976-1982. HOU J X，QIAN G，ZHOU X G.Effect of gas inlet angle and cross-section aspect ratio on Taylor bubble behavior in microchannels［J］.Journal of Chemical Industry and Engineering，2013，64（6）:1976-1982.（in Chinese）

［15］FRIES D M，VON ROHR P R.Impact of inlet design on mass transfer in gas-liquid rectangular microchannels［J］. Microfluidics and Nanofluidics，2009，6（1）:27-35.

［16］SHAO N，GAVRIILIDIS A，ANGELI P.Effect of inlet conditions on Taylor bubble length in microchannels［J］. Heat Transfer Engineering，2011，32（13/14）:1117-1125.

［17］YAN Y，GUO D，WEN S Z.Numerical simulation of junction point pressure during droplet formation in a microfluidicT-junction［J］.ChemicalEngineering Science，2012，84:591-601.

［18］LI X B，LI F C，YANG J C，et al.Study on the mechanismofdropletformationinT-junction microchannel［J］.Chemical Engineering Science，2012，69（1）:340-351.

［19］SIVASAMY J，WONG T N，NGUYEN N T，et al.An investigation on the mechanism of droplet formation in a microfluidicT-junction［J］.Microfluidicsand Nanofluidics，2011，11（1）:1-10.

［20］GARSTECKI P，F(xiàn)UERSTMAN M J，STONE H A，et al. Formation of droplets and bubbles in a microfluidic T-junction—scaling and mechanism of break-up［J］.Lab on a Chip，2006，6（3）:437-446.

（責(zé)任編輯 呂小紅）

Effect of Inlet Angle on Droplet Generation in Microchannel

LI Yifan，XIA Guodong，WANG Jun
（College of Environment and Energy Engineering，Beijing University of Technology，Beijing 100124，China）

To study the effect of microchannel inlet angle on the droplet generation process under the condition of varying fluid physical parameters，a three dimensional numerical simulation in the cross-flow microchannel was proposed by the volume of fluid（VOF）method.The two-phase pressure difference in the droplet generation process，the time period and size of the microdroplet were studied.The results show that the microdroplet size decreases firstly and then increases gradually with the increase of the inlet angle in the squeezing regime.A 90°inlet angle is helpful to generating shortest droplet for high surface tension or low viscosity of continuous phase；A 120°inlet angle is conducive to generating shortest droplet for low surface tension or high viscosity of continuous phase.For constant inlet angle，the microdroplet size rises with the increase of the surface tension coefficient and reduces with the increase of the viscosity of continuous phase.The shortest droplet is generated in different microchannel for varying fluid physical parameters.Considering the inlet angle，surface tension coefficient and viscosity of continuous phase，the new correlation of the dimensionless droplet length for the squeezing regime is proposed which can provide guidance for the design and application of the microdroplet system.

two-phase flow；microdroplet；inlet angle；pressure difference；fluid physical parameter

O 359

A

0254-0037（2016）09-1414-08

10.11936/bjutxb2016010069

2016-01-30

北京市自然科學(xué)基金資助項目（3142004）；北京市優(yōu)秀博士學(xué)位論文指導(dǎo)教師科技項目（20131000503）

李藝凡（1988—），女，博士研究生，主要從事微電子設(shè)備的冷卻及微尺度流動與傳熱、氣液兩相流動與傳熱方面的研究，E-mail:yifan100zixuan@163.com

夏國棟（1965—），男，教授，主要從事微電子設(shè)備的冷卻及微尺度流動與傳熱、氣液兩相流動與傳熱方面的研究，E-mail:xgd@bjut.edu.cn