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

    微通道內(nèi)液-液泰勒流傳熱的計(jì)算流體力學(xué)模擬

    2017-06-09 18:13:50李婷許松林
    化工進(jìn)展 2017年6期
    關(guān)鍵詞:界面

    李婷,許松林

    ?

    微通道內(nèi)液-液泰勒流傳熱的計(jì)算流體力學(xué)模擬

    李婷,許松林

    (天津大學(xué)化工學(xué)院,系統(tǒng)生物工程教育部重點(diǎn)實(shí)驗(yàn)室,天津 300072)

    微通道內(nèi)的液-液兩相流型在低流速時(shí)以泰勒流為主,本文使用計(jì)算流體力學(xué)方法,對微通道內(nèi)液-液泰勒流的傳熱特性進(jìn)行了研究。首先考察了分散相流速、物系和管徑對微通道壁面溫度的影響。結(jié)果表明:分散相流速和物系熱導(dǎo)率增大使得微通道壁面溫度降低,管徑的改變對微通道壁面溫度影響較小。然后針對當(dāng)量直徑為0.5mm的微通道內(nèi)工作介質(zhì)為甲苯和水的兩相泰勒流模型,考察了微通道壁面剪切力、界面渦度和努塞爾數(shù)對壁面和內(nèi)部溫度的影響,并與文獻(xiàn)中氣-液兩相泰勒流的傳熱性質(zhì)作了比較。結(jié)果表明:壁面剪切力和界面渦度對管壁和界面溫度的波動性變化有一定影響,壁面剪切力和界面渦度的波峰往往出現(xiàn)在溫度的波峰附近,并且有一定的時(shí)間滯后性。甲苯-水兩相泰勒流動下的壁面努塞爾數(shù)比氣-液兩相流大得多,液彈單元的平均努塞爾數(shù)是相同條件下單相流體流動的1.3倍。

    微通道;泰勒流;傳熱;計(jì)算流體力學(xué)

    微通道,也稱為微通道換熱器,是一類結(jié)構(gòu)緊湊、輕巧、高效的換熱器,在航空航天、現(xiàn)代醫(yī)療、化學(xué)生物工程等諸多領(lǐng)域都具有廣闊的應(yīng)用前 景[1]。微通道換熱器的當(dāng)量直徑通常在10~1000μm之間,與常規(guī)尺度通道相比,微通道內(nèi)流體間的界面張力和壁面效應(yīng)起主要作用,而重力、浮力等作用很小,甚至有時(shí)可以忽略不計(jì),因而其流體特性與常規(guī)尺度通道內(nèi)的流體具有非常大的區(qū)別[2-4]。在低流速條件下,微通道內(nèi)兩相流動的形式通常為泰勒流,即一連串子彈形狀的分散相被連續(xù)相的液柱隔開,均勻分布在微通道中,而且在彈狀流與壁面之間存在著一層薄薄的液膜[5]。低流速下微通道內(nèi)流體的流動形式通常為層流,與常規(guī)通道中主要利用湍流擾動的方式強(qiáng)化換熱不同,泰勒流型由于存在著周期性的兩相界面,能夠使界面兩側(cè)流體產(chǎn)生擾動,繼而對傳熱傳質(zhì)過程產(chǎn)生強(qiáng)化,因而微通道內(nèi)穩(wěn)定周期性的泰勒流動可以滿足化工以及芯片內(nèi)的高通量散熱要求[6]。

    對于微通道內(nèi)兩相泰勒流動和傳熱已有較多的實(shí)驗(yàn)和模擬研究[7-12],主要針對氣-液兩相體系流動形態(tài)、液膜厚度、通道壓降和孔隙率等方面進(jìn)行分析。GUPTA等[12]通過模擬微通道中氣-液兩相泰勒流的流動和傳熱情況,發(fā)現(xiàn)只有大約0.13%的壁面熱流量從氣泡的突出部分流入連續(xù)相,因而泰勒流中的氣體在熱量傳遞過程中幾乎不起作用。相對于氣-液體系,液-液兩相流由于本身熱導(dǎo)率大,在增強(qiáng)傳熱效果方面有著更高的潛能,對于具有超強(qiáng)導(dǎo)熱性的金屬流體和納米流體,這一優(yōu)勢更為顯著[13]。目前雖有一些關(guān)于液-液兩相泰勒流流動的實(shí)驗(yàn)和模擬研究[14-19],但對于微通道中液-液兩相流動和傳熱特性的分析還不是很清楚。JOVANOVI? 等[14]通過實(shí)驗(yàn)證明了滯留膜模型比移動膜模型更能準(zhǔn)確預(yù)測泰勒流的壓降情況。SALIM等[15]通過實(shí)驗(yàn)研究了微通道中油-水兩相流的壓降、流型和潤濕性,他們在最初充滿水的石英和玻璃這兩種材質(zhì)的微通道中分別觀察到了不同的流型,這表明界面張力對兩相流有著重要影響。FOROUGHI等[16]通過實(shí)驗(yàn)將兩相流流型與壓降作了關(guān)聯(lián)。ASTHANA 等[17]通過實(shí)驗(yàn)對比兩相流和單相流的傳熱效果,發(fā)現(xiàn)泰勒流動下的努塞爾數(shù)是純水的4倍。URBANT等[18]通過數(shù)值模擬礦物油夾帶水滴的兩相流動,發(fā)現(xiàn)增大液滴的尺寸能夠顯著增加傳熱效率和努塞爾數(shù)。FISCHER等[19]將不溶于水的流體作為攜帶液進(jìn)行模擬,考慮馬蘭戈尼效應(yīng)和膠體效應(yīng)的影響,發(fā)現(xiàn)采用懸浮液作為第二相能夠有效提高傳熱效果,并且沒有很大的壓力損失。

    到目前為止,對微通道內(nèi)液-液兩相流的傳熱規(guī)律和機(jī)理研究還有待進(jìn)一步開展。本文作者采用計(jì)算流體力學(xué)方法,構(gòu)建了微通道內(nèi)液-液兩相流的流動和傳熱模型,首先分析了不同因素對微通道壁面溫度分布的影響,然后研究了壁面剪切力和界面渦度對微通道內(nèi)泰勒流傳熱的影響,以及界面處局部努塞爾數(shù)隨管長的變化,對泰勒流的傳熱作用作出了評價(jià)。

    1?模型和方法

    1.1?物理模型

    在以往的研究中,2D和3D網(wǎng)格模型均有使用,其中3D模型主要是用來模擬橫截面為正方形或矩形的微通道內(nèi)泰勒流的流動情況[20-21],而對于截面為圓形的管道中的泰勒流動,通常采用2D模型以減少運(yùn)算時(shí)間[5],采用Axis邊界條件利用二維模擬三維軸對稱圓管中的泰勒流動。

    本研究采用Gambit軟件構(gòu)造二維軸對稱幾何結(jié)構(gòu)的微通道物理模型。如圖1所示,微通道直徑=0.5mm,長度為10。該模型由入口區(qū)、通道區(qū)和出口區(qū)組成。在入口區(qū),甲苯由管道中心半徑為0.3的圓形入口流入,而水以包圍著甲苯入口的環(huán)形入口流入,甲苯和水在微通道內(nèi)遵循特定機(jī)理形成兩相泰勒流,并從微通道右側(cè)出口流出。

    1.2?控制方程

    多相流的模擬可以采用不同模型,如水平集、流體體積函數(shù)、標(biāo)記粒子、離散格子玻爾茲曼、界面追蹤等[22]。本研究中使用的流體體積函數(shù)(volume of fluid,VOF)模型是最常用的隱式自由表面構(gòu)建方法,該模型具有追蹤兩相界面移動變化的功能,并且具有易實(shí)現(xiàn)、復(fù)雜度小、精度高的特點(diǎn)[23]。

    VOF模型關(guān)于多相流模擬計(jì)算的控制方程如式(1)~式(4)。

    連續(xù)性方程

    動量方程

    (2)

    能量守恒方程

    顏色函數(shù)方程

    (4)

    式中,為時(shí)間;為速度向量;為壓力;為內(nèi)能;為流體的焓;為溫度;為顏色函數(shù),在不同的界面捕捉方法中具有不同的意義[12]。

    本研究中所使用的微通道當(dāng)量直徑屬于亞毫米級,通道內(nèi)液-液兩相間的界面張力和壁面吸附力起主要作用,在應(yīng)用VOF模型時(shí),需要定義動量方程中由界面張力和壁面黏附作用產(chǎn)生的源項(xiàng)[24]。Fluent中界面張力模型是BRACKBILL等[25]提出的連續(xù)界面力(continuum surface force,CSF)模型,兩相界面的界面張力影響包括在這個(gè)模型中。使用這個(gè)模型,VOF計(jì)算中附加的界面張力導(dǎo)致了動量方程中的源項(xiàng),即如式(5)。

    SV=(–int)(5)

    式中,SV為界面張力;為界面張力系數(shù);為界面曲率;為狄拉克函數(shù);為徑向位置向量;為單位法向量。

    1.3?網(wǎng)格劃分策略和相關(guān)參數(shù)設(shè)定

    本實(shí)驗(yàn)使用有限體積元軟件Fluent16.2模擬二維軸對稱結(jié)構(gòu)微通道中的非穩(wěn)態(tài)兩相泰勒流。物理模型計(jì)算區(qū)域采用四邊形結(jié)構(gòu)網(wǎng)格,模擬過程采用離散型非穩(wěn)態(tài)求解,動量方程采用二階迎風(fēng)離散法,壓力項(xiàng)關(guān)聯(lián)采用交錯壓力格式(PRESTO?。?,流場采用壓力和速度耦合的PISO方法求解,流體體積函數(shù)方程采用隱式時(shí)間格式離散求解,對壁面附近區(qū)域采用幾何重構(gòu)法處理,時(shí)間步長取1×10–6s,Courant數(shù)為0.25能夠使收斂比較平穩(wěn)的進(jìn)行。設(shè)定甲苯和水兩個(gè)入口皆為速度入口,通道出口為壓力出口,表壓為0。管壁設(shè)為無滑移壁面條件,壁面接觸角為90°。入口溫度為300K,管壁恒定熱通量32000W/m2,對稱軸上所有變量的法向梯度均為0。水相進(jìn)口處液相體積分?jǐn)?shù)wa=1,即該處全為水,甲苯相進(jìn)口處wa=0,即該處全為甲苯。30℃時(shí)甲苯在水中的溶解度僅為0.57g/L,可視為不溶于水,因此甲苯和水可以在微通道內(nèi)形成相界面清晰的兩相流。各相介質(zhì)的物理性質(zhì)如表1 所示[26-27]。

    表1?模擬計(jì)算中各相物理參數(shù)(T=300K)

    在進(jìn)行網(wǎng)格獨(dú)立性分析時(shí),采用了4種網(wǎng)格劃分方式,其中前3種網(wǎng)格劃分為均勻的正方形網(wǎng)格,邊長分別為10μm、5μm、2.5μm,第4種網(wǎng)格在通道主體部分劃分為邊長5μm的正方形網(wǎng)格,對靠近壁面處12μm厚的液膜區(qū)域進(jìn)行網(wǎng)格加密,劃分為5μm×2.4μm的矩形網(wǎng)格。在相同操作條件下使用以上4種網(wǎng)格對甲苯-水兩相流進(jìn)行模擬,結(jié)果發(fā)現(xiàn)網(wǎng)格邊長在2.5μm以上時(shí)無法捕捉到壁面附近的液膜,這是由于泰勒流液膜的厚度非常小,必須保證微通道壁面附近的網(wǎng)格具有足夠高的精度,才能在模擬結(jié)果中捕捉到液膜。將網(wǎng)格3和網(wǎng)格4的模擬結(jié)果與文獻(xiàn)[14]實(shí)驗(yàn)結(jié)果進(jìn)行對比,結(jié)果如圖2所示??梢钥闯?,當(dāng)泰勒流液彈速度變化時(shí),模擬結(jié)果中液彈單元長度、液膜厚度與文獻(xiàn)實(shí)驗(yàn)結(jié)果的變化趨勢相同,差別很小,說明本文模型的合理性。因此在滿足模擬要求的前提下,為了減少計(jì)算時(shí)間,提高計(jì)算效率,采用網(wǎng)格4的劃分方式。

    2 模擬結(jié)果與討論

    2.1?操作條件與結(jié)構(gòu)參數(shù)對泰勒流傳熱的影響

    圖3考察了在其他條件相同的情況下,改變某一操作條件或結(jié)構(gòu)參數(shù)對微通道壁面溫度的影響,圖中均為同一時(shí)刻微通道內(nèi)壁面溫度隨管長的變化。從圖3(a)中可以看出,在同一體系中,隨著分散相速率增大,微通道內(nèi)壁面溫度波動性減小,整體溫度降低,這是因?yàn)殡S著分散相速率的增大,相同時(shí)間內(nèi)微通道中產(chǎn)生的彈狀流增加,兩相交界面面積的增加使得通道內(nèi)流體擾動增強(qiáng),溫度擴(kuò)散加快,于是微通道壁面處的溫度相對降低。從圖3(b)中可以看出,在相同操作條件下,四氯化碳-水體系的整體壁溫高于甲苯-水體系,苯-水體系的壁溫分布與甲苯-水體系近似。這是由于在四氯化碳、甲苯和苯這3種不同的分散相中,甲苯[to=0.151W/(m·K)]與苯[be=0.159W/(m·K)]的熱導(dǎo)率相近,而四氯化碳熱導(dǎo)率[ct=0.0454W/(m·K)]最小,四氯化碳-水體系兩相流的傳熱最慢,壁面上的熱量來不及傳導(dǎo)出去,因而壁面溫度最高。從圖3(c)中可以看出,相同進(jìn)口速率下,隨著管徑變大,生成的泰勒流液彈長度變大,個(gè)數(shù)減小,管壁的溫度變化和分布波動變大,但管徑對管壁整體溫度的影響不大。

    2.2?壁面剪切力與界面渦度對泰勒流傳熱的影響

    為了進(jìn)一步研究微通道內(nèi)液-液兩相泰勒流的傳熱規(guī)律,以下針對當(dāng)量直徑為0.5mm的微通道內(nèi),工作介質(zhì)為甲苯和水,進(jìn)口速度分別為0.3m/s和0.2m/s的泰勒流模型進(jìn)行討論。

    圖4給出了= 0.042s時(shí)刻相界面和微通道中心液體溫度隨管長的變化,圖中也給出了泰勒流液彈的位置和形狀。從圖4中可以看出,泰勒流液彈中部的液膜厚度并不是恒定的,在液彈的尾部附近出現(xiàn)了局部凸起,使得此處的液膜厚度達(dá)到最小值,由于液膜的波動,使得相界面的溫度變化也呈波動形式,但液膜波動與溫度變化并不是完全同步的,相界面上溫度的最低點(diǎn)出現(xiàn)在液彈尾部凸起處略后一點(diǎn)。液彈頭部相界面的溫度高于尾部的溫度,尾部相界面的溫度從局部最低升到局部最高所用距離很短,導(dǎo)致曲線很陡。

    由于甲苯和水的熱導(dǎo)率不同,微通道內(nèi)流體的溫度分布呈周期性波動變化。從圖4中心線處流體的溫度變化曲線可以看出,曲線的波谷處在泰勒流尾部兩相界面的右側(cè),波峰處在泰勒流頭部兩相界面的右側(cè)。對于一個(gè)甲苯液彈,頭部的溫度高于尾部,而且從頭部到尾部的溫度下降速度先增大后減小,液彈中部的溫度變化曲線最陡。對于甲苯液彈和它兩側(cè)的水相,可以看到液彈尾部的兩相界面處溫度變化最劇烈,而頭部的兩相界面處的溫度變化則較為平緩。

    圖5給出了=0.050s時(shí)刻微通道內(nèi)第5個(gè)泰勒流液彈處的局部溫度分布相圖,與GUPTA等[12]在氣-液兩相泰勒流的模擬結(jié)果較為相似。但GUPTA等僅在作為液相的水相中觀察到了渦流內(nèi)循環(huán),而從圖5中可以看出,在甲苯-水兩相泰勒流動中,無論是作為分散相的甲苯,還是作為連續(xù)相的水中,均存在著渦流內(nèi)循環(huán),分別位于壁面附近,液彈的頭部和尾部,因而液-液泰勒流動要比氣-液泰勒流動具有更強(qiáng)的傳熱能力。從圖5中還可以看到,在微通道中,泰勒流與壁面之間的液膜處,尤其是在泰勒流液彈尾部凸起處溫度最高,這是因?yàn)橐耗ず穸鹊臏p少,使得此處熱阻力降低,因而溫度較高[28]。

    圖6給出了= 0.042s時(shí)刻微通道壁面剪切力和壁面溫度隨管長的變化。從圖6中可以看出,泰勒流動下微通道管壁溫度呈周期性的波動上升趨勢,液彈處壁面溫度高于連續(xù)相處的壁面溫度。壁面溫度的波峰位于液彈尾部附近,液彈中部附近的壁面溫度則較為平緩,而文獻(xiàn)[12]中的氣-液泰勒流壁面溫度的波峰出現(xiàn)在氣彈的中部。從圖6中還可以看出,泰勒流動下的微通道壁面剪切力也呈周期性變化。液彈處的壁面剪切力高于連續(xù)相處的壁面剪切力。在每一個(gè)周期變化中有3個(gè)主要波峰,前兩個(gè)波峰位于液彈尾部,對應(yīng)著液膜最薄的位置,第3個(gè)波峰位于液彈中部偏右側(cè),對應(yīng)著液膜變厚的位置。甲苯-水泰勒流的壁面剪切力分布與文獻(xiàn)[12]中氣-液兩相泰勒流較為相似,但前者的壁面剪切力分布較為平均,為20Pa左右,而氣-液泰勒流的壁面剪切力的局部最低值達(dá)到了–80Pa。泰勒流動下微通道壁面剪切力與壁面溫度的周期性變化趨勢較為類似,壁面剪切力較大的位置壁面溫度也較高,說明壁面上的溫度變化與壁面剪切力的作用有關(guān),但這兩者的變化波動并不完全對應(yīng),壁面溫度的變化有一定的時(shí)間滯后性。

    圖7給出了= 0.042s時(shí)刻界面處渦度和溫度隨管長的變化。從圖7中可以看出,泰勒流兩相界面處的渦度和溫度呈周期性變化。渦度變化曲線在液彈頭部和尾部具有兩個(gè)比較尖銳的波峰,而且液彈尾部的渦度更高,達(dá)到24000s–1。除去液體吸熱逐漸升溫的趨勢來看,界面溫度的變化同樣呈周期性變化趨勢,液彈尾部的界面溫度在一個(gè)液彈周期內(nèi)最高。界面溫度的變化和渦度具有很強(qiáng)的關(guān)聯(lián)性,渦度的波峰往往出現(xiàn)在壁面溫度的波峰附近,但界面渦度與壁面溫度的變化不是完全同步的,渦度的最大值出現(xiàn)在溫度曲線波峰的偏右側(cè)一點(diǎn),說明渦度對溫度的影響有一個(gè)延遲現(xiàn)象。

    圖8給出了= 0.050s微通道內(nèi)第5個(gè)泰勒流液彈處的渦度分布相圖。從圖8中可以看出,液膜處的流體旋轉(zhuǎn)強(qiáng)度遠(yuǎn)大于液彈內(nèi)部,且液彈尾部凸起處的流體旋轉(zhuǎn)強(qiáng)度最為劇烈。微通道內(nèi)低速流體的流動形式通常為層流,通過兩相界面處的渦流內(nèi)循環(huán),熱量與分子更快速的傳遞到分散相的主體中,這也正是彈狀流的傳熱和傳質(zhì)效果優(yōu)于平行流的主要原因。

    圖9給出了= 0.042s時(shí)刻界面處局部努塞爾數(shù)和壁面溫度隨管長的變化。局部努塞爾數(shù)Nu定義如式(7)。

    式中,qw表示傳熱量;Tw表示壁面溫度;Tb表示靠近壁面處的液體主體溫度;kwa表示水的熱導(dǎo)率。從圖9中可以看出,Nux曲線的波峰對應(yīng)著液膜最薄的區(qū)域,壁面溫度波峰的偏右側(cè)。當(dāng)微通道內(nèi)的泰勒流液彈發(fā)展完全時(shí),液彈界面處的對流換熱強(qiáng)度是呈周期性穩(wěn)定變化的,液膜最薄處附近的對流換熱強(qiáng)度最高,其次是在液彈的頭部,而且液彈內(nèi)的對流換熱強(qiáng)度大于連續(xù)相部分。文獻(xiàn)[12]中氣-液泰勒流的壁面努塞爾數(shù)變化較為平緩,在200以下,而甲苯-水泰勒流的壁面努塞爾數(shù)在400左右,液彈尾部的局部壁面努塞爾數(shù)最大值超過了1000,這說明液-液兩相泰勒流動的壁面換熱強(qiáng)度要比氣-液兩相泰勒流大得多。

    由于局部努塞爾數(shù)不能很好地體現(xiàn)由于泰勒流液彈和連續(xù)相長度的不同造成的傳熱差異,引入泰勒流液彈平均努塞爾數(shù)這一概念,來表示泰勒流液彈的傳熱效果[12]。泰勒流液彈的平均努塞爾數(shù)定義如式(8)。

    式中,bav表示一個(gè)液彈主體的平均溫度。通過實(shí)驗(yàn)?zāi)M結(jié)果計(jì)算出=0.034s時(shí)刻甲苯-水兩相泰勒流中第4和第5個(gè)液彈的平均努賽爾數(shù)分別為14.3和14.42,是相同條件下單相流的1.3倍,與文獻(xiàn)[12]中2.5倍的結(jié)果相比偏小,可能是流體入口速度、兩相流量比等實(shí)驗(yàn)條件設(shè)置的差異所造成的。

    3 結(jié)論

    本文采用計(jì)算流體力學(xué)方法,針對微通道內(nèi)工作介質(zhì)為甲苯和水的兩相泰勒流的傳熱性質(zhì)進(jìn)行了研究,分析了不同因素對微通道壁面溫度分布的影響,以及壁面剪切力和界面渦度對微通道內(nèi)泰勒流傳熱的影響,發(fā)現(xiàn)壁面剪切力和界面渦度對微通道壁面和界面溫度的波動性變化有一定影響,并且有一定延遲現(xiàn)象。通過計(jì)算界面處局部努塞爾數(shù)隨管長的變化以及泰勒流液彈的平均努賽爾數(shù),對泰勒流的傳熱效果進(jìn)行了評價(jià)。泰勒流流型使得微通道內(nèi)的傳熱作用產(chǎn)生穩(wěn)定的周期性波動,兩相交界面處由于渦流內(nèi)循環(huán)的存在而成為溫度變化最為劇烈的位置,泰勒流動相較于單相流動和平行流更能滿足高通量下的換熱要求。

    符號說明

    C——顏色函數(shù) Ca——毛細(xì)管數(shù),μLUTP/σ cp——比熱容,J/(kg?K) D——通道直徑,m e——內(nèi)能,J/kg FSV——界面張力,N H——焓,J/kg K——熱導(dǎo)率,W/(m?K) Nuav——平均努塞爾數(shù) Nux——局部努塞爾數(shù) n——單位法向量 P——壓力,Pa qw——壁面熱通量,W/m2 R——徑向位置向量,m Re——雷諾數(shù),UTPρL D/μL T——溫度,K Tb——主體溫度,K Tw——壁面溫度,K t——時(shí)間,s UTP——混合物表觀速度,m/s v——速率,m/s x——軸向坐標(biāo),m α——體積分?jǐn)?shù) δ(r)——狄拉克函數(shù),m–1 κ——界面曲率,m–1 μ——動力黏度,Pa?s ρ——密度,kg/m3 σ——界面張力系數(shù),N/m 下角標(biāo) av——平均值 be——苯 ct——四氯化碳 int——界面 to——甲苯 w——壁面 wa——水

    [1] 蘇尚美,張亞男,成方園,等.微通道換熱器的特性分析及其應(yīng)用前景[J].區(qū)域供熱,2007(5):34-38.

    SU S M,ZHANG Y N,CHENG F Y,et al.The feature analysis of the micro channel heat exchanger and its application prospect[J].District Heating,2007(5):34-38.

    [2] TRIPLETT K A,GHIAASIAAN S M,ABDEL-KHALIK S I,et al.Gas–liquid two-phase flow in microchannels. Part I:two-phase flow patterns[J].International Journal of Multiphase Flow,1999,25(3):377-394.

    [3] KASHID M N,AGAR D W.Hydrodynamics of liquid-liquid slug flow capillary microreactor:flow regimes,slug size and pressure drop[J].Chemical Engineering Journal,2007,131(1/2/3):1-13.

    [4] 梁曉光,鄭揚(yáng),許松林.毛細(xì)管內(nèi)氣液兩相流動的CFD模擬[J].高校化學(xué)工程學(xué)報(bào),2010,24(3):370-375.

    LIANG X G,ZHENG Y,XU S L.CFD modeling for gas-liquid two phases Taylor flow in capillary[J].Journal of Chemical Engineering of Chinese Universities,2010,24(3):370-375.

    [5] QIAN D,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.

    [6] 張井志,李蔚.毛細(xì)管內(nèi)氣液Taylor流動的氣泡及阻力特性[J].化工學(xué)報(bào),2015,66(3):942-948.

    ZHANG J Z,LI W.Bubble and frictional characteristics of gas-liquid Taylor flow in capillary tube[J].CIESC Journal,2015,66(3):942-948.

    [7] BRETHERTON F P.The motion of long bubbles in tubes[J]. Journal of Fluid Mechanics,1961,10(2):166-188.

    [8] CHERLO S K R,KARIVETI S,PUSHPAVANAM S.Experimental and numerical investigations of two-phase(liquid-liquid) flow behavior in rectangular microchannels[J].Industrial & Engineering Chemistry Research,2009,49(2):893-899.

    [9] KREUTZER M T,KAPTEIJN F,MOULIJN J A,et al.Multiphase monolith reactors:chemical reaction engineering of segmented flow in microchannels[J].Chemical Engineering Science,2005,60(22):5895-5916.

    [10] LEUNG S S Y,GUPTA R,F(xiàn)LETCHER D F,et al.Effect of flow characteristics on Taylor flow heat transfer[J].Industrial & Engineering Chemistry Research,2011,51(4):2010-2020.

    [11] WALSH P A,WALSH E J,MUZYCHKA Y S.Heat transfer model for gas-liquid slug flows under constant flux[J].International Journal of Heat & Mass Transfer,2010,53(15/16):3193-3201.

    [12] GUPTA R,F(xiàn)LETCHER D F,HAYNES B S.CFD modelling of flow and heat transfer in the Taylor flow regime[J].Chemical Engineering Science,2010,65(6):2094-2107.

    [13] BANDARA T,NGUYEN N T,ROSENGARTEN G.Slug flow heat transfer without phase change in microchannels:a review[J].Chemical Engineering Science,2015,126:283-295.

    [14] JOVANOVI? J,ZHOU W,REBROV E V,et al.Liquid-liquid slug flow:hydrodynamics and pressure drop[J].Chemical Engineering Science,2011,66(1),42-54.

    [15] SALIM A,F(xiàn)OURAR M,PIRONON J,et al.Oil-water two-phase flow in microchannels:flow patterns and pressure drop measurements[J].The Canadian Journal of Chemical Engineering,2008,86(6):978-988.

    [16] FOROUGHI H,KAWAJI M.Viscous oil-water flows in a microchannel initially saturated with oil:flow patterns and pressure drop characteristics[J].International Journal of Multiphase Flow,2011,37(9):1147-1155.

    [17] ASTHANA A,ZINOVIK I,WEINMUELLER C,et al.Significant Nusselt number increase in microchannels with a segmented flow of two immiscible liquids:an experimental study[J].International Journal of Heat & Mass Transfer,2011,54(7):1456-1464.

    [18] URBANT P,LESHANSKY A,HALUPOVICH Y.On the forced convective heat transport in a droplet-laden flow in microchannels[J].Microfluidics & Nanofluidics,2008,4(6):533-542.

    [19] FISCHER M,JURIC D,POULIKAKOS D.Large convective heat transfer enhancement in microchannels with a train of coflowing immiscible or colloidal droplets[J].Journal of Heat Transfer,2010,132(11):184-191.

    [20] LIU D,WANG S.Hydrodynamics of Taylor flow in noncircular capillaries[J]. Chemical Engineering & Processing Process Intensification,2008,47(12):2098-2106.

    [21] ZHANG J Z,F(xiàn)LETCHER D F,LI W.Heat transfer and pressure drop characteristics of gas–liquid Taylor flow in mini ducts of square and rectangular cross-sections[J].International Journal of Heat & Mass Transfer,2016,103:45-56.

    [22] ANNALAND M V S,DIJKHUIZEN W,DEEN N G,et al.Numerical simulation of behavior of gas bubbles using a 3-D front-tracking method[J].AIChE Journal,2005,52(1): 99-110.

    [23] KASHID M N,RENKEN A,KIWI-MINSKER L.CFD modelling of liquid–liquid multiphase microstructured reactor:slug flow generation[J].Chemical Engineering Research & Design,2010,88(3A):362-368.

    [24] 許松林,趙嬋.氣液并流垂直液膜流動的數(shù)值模擬[J].天津大學(xué)學(xué)報(bào)(自然科學(xué)與工程技術(shù)版),2014(12):1039-1046.

    XU S L,ZHAO C.Numerical simulation of-current gas-liquid vertical liquid film flow[J].Journal of Tianjin University (Science and Technology),2014(12):1039-1046.

    [25] BRACKBILL J U,KOTHE D B,ZEMACH C.A continuum method for modeling surface tension[J]. Journal of Computational Physics,1992,100(2):335-354.

    [26] HAYNES W M.CRC Handbook of chemistry and physics[M]. 95th ed.Boca Raton:CRC Press,2014:1071-1115.

    [27] 顧惕人.表面化學(xué)[M].北京:科學(xué)出版社,1994:76.

    GU T R. Surface chemistry[M].Beijing:Science Press,1994:76.

    [28] PATTAMATTA A,SIELAFF A,Stephan P.A numerical study on the hydrodynamic and heat transfer characteristics of oscillating Taylor bubble in a capillary tube[J].Applied Thermal Engineering,2015,89:628-639.

    CFD simulations of heat transfer of liquid-liquid Taylor flow in microchannels

    LI Ting,XU Songlin

    (Key Laboratory of Systems Bioengineering(Ministry of Education),School of Chemical Engineering and Technology,Tianjin University,Tianjin 300072,China)

    Under low flow velocity,the flow pattern in a microchannel is typically the so-called Taylor flow regime.A computational fluid dynamics package FLUENT was adopted for simulation of the heat characteristics of liquid-liquid two-phase Taylor flow in a microchannel. The effects of the dispersed phase velocity,the materials and the diameter of the microchannel on the wall temperature were investigated first. The results indicated that when the velocity of the dispersed phase and the thermal conductivity of the material increased,the wall temperature decreased,and the change of diameter had little effect on the wall temperature. The effect of the wall shear stress,interfacial vorticity and the Nusselt number on the wall and internal temperature were investigated subsequently in view of toluene-water two-phase Taylor flow system in a microchannel with a diameter of 0.5mm. The heat transfer characteristics were compared with the gas-liquid two-phase flow in literature. The results showed that the wall shear stress and interface vorticity had certain effect on the fluctuating changes of wall and interface temperature. The wave crests of wall shear stress and interface vorticity often appeared near the peaks of temperature,and had a certain time lag. The wall Nusselt number of toluene-water two-phase Taylor flow was much greater than that of gas-liquid two-phase flow. The average Nusselt number of a liquid slug unit was 1.3 times of single-phase flow under the same condition.

    microchannels;Taylor flow;heat transfer;computational fluid dynamics(CFD)

    TQ021.1

    A

    1000–6613(2017)06–2078–08

    10.16085/j.issn.1000-6613.2017.06.016

    2016-10-08;

    2017-02-11。

    國家重點(diǎn)研究計(jì)劃項(xiàng)目(2016YFF0203802)。

    李婷(1991—),女,碩士研究生,主要從事計(jì)算流體力學(xué)研究。E-mail:litingtj@foxmail.com。聯(lián)系人:許松林,副研究員,主要從事化工和制藥分離技術(shù)、計(jì)算流體力學(xué)等研究。E-mail:slxu@tju.edu.cn。

    猜你喜歡
    界面
    聲波在海底界面反射系數(shù)仿真計(jì)算分析
    微重力下兩相控溫型儲液器內(nèi)氣液界面仿真分析
    國企黨委前置研究的“四個(gè)界面”
    基于FANUC PICTURE的虛擬軸坐標(biāo)顯示界面開發(fā)方法研究
    西門子Easy Screen對倒棱機(jī)床界面二次開發(fā)
    空間界面
    金秋(2017年4期)2017-06-07 08:22:16
    鐵電隧道結(jié)界面效應(yīng)與界面調(diào)控
    電子顯微打開材料界面世界之門
    人機(jī)交互界面發(fā)展趨勢研究
    手機(jī)界面中圖形符號的發(fā)展趨向
    新聞傳播(2015年11期)2015-07-18 11:15:04
    9热在线视频观看99| 久久国产精品人妻蜜桃| 日本av手机在线免费观看| 最近最新免费中文字幕在线| 久久人人97超碰香蕉20202| 黄片播放在线免费| 国产精品偷伦视频观看了| 久热爱精品视频在线9| 国产97色在线日韩免费| 欧美日韩亚洲综合一区二区三区_| 蜜桃国产av成人99| 国产片内射在线| 十八禁人妻一区二区| 国产一卡二卡三卡精品| 国产一区二区三区视频了| 黄色a级毛片大全视频| 久久久精品国产亚洲av高清涩受| 99re6热这里在线精品视频| 免费在线观看影片大全网站| 国产成人免费无遮挡视频| 国产麻豆69| 久久人妻熟女aⅴ| 巨乳人妻的诱惑在线观看| 午夜激情av网站| av一本久久久久| 99香蕉大伊视频| 深夜精品福利| 日韩制服丝袜自拍偷拍| 热99国产精品久久久久久7| 丰满迷人的少妇在线观看| 国产欧美日韩综合在线一区二区| 在线天堂中文资源库| 日韩制服丝袜自拍偷拍| 一边摸一边做爽爽视频免费| 黑人操中国人逼视频| 熟女少妇亚洲综合色aaa.| 别揉我奶头~嗯~啊~动态视频| 国产av精品麻豆| 国产精品.久久久| 涩涩av久久男人的天堂| 亚洲中文字幕日韩| 咕卡用的链子| 亚洲国产中文字幕在线视频| 成人黄色视频免费在线看| 黑丝袜美女国产一区| 1024香蕉在线观看| 色在线成人网| 日韩免费av在线播放| 国产精品久久久久久精品电影小说| 国产一区二区 视频在线| 99热网站在线观看| 桃红色精品国产亚洲av| 国产午夜精品久久久久久| 国产欧美日韩一区二区三区在线| 老司机在亚洲福利影院| 淫妇啪啪啪对白视频| 午夜久久久在线观看| 日本一区二区免费在线视频| 90打野战视频偷拍视频| 十八禁网站免费在线| 一区福利在线观看| 国产精品影院久久| 色老头精品视频在线观看| 午夜免费鲁丝| 91麻豆av在线| 一区二区av电影网| 人成视频在线观看免费观看| 99精国产麻豆久久婷婷| 亚洲精品美女久久久久99蜜臀| 亚洲精品粉嫩美女一区| 天天添夜夜摸| 亚洲av第一区精品v没综合| 色精品久久人妻99蜜桃| 69精品国产乱码久久久| bbb黄色大片| 日韩大片免费观看网站| 黑人猛操日本美女一级片| 人人澡人人妻人| 91精品国产国语对白视频| 亚洲专区字幕在线| 777久久人妻少妇嫩草av网站| 国产视频一区二区在线看| 午夜精品国产一区二区电影| 国产在线一区二区三区精| 啦啦啦免费观看视频1| 精品乱码久久久久久99久播| 精品国产一区二区三区久久久樱花| 国产精品.久久久| 肉色欧美久久久久久久蜜桃| 欧美在线一区亚洲| 1024视频免费在线观看| 日韩大码丰满熟妇| 精品国产一区二区久久| 国产亚洲欧美精品永久| 国产激情久久老熟女| 久久影院123| 亚洲av片天天在线观看| 亚洲自偷自拍图片 自拍| www日本在线高清视频| 国产精品亚洲一级av第二区| 男人舔女人的私密视频| 亚洲国产毛片av蜜桃av| 国产av一区二区精品久久| 一本一本久久a久久精品综合妖精| 精品人妻熟女毛片av久久网站| 午夜福利在线免费观看网站| 狠狠婷婷综合久久久久久88av| 久久久久国内视频| 极品人妻少妇av视频| 国产精品免费大片| 国产aⅴ精品一区二区三区波| 国精品久久久久久国模美| 精品少妇内射三级| 色婷婷久久久亚洲欧美| 水蜜桃什么品种好| 亚洲精品国产区一区二| 在线观看人妻少妇| 免费少妇av软件| 欧美变态另类bdsm刘玥| 国产三级黄色录像| 亚洲国产精品一区二区三区在线| 电影成人av| tocl精华| 成人国产av品久久久| 777米奇影视久久| 国产精品一区二区在线不卡| 国产高清国产精品国产三级| 夫妻午夜视频| 黄色片一级片一级黄色片| 亚洲人成电影免费在线| 亚洲国产欧美在线一区| av天堂久久9| 中文字幕高清在线视频| 在线十欧美十亚洲十日本专区| 久久久精品区二区三区| 乱人伦中国视频| 久久狼人影院| 美女高潮喷水抽搐中文字幕| 日韩视频一区二区在线观看| 1024视频免费在线观看| 亚洲伊人久久精品综合| 欧美国产精品一级二级三级| 日本av手机在线免费观看| 国产亚洲精品一区二区www | 超碰成人久久| 女人被躁到高潮嗷嗷叫费观| 国产又色又爽无遮挡免费看| 老汉色av国产亚洲站长工具| 高清黄色对白视频在线免费看| 制服诱惑二区| 亚洲国产欧美在线一区| 国产1区2区3区精品| 变态另类成人亚洲欧美熟女 | 日韩欧美一区二区三区在线观看 | 搡老熟女国产l中国老女人| 超碰成人久久| 中文字幕高清在线视频| 欧美日韩av久久| 国产高清视频在线播放一区| 99re在线观看精品视频| 国产欧美日韩精品亚洲av| 亚洲av欧美aⅴ国产| 日本撒尿小便嘘嘘汇集6| 免费av中文字幕在线| 男男h啪啪无遮挡| 亚洲成人免费电影在线观看| 国产一区有黄有色的免费视频| 欧美日本中文国产一区发布| 亚洲黑人精品在线| 亚洲成a人片在线一区二区| 91老司机精品| 欧美日韩av久久| 国产成人系列免费观看| 亚洲熟妇熟女久久| 一边摸一边抽搐一进一小说 | 午夜久久久在线观看| 老司机深夜福利视频在线观看| 老熟妇仑乱视频hdxx| 18禁国产床啪视频网站| 国产成人精品久久二区二区91| 两人在一起打扑克的视频| 午夜免费鲁丝| 亚洲熟女毛片儿| 国产高清激情床上av| 狠狠婷婷综合久久久久久88av| 国产成人一区二区三区免费视频网站| 久久热在线av| 成年版毛片免费区| 十八禁高潮呻吟视频| 91精品三级在线观看| 18在线观看网站| 成年动漫av网址| 国产精品麻豆人妻色哟哟久久| 欧美日韩中文字幕国产精品一区二区三区 | 成人影院久久| 午夜日韩欧美国产| 丝袜在线中文字幕| 国产av一区二区精品久久| 一夜夜www| 操出白浆在线播放| 国产精品一区二区免费欧美| av免费在线观看网站| 啦啦啦免费观看视频1| 男女之事视频高清在线观看| 亚洲国产av新网站| 精品国产乱子伦一区二区三区| 在线观看免费视频日本深夜| 久久狼人影院| 亚洲七黄色美女视频| 性少妇av在线| 色尼玛亚洲综合影院| 国产成人免费无遮挡视频| 亚洲精品在线美女| 国产精品香港三级国产av潘金莲| 99国产综合亚洲精品| 一区二区av电影网| 黄色毛片三级朝国网站| 在线 av 中文字幕| 国产在线观看jvid| 成人亚洲精品一区在线观看| 国产欧美日韩一区二区精品| 色尼玛亚洲综合影院| 老熟妇乱子伦视频在线观看| 国产精品麻豆人妻色哟哟久久| 天天添夜夜摸| 性少妇av在线| 国产aⅴ精品一区二区三区波| 少妇 在线观看| 国产男靠女视频免费网站| 国产男女超爽视频在线观看| 最黄视频免费看| 精品少妇久久久久久888优播| 亚洲精品在线美女| 久久av网站| 国产成人影院久久av| 正在播放国产对白刺激| 女同久久另类99精品国产91| 精品一区二区三卡| videosex国产| bbb黄色大片| 啦啦啦在线免费观看视频4| 亚洲精品国产区一区二| 精品人妻熟女毛片av久久网站| 制服诱惑二区| 精品一区二区三区av网在线观看 | kizo精华| 国产成人精品久久二区二区91| 精品久久蜜臀av无| 热99国产精品久久久久久7| 欧美在线黄色| 国产福利在线免费观看视频| 老司机午夜福利在线观看视频 | 久久人人97超碰香蕉20202| 岛国毛片在线播放| 男人操女人黄网站| 国产在线精品亚洲第一网站| 欧美日韩黄片免| 精品少妇黑人巨大在线播放| 欧美乱妇无乱码| 天堂8中文在线网| 动漫黄色视频在线观看| 日韩欧美三级三区| 嫁个100分男人电影在线观看| 黑人猛操日本美女一级片| 亚洲五月色婷婷综合| 亚洲成av片中文字幕在线观看| 超色免费av| 精品一区二区三区av网在线观看 | 成年人免费黄色播放视频| av福利片在线| 国产成人精品久久二区二区免费| 另类精品久久| 欧美日本中文国产一区发布| 在线观看免费午夜福利视频| 真人做人爱边吃奶动态| 老熟女久久久| 男女之事视频高清在线观看| 亚洲精品国产一区二区精华液| 亚洲精品一卡2卡三卡4卡5卡| a级片在线免费高清观看视频| 高清毛片免费观看视频网站 | 日本五十路高清| 少妇裸体淫交视频免费看高清 | 久久人妻av系列| 精品一区二区三区视频在线观看免费 | 久久久水蜜桃国产精品网| bbb黄色大片| 老司机深夜福利视频在线观看| 日本欧美视频一区| 亚洲一区二区三区欧美精品| 久久国产精品人妻蜜桃| 欧美 日韩 精品 国产| 久久精品人人爽人人爽视色| bbb黄色大片| 无限看片的www在线观看| 啦啦啦免费观看视频1| 最近最新免费中文字幕在线| 国产高清国产精品国产三级| 超碰97精品在线观看| 国产欧美日韩精品亚洲av| 性少妇av在线| 老司机亚洲免费影院| 久久久精品免费免费高清| 黄片大片在线免费观看| 一级毛片精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲熟妇熟女久久| 国产麻豆69| 免费观看人在逋| 法律面前人人平等表现在哪些方面| 黄色成人免费大全| 99国产精品免费福利视频| 日韩有码中文字幕| 精品卡一卡二卡四卡免费| 欧美精品高潮呻吟av久久| 精品福利永久在线观看| 国产人伦9x9x在线观看| 欧美 日韩 精品 国产| 美女高潮到喷水免费观看| 亚洲国产av影院在线观看| 99国产精品99久久久久| 好男人电影高清在线观看| 中文字幕另类日韩欧美亚洲嫩草| 成人三级做爰电影| 女人爽到高潮嗷嗷叫在线视频| 精品人妻1区二区| 男女下面插进去视频免费观看| 搡老乐熟女国产| 午夜视频精品福利| 国产精品.久久久| 午夜精品久久久久久毛片777| 亚洲精品久久午夜乱码| 亚洲全国av大片| 日韩视频一区二区在线观看| 在线av久久热| 国产精品1区2区在线观看. | 老熟妇乱子伦视频在线观看| 日韩欧美三级三区| 国产熟女午夜一区二区三区| 国产1区2区3区精品| 国产在线一区二区三区精| 丝袜美腿诱惑在线| 亚洲精品国产一区二区精华液| 亚洲自偷自拍图片 自拍| 一区二区日韩欧美中文字幕| 又大又爽又粗| 人人妻人人添人人爽欧美一区卜| 亚洲全国av大片| 人人妻人人澡人人爽人人夜夜| 一本色道久久久久久精品综合| 麻豆成人av在线观看| 十八禁网站网址无遮挡| 久久这里只有精品19| 久久久久网色| 亚洲av美国av| 精品久久久久久电影网| 亚洲va日本ⅴa欧美va伊人久久| 国产主播在线观看一区二区| 日本av免费视频播放| 中文字幕高清在线视频| 日本五十路高清| 亚洲国产精品一区二区三区在线| 我的亚洲天堂| 成年版毛片免费区| h视频一区二区三区| 国产黄频视频在线观看| 国产1区2区3区精品| 男女无遮挡免费网站观看| a在线观看视频网站| 免费一级毛片在线播放高清视频 | av欧美777| 日本一区二区免费在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产a三级三级三级| 美女高潮到喷水免费观看| 十分钟在线观看高清视频www| 又紧又爽又黄一区二区| 亚洲精品国产精品久久久不卡| 色综合欧美亚洲国产小说| 超碰97精品在线观看| 国产有黄有色有爽视频| 亚洲欧洲日产国产| 另类精品久久| 免费高清在线观看日韩| av网站在线播放免费| 999精品在线视频| 人人妻人人澡人人爽人人夜夜| 人妻久久中文字幕网| 美女视频免费永久观看网站| 99在线人妻在线中文字幕 | 菩萨蛮人人尽说江南好唐韦庄| 一本大道久久a久久精品| 精品一品国产午夜福利视频| 一区二区三区精品91| 99国产精品99久久久久| 新久久久久国产一级毛片| 啪啪无遮挡十八禁网站| 50天的宝宝边吃奶边哭怎么回事| 涩涩av久久男人的天堂| 满18在线观看网站| 我的亚洲天堂| 欧美精品高潮呻吟av久久| 精品久久久久久久毛片微露脸| 亚洲五月婷婷丁香| 国产精品av久久久久免费| 久久久欧美国产精品| 侵犯人妻中文字幕一二三四区| 国产1区2区3区精品| 人人澡人人妻人| 美女国产高潮福利片在线看| 少妇被粗大的猛进出69影院| 在线观看免费日韩欧美大片| 91九色精品人成在线观看| 欧美av亚洲av综合av国产av| 亚洲av国产av综合av卡| 精品国内亚洲2022精品成人 | 国产成+人综合+亚洲专区| 热99久久久久精品小说推荐| 国产熟女午夜一区二区三区| 欧美乱码精品一区二区三区| 国产av一区二区精品久久| 国产欧美日韩一区二区精品| 中亚洲国语对白在线视频| 久久热在线av| 精品亚洲成国产av| 丝袜人妻中文字幕| 啦啦啦在线免费观看视频4| 久久精品国产亚洲av香蕉五月 | 一本—道久久a久久精品蜜桃钙片| 女同久久另类99精品国产91| 久久中文字幕一级| 亚洲精品久久午夜乱码| 久久精品国产亚洲av香蕉五月 | 无人区码免费观看不卡 | 色精品久久人妻99蜜桃| 丝袜美足系列| 他把我摸到了高潮在线观看 | 亚洲伊人色综图| 久久久久国内视频| 国产在视频线精品| 亚洲精品美女久久久久99蜜臀| 欧美另类亚洲清纯唯美| 三上悠亚av全集在线观看| 男女下面插进去视频免费观看| 9191精品国产免费久久| 久久久国产成人免费| 亚洲国产中文字幕在线视频| 精品人妻1区二区| 在线观看免费日韩欧美大片| 又黄又粗又硬又大视频| 十八禁网站免费在线| 一级毛片电影观看| 免费看a级黄色片| 久9热在线精品视频| 黄色成人免费大全| 久久人妻福利社区极品人妻图片| 久久 成人 亚洲| 91麻豆精品激情在线观看国产 | 午夜福利视频在线观看免费| 波多野结衣一区麻豆| 欧美日韩av久久| 欧美日韩中文字幕国产精品一区二区三区 | 日韩三级视频一区二区三区| 精品午夜福利视频在线观看一区 | 91成年电影在线观看| 国产精品 欧美亚洲| 黑人巨大精品欧美一区二区蜜桃| 超色免费av| 国产精品免费视频内射| 黑丝袜美女国产一区| 亚洲欧美日韩高清在线视频 | 亚洲国产欧美在线一区| 成人国产一区最新在线观看| 如日韩欧美国产精品一区二区三区| 9热在线视频观看99| 极品教师在线免费播放| 一区二区三区国产精品乱码| 亚洲av美国av| 国产精品免费大片| 亚洲美女黄片视频| 最黄视频免费看| 一区二区三区激情视频| 欧美精品人与动牲交sv欧美| 午夜91福利影院| 黄色视频不卡| 黑人欧美特级aaaaaa片| 久久狼人影院| 女人精品久久久久毛片| 嫩草影视91久久| 波多野结衣一区麻豆| 热re99久久精品国产66热6| 亚洲欧美一区二区三区黑人| 在线观看免费高清a一片| 中文字幕人妻丝袜一区二区| 一级毛片电影观看| 丝袜喷水一区| 如日韩欧美国产精品一区二区三区| 欧美日韩精品网址| 91av网站免费观看| 国产男女内射视频| 国产无遮挡羞羞视频在线观看| 男女免费视频国产| 制服人妻中文乱码| 在线观看人妻少妇| 日本wwww免费看| 日韩制服丝袜自拍偷拍| 欧美成狂野欧美在线观看| 亚洲午夜精品一区,二区,三区| 国产精品一区二区在线观看99| 久久99热这里只频精品6学生| 久久久欧美国产精品| 成人三级做爰电影| 国产精品久久久av美女十八| 精品国产一区二区久久| 纯流量卡能插随身wifi吗| 香蕉国产在线看| 建设人人有责人人尽责人人享有的| 又大又爽又粗| 老司机在亚洲福利影院| 18禁美女被吸乳视频| 咕卡用的链子| 女性生殖器流出的白浆| 丁香六月欧美| 国产在线免费精品| 黄色怎么调成土黄色| 黄色视频不卡| 淫妇啪啪啪对白视频| 深夜精品福利| 国产精品 欧美亚洲| h视频一区二区三区| √禁漫天堂资源中文www| av欧美777| 欧美成人免费av一区二区三区 | 一级a爱视频在线免费观看| 高清黄色对白视频在线免费看| 国产精品熟女久久久久浪| 18禁美女被吸乳视频| 免费日韩欧美在线观看| 露出奶头的视频| 国产av国产精品国产| 国产精品.久久久| 久久 成人 亚洲| 亚洲精品美女久久av网站| 99精品在免费线老司机午夜| 男女床上黄色一级片免费看| 国产aⅴ精品一区二区三区波| a级片在线免费高清观看视频| 亚洲欧洲精品一区二区精品久久久| 少妇精品久久久久久久| 99热国产这里只有精品6| 久久久久久免费高清国产稀缺| 丝袜美腿诱惑在线| 真人做人爱边吃奶动态| 他把我摸到了高潮在线观看 | 中文字幕另类日韩欧美亚洲嫩草| 欧美大码av| 国产日韩欧美在线精品| 女警被强在线播放| 黄色视频在线播放观看不卡| 叶爱在线成人免费视频播放| 大片免费播放器 马上看| 国产欧美日韩精品亚洲av| 两个人看的免费小视频| 国产精品久久久人人做人人爽| 人妻久久中文字幕网| 午夜福利视频在线观看免费| 久久这里只有精品19| 91字幕亚洲| 国产精品成人在线| 9热在线视频观看99| 男女免费视频国产| 日韩免费高清中文字幕av| 亚洲av电影在线进入| 午夜日韩欧美国产| 五月开心婷婷网| 在线观看免费高清a一片| 人人妻人人爽人人添夜夜欢视频| 欧美国产精品一级二级三级| 少妇被粗大的猛进出69影院| 一级片'在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 午夜老司机福利片| 国产高清视频在线播放一区| 在线观看免费午夜福利视频| 一级a爱视频在线免费观看| 99国产精品一区二区蜜桃av | 欧美一级毛片孕妇| 久久国产精品男人的天堂亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 一边摸一边做爽爽视频免费| 国产成人影院久久av| 亚洲午夜精品一区,二区,三区| 中文字幕人妻丝袜制服| 大片免费播放器 马上看| 欧美激情极品国产一区二区三区| 欧美大码av| 1024香蕉在线观看| 在线播放国产精品三级| 日韩欧美一区二区三区在线观看 | 黄片小视频在线播放| 久久人妻av系列| 国产一卡二卡三卡精品| 正在播放国产对白刺激| kizo精华| 欧美久久黑人一区二区| 99精品欧美一区二区三区四区| 淫妇啪啪啪对白视频| 国产精品.久久久| 高清欧美精品videossex| 中国美女看黄片| 一边摸一边做爽爽视频免费| 波多野结衣一区麻豆| 亚洲伊人色综图| 欧美中文综合在线视频| 国产精品免费一区二区三区在线 |