• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于微流控技術(shù)的微液滴融合研究進(jìn)展

    2015-01-04 07:25:56申峰李易劉趙淼曹刃拓王貴人北京工業(yè)大學(xué)機(jī)械工程與應(yīng)用電子技術(shù)學(xué)院北京004DepartmentofMechanicalEngineeringBiomedicalEngineeringProgramUniversityofSouthCarolinaColumbiaSC908USA
    分析化學(xué) 2015年12期
    關(guān)鍵詞:排液微流液滴

    申峰李 易劉趙淼曹刃拓王貴人(北京工業(yè)大學(xué)機(jī)械工程與應(yīng)用電子技術(shù)學(xué)院,北京004)(Department of Mechanical Engineering&Biomedical Engineering Program, University of South Carolina,Columbia,SC,908,USA)

    評述與進(jìn)展

    基于微流控技術(shù)的微液滴融合研究進(jìn)展

    申峰1,2李 易1劉趙淼*E-mail:lzm@bjut.edu.cn1曹刃拓1王貴人21(北京工業(yè)大學(xué)機(jī)械工程與應(yīng)用電子技術(shù)學(xué)院,北京100124)
    2(Department of Mechanical Engineering&Biomedical Engineering Program, University of South Carolina,Columbia,SC,29208,USA)

    隨著芯片實(shí)驗(Lab-on-a-chip)研究的不斷發(fā)展,利用微流控技術(shù)(Microfluidics)對微通道中的液滴進(jìn)行精準(zhǔn)操控,引起了越來越多的研究者關(guān)注。本文介紹了微通道液滴基本理論及融合機(jī)理,通過改變微通道幾何結(jié)構(gòu)和表面特性等影響因素實(shí)現(xiàn)微液滴融合的被動融合方式,通過施加外電場、磁場、溫度場、表面聲波和激光聚焦等作用實(shí)現(xiàn)微液滴融合的主動融合方式,以及液滴融合過程中的流體動力學(xué)特性。最后對微流控芯片液滴融合技術(shù)研究的發(fā)展前景進(jìn)行了展望,為微流控器件設(shè)計、微液滴控制技術(shù)的廣泛利用提供有價值的參考。

    微流控技術(shù);液滴融合;微通道;芯片實(shí)驗室;評述

    1 引 言

    微流控芯片(Microfluidic chip)也稱作芯片實(shí)驗室(Lab-on-a-chip),是近十幾年來發(fā)展迅速的多學(xué)科交叉研究領(lǐng)域[1],它是一種以在微米尺度空間對流體進(jìn)行操控為主要特征的科學(xué)技術(shù)[2]。微流控芯片利用微電子機(jī)械系統(tǒng)(MEMS)技術(shù)在硅、玻璃或高分子材料等基片上制作微通道、微泵、微閥、微反應(yīng)器等功能單元,并施加電場或壓力等各種外力場驅(qū)動,對芯片中的微量流體進(jìn)行控制和處理,充分利用微通道中流體的物理特性,可以實(shí)現(xiàn)對微量流體或樣品的精確操縱、處理與控制的功能[3],可以將生物、化學(xué)等實(shí)驗室的各種操作功能,如樣品制備、反應(yīng)、分離、檢測及細(xì)胞培養(yǎng)、分選、裂解等,集成到一塊幾平方厘米(甚至更?。┑男酒蟍4]。相比于傳統(tǒng)方法,微流控芯片具有成本低、樣品消耗少、易自動化、靈活快速、易微型化等優(yōu)勢,在疾病診斷、細(xì)胞研究、藥物篩選、環(huán)境監(jiān)測、綠色能源、材料合成等領(lǐng)域展示出廣闊的應(yīng)用前景[4]。因此,Nature雜志評論認(rèn)為芯片實(shí)驗室可能成為“本世紀(jì)的技術(shù)”[5~7]。

    微流控芯片液滴是一種新的操縱微小體積液體的技術(shù)。微液滴操控技術(shù)一般是指將兩種不相融的流體,如水和油,其中一種作為連續(xù)相,另一種作為分散相微液滴體積單元(10-15~10-9L),并利用微流控技術(shù),控制兩相流動行為,以實(shí)現(xiàn)微液滴的生成、輸運(yùn)、分裂、融合、分選、定位和捕獲等控制功能[2,4]。微流控芯片液滴技術(shù)已成為微流控芯片研究的一個重要發(fā)展方向[5~8]。微液滴具有高比表面積、傳質(zhì)傳熱時間短和擴(kuò)散距離短等特點(diǎn),在微液滴中進(jìn)行反應(yīng)的時間比常規(guī)體系短,因此,微流控芯片中的微液滴具有微反應(yīng)器的功能[9~13]。此外,微液滴具有較高的單分散性、微小的體積,以及不受限的液滴生成數(shù)量,為化學(xué)和生物研究提供了一個獨(dú)特的微型化平臺,而且微液滴體積及運(yùn)動行為可以精確操控,并具有可重現(xiàn)性。由于液滴界面張力的作用,微液滴內(nèi)部可以產(chǎn)生與外界連續(xù)相流體不同的物理和化學(xué)反應(yīng)環(huán)境,同時,微液滴內(nèi)部流體可以與微通道壁面隔離,以減小通道表面性質(zhì)對內(nèi)部反應(yīng)的影響[8]。由于微液滴平臺的這些獨(dú)特優(yōu)勢,其應(yīng)用領(lǐng)域正在不斷擴(kuò)展,包括細(xì)胞和亞細(xì)胞結(jié)構(gòu)分析[14~17]、蛋白質(zhì)結(jié)晶[18,19]、納米粒子的合成[20,21]和納米分子自組裝[22,23]以及高通量注射等[24,25]。

    基于微液滴的生物化學(xué)反應(yīng),如微納顆粒的材料合成過程中,需要將兩種不同的液滴融合在一起,以便獲得較好的混合反應(yīng)效果,因此,液滴的可控融合(Coalescence)已成為一項重要的微流控技術(shù),對于在液滴內(nèi)實(shí)現(xiàn)多步反應(yīng)十分重要[8~13]。微通道液滴融合是指將兩個或多個微液滴接觸并融合在一起。融合過程中需要克服液滴之間流體薄膜的表面張力作用,并使得液滴界面失穩(wěn)才能實(shí)現(xiàn)融合。由于微通道上下壁面的限制作用,液滴之間的連續(xù)相流體需要從有限的空間排出,增加了融合的難度。同時,液滴界面通常不再是球形,液滴融合過程中毛細(xì)力-黏性力-慣性力的平衡方式與宏觀尺度下的液滴融合過程不同,現(xiàn)有宏觀尺度的液滴融合理論難以對其進(jìn)行定量的描述[2]。此外,由于液滴界面的可變形性,導(dǎo)致了界面張力的不均勻變化和液滴的非線性問題等一系列新的流體力學(xué)問題[8,26],如微通道入口壓力的微小變化可能引發(fā)流動由液滴向射流的轉(zhuǎn)變,導(dǎo)致非線性不穩(wěn)定流動的產(chǎn)生[2,8]。

    目前,微液滴融合主要分為被動融合和主動融合兩類[2,8,9]。被動融合主要利用液滴及連續(xù)相的流動特性,通過改變微通道幾何結(jié)構(gòu)、增加特殊的微結(jié)構(gòu)、調(diào)節(jié)液滴生成速率和改變通道表面浸潤性等方法,使上游液滴流速變慢(或暫時被捕獲在某特定位置),使得下游液滴能追上上游液滴,并在特定位置處發(fā)生接觸并融合。在兩液滴相互靠近的過程中,液滴之間的連續(xù)相流體逐漸被排空;在接觸過程中,由于壓力或者剪切力作用,液滴界面變得不穩(wěn)定并破裂,誘導(dǎo)發(fā)生兩液滴融合[2]。被動融合方式不需要主動控制或電子設(shè)備,比主動控制的微流控系統(tǒng)更簡單,但是需要對微通道結(jié)構(gòu)進(jìn)行更科學(xué)合理的設(shè)計,必須充分掌握微通道兩相流動特性。被動融合方式的一個主要優(yōu)勢是對液滴內(nèi)部物質(zhì)的干擾較小,特別是進(jìn)行細(xì)胞等生化分析時非常重要。此外,液滴被動融合受到流體表面張力波動等自然現(xiàn)象的限制,其融合效率一般比主動融合速率低[27]。

    主動融合是指通過施加外加電場[28]、磁場[29]、溫度場[30]、表面聲波[31,32]和激光聚焦[33]等作用,誘導(dǎo)液滴界面失穩(wěn),發(fā)生融合。相比被動融合,主動融合方式本質(zhì)上更加復(fù)雜,因為需要在微通道上制作微電極,并對電信號進(jìn)行精確的時間控制。此外,微電極可能會對液滴內(nèi)部物質(zhì)造成干擾或污染,比如液滴內(nèi)部物質(zhì)可能沉積在微電極上,電信號可能對生物分子(DNA或蛋白質(zhì))的生物相融性造成影響。主動融合方式可以利用電極加速液滴表面張力的不穩(wěn)定性促進(jìn)融合發(fā)生,具有更高的效率[29]。

    本文首先闡述了微通道內(nèi)液滴融合的基本理論及其影響因素,然后詳細(xì)綜述了目前國內(nèi)外研究中通過控制微通道內(nèi)液滴融合影響因素,以實(shí)現(xiàn)微液滴融合的被動融合和主動融合兩種融合方式,主要包括控制通道幾何結(jié)構(gòu)、表面活性劑和外加電場、磁場、溫度場、表面聲波和激光聚焦等方法,并闡述了微通道內(nèi)液滴融合的流體動力學(xué)特性,最后對微通道內(nèi)液滴融合技術(shù)的發(fā)展進(jìn)行了展望。

    2 微液滴基礎(chǔ)理論

    從微流體力學(xué)角度來看,基于微液滴的生成、輸運(yùn)、反應(yīng)及合成等一系列控制功能都是以多相流形式進(jìn)行的[2]。多相流的根本特點(diǎn)是流動中同時存在著被相界面明顯分開的兩種或者多種物質(zhì)組分,界面的存在將非線性效應(yīng)引入到微流動中,微液滴流動行為變得復(fù)雜[2,8,26]。同時,由于微通道壁面的限制作用,液滴的比表面積非常大,壁面浸潤性引起的邊界效應(yīng)對微液滴流動行為也有重要影響,微液滴的演化規(guī)律與宏觀流動下有顯著不同[2]。深入研究微通道中液滴流動行為問題,將推進(jìn)微流控芯片結(jié)構(gòu)的優(yōu)化設(shè)計和微液滴控制技術(shù)的創(chuàng)新,具有重要理論指導(dǎo)和實(shí)際應(yīng)用價值[2,3,8,26,34~37]。

    2.1 表面張力

    由于液滴界面的存在,必須考慮界面張力的作用。界面張力σ是由于液體交界面處分子引力不均衡而產(chǎn)生的沿著界面的張力,定義為單位長度上的作用力(N/m),或者為單位面積上的能量(J/m2)[2,8]。當(dāng)界面張力不平衡時,可以引發(fā)從低界面張力處到高界面張力處的流動(Marangoni對流)。由溫度梯度或表面活性劑分子改變引起的界面張力驅(qū)動的流動分別稱為熱毛細(xì)對流或溶質(zhì)毛細(xì)對流[8]。同時,界面張力作用使得液滴界面面積趨于最小化,以減小表面能量。孤立的液滴或氣泡其表面是球形時表面積最小,然而,微通道中液滴還受到壁面限制作用,其界面往往不能收縮為球形,而是變?yōu)閺澢?。界面的曲率會?dǎo)致內(nèi)外產(chǎn)生壓力跳躍,即拉普拉斯壓力(Laplace pressure):

    其中R1和R2為液滴界面不同位置處的曲率半徑。半徑越小,ΔP越大。由于液滴界面不同空間位置的曲率半徑不同,液滴內(nèi)部的壓力分布也不同,其對液滴內(nèi)外流場的流動演變過程起重要作用[2]。

    2.2 無量綱參數(shù)

    微液滴多相流的流動行為依賴于各種作用力,包括粘性力、慣性力、重力、表面張力、剪切力、分離壓力等,它們對流動的相對影響可以用以下無量綱參數(shù)表述[2],即:

    其中,ρ,μ,U和l分別為流體密度、粘度、特征速度及特征長度,Re是慣性力(ρU2/l)與粘性力(ρU/l2)的比值,一般小于10,流動表現(xiàn)為層流。多數(shù)情況下慣性力可以忽略,然而當(dāng)微液滴速度很快或者液滴破裂瞬間,慣性力也有一定的影響作用。

    We表示慣性力與界面張力的比值。在微通道中,We一般較宏觀尺度時小很多,界面張力對流動的影響明顯。

    其中,g、Δρ分別為重力加速度和不同流體的密度差。Bo表示重力作用與界面張力之比,一般Bo遠(yuǎn)小于1,即重力的影響相對表面張力較弱,通??梢院雎灾亓Φ淖饔?。

    Ca表示粘性應(yīng)力與界面張力的比值,其中,μc是兩相流體中連續(xù)相粘度,U一般是指液滴融合處的特征速度,σ是連續(xù)相和分散相的界面張力。若Ca較小時,界面張力起主導(dǎo)作用,界面較穩(wěn)定,液滴的端面易產(chǎn)生球形;反之,粘性力起主導(dǎo)液滴流動行為,液滴界面易產(chǎn)生較大的變形或者其形狀變得不對稱。

    從上述無量綱參數(shù)的定義可以發(fā)現(xiàn),微液滴流動的Re數(shù)、We數(shù)和Bo數(shù)都較小,具有粘性層流的特點(diǎn),多數(shù)情況下粘性應(yīng)力和表面張力占主導(dǎo)作用(Ca)。此外,微通道壁面對液滴界面的限制作用、通道壁面浸潤性、微通道幾何尺寸的比值及液滴大小與通道尺寸的關(guān)系等因素都對微液滴多相流有重要影響[8,26]。對于外加力場(如電場、光熱、聲波等)的情況,則需要表征外加力與表面張力的相對重要性,并描述其作用機(jī)制[26]。

    2.3 表面活性劑作用

    為了有利于微液滴的生成和輸運(yùn),通常在微通道連續(xù)相中添加表面活性劑以降低界面張力。表面活性劑可使液滴界面變得相對穩(wěn)定,液滴間不容易融合,從而使液滴成為一種微反應(yīng)器。但是,由于流動和分子擴(kuò)散等因素的影響,使得表面活性劑濃度的空間分布變得不均勻,導(dǎo)致界面張力的空間變化,進(jìn)而影響界面流動。表面活性劑使液滴處于亞穩(wěn)狀態(tài)[38],當(dāng)兩個液滴相互靠近時,如圖1a所示,有兩種作用機(jī)制[26,37]:一是表面活性劑分子的排斥作用使相互靠近的界面更加穩(wěn)定;二是連續(xù)相在排出過程中,流動導(dǎo)致界面表面活性劑濃度發(fā)生變化,改變了界面各點(diǎn)的界面張力σ,產(chǎn)生了Marangoni應(yīng)力[8,38]:

    其中,μd和μc分別為離散相和連續(xù)相的黏度,μ‖為切向速度,r為半徑方向,▽‖為界面某點(diǎn)處界面張力在切線方向的導(dǎo)數(shù)。Marangoni應(yīng)力作用下,產(chǎn)生沿著界面向液膜中心的Marangoni對流,如圖1b所示,阻礙了連續(xù)相的排出,延長了排液時間,使得微液滴更加穩(wěn)定[38]。

    2.4 微液滴融合機(jī)理

    微液滴融合基本包括4個過程:(1)液滴的捕捉或定位;(2)兩液滴的相互靠近、碰撞及變形過程;(3)液滴間連續(xù)相液膜排液過程;(4)液滴界面薄膜破裂及融合發(fā)生過程[39]。目前多數(shù)研究著重于微液滴的生成過程,而對于微液滴融合機(jī)理尚需深入研究。影響液滴融合的因素復(fù)雜,包括:通道幾何結(jié)構(gòu)、入口流速、表面活性劑、壁面性質(zhì)和外加力場等。目前,外加力場作用對液滴的作用強(qiáng)度尚無法進(jìn)行定量化的測量,無法準(zhǔn)確表征外場力與表面張力的相對重要性。此外,由于外加力場作用下產(chǎn)生的液滴界面的變形及瞬態(tài)非線性振動等問題,一直是學(xué)術(shù)研究的難題,目前外加力場作用下的液滴融合尚無精確的理論分析,多數(shù)依靠實(shí)驗經(jīng)驗進(jìn)行。由于篇幅限制,本文未介紹這方面內(nèi)容。

    圖1 微液滴排液過程示意圖[26,37]Fig.1 Schematic representation of film drainage of continuous phase between drops[26,37](a)無表面活性劑排液過程;(b)表面活性劑濃度梯度引起的Marangoni應(yīng)力(紅色箭頭)阻礙了排液。(a)Absence of surfactant;(b)Marangoni stress(red arrows)induced by the gradient of surface density of the surfactant counteracts the film drainage.

    目前,許多學(xué)者用液膜排液理論(Film drainage model)描述兩液滴的被動融合過程[39]。當(dāng)下游液滴與上游液滴接觸后,可能依附在一起,由于中間存在液膜而無法相互融合。兩液滴接觸后,具有相互擠壓或者相對滑移的趨勢,可能發(fā)生旋轉(zhuǎn),導(dǎo)致兩液滴中心點(diǎn)連線不再與通道軸線平行。同時,液滴間的連續(xù)相液膜開始排液,而液膜之間多余的毛細(xì)力同時驅(qū)動其排液[40],排液速度由碰撞速度和液膜粘度等因素決定,即毛細(xì)數(shù)Ca。臨界毛細(xì)數(shù)是判斷融合是否發(fā)生的關(guān)鍵參數(shù)[3]。研究發(fā)現(xiàn),液滴融合多數(shù)發(fā)生在低速率碰撞時,臨界毛細(xì)數(shù)Cac≈10-2,低于臨界值時更容易發(fā)生融合[3],如圖2a所示。當(dāng)毛細(xì)數(shù)大于臨界值時,液滴可能會滑移過先前的液滴不發(fā)生融合[3],如圖2b所示。

    圖2 毛細(xì)數(shù)對液滴融合的影響[3]Fig.2 Influence of capillary number(Ca)on drop coalescence[3](a)毛細(xì)數(shù)低于臨界值,液滴發(fā)生融合;(b)毛細(xì)數(shù)大于臨界值,液滴發(fā)生滑移而不融合。(a)Ca<Cac,coalescence;(b)Ca>Cac,slipping.

    如果液滴保持足夠長時間的接觸,且液膜變得越來越薄,當(dāng)液膜足夠薄時,范德華力和其它分子間作用力開始起主導(dǎo)作用,液膜破裂且液滴融合發(fā)生[3]。通常液膜排液效率隨毛細(xì)數(shù)的增大而增大。另外,液滴碰撞角度對排液時間也有一定影響[42]。Klaseboer等[43]給出了液滴正面碰撞時液膜排液時間的表達(dá)式:

    此外,Hu等[44]發(fā)現(xiàn),兩液滴正面碰撞時液滴融合的臨界毛細(xì)數(shù)依賴于液滴的半徑r和兩相流體的粘度比λ=μd/μc,可以用以下公式表示:

    目前,大量實(shí)驗及數(shù)值模擬研究表明,液膜排液時間、融合的臨界毛細(xì)數(shù)及液滴融合角度還受到液滴尺寸[45,46]、粘度比[45,47]、液滴的整體變形[48,49]和碰撞的偏移[42]等因素影響。同時,壁面親疏水性對液膜排液也有重要影響。Li等[50]研究了親水微通道液滴和液柱間的融合,發(fā)現(xiàn)更好的親水表面、更小的通道尺寸、更短的液滴到入口的距離、更大的液滴尺寸會誘發(fā)較大的毛細(xì)壓力。融合過程的其它形式也得到了研究者的關(guān)注,例如簡單的剪切流引發(fā)液滴融合[51,52],大量乳液在特殊結(jié)構(gòu)中通過平滑表面促進(jìn)液滴融合等[53,54]。

    3 微液滴融合的控制方法

    3.1 控制微通道幾何結(jié)構(gòu)實(shí)現(xiàn)微液滴融合的方法

    3.1.1 利用流動阻礙結(jié)構(gòu)實(shí)現(xiàn)液滴融合的方法通過微加工方法在通道內(nèi)增加流動阻礙結(jié)構(gòu)是最常見的實(shí)現(xiàn)微液滴融合的控制方法。Xu等[55]報道了壁面突起、圓形擴(kuò)張截面和橢圓形擴(kuò)張截面3種不同的流動阻礙結(jié)構(gòu)。液滴融合的過程大致可以描述成:增加阻礙結(jié)構(gòu),使先到達(dá)的液滴停留,捕捉前面液滴,待后面的液滴到達(dá)時與其發(fā)生融合,之后共同運(yùn)動。

    為了實(shí)現(xiàn)前面液滴的流動減速,Tan等[56]設(shè)計并實(shí)驗研究了3種液滴融合的通道結(jié)構(gòu),包括:矩形擴(kuò)張結(jié)構(gòu),錐形擴(kuò)張結(jié)構(gòu)和流動聚焦結(jié)構(gòu),如圖3所示。矩形擴(kuò)張結(jié)構(gòu)的作用是可以使液滴在擴(kuò)張段減速,并有利于液滴間連續(xù)相液體的排出,減小液滴間距,使兩液滴碰撞并而發(fā)生融合。擴(kuò)張段長度與寬度以及兩相流速度對液滴是否融合有重要影響。錐形擴(kuò)張結(jié)構(gòu)通道中,由于離散相液滴的速度急劇減小,以及連續(xù)相在錐形擴(kuò)張中的連續(xù)排液,可以造成多個液滴的不斷融合,融合液滴的數(shù)量不易控制。與前兩種結(jié)構(gòu)不同,十字形的流動整流結(jié)構(gòu)可以通過控制下游上下兩通道的流速,以控制兩液滴間連續(xù)相液體的排液速度,從而實(shí)現(xiàn)液滴在十字交叉處的融合控制[56]。

    圖3 通過流動控制實(shí)現(xiàn)液滴融合的不同微通道結(jié)構(gòu)[56]Fig.3 Three differentmicrochannels for droplet fusion[56]

    運(yùn)用相似的流動整流方法,Gunes等[57]在微通道上、下兩側(cè)增加了許多垂直方向的流動可控的微通道,如圖4a所示。通過控制垂直方向微通道內(nèi)流體的流動,以加速兩液滴間的排液速度,實(shí)現(xiàn)液滴的可控融合,并可以實(shí)現(xiàn)多液滴融合。此結(jié)構(gòu)的優(yōu)點(diǎn)為:液滴之間的碰撞與分離可以獨(dú)立控制,不會產(chǎn)生明顯的流動干擾。Tan等[58]研究了微整流通道寬度對通道交叉處液滴融合的影響,通過增加整流通道的寬度,可以實(shí)現(xiàn)2~6個多液滴融合,如圖4b所示為2~3個不同液滴的融合過程。

    圖4 利用微整流通道控制的液滴融合(a)并列多個微整流通道控制[57];(b)增加整流通道寬度[58]。Fig.4 Controll of droplet fusion using(a)listof smallmediated flow channel[57]and(b)Widenedmediated flow channel[58].

    Niu等[59]在通道擴(kuò)大段設(shè)計加工了兩排微柱作為液滴流動阻礙結(jié)構(gòu)。由于微柱對液滴的擠壓作用,液滴被捕捉在微柱腔內(nèi),直到與后來液滴發(fā)生擠壓融合。當(dāng)連續(xù)相液體壓力大于液滴的界面張力時,融合后的液滴離開微柱腔。

    在微通道中增加一段漸寬的通道已經(jīng)成為液滴融合廣泛采用的一種結(jié)構(gòu)[56,58,60]。通道變寬處由于通道寬度增大,連續(xù)相速度降低,致使前面液滴被減速,并與后面液滴不斷靠近直至發(fā)生融合。然而,與普遍認(rèn)為的液滴間由于擠壓融合的現(xiàn)象不同,Bremond等[61]實(shí)驗發(fā)現(xiàn),在擴(kuò)展通道中液滴并不是第一次相遇時融合,而是在分離的過程中當(dāng)前一個液滴進(jìn)入到窄通道時與后一個液滴發(fā)生融合,如圖5a所示。這種分離引起表面張力的瞬時不穩(wěn)定,并產(chǎn)生了液滴表面的突起(圖5b),進(jìn)而打破了液滴的界面張力,使得這兩個液滴能夠迅速靠近并發(fā)生融合。此后,Lai等[62]研究了此現(xiàn)象的產(chǎn)生機(jī)理,結(jié)果表明,這種突起使得局部表面積迅速增加,從而導(dǎo)致液滴表面某些位置的不穩(wěn)定性,進(jìn)而使液滴發(fā)生融合。

    圖5 液滴在擴(kuò)展通道分離時流動擾動引起的液滴(a)融合過程和(b)融合機(jī)理[61]Fig.5 Droplet in separation caused by flow destabilization of extend the channel(a)fusion process and(b)fusionmechanism[61]

    3.1.2 利用交叉結(jié)構(gòu)實(shí)現(xiàn)液滴融合的方法除了流動阻礙結(jié)構(gòu),微通道交叉結(jié)構(gòu)也是一種常見的能夠有效促進(jìn)液滴融合的方法。其融合的基本原理是保證兩相向運(yùn)動的微液滴在交叉處相遇,由于擠壓作用實(shí)現(xiàn)排液,并使液膜破裂,從而實(shí)現(xiàn)融合。相對于同方向運(yùn)動的液滴融合,相向碰撞的液滴融合方法的效率較高,且不會產(chǎn)生多液滴融合的現(xiàn)象,液滴融合的可控性強(qiáng)。但是,為保證液滴相向流動時能夠同時到達(dá)交匯處,必須對上游的微液滴運(yùn)動進(jìn)行精確的同步控制,增加了實(shí)際操作的困難。

    微通道交叉結(jié)構(gòu)一般包括T型和Y型兩種。液滴融合過程中,液體粘度和兩相流速度對是否融合有重要影響,一般增大連續(xù)相液體粘度與兩相流速度會阻礙液滴融合。同時,根據(jù)液膜排液理論,兩液滴表面接觸時間與排液時間也是判斷液滴是否融合的關(guān)鍵因素。其中,液滴在通道內(nèi)的流速與總流量、液滴尺寸、微通道寬高比等對排液效率有重要影響。Yang等[63]實(shí)驗研究了T型交叉微通道內(nèi)氣泡的融合,觀察到了微氣泡碰撞后會產(chǎn)生的3種融合現(xiàn)象,即:絕對融合(Absolute coalescence)、概率性融合(Probabilistic coalescence)和不融合(Non-coalescence),并對融合概率的影響因素進(jìn)行了分析。

    Christopher等[41]利用T型交叉微通道研究了液滴尺寸及流動速度對融合過程的影響,發(fā)現(xiàn)在通道的T型交叉處碰撞后會發(fā)生6種現(xiàn)象,即:除融合和滑移兩種現(xiàn)象(圖2)外,還會發(fā)生分裂、多次分裂、后融合和后滑移4種現(xiàn)象,如圖6所示。研究發(fā)現(xiàn),液滴尺寸及流動速度對融合過程有重要影響,碰撞速度較低時容易發(fā)生融合,小液滴流動速度較快時不發(fā)生融合,流速快的大液滴相互擠壓會分裂成多個液滴。只有對液滴到達(dá)交叉處的時間進(jìn)行精確的同步控制,使兩微液滴到達(dá)的時間差遠(yuǎn)低于在交叉處停留的時間,才可使兩液滴發(fā)生融合。

    圖6 液滴在通道交叉處碰撞的(a)分裂;(b)多次分裂;(c)后融合;(d)后滑移4種流動現(xiàn)象[41]Fig.6 Four different forms of droplet collision atmicrofluidic junctions(a)splitting,(b)multiple splitting, (c)late coalescence and(d)late splitting[41]

    Liu等[64]研究了Y型交叉微通道不同夾角30°,60°,90°,120°,150°和180°對液滴融合的影響,發(fā)現(xiàn)液滴融合的臨界毛細(xì)數(shù)隨通道夾角的增大而減小,增大兩相流速和通道夾角可以減小液膜排液時間。Wang等[65]設(shè)計了夾角為60°,120°和180°的3種Y型交叉微通道,觀測了微液滴在交叉處接觸、擠壓、融合或分離的過程,分析了液滴尺寸、碰撞角度、流體物理性質(zhì)對液膜排液時間、液滴接觸時間以及融合臨界毛細(xì)數(shù)的影響,實(shí)驗結(jié)果表明,減小液滴碰撞角度能夠更好的促進(jìn)液滴融合,并提出了預(yù)測融合的臨界毛細(xì)數(shù)公式。

    上述兩種交叉通道都需要對液滴到達(dá)交叉處的時間進(jìn)行精確的同步控制,增加了實(shí)驗的難度,為此,Wang等[66]設(shè)計了一種十字交叉型通道,解決了微液滴同步問題。與T型和Y型交叉結(jié)構(gòu)相比,十字交叉型通道能提供上下對稱相對穩(wěn)定的流場,使先到達(dá)交叉處的液滴停留下來,直到對面的液滴運(yùn)動到交叉處并與其發(fā)生接觸,允許液滴到達(dá)交叉的時間有一定偏差。

    Mazutis等[67]改變T型交叉微通道流動方式,讓兩種不同液滴在交叉處相遇形成配對,并在下游較寬的融合通道中共同運(yùn)動,由流動致使液滴發(fā)生融合,其過程如圖7a所示。由于液滴表面活性劑分子作用,兩液滴表面接觸,但并不融合。在融合通道中,由于流場作用使直徑小的乳液滴繞大液滴發(fā)生旋轉(zhuǎn),旋轉(zhuǎn)一定角度后,小液滴被拉長變形,兩液滴中心點(diǎn)間距增大到臨界值時,液滴表面張力變得不平衡,使兩液滴發(fā)生融合,如圖7b所示。此方法不但實(shí)現(xiàn)了單通道液滴旋轉(zhuǎn)融合,還可以實(shí)現(xiàn)液滴與乳液數(shù)量比為1∶1,1∶2和1∶3的融合,如圖7c所示。

    圖7 (a)微流控芯片示意圖;(b)流動致液滴旋轉(zhuǎn)融合;(c)液滴與乳液數(shù)量比為1∶2的融合過程[67]Fig.7 (a)Schematic of the microfluidic chip;(b)Flow-induced droplet rotation and coalescence;(c)Process of 1∶2 coalescence[67]

    3.1.3 利用微刀片結(jié)構(gòu)實(shí)現(xiàn)液滴融合的方法Deng等[68]提出一種新穎的利用微刀片結(jié)構(gòu)實(shí)現(xiàn)液滴融合的方法,通過在兩液滴交匯處加入微刀片來促進(jìn)液滴融合。當(dāng)液滴流過微刀片時,液滴表面液膜平衡狀態(tài)被刀片破壞,使表面活性劑發(fā)生臨時的局部分散,通過“傷口”處誘導(dǎo)液滴融合,并分析了微刀片表面親疏水性和尖端形狀對融合的影響,如圖8所示。結(jié)果表明,親水的帶尖端的刀片更有利于破壞液滴界面穩(wěn)定性,誘發(fā)液滴融合。

    圖8 利用微刀片結(jié)構(gòu)誘導(dǎo)液滴融合過程[68]Fig.8 Droplets fusion induced bymicro-lancets[68]

    3.2 通過表面活性劑實(shí)現(xiàn)微液滴融合的方法

    除上述運(yùn)用通道幾何結(jié)構(gòu)和流動特性促進(jìn)液滴融合方法外,在微通道連續(xù)相中添加表面活性劑以降低界面張力,成為促進(jìn)液滴融合的另一種方法。Baret等[37]綜述了表面活性劑對微流控芯片液滴控制的重要影響,闡述了其在生物系統(tǒng)相融方面和液滴間的分子交換領(lǐng)域中破壞界面穩(wěn)定性的機(jī)理,指出表面活性劑可以溶解在油分子中,通過油分子的擴(kuò)散與另一個液滴融合。Mazutis等[69]研究了一種簡單可靠的裝置,可以在KHz的頻率下達(dá)到液滴對在微流控系統(tǒng)中一對一的融合。這種技術(shù)需要在穩(wěn)定的液滴表面加入活性劑,連接到合適的幾何結(jié)構(gòu)中實(shí)現(xiàn)液滴融合。Fryd等[70]研究了油包水的乳狀液,應(yīng)用一種離子表面活性劑保持液滴穩(wěn)定,應(yīng)用另一種帶相反電荷的離子表面活性劑溶液與其發(fā)生融合。然而,表面活性劑的使用具有一定的雙重性。Liu等[71]研究發(fā)現(xiàn),隨著表面活性劑濃度的增加,液滴尺寸逐漸減小并導(dǎo)致液滴不易發(fā)生融合。此外,Shestopalov等[21]研究發(fā)現(xiàn),當(dāng)表面活性劑重量增加到5% 時,額外的試劑會導(dǎo)致液滴分離,不發(fā)生融合。

    Tullis等[72]提出了一種微通道內(nèi)平行陣列中指定液滴對融合的有效方法。通過減小外部流動油相表面活性劑的含量移走液滴表面的活性劑來達(dá)到液滴成對的融合。研究發(fā)現(xiàn),油相表面不加活性劑可以縮短液滴融合時間。然而,不使用表面活性劑也有其缺點(diǎn),例如發(fā)生意外的融合,并且融合后進(jìn)一步操控液滴受到限制。

    還有一些學(xué)者研究了帶有表面活性劑液滴的被動融合。Mazutis等[69]設(shè)計了一種尺寸明顯不對稱,在表面活性劑作用下實(shí)現(xiàn)液滴融合的微通道,發(fā)現(xiàn)液滴并不是按照預(yù)定的方式進(jìn)行融合。因此,較好的融合方式應(yīng)先使用表面活性劑使液滴穩(wěn)定,后再借用外力使液滴融合。

    3.3 利用外場作用的液滴主動融合方式

    由于通道尺寸微小,液滴運(yùn)動規(guī)律不易精確控制,要實(shí)現(xiàn)液滴捕獲使其融合進(jìn)行反應(yīng),可能需要施加外力作用。利用電場進(jìn)行融合(Electro-coalescence)是目前最常用的主動融合方式[73~76],其原理為通過電誘導(dǎo)使液滴帶相反電荷從而實(shí)現(xiàn)液滴對融合[73]。液滴可以在靜電力作用下相互靠近并融合,或者電場產(chǎn)生的麥克斯韋應(yīng)力使液滴界面變形并破裂,促進(jìn)液滴融合。電場的作用條件可以很寬泛,電壓為1~103V,頻率從DC到數(shù)kHz,均能較好地控制液滴融合。此外,利用電場進(jìn)行液滴融合時,可以將電極嵌入到微通道中或者是距離通道幾mm遠(yuǎn)的位置,電場方向可以平行或垂直于液滴接觸面。

    Wang等[74]利用外加電場來捕獲液滴,從而實(shí)現(xiàn)液滴對的融合,微流控芯片結(jié)構(gòu)如圖9a所示。研究發(fā)現(xiàn),不加電場時液滴會沿著主通道流過,不會進(jìn)入外凹槽。加入電場時主通道的液滴會進(jìn)入微凹槽被捕獲,并與下一個液滴融合,再去掉電場時,液滴又沿主通道運(yùn)動,如圖9b所示。

    圖9 外加電場實(shí)現(xiàn)液滴捕捉融合:(a)微流控芯片結(jié)構(gòu)示意圖;(b)液滴捕捉融合過程[74]Fig.9 Droplet trapping and fusing using electric field:(a)Microfluidic platform diagram;(b)Droplet trapping and fusing process[74]

    Hao等[75]設(shè)計了一種可控的嵌入式局部電場,控制液滴在微通道內(nèi)的融合。該系統(tǒng)利用電潤濕原理,結(jié)合流體動力學(xué)控制,可以單獨(dú)改變液滴的尺寸和生成頻率。通過幾百次的觀察,發(fā)現(xiàn)按需求融合的效率達(dá)到98%。該系統(tǒng)采用了一種交叉電場促進(jìn)液滴融合,關(guān)掉這種電場,此時的表面活性劑穩(wěn)定,即便它們之間相互擠壓也不發(fā)生融合。需要融合時,打開電場開關(guān),液滴表面的不穩(wěn)定導(dǎo)致液滴融合。

    Schoeman等[76]設(shè)計了一種實(shí)現(xiàn)細(xì)胞液滴生成與融合的微通道結(jié)構(gòu),通過雙T通道形成包裹單細(xì)胞的液滴,通過調(diào)節(jié)兩相流速等參數(shù),使兩個包裹不同細(xì)胞的液滴同時到達(dá)Y型微通道交叉處實(shí)現(xiàn)液滴配對,再利用微電極施加電場作用使液滴融合,如圖10所示。

    圖10 單細(xì)胞液滴生成并配對以及融合過程[76]Fig.10 Process of single-cell encapsulation,droplet pairing and fusion[76]

    除了外加電場實(shí)現(xiàn)微液滴融合方法外,其它主動控制的液滴融合方法還包括:磁場[29]、溫度場[30]、表面聲波[31,32]和激光聚焦[33]等。其中,表面聲波不但可以進(jìn)行微液滴操控[26,77],還可以控制微液滴的融合。Fu等[78]在一種壓電基板上應(yīng)用微電子技術(shù)安裝了一個交互傳感器和一個反射器,彼此間由表面聲波互相干擾,實(shí)現(xiàn)微通道內(nèi)液滴融合。研究表明,液滴的大小、距離和信號強(qiáng)度能夠顯著的影響液滴融合?;诼暡ㄔ趦上嗔黧w中的傳播阻力不同會產(chǎn)生聲散射作用力的原理,Muhsincan[79]等制作了表面聲波裝置,如圖11a所示。利用表面聲波將液滴捕捉在矩形微凹槽內(nèi),等待下一個液滴到來發(fā)生融合,可以實(shí)現(xiàn)2~3個液滴的融合。當(dāng)融合后的液滴體積超過一定臨界值時,在連續(xù)相流體作用力下從凹槽流出,如圖11b所示。

    圖11 表面聲波液滴融合:(a)微流控芯片結(jié)構(gòu)示意圖;(b)液滴捕捉及融合過程[79]Fig.11 Dropletmerging using surface acoustic waves:(a)Microfluidic chip design;(b)Droplet trapping and merging[79]

    利用電場進(jìn)行融合一般需要相對較高的電壓,而熱融合的提出克服了被動融合與電場融合的弊端。文獻(xiàn)[80,81]研究了液滴在融合室內(nèi)加熱時的融合過程,實(shí)驗表明流速對臨界溫度具有顯著影響。Peng等[82]研究了在加熱的微室中,由于兩個相鄰液滴間間隙減小導(dǎo)致液滴融合的物理過程,并通過數(shù)學(xué)公式的推導(dǎo),物理方程的求解進(jìn)行了分析。液滴融合還可以通過引入激光束來實(shí)現(xiàn),一些學(xué)者通過渦的捕捉和激光誘導(dǎo)的熱梯度實(shí)現(xiàn)了光作用下的微液滴融合[83,84],利用激光的熱效應(yīng)對兩液滴接觸處進(jìn)行加熱,一方面排出了液滴界面處的表面活性劑分子,一方面誘導(dǎo)產(chǎn)生復(fù)雜的三維流動,這兩方面作用都可以有效地排除液滴間的流體薄層,從而促進(jìn)液滴融合[8]。Jin等[85]提出了一種高精度的光控法來誘導(dǎo)液滴融合。在融合前,利用光的分散性和梯度力擠壓捕獲液滴,并使用分析模型進(jìn)行實(shí)驗,實(shí)驗數(shù)據(jù)和光對液滴行為的影響是一致的。

    3.4 液滴融合過程中的流場特性研究

    利用流體力學(xué)實(shí)驗方法和顯微流場測速技術(shù)(Micro-PIV),開展液滴融合過程流場特性研究,深入揭示液滴融合機(jī)理,成為目前微液滴控制研究重要內(nèi)容。Jin等[86]研究了液滴在直通道和收斂通道中的融合過程。觀測結(jié)果表明,液滴融合要通過接觸、排液、表面融合、滲透或者包裹等步驟,主要依賴于通道的幾何尺寸和各種力的作用。在直通道中,后面的液滴滲透入前面的液滴瞬間合并,如圖12a所示。而在收斂通道中,后面的液滴包圍住前面的液滴時有強(qiáng)烈的渦旋運(yùn)動發(fā)生,如圖12b所示。

    Wang等[66]研究了十字通道交叉處液滴擠壓和旋轉(zhuǎn)兩種狀態(tài)下的液滴內(nèi)部速度矢量場,并對其所受到的粘性剪應(yīng)力進(jìn)行分析。如圖13a所示,兩液滴接觸處邊緣的速度矢量沿著垂直方向,加之剪應(yīng)力作用使液滴間相互擠壓(圖13b),有利于液膜排液,直至融合。圖13c與13a速度矢量場完全不同,且剪應(yīng)力是非對稱的,液滴的移動對液膜排液并無幫助,因此液滴旋轉(zhuǎn)后并未發(fā)生融合[66]。

    Liu等[64]研究了30°的Y型通道中液滴在交叉處發(fā)生碰撞融合和分離融合過程的速度矢量場。研究發(fā)現(xiàn),液滴前部碰撞接觸處的相對剪切運(yùn)動導(dǎo)致液滴內(nèi)部渦的產(chǎn)生,速度最小的部分出現(xiàn)在液滴與通道壁面接觸的位置。由于融合過程中后一個液滴對前一個液滴的擠壓作用,導(dǎo)致前一個液滴在靠近接觸區(qū)的部分速度顯著增加。分離融合過程中,兩液滴相對運(yùn)動,由于液滴表面的突起作用發(fā)生急速融合,其中突起的相接觸附近速度顯著增加。

    Yan等[87]將表面張力引入到納維-斯托克斯方程(N-S方程),數(shù)值模擬研究了液滴在收斂微通道內(nèi)融合的流動特性。液滴的初始速度決定了其在收斂微通道內(nèi)的流動情況,并首次提出了描述收斂通道幾何特性的微通道尺寸參數(shù)“Da”:

    圖12 微通道中(a)直通道和(b)收斂通道液滴融合瞬間的速度矢量場及各種力的作用示意圖[84]Fig.12 Velocity vector and schematics of forces acting on the two droplets at the instant of coalescences in differentmicrochannels:(a)Straight channel and(b)Convergent channel[84]

    圖13 不同狀態(tài)下的液滴速度矢量場Fig.13 Velocity fields of droplet in two different conditions(a)液滴擠壓;(b)粘性剪應(yīng)力示意圖;(c)液滴旋轉(zhuǎn);(d)剪應(yīng)力導(dǎo)致液滴旋轉(zhuǎn)示意圖[66]。(a)Droplets squeezing;(b)Schematic diagram of the viscous shearing force;(c)Droplets rotating;(d)Schematic diagram of droplet rotation driven by shearing forces[66].

    其中,α為通道兩端寬度比率,β為液滴直徑與間距比,γ為通道長度與液滴間距比,根據(jù)數(shù)值模擬結(jié)果歸納了不同尺寸參數(shù)Da與液滴初始速度情況下可能產(chǎn)生的3種融合情況,即完全融合、部分融合和不融合。研究結(jié)果表明,液滴融合位置和通道寬度比呈非線性關(guān)系,增大通道入口寬度、通道長度、液滴直徑有利于促進(jìn)液滴融合。

    4 結(jié)論與展望

    本文綜述了微通道液滴基礎(chǔ)理論及其影響因素,詳細(xì)介紹了被動融合和主動融合兩種微液滴融合方式的研究進(jìn)展[2,8,9],闡述了液滴融合過程中的流體力學(xué)特性。目前,液滴融合技術(shù)的研究已取得了一定的進(jìn)展,同時,對影響液滴融合的因素也進(jìn)行了研究。然而,液滴融合技術(shù)仍面臨著許多挑戰(zhàn):

    微液滴的精確控制方面,特別是對多液滴融合的自動化控制尚存在一定的不確定性。目前研究多集中于兩個液滴融合,多液滴融合技術(shù)研究較少。精確可控的多液滴融合技術(shù)將成為研究的重要方向。

    目前,被動融合方式主要依靠液滴的流動特性實(shí)現(xiàn)融合,融合效率較低。多通道、并行化的高通量液滴融合技術(shù)將成為下一步實(shí)際應(yīng)用研究的重點(diǎn)內(nèi)容。

    主動融合方式需要利用外場作用,需要外部設(shè)備并加工微電極等,增加了微流控芯片系統(tǒng)集成化的難度,限制了大規(guī)模集成化發(fā)展。此外,主動控制方式需要對外場工作條件進(jìn)行優(yōu)化調(diào)節(jié),以實(shí)現(xiàn)融合并減小對生物細(xì)胞的影響,也是研究中的重要內(nèi)容。

    微液滴融合技術(shù)成為了生物細(xì)胞分析、多重乳液的精確控制等領(lǐng)域研究的新手段,其應(yīng)用領(lǐng)域正在不斷擴(kuò)大,但實(shí)際應(yīng)用中尚需克服許多具體難題,將繼續(xù)成為未來研究的重要內(nèi)容。

    此外,液滴操控過程屬于微尺度兩相流動的研究范疇,其中的影響因素復(fù)雜、難以定量測量及控制,特別是由于微通道中液滴界面的可變形性導(dǎo)致的界面張力的非均勻變化和液滴的非線性問題引發(fā)的流體力學(xué)問題,尚需進(jìn)行深入研究。主動融合過程中,外場作用更加復(fù)雜,尚無法進(jìn)行定量化的研究,其融合機(jī)理尚不清楚。

    總之,液滴控制技術(shù)已成為微流控芯片研究的一個重要發(fā)展方向[5~8],微液滴融合已成為其重要的研究內(nèi)容之一[8~13]。微液滴反應(yīng)器具有樣品消耗少、混合速度快、通量高和反應(yīng)條件穩(wěn)定等優(yōu)點(diǎn),利用多學(xué)科交叉研究優(yōu)勢,可實(shí)現(xiàn)微液滴融合在多個領(lǐng)域的應(yīng)用,其應(yīng)用領(lǐng)域也在不斷擴(kuò)展,包括細(xì)胞和亞細(xì)胞結(jié)構(gòu)分析[14~17]、蛋白質(zhì)結(jié)晶[18,19]、納米粒子的合成[20,21]和納米分子自組裝[22,23]以及高通量注射等[24,25]。深入揭示融合過程中的流體力學(xué)問題,特別是主動融合過程中外場力的影響,尚需流體力學(xué)研究者的繼續(xù)努力。微液滴融理論研究尚需繼續(xù)深入,為高精確的液滴控制技術(shù)提供理論指導(dǎo)。這一多學(xué)科交叉領(lǐng)域的發(fā)展,需要微電子機(jī)械學(xué)家、化學(xué)家、生物學(xué)家以及研究微流控的物理學(xué)家的努力與合作。

    1 FANG Zhao-Lun.Microfluidic Chip.Beijing:Science Press,2003:23-78

    方肇倫.微流控分析芯片.北京:科學(xué)出版社,2003:23-78

    2 LIZhan-Hua,WU Jian-Kang,HU Guo-Qing,HU Guo-Hui.Fluid Flow in Microfluidic Chips.Beijing:Science Press, 2012:121-138

    李戰(zhàn)華,吳健康,胡國慶,胡國輝.微流控芯片中的流體流動.北京:科學(xué)出版社,2012:121-138

    3 Stone H A,Stroock A D,Ajdari A.Annu.Rev.Fluid.Mech.,2004,36(1):381-411

    4 LIN Bing-Cheng,QIN Jian-Hua.Graphic Laboratory on a Microfluidic Chip.Beijing:Science Press,2008:2-86

    林炳承,秦建華.圖解微流控芯片實(shí)驗室.北京:科學(xué)出版社,2008:2-86

    5 Whitesides G M.Nature,2006,442(7101):368-373

    6 Joanicot M,Ajdari A.Science,2005,309(5736):887-888

    7 Utada A S,Lorenceau E,Link D R,Kaplan PD,Stone H A,Weitz D A.Science,2005,308(5721):537-541

    8 Charles N B,Francois G,Remi D.Lab Chip,2010,10(16):2032-2045

    9 CHEN Jiu-Sheng,JIANG Jia-Huan.Chinese J.Anal.Chem.,2012,40(8):1293-1300

    陳九生,蔣稼歡.分析化學(xué),2012,40(8):1293-1300

    10 ZHANG Kai,HU Ping,LIANG Qiong-Lin,LUO Guo-An.Chinese J.Anal.Chem.,2008,36(4):556-562

    張凱,胡坪,梁瓊麟,羅國安.分析化學(xué),2008,36(4):556-562

    11 SONGWen-Bin,DONG Chao-Qing,REN Ji-Cun.J.Anal.Sci.,2011,27(1):106-112

    宋文斌,董朝青,任吉存.分析科學(xué)學(xué)報,2011,27(1):106-112

    12 DENG Nan-Nan,WANG Wei,JU Xiao-Jie,XIE Rui,LIU Zhuang,CHU Liang-Yin.Scientia Sinica Chimica,2015, 45(1):7-15

    鄧楠楠,汪偉,巨曉潔,謝銳,劉壯,褚良銀.中國科學(xué):化學(xué),2015,45(1):7-15

    13 XIAO Zhi-Liang,ZHANG Bo.Chinese Journal ofChromatography,2011,29(10):949-956

    肖志良,張博.色譜,2011,29(10):949-956

    14 He M Y,Edgar JS,Jeffries G D M,Lorenz R M,Shelby JP,Chiu D T.Anal.Chem.,2005,77(6):1539-1544

    15 Demirci U,Montesano G.Lab Chip,2007,7(9):1139-1145

    16 Funfak A,Brosing A,Brand M,Kohler JM.Lab Chip,2007,7(9):1132-1138

    17 ShiW,Qin J,Ye N,Lin B.Lab Chip,2008,8(9):1432-1435

    18 Hansen C L,Classen S,Berger JM,Quake SR.J.Am.Chem.Soc.,2006,128(10):3142-3143

    19 Li L,Mustafi D,Fu Q,Tereshko V,Chen D L,Tice JD,Ismagilov R F.Pro.Natl.Acad.Sci.USA,2006,103(51):19243-19248

    20 Hung L H,Choi K M,TsengW Y,Tan Y C,Shea K J,Lee A P.Lab Chip,2006,62(2):174-178

    21 Shestopalov I,Tice JD,Ismagilov R F.Lab Chip,2004,4(4):316-321

    22 Steinbacher JL,Moy RW Y,Price K E,Cummings M A,Roychowdhury C,Buffy JJ,OlbrichtW L,HaafM,McQuade D T.J.Am.Chem.Soc.,2006,128(29):9442-9447

    23 Schabas G,Wang CW,Oskooei A,Yusuf H,Moffitt M G,Sinton D.Langmuir,2008,24(19):10596-10603

    24 Niu X,Gielen F,Edel JB,deMello A J.Nat.Chem.,2011,3(6):437-442

    25 Brouzes E,Medkova M,Savenelli N,Marran D,Twardowski M,Hutchison JB,Rothberg JM,Link D R,Perrimon N, Samuels M L.Proc.Natl.Acad.Sci.USA,2009,106(34),14195-14200

    26 CHEN Xiao-Dong,HU Guo-Qing.Advances in Mechanics,2015,45:201503

    陳曉東,胡國慶.力學(xué)進(jìn)展,2015,45:201503

    27 Simon M G,Lee A P.Microdroplet Technology.New York:Springer,2012:23-50

    28 ZagnoniM,Cooper JM.Lab Chip,2009,9(18):2652-2658

    29 García A A,Egatz-Gómez A,Lindsay SA,Dominguez-Garcia P,Melle S.J.Magn.Magn.Mater.,2007,311(1):238-243

    30 K?hler JM,Henkel T,Grodrian A,Kirner T,Roth M,Martin K,Metze J.Chem.Eng.J.,2004,101(1-3):201-216

    31 Lin SC S,Mao X,Huang T J.Lab Chip,2012,12(16):2766-2770

    32 ZHANG Yue,GAO Ting,HU Chu,HUANG Chang,YU Yi-Qing,ZHANG An-Liang.Piezoelectrics&Acoustooptics, 2012,34(5):676-679

    張悅,高挺,胡楚,黃昶,尉一卿,章安良.壓電與聲光,2012,34(5):676-679

    33 Baroud C N,de Saint Vincent M R,Delville JP.Lab Chip,2007,7(8):1029-1033

    34 DONG Li-Chun,WU Ji-Zhou,REN Gui-Xiang,TAN Shi-Yu,HUANG Da-Fu,LIWen-Ping,HUANG Shao-Jian.Journal ofChongqing University,2011,34(1):134-139

    董立春,吳紀(jì)周,任桂香,譚世語,黃大富,李文平,黃少建.重慶大學(xué)學(xué)報,2011,34(1):134-139

    35 Wu P,Luo Z F,Liu Z F,Li ZD,Chen C,Feng L L,He LQ.Chinese J.Chem.Eng.,2015,23(1):7-14

    36 WANG Jia-Nan,WANG Jia-Jun,FENG Lian-Fang,GU Xue-Ping.J.Chem.Eng.Chin.Univ.,2014,28(2):218-222

    王佳男,王嘉駿,馮連芳,顧雪萍.高?;瘜W(xué)工程學(xué)報,2014,28(2):218-222

    37 Baret JC.Lab Chip,2012,12(3):422-433

    38 Dai B,Leal L.Phys.Fluids,2008,20(4):1569-1579

    39 Gunes D Z,Clain X,Breton O,Mayor G,Burbidge A S.J.Colloid Interface Sci.,2010,343(1):79-86

    40 Chesters A K.Chem.Eng.Res.Des.,1991,69(4):259-270

    41 Christopher G F,Bergstein J,End N B,Poon M,Nguyen C,Anna S L.Lab Chip,2009,9:1102-1109

    42 Borrell M,Yoon Y,Leal L G.Phys.Fluids,2004,16(11):3945-3954

    43 Klaseboer E,Chevaillier JP,Gourdon C,Masbernaty O.J.Colloid Interface Sci.,2000,229(1):472-485

    44 Hu Y T,Pine D J,Leal L G.Phys.Fluids,2000,12(3):484-489

    45 Baldessari F,Homsy G M,Leal LG.J.Colloid Interface Sci.,2007,307(1):188-202

    46 Yang H,Park C C,Hu Y T,Leal L G.Phys.Fluids,2001,13(5):1087-1106

    47 Yoon Y,Borrell M,Park C C,Leal L G.J.Fluid Mech.,2005,525(2):355-379

    48 Baldessari F,Leal L G.Phys.Fluids,2006,18(1):302-309

    49 Rother M A,Davis R H.Phys.Fluids,2001,13(5):1178-1190

    50 Li SZ,Chen R,Wang H,Liao Q,Zhu X,Wang Z B.Appl.Therm.Eng.,2014,64(1-2):129-138

    51 Burkhart B E,Gopalkrishnan P V,Hudson S D,Jamieson A M,Rother M A,Davis R H.Phys.Rev.Lett.,2001, 87(9):098304

    52 Guido S,Simeone M.J.Fluid Mech.,1998,357(2):1-20

    53 Aryafar H,Kavehpour H P.Phys.Fluids,2006,18(7):3453-3475

    54 Mohamed-Kassim Z,Longmire E K.Phys.Fluids,2003,15(11):3263-3273

    55 Xu B,Nguyen N T,Wong TN.Micro and Nanosystems,2011,3:131-136

    56 Tan Y C,Fisher JS,Lee A L,Cristini V,Lee A P.Lab Chip,2004,4(4):292-298

    57 Gunes D Z,Bercy M,Watzke B,Breton O,Burbidge A S.Soft Matter,2013,9(31):7526-7537

    58 Tan Y C,Ho Y L,Lee A P.Microfluid.Nanofluid.,2007,3(4):495-499

    59 Niu X,Gulati S,Edel JB,Demello A J.Lab Chip,2008,8(11):1837-1841

    60 Liu K,Ding H J,Chen Y,Zhao X Z.Microfluid.Nanofluid.,2007,3(2):239-243

    61 Bremond N,Thiam A R,Bibette J.Phys.Rev.Lett.,2008,100(2):024501

    62 Lai A,Bremond N,Stone H A.J.Fluid Mech.,2009,632(10):97-107

    63 Yang L,Wang K,Tan J,Lu Y C,Luo G S.Microfluid.Nanofluid.,2012,12(5):715-722

    64 Liu ZM,Cao R T,Pang Y,Shen F.Exp.Fluids,2015,56(2):1-4

    65 Wang K,Lu Y,Yang L,Luo G S.AICHE Journal,2013,59(2):643-649

    66 Wang K,Lu Y,Tostado C P,Yang L,Luo G S.Chem.Eng.J.,2013,227(7):90-96

    67 Mazutis L,Griffiths A D.Lab Chip,2012,12(10):1800-1806

    68 Deng N N,Sun SX,Wang W,Ju X J,Xie R,Chu L.Lab Chip,2013,13(18):3653-3657

    69 Mazutis L,Baret JC,Griffiths A D.Lab Chip,2009,9(18):2665-2672

    70 Fryd M M,Mason TG.SoftMatter,2014,10(26):4662-4673

    71 Liu K,Qin J.Nanotechnology,2013,24(12):1079-1094

    72 Tullis J,Park C L,Abbyad P.Lab Chip,2014,14(17):3285-3289

    73 Sarrazin F,Prat L,Miceli N D,Cristobal G,Link D R,Weitz D A.Chem.Eng.Sci.,2007,62(4):1042-1048

    74 WangW.Lab Chip,2009,9(11):1504-1506

    75 Hao G,Duits M H G,Mugele F.Int.J.Mol.Sci.,2011,12(4):2572-2597

    76 Schoeman R M,Kemna EW M,Wolbers F,Albert V D B.Electrophoresis,2014,35(2-3):385-392

    77 ZHANG An-Liang,XIA Xing-Hua.Chinese J.Anal.Chem.,2011,39(5):765-769

    章安良,夏興華.分析化學(xué),2011,39(5):765-769

    78 Fu X T,Zha Y,Zhang A L.TELKOMNIKA,2013,11(7):3767-3773

    79 Sesen M,Alan T,Neild A.Lab Chip,2014,14(17):3325-3333

    80 Luong T D,Nguyen N T,Sposito A.Appl.Phys.Lett.,2012,100(25):254105

    81 Xu B,Nguyen N T,Wong T N.Biomicrofluidics,2012,6(1):263-268

    82 Peng C H,Nguyen N T.Phys.Fluids,2013,25(8):840-840

    83 Kotz K T,Gu Y,Gw.F.J.Am.Chem.Soc.,2005,127(16):5736-5737

    84 Lorenz R M,Edgar JS,Jeffries G D M,Zhao Y,McGloin D,Chiu D T.Anal.Chem.,2007,79(1):224-228

    85 Jin H J,Lee K H,Destgeer G,Kang S L,Cho H,Ha B H,Sung H J.Microfluid.Nanofluid.,2014,18(11):1522-1529

    86 Jin B J,Yoo JY.Exp.Fluids,2012,52(1):235-245

    87 Yan Y,Guo D,Luo J,Wen S.Biochip Journal,2013,7(4):325-334

    (Received 22 June 2015;accepted 16 August2015)

    This work was supported by the Natural Science Foundation of Beijing,China(No.7152012)

    Advances in M icro-Droplets Coalescence Using M icrofluidics

    SHEN Feng1,2,LIYi1,LIU Zhao-Miao*1,CAO Ren-Tuo1,WANG Gui-Ren21(College ofMechanical Engineering and Applied Electronics Technology, Beijing University of Technology,Beijing 100124,China)
    2(Department ofMechanical Engineering&Biomedical Engineering Program, University ofSouth Carolina,Columbia,SC,29208,USA)

    Recently,with the development of lab-on-a-chip researches,the use ofmicrofluidics to precisely control the droplet behavior in microchannel has received more and more research attention.This article introduces the basic theory and mechanism of the droplets coalescence,and then carefully reviews the passive methods of droplets coalescence used in current researches by changing microchannel geometry and adding surfactants.The active methods to trigger droplet fusion such as applied electric field,magnetic fields, temperature gradients,surface acoustic wave and focused lasers are also introduced briefly.We also introduce the fluid dynamics of the droplets coalescence and forecast its progress at last to provide useful guidance for microfluidic devices design and wide use of the drop-based microfluidics.

    Microfluidics;Droplets coalescence;Microchannel;Lab-on-a-chip;Review

    10.11895/j.issn.0253-3820.100509

    2015-06-22收稿;2015-08-16接受

    本文系北京市自然科學(xué)基金(No.7152012)、國家留學(xué)基金(No.201406545031)及北京工業(yè)大學(xué)“日新人才”培養(yǎng)計劃項目(No.2015-RX-L02)資助

    猜你喜歡
    排液微流液滴
    液滴間相互碰撞融合與破碎的實(shí)驗研究
    噴淋液滴在空氣環(huán)境下的運(yùn)動特性
    絮凝菌處理頁巖氣壓裂返排液的響應(yīng)面優(yōu)化
    微流控法制備P(NIPA-co-MAA)水凝膠微球及其性能表征
    民用飛機(jī)輔助動力裝置艙排液裝置設(shè)計方法研究
    排液法測物體重力
    微流控芯片在食品安全分析中的應(yīng)用進(jìn)展
    微流控SERS芯片的設(shè)計制備及其在細(xì)菌檢測中的應(yīng)用
    紙芯片微流控技術(shù)的發(fā)展及應(yīng)用
    排液采氣技術(shù)在凝析氣田開發(fā)中的應(yīng)用
    化工管理(2015年5期)2015-12-22 08:25:16
    欧美日韩黄片免| 亚洲精品久久国产高清桃花| 最好的美女福利视频网| 精品欧美国产一区二区三| 99国产精品一区二区蜜桃av| a级毛片免费高清观看在线播放| 亚洲无线在线观看| 免费搜索国产男女视频| 久久久久九九精品影院| 成人无遮挡网站| 桃红色精品国产亚洲av| 久久九九热精品免费| 中文字幕精品亚洲无线码一区| 日韩亚洲欧美综合| 亚洲人成伊人成综合网2020| 91麻豆av在线| 日本黄色片子视频| 国产淫片久久久久久久久 | 午夜老司机福利剧场| 久久久久九九精品影院| av女优亚洲男人天堂| 麻豆av噜噜一区二区三区| 麻豆国产97在线/欧美| 亚洲激情在线av| 久久久久久国产a免费观看| 99国产综合亚洲精品| 18美女黄网站色大片免费观看| 亚洲精华国产精华精| 在现免费观看毛片| 日本在线视频免费播放| 两性午夜刺激爽爽歪歪视频在线观看| 国产一级毛片七仙女欲春2| 精品久久久久久,| 亚洲自偷自拍三级| 91麻豆精品激情在线观看国产| 日韩中字成人| 午夜免费成人在线视频| av国产免费在线观看| 国产三级在线视频| av中文乱码字幕在线| 国产乱人视频| 人人妻,人人澡人人爽秒播| 无遮挡黄片免费观看| 国产主播在线观看一区二区| 欧美在线黄色| 日本精品一区二区三区蜜桃| 亚洲av一区综合| 18禁裸乳无遮挡免费网站照片| 俺也久久电影网| 国产高清有码在线观看视频| 俄罗斯特黄特色一大片| 欧美一区二区亚洲| 美女xxoo啪啪120秒动态图 | 人妻丰满熟妇av一区二区三区| 一本精品99久久精品77| 男女那种视频在线观看| 夜夜躁狠狠躁天天躁| 亚洲第一欧美日韩一区二区三区| 老司机午夜福利在线观看视频| 日本三级黄在线观看| 99热精品在线国产| 18+在线观看网站| 欧美绝顶高潮抽搐喷水| 国产精品98久久久久久宅男小说| 亚洲欧美日韩东京热| 网址你懂的国产日韩在线| 人妻制服诱惑在线中文字幕| 亚洲一区二区三区不卡视频| 天天一区二区日本电影三级| 国产精品一区二区三区四区免费观看 | 成人欧美大片| 在线国产一区二区在线| 人人妻人人看人人澡| 免费看光身美女| 国产视频一区二区在线看| 国产乱人伦免费视频| 国产老妇女一区| 12—13女人毛片做爰片一| 五月伊人婷婷丁香| 五月玫瑰六月丁香| 久久精品综合一区二区三区| 一级毛片久久久久久久久女| 亚洲欧美日韩卡通动漫| 精品一区二区三区视频在线| 18禁黄网站禁片午夜丰满| 99久国产av精品| 亚洲国产精品久久男人天堂| 99精品在免费线老司机午夜| 亚洲,欧美,日韩| 免费在线观看影片大全网站| 97碰自拍视频| av欧美777| 九色成人免费人妻av| 18禁在线播放成人免费| 日本一二三区视频观看| 性欧美人与动物交配| 深夜a级毛片| 国产精品久久久久久精品电影| 欧美日韩乱码在线| 欧美黄色淫秽网站| 欧美不卡视频在线免费观看| 亚洲,欧美精品.| 亚洲欧美日韩高清专用| 小蜜桃在线观看免费完整版高清| 最好的美女福利视频网| av视频在线观看入口| 免费黄网站久久成人精品 | 美女cb高潮喷水在线观看| 精品国内亚洲2022精品成人| 午夜福利在线在线| 岛国在线免费视频观看| 级片在线观看| 悠悠久久av| 噜噜噜噜噜久久久久久91| 久久99热6这里只有精品| 特级一级黄色大片| av欧美777| 午夜免费激情av| 怎么达到女性高潮| 午夜福利在线在线| 偷拍熟女少妇极品色| 91久久精品国产一区二区成人| 亚洲无线在线观看| 性色av乱码一区二区三区2| 女人被狂操c到高潮| 在线观看舔阴道视频| 亚洲 国产 在线| 国产免费av片在线观看野外av| 欧美性猛交╳xxx乱大交人| 1024手机看黄色片| 美女被艹到高潮喷水动态| 悠悠久久av| 精品久久久久久久久av| 国产精品女同一区二区软件 | 亚洲内射少妇av| 国内毛片毛片毛片毛片毛片| 免费人成视频x8x8入口观看| 可以在线观看的亚洲视频| 亚洲美女搞黄在线观看 | 丝袜美腿在线中文| 日韩高清综合在线| 蜜桃亚洲精品一区二区三区| 9191精品国产免费久久| 国产乱人视频| 草草在线视频免费看| 嫁个100分男人电影在线观看| 脱女人内裤的视频| 97人妻精品一区二区三区麻豆| 免费搜索国产男女视频| 国产精品电影一区二区三区| 久久国产乱子伦精品免费另类| 男女视频在线观看网站免费| 精品国内亚洲2022精品成人| 精华霜和精华液先用哪个| 高清毛片免费观看视频网站| 中文字幕av在线有码专区| 搡老岳熟女国产| 国模一区二区三区四区视频| 日韩欧美在线乱码| 久久精品夜夜夜夜夜久久蜜豆| 成人午夜高清在线视频| 男人和女人高潮做爰伦理| 麻豆av噜噜一区二区三区| 内地一区二区视频在线| 久久伊人香网站| 十八禁国产超污无遮挡网站| 国产色婷婷99| 成人高潮视频无遮挡免费网站| 日本黄大片高清| 亚洲 国产 在线| 丰满的人妻完整版| 乱人视频在线观看| 成年人黄色毛片网站| 亚洲人成网站在线播| 真人做人爱边吃奶动态| 日韩国内少妇激情av| 国产精品亚洲一级av第二区| 成人毛片a级毛片在线播放| 亚洲精品成人久久久久久| 国产精品一区二区三区四区久久| 欧美激情在线99| 一进一出抽搐动态| 亚洲精品乱码久久久v下载方式| 波多野结衣巨乳人妻| 亚洲色图av天堂| 亚洲精品在线观看二区| 国产一区二区亚洲精品在线观看| 麻豆av噜噜一区二区三区| 精品人妻视频免费看| 一级毛片久久久久久久久女| 校园春色视频在线观看| 日本熟妇午夜| 一a级毛片在线观看| 草草在线视频免费看| 久久人人爽人人爽人人片va | 男插女下体视频免费在线播放| 久久久久国内视频| 麻豆国产av国片精品| 国产亚洲欧美98| 亚洲成av人片在线播放无| 精品人妻1区二区| 美女xxoo啪啪120秒动态图 | 人妻夜夜爽99麻豆av| 草草在线视频免费看| 午夜福利视频1000在线观看| 精品久久久久久久久久久久久| 波野结衣二区三区在线| 美女免费视频网站| 日本a在线网址| 亚洲av不卡在线观看| 在线十欧美十亚洲十日本专区| 午夜影院日韩av| 变态另类丝袜制服| 日本一二三区视频观看| av在线蜜桃| 人人妻,人人澡人人爽秒播| 亚洲18禁久久av| 色精品久久人妻99蜜桃| 国产单亲对白刺激| 可以在线观看毛片的网站| 一区二区三区激情视频| eeuss影院久久| 欧美成人一区二区免费高清观看| 久久久久久久亚洲中文字幕 | 免费看美女性在线毛片视频| 精品午夜福利视频在线观看一区| 日韩国内少妇激情av| av中文乱码字幕在线| 淫秽高清视频在线观看| 精品午夜福利视频在线观看一区| 亚洲av五月六月丁香网| 久久婷婷人人爽人人干人人爱| 99久久无色码亚洲精品果冻| 国产视频内射| 最近最新中文字幕大全电影3| 国产精品99久久久久久久久| 好男人电影高清在线观看| 精品乱码久久久久久99久播| 老鸭窝网址在线观看| 国产精品伦人一区二区| 91九色精品人成在线观看| 久久国产精品影院| 精品国产三级普通话版| 嫩草影院新地址| 欧美日韩亚洲国产一区二区在线观看| 国产精品国产高清国产av| 精品国产三级普通话版| 97热精品久久久久久| 久久99热6这里只有精品| 亚洲三级黄色毛片| 精品一区二区三区视频在线观看免费| 88av欧美| 男人和女人高潮做爰伦理| 高潮久久久久久久久久久不卡| 亚洲 欧美 日韩 在线 免费| 91麻豆精品激情在线观看国产| 午夜老司机福利剧场| 色尼玛亚洲综合影院| 亚洲av五月六月丁香网| 精品久久久久久久久亚洲 | av在线观看视频网站免费| 欧美黑人欧美精品刺激| 亚洲精华国产精华精| 精品久久久久久成人av| 欧美+日韩+精品| 成人av一区二区三区在线看| 91麻豆精品激情在线观看国产| 最近中文字幕高清免费大全6 | 一二三四社区在线视频社区8| 久久久久久久精品吃奶| 校园春色视频在线观看| 亚洲成人免费电影在线观看| 国产av不卡久久| aaaaa片日本免费| 欧美激情久久久久久爽电影| 久久精品国产清高在天天线| 久久久色成人| 日日摸夜夜添夜夜添av毛片 | 搡老岳熟女国产| 极品教师在线免费播放| 午夜激情欧美在线| 中文资源天堂在线| 一a级毛片在线观看| 夜夜爽天天搞| 久久精品国产99精品国产亚洲性色| 国产亚洲精品综合一区在线观看| 国产一级毛片七仙女欲春2| 在线免费观看的www视频| 99国产精品一区二区蜜桃av| 首页视频小说图片口味搜索| 亚洲一区二区三区不卡视频| 亚洲美女黄片视频| 欧美性猛交黑人性爽| 久久性视频一级片| av在线老鸭窝| 久久香蕉精品热| 一进一出好大好爽视频| 久久午夜福利片| 久久精品国产亚洲av天美| 国产乱人伦免费视频| 美女大奶头视频| 国产精品精品国产色婷婷| 免费在线观看成人毛片| 国产av麻豆久久久久久久| 国产高潮美女av| 免费看a级黄色片| 久久中文看片网| 最近视频中文字幕2019在线8| 18美女黄网站色大片免费观看| 一本精品99久久精品77| 九九在线视频观看精品| 97超视频在线观看视频| 久久久色成人| 观看免费一级毛片| 国内精品美女久久久久久| 国产精品,欧美在线| 久久午夜亚洲精品久久| 蜜桃久久精品国产亚洲av| 国内精品一区二区在线观看| 村上凉子中文字幕在线| 一个人免费在线观看的高清视频| 国产精品美女特级片免费视频播放器| 亚洲在线自拍视频| 日本黄大片高清| 99在线人妻在线中文字幕| 欧美色视频一区免费| 免费在线观看日本一区| 亚洲 国产 在线| 中文字幕熟女人妻在线| 简卡轻食公司| 看免费av毛片| 男人舔奶头视频| 偷拍熟女少妇极品色| 内地一区二区视频在线| .国产精品久久| 亚洲精品日韩av片在线观看| 性色av乱码一区二区三区2| www.熟女人妻精品国产| 熟女电影av网| 变态另类丝袜制服| 色av中文字幕| 99视频精品全部免费 在线| 欧美黄色片欧美黄色片| x7x7x7水蜜桃| 欧美最新免费一区二区三区 | a在线观看视频网站| 小蜜桃在线观看免费完整版高清| 久久国产乱子免费精品| 一区二区三区免费毛片| 中文字幕久久专区| 又粗又爽又猛毛片免费看| 欧美性感艳星| 级片在线观看| 亚洲精品在线观看二区| www.www免费av| 国产 一区 欧美 日韩| 桃色一区二区三区在线观看| 黄色配什么色好看| 成年女人永久免费观看视频| 欧美另类亚洲清纯唯美| 国产免费一级a男人的天堂| 91在线精品国自产拍蜜月| 免费av观看视频| 国产国拍精品亚洲av在线观看| 直男gayav资源| www日本黄色视频网| 午夜精品在线福利| 国产毛片a区久久久久| 波多野结衣高清作品| 欧美另类亚洲清纯唯美| 人人妻,人人澡人人爽秒播| 欧美成人一区二区免费高清观看| 在线国产一区二区在线| 男女做爰动态图高潮gif福利片| 久久性视频一级片| 9191精品国产免费久久| 黄色女人牲交| 一级黄片播放器| .国产精品久久| 成人特级av手机在线观看| 麻豆av噜噜一区二区三区| 国产爱豆传媒在线观看| a在线观看视频网站| 国产欧美日韩精品一区二区| 日本与韩国留学比较| 露出奶头的视频| 桃色一区二区三区在线观看| 中文亚洲av片在线观看爽| 18禁在线播放成人免费| 天堂影院成人在线观看| 国产高潮美女av| 色综合站精品国产| 免费看日本二区| 国产一区二区三区视频了| 国产精品爽爽va在线观看网站| 成人亚洲精品av一区二区| 亚洲 国产 在线| 十八禁人妻一区二区| 国产成+人综合+亚洲专区| 男女床上黄色一级片免费看| 国产精品亚洲av一区麻豆| 亚洲专区中文字幕在线| 俺也久久电影网| .国产精品久久| 成人午夜高清在线视频| 亚洲无线观看免费| 亚洲av日韩精品久久久久久密| 欧美一级a爱片免费观看看| 色综合亚洲欧美另类图片| 国产精品亚洲av一区麻豆| 中文亚洲av片在线观看爽| 日本免费a在线| 黄色丝袜av网址大全| 精品一区二区三区视频在线观看免费| 日本免费a在线| 中文在线观看免费www的网站| 久久国产乱子伦精品免费另类| 久久九九热精品免费| 少妇高潮的动态图| 亚洲无线在线观看| 欧美成人免费av一区二区三区| 亚洲专区中文字幕在线| 最近最新免费中文字幕在线| a级一级毛片免费在线观看| 免费av毛片视频| 夜夜夜夜夜久久久久| 国产乱人伦免费视频| 久久久久久大精品| 日韩欧美在线二视频| 午夜影院日韩av| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩无卡精品| 久久精品国产亚洲av天美| 久久中文看片网| 欧美最黄视频在线播放免费| or卡值多少钱| 激情在线观看视频在线高清| 免费观看的影片在线观看| 国产精品久久久久久久电影| 成人av在线播放网站| 国产av在哪里看| 特大巨黑吊av在线直播| 色吧在线观看| 赤兔流量卡办理| 一个人免费在线观看电影| 看十八女毛片水多多多| 日韩亚洲欧美综合| 最好的美女福利视频网| 亚洲美女搞黄在线观看 | 久久精品91蜜桃| 精品人妻熟女av久视频| 亚洲国产精品sss在线观看| 国产成人啪精品午夜网站| 日韩欧美三级三区| 亚洲中文日韩欧美视频| 亚洲av熟女| 亚洲成人免费电影在线观看| 午夜激情欧美在线| 久久精品国产自在天天线| 国产aⅴ精品一区二区三区波| 99在线人妻在线中文字幕| 国产熟女xx| 99久久无色码亚洲精品果冻| 日本黄大片高清| 国内精品久久久久久久电影| 最新中文字幕久久久久| 又黄又爽又刺激的免费视频.| 国产精品日韩av在线免费观看| 特级一级黄色大片| 十八禁人妻一区二区| 淫秽高清视频在线观看| 日韩国内少妇激情av| 一区二区三区免费毛片| 我的老师免费观看完整版| 99久久精品国产亚洲精品| www日本黄色视频网| 老女人水多毛片| 亚洲国产欧美人成| 亚洲欧美日韩东京热| 高清毛片免费观看视频网站| 此物有八面人人有两片| 亚洲久久久久久中文字幕| 欧美黑人欧美精品刺激| 亚洲无线观看免费| 97热精品久久久久久| 久久久久久九九精品二区国产| 久久久国产成人精品二区| 99久久精品热视频| 亚洲av中文字字幕乱码综合| 最近中文字幕高清免费大全6 | 色5月婷婷丁香| 黄色视频,在线免费观看| 欧美+日韩+精品| 露出奶头的视频| 好看av亚洲va欧美ⅴa在| 一级毛片久久久久久久久女| 美女被艹到高潮喷水动态| 五月玫瑰六月丁香| 天堂网av新在线| 亚洲,欧美精品.| 在线播放国产精品三级| 黄色女人牲交| 一个人观看的视频www高清免费观看| 久久久久免费精品人妻一区二区| 国产三级黄色录像| 欧美潮喷喷水| 国产一区二区在线av高清观看| av专区在线播放| 国产一区二区激情短视频| a级一级毛片免费在线观看| 国产精品日韩av在线免费观看| 脱女人内裤的视频| 欧美精品国产亚洲| 国产探花在线观看一区二区| 黄片小视频在线播放| 亚洲人成网站高清观看| 亚洲自偷自拍三级| 中文在线观看免费www的网站| 黄色女人牲交| 中文在线观看免费www的网站| 精品乱码久久久久久99久播| 怎么达到女性高潮| 变态另类丝袜制服| 小蜜桃在线观看免费完整版高清| 亚洲狠狠婷婷综合久久图片| 好看av亚洲va欧美ⅴa在| 一个人免费在线观看电影| 国产精品乱码一区二三区的特点| 一区二区三区四区激情视频 | 大型黄色视频在线免费观看| 亚洲人与动物交配视频| 精品人妻1区二区| 欧美丝袜亚洲另类 | 一级a爱片免费观看的视频| 精品午夜福利视频在线观看一区| 3wmmmm亚洲av在线观看| 亚洲人成网站在线播| 久久99热这里只有精品18| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美三级三区| 精品国内亚洲2022精品成人| 成年女人永久免费观看视频| 免费av观看视频| 制服丝袜大香蕉在线| 老司机深夜福利视频在线观看| 一区二区三区四区激情视频 | 十八禁人妻一区二区| 国产男靠女视频免费网站| 久久久久精品国产欧美久久久| 噜噜噜噜噜久久久久久91| 欧美日韩综合久久久久久 | 成年女人毛片免费观看观看9| 午夜福利欧美成人| 亚洲天堂国产精品一区在线| 小说图片视频综合网站| 国产三级黄色录像| 国产伦一二天堂av在线观看| 日韩av在线大香蕉| 精品一区二区三区视频在线观看免费| bbb黄色大片| 波多野结衣高清无吗| 久久久久性生活片| 美女免费视频网站| 中文字幕人妻熟人妻熟丝袜美| 亚洲18禁久久av| 午夜日韩欧美国产| 一区二区三区四区激情视频 | 亚洲成人中文字幕在线播放| www.熟女人妻精品国产| 麻豆国产97在线/欧美| 天天一区二区日本电影三级| 国产精华一区二区三区| av中文乱码字幕在线| 久久伊人香网站| 亚洲熟妇中文字幕五十中出| 亚洲无线观看免费| 亚洲欧美日韩卡通动漫| 亚洲电影在线观看av| 久99久视频精品免费| 亚洲国产精品sss在线观看| 久久草成人影院| 国产精品免费一区二区三区在线| a级一级毛片免费在线观看| 51午夜福利影视在线观看| 日韩中文字幕欧美一区二区| 一进一出抽搐动态| 午夜福利欧美成人| 啦啦啦观看免费观看视频高清| 我要看日韩黄色一级片| 一个人看的www免费观看视频| 久久久久九九精品影院| 亚洲精品在线美女| 夜夜夜夜夜久久久久| 高清日韩中文字幕在线| 在线国产一区二区在线| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩综合久久久久久 | 久久天躁狠狠躁夜夜2o2o| 亚洲一区二区三区色噜噜| 国产精品亚洲av一区麻豆| 日韩欧美一区二区三区在线观看| 怎么达到女性高潮| 欧美性感艳星| 波野结衣二区三区在线| 久久性视频一级片| 亚洲精品在线观看二区| 欧美极品一区二区三区四区| 丰满的人妻完整版| 日日夜夜操网爽| 一区二区三区四区激情视频 | 国产蜜桃级精品一区二区三区| 高潮久久久久久久久久久不卡| 亚洲最大成人av| 欧美不卡视频在线免费观看| 久久九九热精品免费|