• <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
    国产亚洲精品av在线| 亚洲中文日韩欧美视频| 日韩欧美在线乱码| 大型黄色视频在线免费观看| 国产av不卡久久| 1024视频免费在线观看| 级片在线观看| 黑人巨大精品欧美一区二区mp4| 啦啦啦免费观看视频1| 最近视频中文字幕2019在线8| 成人手机av| 日韩欧美在线二视频| 超碰成人久久| 国产野战对白在线观看| 国内毛片毛片毛片毛片毛片| 欧美丝袜亚洲另类 | 男女那种视频在线观看| 国产视频一区二区在线看| 男女下面进入的视频免费午夜| 成在线人永久免费视频| aaaaa片日本免费| 欧美最黄视频在线播放免费| 天天躁夜夜躁狠狠躁躁| 听说在线观看完整版免费高清| 在线观看日韩欧美| 99国产精品99久久久久| 欧美日本亚洲视频在线播放| 精品免费久久久久久久清纯| 成人亚洲精品av一区二区| 无限看片的www在线观看| 国产三级中文精品| 国内毛片毛片毛片毛片毛片| av片东京热男人的天堂| 日韩国内少妇激情av| 亚洲真实伦在线观看| 亚洲av日韩精品久久久久久密| 毛片女人毛片| 国产精品电影一区二区三区| 国产高清视频在线播放一区| 亚洲无线在线观看| 日本一区二区免费在线视频| 日韩欧美在线乱码| 亚洲成人中文字幕在线播放| 亚洲黑人精品在线| 国产av一区在线观看免费| 欧美日韩精品网址| 岛国视频午夜一区免费看| 99riav亚洲国产免费| 丰满人妻熟妇乱又伦精品不卡| 亚洲色图av天堂| 国产一区二区在线观看日韩 | 亚洲成人精品中文字幕电影| 亚洲aⅴ乱码一区二区在线播放 | 长腿黑丝高跟| 在线观看免费午夜福利视频| 窝窝影院91人妻| 久久久久久人人人人人| 舔av片在线| 很黄的视频免费| 亚洲国产欧美网| 伊人久久大香线蕉亚洲五| 99精品欧美一区二区三区四区| 青草久久国产| 天堂av国产一区二区熟女人妻 | 又大又爽又粗| 怎么达到女性高潮| 国产精品影院久久| 亚洲七黄色美女视频| 中出人妻视频一区二区| 国产精品免费视频内射| 亚洲精品色激情综合| 欧美乱妇无乱码| 精品国产乱子伦一区二区三区| 日本黄大片高清| www.999成人在线观看| 欧美成狂野欧美在线观看| 黄色a级毛片大全视频| 欧美日韩黄片免| 天堂动漫精品| 国产视频内射| 日韩国内少妇激情av| 中文字幕人成人乱码亚洲影| 这个男人来自地球电影免费观看| 脱女人内裤的视频| 亚洲自拍偷在线| 男女之事视频高清在线观看| 国产片内射在线| 午夜两性在线视频| 亚洲中文av在线| 99热6这里只有精品| 中文字幕av在线有码专区| 精品人妻1区二区| 国内毛片毛片毛片毛片毛片| 亚洲精品国产一区二区精华液| 婷婷精品国产亚洲av| 色综合亚洲欧美另类图片| 桃红色精品国产亚洲av| 人人妻人人看人人澡| www.自偷自拍.com| 日韩av在线大香蕉| 法律面前人人平等表现在哪些方面| 亚洲精品一卡2卡三卡4卡5卡| 久久婷婷人人爽人人干人人爱| www.www免费av| 美女 人体艺术 gogo| 人人妻人人看人人澡| 国产亚洲精品第一综合不卡| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲精品久久久久5区| 久久久久国产一级毛片高清牌| 亚洲国产精品合色在线| 深夜精品福利| 亚洲欧美激情综合另类| 欧美不卡视频在线免费观看 | 亚洲精品av麻豆狂野| 好男人电影高清在线观看| 极品教师在线免费播放| 99热6这里只有精品| 亚洲自偷自拍图片 自拍| 国产不卡一卡二| 1024视频免费在线观看| 中文字幕av在线有码专区| 久久人人精品亚洲av| 岛国在线免费视频观看| 久久久精品大字幕| 精品人妻1区二区| 露出奶头的视频| 精品无人区乱码1区二区| 欧美 亚洲 国产 日韩一| 亚洲国产看品久久| 久久精品成人免费网站| 日日夜夜操网爽| 亚洲精品在线观看二区| 啪啪无遮挡十八禁网站| 丰满人妻熟妇乱又伦精品不卡| 国产视频一区二区在线看| 丰满人妻一区二区三区视频av | 黄色视频不卡| 成年女人毛片免费观看观看9| 欧美精品亚洲一区二区| 国产精品日韩av在线免费观看| 日本撒尿小便嘘嘘汇集6| 少妇熟女aⅴ在线视频| 国产欧美日韩一区二区三| 久久午夜亚洲精品久久| 两性午夜刺激爽爽歪歪视频在线观看 | av在线天堂中文字幕| 精品午夜福利视频在线观看一区| 两个人的视频大全免费| 大型av网站在线播放| 欧美大码av| 大型黄色视频在线免费观看| 精品久久久久久久末码| 9191精品国产免费久久| 一进一出好大好爽视频| 最新在线观看一区二区三区| 亚洲中文字幕日韩| 在线观看美女被高潮喷水网站 | 日韩欧美精品v在线| 母亲3免费完整高清在线观看| 观看免费一级毛片| 亚洲美女黄片视频| 1024手机看黄色片| 色老头精品视频在线观看| 在线免费观看的www视频| 无人区码免费观看不卡| 日本成人三级电影网站| or卡值多少钱| 久久国产精品人妻蜜桃| 亚洲国产精品999在线| 日本 欧美在线| 成人一区二区视频在线观看| 久久精品夜夜夜夜夜久久蜜豆 | 成年版毛片免费区| 国产一区二区激情短视频| 免费看日本二区| 麻豆成人av在线观看| 床上黄色一级片| 国产单亲对白刺激| 亚洲精品中文字幕在线视频| 两个人视频免费观看高清| 九色国产91popny在线| 波多野结衣高清无吗| 岛国在线免费视频观看| 99国产精品99久久久久| 免费电影在线观看免费观看| 精品欧美一区二区三区在线| 国产爱豆传媒在线观看 | 亚洲国产高清在线一区二区三| 久久久久久大精品| 成人18禁高潮啪啪吃奶动态图| 国产精品影院久久| 黄片大片在线免费观看| 成人高潮视频无遮挡免费网站| 国产91精品成人一区二区三区| 99久久国产精品久久久| 全区人妻精品视频| 狠狠狠狠99中文字幕| 国产精品一区二区免费欧美| 天天躁狠狠躁夜夜躁狠狠躁| 日韩高清综合在线| 欧美 亚洲 国产 日韩一| 亚洲av美国av| 亚洲男人的天堂狠狠| 欧美一区二区精品小视频在线| 欧美成人性av电影在线观看| 99国产精品99久久久久| 岛国视频午夜一区免费看| 成年免费大片在线观看| 黄色 视频免费看| 免费观看人在逋| 两个人看的免费小视频| 美女午夜性视频免费| 给我免费播放毛片高清在线观看| 国产精品国产高清国产av| 亚洲18禁久久av| 性欧美人与动物交配| 亚洲精品在线美女| 亚洲五月天丁香| 亚洲中文日韩欧美视频| 亚洲18禁久久av| 91麻豆av在线| 亚洲全国av大片| av超薄肉色丝袜交足视频| 99国产精品99久久久久| 99精品欧美一区二区三区四区| 国产单亲对白刺激| 久久欧美精品欧美久久欧美| 成年人黄色毛片网站| 99国产精品一区二区三区| 亚洲精品美女久久久久99蜜臀| 欧美中文日本在线观看视频| 999久久久国产精品视频| 亚洲国产精品成人综合色| 国产激情欧美一区二区| 18禁国产床啪视频网站| 欧美极品一区二区三区四区| 国产三级中文精品| 麻豆成人午夜福利视频| 亚洲熟女毛片儿| 在线观看美女被高潮喷水网站 | 国产精品久久久人人做人人爽| 99久久国产精品久久久| 老汉色∧v一级毛片| 人成视频在线观看免费观看| 国产亚洲精品一区二区www| 精品乱码久久久久久99久播| 法律面前人人平等表现在哪些方面| 国产三级中文精品| 99精品久久久久人妻精品| 欧美激情久久久久久爽电影| 三级国产精品欧美在线观看 | 亚洲avbb在线观看| 亚洲精品一区av在线观看| 毛片女人毛片| 国产亚洲av高清不卡| 日本免费a在线| 日日夜夜操网爽| 精品国产美女av久久久久小说| 两人在一起打扑克的视频| 国产99白浆流出| 国产黄色小视频在线观看| 欧美最黄视频在线播放免费| 日韩 欧美 亚洲 中文字幕| 成人av在线播放网站| 精品福利观看| 欧美成人性av电影在线观看| 免费在线观看日本一区| 国产三级中文精品| 日韩欧美三级三区| 国产亚洲精品第一综合不卡| 中文字幕熟女人妻在线| 亚洲成人久久爱视频| 两性午夜刺激爽爽歪歪视频在线观看 | 19禁男女啪啪无遮挡网站| 免费在线观看黄色视频的| 韩国av一区二区三区四区| 999精品在线视频| 久久久久久久午夜电影| 午夜福利成人在线免费观看| 久久久久久人人人人人| 91av网站免费观看| 最近最新中文字幕大全电影3| 欧美又色又爽又黄视频| 麻豆av在线久日| 精品久久久久久久人妻蜜臀av| 亚洲国产精品999在线| 性色av乱码一区二区三区2| 91大片在线观看| 欧美日本亚洲视频在线播放| 国产伦人伦偷精品视频| 亚洲欧美精品综合久久99| 国产乱人伦免费视频| 男女做爰动态图高潮gif福利片| 久久精品成人免费网站| 亚洲人成网站高清观看| 一区二区三区高清视频在线| 精品久久久久久久毛片微露脸| 成年免费大片在线观看| 成人国产综合亚洲| 免费电影在线观看免费观看| 亚洲人成电影免费在线| 搡老熟女国产l中国老女人| 一a级毛片在线观看| 午夜免费激情av| 亚洲av中文字字幕乱码综合| 国产亚洲精品第一综合不卡| 日本 欧美在线| 欧美一级a爱片免费观看看 | 日韩有码中文字幕| 久久婷婷成人综合色麻豆| 日韩大码丰满熟妇| 精品少妇一区二区三区视频日本电影| 麻豆av在线久日| 又黄又爽又免费观看的视频| 两个人视频免费观看高清| 国产熟女xx| 国产一区二区激情短视频| 久久热在线av| 日韩欧美三级三区| 亚洲av成人精品一区久久| 日韩欧美一区二区三区在线观看| 亚洲男人的天堂狠狠| 成人三级黄色视频| 久久草成人影院| av中文乱码字幕在线| 日本一区二区免费在线视频| 欧美+亚洲+日韩+国产| 别揉我奶头~嗯~啊~动态视频| 欧美成人午夜精品| 在线观看午夜福利视频| 欧美成人午夜精品| 禁无遮挡网站| 999精品在线视频| 在线观看免费日韩欧美大片| 亚洲电影在线观看av| 久久精品成人免费网站| 熟女电影av网| 午夜福利在线观看吧| 精品福利观看| 1024视频免费在线观看| 欧美成人午夜精品| 成人手机av| 一区二区三区激情视频| 亚洲成av人片免费观看| 悠悠久久av| 最近在线观看免费完整版| 国产一级毛片七仙女欲春2| av超薄肉色丝袜交足视频| 国产精品国产高清国产av| а√天堂www在线а√下载| 欧美高清成人免费视频www| 亚洲国产欧美人成| 中文字幕av在线有码专区| 亚洲国产欧美人成| 亚洲av片天天在线观看| 俄罗斯特黄特色一大片| 欧美 亚洲 国产 日韩一| 欧美日韩乱码在线| 久久久久国产精品人妻aⅴ院| 国产伦人伦偷精品视频| 又黄又爽又免费观看的视频| 亚洲五月天丁香| 国产精品99久久99久久久不卡| 国产成人精品久久二区二区91| 亚洲最大成人中文| 99国产精品一区二区蜜桃av| 欧美精品亚洲一区二区| 久久久久九九精品影院| 全区人妻精品视频| 久9热在线精品视频| 999精品在线视频| 国产欧美日韩一区二区精品| 又粗又爽又猛毛片免费看| 18禁黄网站禁片午夜丰满| 午夜福利高清视频| 国产成人av教育| 成人国产一区最新在线观看| 两人在一起打扑克的视频| 2021天堂中文幕一二区在线观| 国产区一区二久久| av国产免费在线观看| 欧美成人免费av一区二区三区| 国产一区二区在线观看日韩 | 夜夜爽天天搞| 小说图片视频综合网站| 久久九九热精品免费| 亚洲中文字幕一区二区三区有码在线看 | 久久香蕉激情| 日韩精品青青久久久久久| 美女 人体艺术 gogo| 午夜福利在线在线| 中出人妻视频一区二区| 久久精品国产亚洲av香蕉五月| 国产片内射在线| 亚洲成人国产一区在线观看| 国内精品久久久久精免费| 天堂√8在线中文| 亚洲精品av麻豆狂野| 久久久久久久午夜电影| 国产伦在线观看视频一区| 性色av乱码一区二区三区2| 欧美日韩瑟瑟在线播放| 午夜福利在线观看吧| 久久香蕉精品热| 男插女下体视频免费在线播放| 欧美日韩黄片免| 亚洲一卡2卡3卡4卡5卡精品中文| 久久欧美精品欧美久久欧美| 老司机午夜福利在线观看视频| 啦啦啦免费观看视频1| 欧美中文综合在线视频| 亚洲av美国av| 波多野结衣巨乳人妻| 成在线人永久免费视频| 又黄又爽又免费观看的视频| 99riav亚洲国产免费| 亚洲中文字幕一区二区三区有码在线看 | 久久久精品国产亚洲av高清涩受| 欧美日本亚洲视频在线播放| 美女扒开内裤让男人捅视频| 宅男免费午夜| 可以在线观看毛片的网站| 国产精品亚洲av一区麻豆| 亚洲国产精品999在线| 熟女少妇亚洲综合色aaa.| 国产高清激情床上av| 级片在线观看| 黄色视频不卡| 亚洲专区字幕在线| 亚洲国产精品成人综合色| 女同久久另类99精品国产91| 久久九九热精品免费| 国产精品98久久久久久宅男小说| 国产精华一区二区三区| 伊人久久大香线蕉亚洲五| 欧美丝袜亚洲另类 | 国产成人精品久久二区二区91| 国产精品一区二区精品视频观看| 两性夫妻黄色片| 999久久久国产精品视频| 国产亚洲欧美98| 国产一区二区在线观看日韩 | 日本一二三区视频观看| 欧美一区二区精品小视频在线| 很黄的视频免费| av国产免费在线观看| 国内久久婷婷六月综合欲色啪| 国产成年人精品一区二区| 日韩精品青青久久久久久| 黄色片一级片一级黄色片| 国产亚洲精品第一综合不卡| 日韩欧美在线二视频| 88av欧美| 亚洲国产日韩欧美精品在线观看 | 亚洲av日韩精品久久久久久密| 国产伦在线观看视频一区| 久久久久久久精品吃奶| 日本黄色视频三级网站网址| 狂野欧美白嫩少妇大欣赏| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久久久久久久久| 国产精品一区二区三区四区久久| 在线观看舔阴道视频| 亚洲精品中文字幕一二三四区| 在线观看66精品国产| 麻豆成人午夜福利视频| 欧美av亚洲av综合av国产av| 制服人妻中文乱码| 亚洲精品国产精品久久久不卡| 国产成+人综合+亚洲专区| 老鸭窝网址在线观看| 全区人妻精品视频| 免费在线观看亚洲国产| 久久精品aⅴ一区二区三区四区| 亚洲欧洲精品一区二区精品久久久| 中文字幕人妻丝袜一区二区| 91九色精品人成在线观看| 成人午夜高清在线视频| 动漫黄色视频在线观看| 老司机午夜十八禁免费视频| 久久精品国产亚洲av香蕉五月| 怎么达到女性高潮| 午夜亚洲福利在线播放| 亚洲成人免费电影在线观看| 久久久精品国产亚洲av高清涩受| 男女下面进入的视频免费午夜| 国产精品一区二区免费欧美| 99国产综合亚洲精品| 午夜激情av网站| 亚洲欧美日韩高清专用| 神马国产精品三级电影在线观看 | 日本精品一区二区三区蜜桃| 国产亚洲精品一区二区www| 精品久久久久久久末码| 亚洲 国产 在线| 欧美极品一区二区三区四区| 长腿黑丝高跟| 日本熟妇午夜| 欧美成狂野欧美在线观看| 日本五十路高清| 国产久久久一区二区三区| 变态另类丝袜制服| 久久久久国产精品人妻aⅴ院| 成人手机av| 人人妻人人澡欧美一区二区| 国产片内射在线| videosex国产| 亚洲色图av天堂| 成人av在线播放网站| 一级a爱片免费观看的视频| 国产高清videossex| 久久久久久大精品| 亚洲欧美精品综合一区二区三区| 别揉我奶头~嗯~啊~动态视频| 国产成人系列免费观看| 91老司机精品| 黄频高清免费视频| 在线a可以看的网站| 18美女黄网站色大片免费观看| 91在线观看av| 一边摸一边做爽爽视频免费| 色播亚洲综合网| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩乱码在线| 日韩国内少妇激情av| 岛国在线免费视频观看| 亚洲中文日韩欧美视频| 亚洲精品粉嫩美女一区| 人人妻人人看人人澡| 日本 欧美在线| 好男人电影高清在线观看| 欧美中文综合在线视频| 琪琪午夜伦伦电影理论片6080| 少妇熟女aⅴ在线视频| 最近在线观看免费完整版| 巨乳人妻的诱惑在线观看| 精品国产乱子伦一区二区三区| 亚洲精品中文字幕一二三四区| 色哟哟哟哟哟哟| 亚洲欧美精品综合久久99| 久久精品影院6| 国内精品一区二区在线观看| 国产熟女午夜一区二区三区| 中文字幕人成人乱码亚洲影| 亚洲精品av麻豆狂野| 99久久无色码亚洲精品果冻| 男人舔女人的私密视频| svipshipincom国产片| a级毛片a级免费在线| 女同久久另类99精品国产91| 最新在线观看一区二区三区| 欧美黑人欧美精品刺激| 欧美日韩瑟瑟在线播放| 日本一区二区免费在线视频| 他把我摸到了高潮在线观看| 香蕉久久夜色| 香蕉国产在线看| 巨乳人妻的诱惑在线观看| 99精品在免费线老司机午夜| 啪啪无遮挡十八禁网站| 欧美又色又爽又黄视频| 观看免费一级毛片| 一夜夜www| 一进一出抽搐gif免费好疼| www.精华液| 一个人免费在线观看的高清视频| 亚洲熟女毛片儿| 国产单亲对白刺激| 19禁男女啪啪无遮挡网站| 国产精品98久久久久久宅男小说| aaaaa片日本免费| 免费在线观看完整版高清| 国产免费男女视频| 99精品久久久久人妻精品| 国产高清激情床上av| netflix在线观看网站| 国产亚洲av嫩草精品影院| 婷婷精品国产亚洲av在线| 观看免费一级毛片| а√天堂www在线а√下载| 国内少妇人妻偷人精品xxx网站 | 久热爱精品视频在线9| 97超级碰碰碰精品色视频在线观看| 亚洲av美国av| 一本大道久久a久久精品| 一级作爱视频免费观看| 国产精品久久久人人做人人爽| 女人爽到高潮嗷嗷叫在线视频| 制服诱惑二区| 久久精品aⅴ一区二区三区四区| 亚洲成av人片在线播放无| 香蕉国产在线看| 一区福利在线观看| 精品高清国产在线一区| 在线观看免费日韩欧美大片| 国产又色又爽无遮挡免费看| 亚洲av成人一区二区三| 青草久久国产| 一区福利在线观看| 欧美乱码精品一区二区三区| 国产欧美日韩一区二区精品| 久久99热这里只有精品18| 欧美乱色亚洲激情| 亚洲人成电影免费在线| 1024香蕉在线观看| www日本在线高清视频| 国产成人一区二区三区免费视频网站| 宅男免费午夜| 亚洲一区中文字幕在线| 精品无人区乱码1区二区| 国产精品一区二区三区四区久久| 777久久人妻少妇嫩草av网站| 亚洲一区高清亚洲精品| 国产一区二区在线av高清观看|