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

    不同截面疏水性微肋陣內(nèi)減阻特性

    2016-07-04 03:43:28姜桂林管寧張承武劉志剛山東省科學(xué)院能源研究所山東濟(jì)南250014
    化工學(xué)報(bào) 2016年4期
    關(guān)鍵詞:壓力降

    姜桂林,管寧,張承武,劉志剛(山東省科學(xué)院能源研究所,山東 濟(jì)南 250014)

    ?

    不同截面疏水性微肋陣內(nèi)減阻特性

    姜桂林,管寧,張承武,劉志剛
    (山東省科學(xué)院能源研究所,山東 濟(jì)南 250014)

    摘要:采用疏水液處理紫銅微肋陣表面,分別在截面為圓形、菱形和橢圓形微肋陣實(shí)驗(yàn)段內(nèi)形成接觸角分別為99.5°、119.5°和151.5°(水為工質(zhì))的疏水性表面,實(shí)驗(yàn)測(cè)試各個(gè)工況下流道內(nèi)流動(dòng)阻力和壓力降,分析了不同截面形狀對(duì)疏水性微肋陣內(nèi)減阻特性的影響規(guī)律。結(jié)果表明,當(dāng)接觸角增大時(shí),壓力降變化率在微肋陣內(nèi)的變化規(guī)律隨截面形狀的改變而發(fā)生變化;同一接觸角下,橢圓形微肋陣內(nèi)壓力降變化率隨流量增加而逐漸減小,而菱形和圓形微肋陣則先減小后保持常數(shù)。相同Reynolds數(shù)(Re)下,3種截面實(shí)驗(yàn)段中減阻率均隨接觸角的增大而增加。接觸角相同時(shí),橢圓形微肋陣內(nèi)阻力系數(shù)變化因子隨Re的增大而逐漸減??;菱形和圓形則先減小后保持常數(shù),接觸角為151.5°時(shí)最小阻力系數(shù)變化因子分別為50.81%和58.68%。

    關(guān)鍵詞:疏水性微肋陣;接觸角;摩擦阻力系數(shù);壓力降;減阻率

    2015-06-10收到初稿,2015-08-04收到修改稿。

    聯(lián)系人:劉志剛。第一作者:姜桂林(1981—),男,碩士,助理研究員。

    Received date: 2015-06-10.

    Foundation item: supported by the National Natural Science Foundation of China (51176105),the Natural Science Foundation of Shandong Province (ZR2012EEQ015) and the Shandong Province Science and Technology Development Plans (2014GGX104008).

    引 言

    近年來(lái),隨著電子科技的發(fā)展,越來(lái)越多的高集成度電子產(chǎn)品滲入到各行各業(yè)的應(yīng)用領(lǐng)域,在給社會(huì)發(fā)展及日常生活帶來(lái)日新月異變革的同時(shí),電子芯片集成度的提高和功率呈幾何級(jí)數(shù)的增長(zhǎng),微肋片結(jié)構(gòu)為主的換熱器因其高效的散熱效率而成為傳熱領(lǐng)域熱點(diǎn)之一[1-4],然而通道內(nèi)顯著的表面效應(yīng)和極高的面體比使流動(dòng)阻力過(guò)高,限制了微肋陣換熱結(jié)構(gòu)的應(yīng)用。降低內(nèi)壁表面能可減小微/納米通道內(nèi)表面附近流體分子的平均密度,從而形成一個(gè)低密度層,使流體流過(guò)通道表面時(shí)的流動(dòng)滑移長(zhǎng)度明顯增加[5-8],因此可有效降低微納米通道內(nèi)的流動(dòng)阻力。

    迄今,國(guó)內(nèi)外學(xué)者已經(jīng)對(duì)微/小肋群的流動(dòng)傳熱和常規(guī)尺寸通道疏水性表面的流動(dòng)傳熱做了一定的研究。國(guó)際上Peles等[9]及Kosar等[10-11]最早開(kāi)始微肋陣的相關(guān)研究,其研究分別從肋片形狀、尺寸以及排列方式對(duì)于流動(dòng)特性和對(duì)流換熱特性的影響進(jìn)行了考察,并首先指出了端壁面效應(yīng)這一影響微肋陣流動(dòng)和傳熱的重要因素;而Mustafa等[12]以及Chang等[13]則針對(duì)端壁面效應(yīng)對(duì)于微肋陣內(nèi)熱流體力學(xué)特性的具體影響進(jìn)行了數(shù)值模擬研究;劉志剛等[14]實(shí)驗(yàn)研究了叉排微柱群內(nèi)頂部縫隙對(duì)傳熱效率影響,研究顯示頂部縫隙隨Re的增加對(duì)流動(dòng)和換熱的影響越來(lái)越明顯。Selvarasu等[15]對(duì)微通道針肋密度進(jìn)行研究,表明層流狀態(tài)下較低的針肋密度具有更好的壓降特征和換熱性能,而針肋密度的增大會(huì)大幅度增加壓力損失并阻礙換熱能力的提高。張承武等[16]實(shí)驗(yàn)研究了不同截面形狀微柱群內(nèi)流動(dòng)阻力特征,發(fā)現(xiàn)流量相同時(shí)圓形截面比橢圓和菱形截面流動(dòng)壓降大;且隨Re增加,橢圓和菱形柱群內(nèi)的f小于圓形微柱群。董麗寧等[17]數(shù)值研究了微針肋陣列通道針肋幾何結(jié)構(gòu)對(duì)層流流動(dòng)的摩擦系數(shù)的影響,擬合了孔隙率、彎曲度、寬高比和Re的摩擦系數(shù)關(guān)聯(lián)式。

    眾多針對(duì)微肋陣的研究表明[18],盡管微肋陣具有卓越的換熱性能,但較高的流動(dòng)阻力始終是其發(fā)展的瓶頸,而表面超疏水處理在微通道減阻中的成功應(yīng)用無(wú)疑為微肋陣性能的改進(jìn)提供了一個(gè)新的思路。研究表明,對(duì)微通道進(jìn)行疏水性處理能使流動(dòng)阻力降低25%~40%[19-21]。Li等[22]數(shù)值分析了疏水表面固體區(qū)域所占平面面積降低到0.4時(shí),微通道內(nèi)流動(dòng)阻力可降低40%?;羲乇蟮萚23]采用Qian等[24]的化學(xué)蝕刻與表面修飾的方法在內(nèi)徑為0.6 mm鋁制微通道內(nèi)制作了超疏水表面,并將流動(dòng)壓力降與超親水微通道進(jìn)行了對(duì)比,發(fā)現(xiàn)流動(dòng)壓力降最高可達(dá)25%。姜桂林等[25]實(shí)驗(yàn)研究了不同接觸角疏水表面對(duì)圓形微肋陣內(nèi)部的流動(dòng)特征,發(fā)現(xiàn)相同Re下接觸角越大,流動(dòng)減阻效果越好,且隨著Re增加,每個(gè)接觸角狀態(tài)下流動(dòng)減阻效果均下降,而且在超疏水處理后流動(dòng)減阻率可達(dá)58.4%。

    目前針對(duì)微肋陣內(nèi)部超疏水表面減阻的研究相對(duì)較少,尤其不同截面微肋陣內(nèi)部超疏水表面的流動(dòng)減阻研究還鮮有公開(kāi)報(bào)道?;谶@一研究現(xiàn)狀,將不同接觸角疏水處理液固載于不同截面紫銅微針肋內(nèi)部,分別在截面為圓形、橢圓形和菱形微肋陣內(nèi)部逐一形成不同接觸角疏水表面。以去離子水為流動(dòng)工質(zhì),實(shí)驗(yàn)研究疏水處理前后在不同Re下微針肋截面形狀以及不同接觸角對(duì)f的影響,并在此基礎(chǔ)上對(duì)疏水性減阻機(jī)制進(jìn)行系統(tǒng)的理論分析。

    1 實(shí)驗(yàn)裝置及測(cè)試原理

    1.1實(shí)驗(yàn)段的制作

    微肋陣使用以紫銅為基底,采用機(jī)械加工制作,實(shí)驗(yàn)段的結(jié)構(gòu)及幾何參數(shù)示意圖如圖1所示,實(shí)物照片如圖2所示。

    以改性有機(jī)硅為基質(zhì),加入2%全氟辛基氟硅烷和微納米粒子后合成超疏水處理液。采用上述處理液對(duì)微針肋截面為圓形、橢圓形和菱形的紫銅微肋陣內(nèi)部進(jìn)行噴涂處理后,將實(shí)驗(yàn)段置于80℃真空環(huán)境下固化形成超疏水性表面涂層。為保證各微柱表面涂層的均勻,利用高壓噴筆霧化疏水處理液,采用噴筆傾斜、實(shí)驗(yàn)段勻速旋轉(zhuǎn)的方法,將疏水處理液噴涂于實(shí)驗(yàn)段表面,所制作實(shí)驗(yàn)段接觸角經(jīng)多點(diǎn)測(cè)試。取多點(diǎn)平均值分,紫銅光面與涂覆不同接觸角疏水涂層的實(shí)驗(yàn)段接觸角采用接觸角測(cè)試儀(上海中晨JC-2000),水滴體積為8 ml,測(cè)試如圖3所示。

    圖1 實(shí)驗(yàn)段示意圖Fig.1 Sketch picture of test section

    涂層固化后表觀接觸角為分別為99.5°、119.5° 和151.5°,各點(diǎn)測(cè)試值之間誤差在±1°以?xún)?nèi)。

    由于涂層固化后具有一定厚度,因此噴涂涂層后微肋陣外徑尺寸有所變化。為保證噴涂后實(shí)驗(yàn)段中各微柱直徑的均一,進(jìn)行了多次實(shí)驗(yàn)段制備操作,并將獲得的實(shí)驗(yàn)段進(jìn)行SEM電鏡掃描,選取各微柱直徑差在±1 μm以?xún)?nèi)的微肋陣作為最終實(shí)驗(yàn)段。未經(jīng)疏水性處理的實(shí)驗(yàn)段原始尺寸列于表1,疏水處理后的尺寸列于表2;實(shí)驗(yàn)段外形尺寸的電鏡掃描測(cè)量照片如圖4所示。

    1.2實(shí)驗(yàn)裝置及誤差分析

    實(shí)驗(yàn)裝置如圖5所示。

    表1 實(shí)驗(yàn)段幾何尺寸Table 1 Geometrical size of test section

    圖2 不同形狀實(shí)驗(yàn)段照片F(xiàn)ig.2 Photos of test sections with different cross section shapes

    圖3 紫銅光面與不同涂層接觸角測(cè)試Fig.3 Contact angle measurements on ordinary copper surface and hydrophobic surfaces

    實(shí)驗(yàn)臺(tái)采用12 MPa的高壓氮?dú)庾鳛閴毫υ?。氮?dú)馄垦b有氧氣減壓閥,可粗調(diào)實(shí)驗(yàn)所需壓力;后接精密減壓閥,實(shí)現(xiàn)實(shí)驗(yàn)所需壓力的準(zhǔn)確調(diào)節(jié)(精確到100 Pa)。儲(chǔ)氣罐能夠使氮?dú)饬鲃?dòng)更加平穩(wěn);儲(chǔ)液罐(容積為20 L)內(nèi)裝有實(shí)驗(yàn)用去離子水,儲(chǔ)液罐后接裝有濾膜為20 mm的過(guò)濾器,過(guò)濾器與實(shí)驗(yàn)段之間安裝微流量計(jì)以測(cè)量通過(guò)實(shí)驗(yàn)段工質(zhì)的流量,實(shí)驗(yàn)段兩端裝有壓力傳感器與T型熱電偶,以測(cè)量進(jìn)出口壓力pin、pout和進(jìn)出口溫度Tin、Tout。所有流量、溫度和壓力參數(shù)直接被數(shù)據(jù)采集儀(Agilent34970A)采集;各實(shí)驗(yàn)組件之間用耐高壓的聚氨酯軟管連接,接頭采用硬密封。

    圖4 不同接觸角各形狀實(shí)驗(yàn)段尺寸掃描Fig.4 Shape experimental section size scan of different contact angle

    表2 疏水處理后實(shí)驗(yàn)段幾何尺寸Table 2 Geometrical size of test section after hydrophobic treated

    圖5 實(shí)驗(yàn)裝置簡(jiǎn)圖Fig.5 Schematic diagram of experimental system

    微肋陣通道通過(guò)快速凝固膠與兩個(gè)傳感器的連接頭連接,出口工質(zhì)直接排入廢液罐。該試驗(yàn)臺(tái)流量是利用微流量計(jì)來(lái)測(cè)量的。當(dāng)微管進(jìn)出口端的溫度與壓力不隨液體流動(dòng)而變化時(shí),采集微流量計(jì)流量信號(hào),通過(guò)數(shù)模轉(zhuǎn)換得到流量值。實(shí)驗(yàn)中采用微流量計(jì)測(cè)量流量,在實(shí)驗(yàn)過(guò)程中,先將壓力調(diào)整到所需壓力值后,待壓力值與溫度值穩(wěn)定后,開(kāi)始實(shí)驗(yàn)。

    實(shí)驗(yàn)裝置中為流量計(jì)為EH8301A(讀數(shù)精度0.01%),壓力傳感器精度為0.1%,T型熱電偶100℃內(nèi)測(cè)量精度為±0.15℃。微肋陣幾何尺寸和微柱外徑用電子顯微鏡測(cè)量,其分辨率高于0.1 mm,測(cè)量誤差為±5 mm。其他誤差按照文獻(xiàn)[26]所提供的方法計(jì)算,結(jié)果列于表3中。

    表3 實(shí)驗(yàn)誤差Table 3 Experimental uncertainties

    2 實(shí)驗(yàn)結(jié)果和分析

    管內(nèi)流動(dòng)摩擦阻力系數(shù)可由式(1)進(jìn)行計(jì)算

    其中

    根據(jù)選取的特征尺寸的不同,Reynolds數(shù)Re有如下兩種計(jì)算公式

    票陣通道和不同接觸角涂層處理后通道內(nèi)Re及相應(yīng)f。

    2.1不同接觸角微肋陣內(nèi)的流動(dòng)特性

    為了考察接觸角對(duì)于微肋陣內(nèi)壓力降的影響規(guī)律,圖6(a)~(c)分別給出了當(dāng)接觸角θ=83°~151.5°時(shí)橢圓、菱形和圓形微肋陣內(nèi)壓力降隨流量的變化曲線。在圖6(a)中,當(dāng)接觸角θ=83°,即未經(jīng)疏水性處理時(shí)去離子水流過(guò)橢圓形微肋陣內(nèi)的壓力降隨接觸角的增加而逐漸降低,其中接觸角θ= 99.5°和119.5°時(shí)微肋陣內(nèi)的壓力降差別較小,而當(dāng)接觸角增加至θ=151.5°時(shí)橢圓形微肋陣內(nèi)壓力降與較小接觸角工況下相比有明顯下降。這是由于本研究中微肋陣表面的疏水性微納米涂層降低了液態(tài)水流過(guò)微肋陣時(shí)的流動(dòng)阻力。

    圖6 不同接觸角微肋陣內(nèi)壓力降隨流量的變化Fig.6 Profiles of pressure drops vs. volume flow rate in micro pin fins with different contact angles

    隨著微肋陣表面接觸角的逐漸增大,微肋陣表面納米粒子涂層在電鏡掃描下呈現(xiàn)出如圖7(a)所示的微結(jié)構(gòu),眾多納米尺度的微觀凸起使得微肋陣表面形成了類(lèi)似于Ou等[19]模型的結(jié)構(gòu),如圖7(b)所示。

    圖7 疏水涂層表面圖及疏水原理圖Fig.7 Schematic diagram of model for super-hydrophobic drag reduction

    由于納米凸起形成的界面上的表面張力作用,水不易進(jìn)入疏水涂層表面相鄰微納米結(jié)構(gòu)內(nèi)部,從而增加了水和空氣的接觸面積,在微肋陣外表面形成滑移效應(yīng),因此對(duì)于具有相同特征尺寸的微肋陣而言,該處理大大降低了流體流過(guò)微肋陣時(shí)的壓力降。盡管采用的納米涂層使得微肋陣的特征尺度隨著接觸角的增加而有所增大,但對(duì)于橢圓形微肋陣而言,特征尺度的增加對(duì)于壓力降的影響要小于疏水性微納米涂層,出現(xiàn)了如圖6(a)所示的變化。而對(duì)于菱形微肋陣,如圖6(b)所示,當(dāng)接觸角從θ=83°增大至θ=119.5°時(shí)微肋陣內(nèi)壓力降明顯降低,然而當(dāng)接觸角繼續(xù)增大至θ=151.5°時(shí),在流量GV<150 ml·min-1時(shí),θ=151.5°的菱形微肋陣內(nèi)的壓力降與θ=99.5°和119.5°時(shí)相比較小。然而隨著流量的增加,θ=151.5°和θ=119.5°的菱形微肋陣內(nèi)的壓力降數(shù)值相差不大,均低于θ=99.5°的微肋陣;對(duì)于圓形微肋陣,如圖6(c)所示,當(dāng)接觸角從θ=83°增加至θ=119.5°時(shí),圓形微肋陣內(nèi)的壓力降呈現(xiàn)下降趨勢(shì),然而當(dāng)表面接觸角繼續(xù)增加至θ=151.5°時(shí),在較低流量下圓形微肋陣內(nèi)壓力降與θ=119.5°時(shí)較為接近,而當(dāng)GV>80 ml·min?1時(shí),不同于橢圓形和菱形微肋陣,接觸角為θ=151.5°的圓形微肋陣內(nèi)的壓力降與θ=119.5°時(shí)相比有所增加。橢圓形、菱形及圓形微肋陣內(nèi)的壓力降隨著接觸角由小增大的過(guò)程出現(xiàn)了上述不同現(xiàn)象,其主要原因有兩方面:一是疏水性微納米涂層對(duì)不同形狀微肋陣的減阻機(jī)理有所不同;二是接觸角變化時(shí)微肋陣特征尺寸發(fā)生了改變。流體繞流微肋片的過(guò)程中的阻力損失由摩擦阻力損失和壓差阻力損失(也稱(chēng)為形狀阻力或漩渦阻力)兩部分構(gòu)成,在低流速下和具有較好流線型肋片的繞流過(guò)程中,繞流邊界層分離較晚,微肋片背風(fēng)區(qū)內(nèi)尚未形成旋渦或者旋渦很小,因此摩擦阻力損失在微肋陣壓力降中占主導(dǎo)地位。然而隨著流體流速的增加,微肋片尾流區(qū)面積逐漸增大,壓差阻力在壓力降中所占比重逐漸增加。尤其是對(duì)于流線型較差的微肋陣,例如圓形微肋陣,流阻中壓差阻力所占比重要明顯高于橢圓和菱形微肋片。表面疏水性微納米涂層能夠降低邊界層內(nèi)的剪切應(yīng)力從而使壓差阻力降低。然而如前所述,不同接觸角的納米涂層的厚度有所不同,隨著接觸角的增加涂層厚度明顯增大,圓形微肋片的接觸角從θ=83°增加至θ=99.5°時(shí),其特征尺寸增加了7.5%,而當(dāng)表面接觸角增加至θ=151.5°時(shí)其特征尺度與θ=83°相比增加了37.75%,與θ=119.5°時(shí)相比增加了12.68%。這不僅使得摩擦阻力損失增大,對(duì)于圓形微肋片而言相同流量下較高的繞流流速也導(dǎo)致壓差阻力迅速增加。盡管表面微納米疏水層仍具有減阻效果,但上述因素的綜合效果使得圓形微肋陣當(dāng)接觸角由θ=119.5°增加至θ=151.5°時(shí)壓力降略有升高,壓差阻力比重略小于圓形微肋陣的菱形微肋陣當(dāng)θ=119.5°和θ=151.5°時(shí)壓力降變化不大,而橢圓形微肋陣內(nèi)的壓力降則由于其良好的流線外形而使得摩擦阻力損失始終占主導(dǎo)地位,其壓力降在研究范圍內(nèi)隨接觸角的增大而始終降低。

    圖8 不同接觸角微肋陣阻力系數(shù)隨Re的變化曲線Fig.8 Profiles of f vs Re in micro pin fins with different contact angles

    圖8給出了橢圓、菱形和圓形微肋陣內(nèi)當(dāng)表面接觸角分別為83°、99.5°、119.5°以及151.5°時(shí)流動(dòng)阻力系數(shù)隨Re的變化曲線。

    由圖可以看出,對(duì)于3種截面形狀的微肋陣,其流動(dòng)阻力系數(shù)均隨著接觸角的增大而減小,且在較低Re降低幅度更加明顯。如圖7所示,隨著微肋陣表面接觸角的增加,壁面上微納米凸起的尺寸明顯減小而分布密度增加,因此當(dāng)液態(tài)水流過(guò)疏水性表面時(shí),壁面-水之間的接觸面積明顯減少,從而增加空氣-水接觸面上的流動(dòng)滑移,降低了邊界層內(nèi)的速度梯度,進(jìn)而使得摩擦阻力明顯降低。這一減阻效果在低Re下更加明顯,一方面流動(dòng)壓力的增加,另一方面是由于邊界層的分離造成的。設(shè)相鄰兩個(gè)納米突起的平均距離為W,則空氣水兩相界面所承受的壓力差為

    當(dāng)壓力差Dp相對(duì)較小時(shí),超疏水表面的表面張力使得水很難進(jìn)入到超疏水表面相鄰微納米尺度的凹槽內(nèi)部,較大比例的空氣-水接觸面積使得邊界層內(nèi)速度梯度明顯減小,從而降低了流動(dòng)阻力;然而當(dāng)壓力差增大到臨界狀態(tài)即f 等于超疏水表面表觀接觸角q 時(shí),此時(shí)

    如果壓力差繼續(xù)增大,水將壓迫或擠出突起之間的空氣,去離子水滴與基體接觸面將顯著增加,從而降低了壁面上的速度滑移,削弱了減阻效果。

    另外,對(duì)比圖8(a)~(c)還可以看出,盡管橢圓、菱形和圓形微肋陣內(nèi)不同接觸時(shí)阻力系數(shù)之間的差距均隨Re的增加而減小,然而在較高Re下圓形微肋陣各接觸角下阻力系數(shù)之差要明顯高于橢圓形微肋陣,菱形則介于兩者之間。這一現(xiàn)象可能是由于邊界層的分離所造成的。在橢圓形微肋陣內(nèi),由于其較好的流線外形,相同接觸角及Re下邊界層分離較晚,摩擦阻力始終對(duì)阻力系數(shù)具有主導(dǎo)影響,因此表面疏水性帶來(lái)的減阻效果受流速增加的影響最為明顯;而對(duì)于圓形微肋片,與橢圓形及菱形微肋片相比其邊界層更早發(fā)生分離而形成尾流區(qū),眾所周知尾流區(qū)壁面附近壓力較低,這在一定程度上降低了高Re下流動(dòng)壓力對(duì)與疏水性壁面減阻效果的影響,從而出現(xiàn)了如圖8所示的現(xiàn)象。

    2.2不同截面疏水性微肋陣內(nèi)的減阻特性

    為了定量分析接觸角對(duì)于微肋陣壓力降的影響,本文定義壓力降變化率和阻力系數(shù)變化率如下所示

    圖9(a)~(c)給出了橢圓、菱形和圓形截面微肋陣內(nèi)當(dāng)接觸角分別為θ=99.5°、119.5°以及151.5°時(shí)壓力降變化率隨流量的變化曲線。如前所述,對(duì)微肋陣表面進(jìn)行疏水性處理后形成的納米結(jié)構(gòu)增加了表面的流動(dòng)滑移,降低了摩擦阻力,因此對(duì)于具有較好流線外形的橢圓形微肋陣,其壓力降變化率隨著表面接觸角的增加而增大,特別是接觸角為θ=151.5°的微肋陣,其dp_coe的值明顯高于其他接觸角,θ=119.5°的橢圓形微肋陣內(nèi)的dp_coe值略高于θ=99.5°時(shí),如圖9(a)所示。同時(shí),3種不同表面接觸角的橢圓形微肋陣內(nèi)的dp_coe值均隨著流量的增加而逐漸降低,如圖9(a)所示,當(dāng)流量20 ml·min-1增至220 ml·min?1時(shí),接觸角為θ=151.5°的橢圓形微肋陣內(nèi)的dp_coe值降低了41.51%,這是由于Re增大過(guò)程中慣性力和表面張力的相對(duì)大小發(fā)生了變化。疏水涂層形成雙重結(jié)構(gòu)的表面張力恒定,當(dāng)微肋陣內(nèi)來(lái)自工質(zhì)內(nèi)部壓力不大時(shí),工質(zhì)與基底凹槽部分接觸面積相對(duì)較小。隨著實(shí)驗(yàn)壓力不斷增加,工質(zhì)內(nèi)壓力克服表面張力不斷滲入至疏水表面凹槽,液固接觸面積逐漸增加,導(dǎo)致減阻效果逐漸降低。如前所述,受壓差阻力及特征尺度變化的影響,菱形微肋陣中壓力降變化率盡管也隨接觸角的增加而增大,但θ=119.5°時(shí)壓力降變化率與θ=151.5°非常接近;而對(duì)于圓形微肋陣,接觸角θ=151.5°時(shí)的壓力降變化率要略低于θ=119.5°,如圖9(b)、(c)所示。另外,對(duì)比圖9(a)~(c)還可以發(fā)現(xiàn),橢圓形微肋陣中壓力降變化率隨流量增加始終下降;菱形和圓形微肋陣中壓力降變化率先減小后保持常數(shù),特別是在大接觸角下更加明顯。對(duì)于菱形和圓形微肋陣,當(dāng)流量高于200 ml·min?1時(shí)其最高壓力降變化率分別為22.92%和28.35%,分別出現(xiàn)于接觸角為151.5°和119.5°時(shí);出現(xiàn)這一現(xiàn)象的原因可能是由于在較大流量下菱形和圓形微肋陣內(nèi)的壓力降中壓差阻力占主導(dǎo)地位,而不同于摩擦阻力,壓差阻力受上述因素的影響較小,同時(shí)肋片在尾流區(qū)內(nèi)壁面附近壓力大致接近于分離點(diǎn)的壓力,壓力梯度和速度梯度較小,因此也降低了流速增加時(shí)慣性力變化對(duì)于疏水性微肋陣的摩擦阻力的減阻效果的影響,從而出現(xiàn)了如圖9(b)、(c)所示的現(xiàn)象。

    圖9 不同接觸角微肋陣內(nèi)壓力降變化率與流量的關(guān)系Fig.9 Relationships between coefficient of pressure drop reduction with volume flow rate in micro pin fins with contact angles

    圖10給出了不同接觸角下橢圓、菱形和圓形微肋陣內(nèi)減阻率隨Re的變化關(guān)系。

    由圖可以看出,由于疏水性表面上微納米凸起內(nèi)空氣-水兩相接觸面的存在,使得3種微肋陣內(nèi)阻力系數(shù)均隨接觸角的增加而增大。且對(duì)于橢圓形和菱形微肋陣,當(dāng)接觸角分別為99.5°、119.5°以及151.5°時(shí)減阻率的增加幅度較為均勻,而圓形微肋陣內(nèi)當(dāng)接觸角為119.5°以及151.5°時(shí)的減阻率較為接近,均遠(yuǎn)高于接觸角為99.5°時(shí),如圖10(c)所示。分別對(duì)比圖10(a)與圖9(a),圖10(b)與圖9(b),以及圖10(c)與圖9(c)可以發(fā)現(xiàn),對(duì)于橢圓形和菱形微肋陣,接觸角為151.5°時(shí)的微肋陣在節(jié)省泵功和降低阻力系數(shù)上均具有最佳減阻效果,而對(duì)于圓形微肋陣,綜合考慮泵功和阻力系數(shù)的最佳接觸角則在119.5°~151.5°之間,如前所示,這一現(xiàn)象是由于接觸角增大對(duì)微肋陣特征尺寸和阻力構(gòu)成的綜合影響造成的。

    圖10 不同接觸角微肋陣內(nèi)減阻率隨Re的變化曲線Fig.10 Comparisons of df_coe with Re in hydrophobic micro pin fins with different cross section shapes

    另外,對(duì)比圖10(a)~(c)還可以發(fā)現(xiàn),在Re=200~1200的范圍內(nèi),各接觸角下橢圓形微肋陣的減阻率始終隨Re的增加而逐漸減小,而接觸角為119.5°和151.5°的菱形微肋陣內(nèi)當(dāng)Re>400后減阻率變得很小,最小值為151.5°時(shí)的50.81%;圓形微肋陣則當(dāng)Re>400后減阻率基本為常數(shù),最小值出現(xiàn)在Re=600時(shí),減阻率為58.68%。作者認(rèn)為這可能是由于尾流區(qū)的流動(dòng)轉(zhuǎn)捩造成的。與橢圓形微肋片相比,菱形和圓形微肋片繞流尾流區(qū)更早發(fā)生轉(zhuǎn)捩,特別是圓形微肋陣,在Re<400時(shí)尾流區(qū)內(nèi)流動(dòng)已經(jīng)進(jìn)入過(guò)渡區(qū),使得尾流區(qū)對(duì)阻力系數(shù)的影響明顯增加。由于尾流區(qū)平均壓力較低,這不但大大降低了較高Re下流體壓力的增加對(duì)表面疏水性帶來(lái)的減阻效果的影響,而且有利于降低尾流區(qū)的漩渦阻力,因此出現(xiàn)了如圖10(b)、(c)所示的現(xiàn)象。

    3 結(jié) 論

    通過(guò)在微肋陣表面固載微納米粒子涂層的方法獲得具有不同接觸角的微肋陣,并建立實(shí)驗(yàn)系統(tǒng)測(cè)試了當(dāng)接觸角分別為83°、99.5°、119.5°以及151.5°時(shí),橢圓形、菱形以及圓形微肋陣內(nèi)的流動(dòng)及減阻特性,定量分析了不同接觸角以及截面形狀對(duì)于微肋陣內(nèi)減阻特性的影響規(guī)律,得到如下結(jié)論。

    (1)當(dāng)接觸角由83°增加至151.5°時(shí),橢圓和菱形微肋陣內(nèi)壓力降變化率均隨接觸角的增加而增大;而對(duì)于圓形微肋陣,當(dāng)接觸角由83°增加至119.5°時(shí)壓力降變化率隨接觸角增加而增大,然而當(dāng)接觸角從119.5°增加至151.5°時(shí)壓力降變化率有所減小。

    (2)橢圓形微肋陣內(nèi)壓力降變化率隨流量的增加而逐漸減小,菱形和圓形微肋陣則先減小后保持常數(shù);對(duì)于菱形和圓形微肋陣,當(dāng)流量高于200 ml·min?1時(shí)其最高壓力降變化率分別為22.92%和28.35%,分別出現(xiàn)于接觸角為151.5°和119.5°時(shí)。

    (3)相同Re下,各實(shí)驗(yàn)段內(nèi)減阻率均隨接觸角的增加而逐漸增大;同一接觸角下橢圓形微肋陣內(nèi)減阻率隨Re的增大而逐漸減小,菱形和圓形則先減小后保持常數(shù),當(dāng)接觸角為151.5°時(shí)其減阻率分別為50.81%和58.68%。

    符號(hào)說(shuō)明

    A——微肋片橫截面積,m2

    d——微柱直徑,m

    df_coe——減阻率

    dp_coe——壓力降變化率

    f——摩擦因子

    G——體積流量,ml·min?1

    H ——微凸起平均高度,m

    L ——通道長(zhǎng)度,m

    n ——流速最大時(shí)微柱個(gè)數(shù)

    Dp ——實(shí)驗(yàn)段兩端壓差,Pa

    Re ——Reynolds數(shù)

    SD——前后排微肋片斜向中心距,m

    SL——前后排微肋片縱向中心距,m

    ST——橫向微肋片中心距,m

    T ——工質(zhì)溫度,℃

    umax——最小通道處流速,m·s?1

    W ——微凸起平均間距,m

    w ——通道寬度,m

    θ ——疏水表面表觀接觸角,(°)

    m ——?jiǎng)恿︷ざ龋琋·s·m-2

    ρ ——工質(zhì)密度,kg·m?3

    σ ——表面張力,mN·m?1

    f ——?dú)庖罕砻娼佑|角,(°)

    下角標(biāo)

    air ——空氣

    hy ——經(jīng)疏水處理

    in ——進(jìn)口

    max ——最大值

    no-hy ——未經(jīng)疏水處理

    out ——出口

    V ——體積

    water ——水

    References

    [1] KOSAR A,MISHRA C,PELES Y. Laminar flow across a bank of low aspect ratio micro pin fins [J]. Journal of Fluids Engineering,2005,127(3): 419-430.

    [2] MOORES K A,KIM J,JOSHI Y K. Heat transfer and fluid flow in shrouded pin fin arrays with and without tip clearance [J]. Int. J.Heat Mass Transfer,2009,52: 5978-5989.

    [3] NDAO S,PELES Y,JENSEN M K. Effects of pin shape and configuration on the single-phase heat transfer characteristics of jet impingement on micro pin fins [J]. Int. J. Heat Mass Transfer,2014,70: 856-863.

    [4] SHENOY S,TULLIUS J F,BAYAZITOGLU Y. Minichannels with carbon nanotube structured surfaces for cooling applications [J]. Int. J. Heat Mass Transfer,2011,54(25/26): 5379-5385.

    [5] BAUDRY J,CHARIAIX E,TONCK A,et al. Experimental evidence for a large slip effect at a nonwetting fluid-solid interface [J]. Langmuir,2001,17(17): 5232-5236.

    [6] CHOI C H,JOHAN K,WESTIN A,et al. Apparent slip flows in hydrophilic and hydrophobic microchannels [J]. Physics of Fluids,2003,15(10): 2897-2902.

    [7] WANG C Y,NG C. Natural convection in a vertical slit microchannel with superhydrophobic slip and temperature jump [J]. J. Heat Transfer,2014,136: 034502.1-034502.6.

    [8] MOHAMMAD N N,SETAREH S,ALIREZA M. Drag reduction in a turbulent channel flow with hydrophobic wall [J]. Journal of Hydrodynamics,2012,24(3): 458-466.

    [9] PELES Y,KOSAR A,MISHRA C,et al. Forced convective heat transfer across a pin fin micro heat sink [J]. Int. J. Heat Mass Transf.,2005,48: 3615-3627.

    [10] KOSAR A,MISHRA C,PELES Y. Laminar flow across a bank of low aspect ratio micro pin fins [J]. J. Fluids Eng.,2005,127(3): 419-430.

    [11] KOSAR A,PELES Y. Thermal-hydraulic performance of memes-based pin fin heat sink [J]. J. Heat Transf.,2006,128: 121-131.

    [12] MUSTAFA K,MEHMED R O,KOSAR A. Parametric study on the effect of end walls on heat transfer and fluid flow across a micro pin-fin [J]. Int. J. Therm. Sci.,2011,50(6): 1073-1084.

    [13] CHANG S W,HU Y-W. Endwall thermal of performances of radially rotating rectangular channel with pin-fins on skewed rib lands [J]. Int. J. Heat Mass Transf.,2014,69: 173-190.

    [14] 劉志剛,張承武,管寧. 叉排微柱群內(nèi)頂部縫隙對(duì)傳熱效率的影響[J]. 化工學(xué)報(bào),2012,63(4):1025-1031. DOI:10.3969/j.issn.0438-1157. 2012.04.005. LIU Z G,ZHANG C W,GUAN N. Influence of tip clearance on heat transfer efficiency in staggered micro-cylinders-group [J].CIESC Journal,2012,63(4):1025-1031. DOI:10.3969/j.issn.0438-1157. 2012.04.005.

    [15] SELVARASU N K C,TAFTI D K,BLACKWELL N E. Effect of pin density on heat-mass transfer and fluid flow at low Reynolds numbers in minichannels [J]. J. Heat Transfer,2010,132: 061702.1-061702.8.

    [16] 張承武,蒲龍梅,姜桂林,等. 不同截面形狀微肋片內(nèi)流動(dòng)阻力特征[J]. 化工學(xué)報(bào),2014,65(6): 2042-2048. DOI: 10.3969/j.issn.0438-1157. 2014.06.013. ZHANG C W,PU L M,JIANG G L,et al. Resistance characteristics of micro pin fins with different cross-section shapes [J]. CIESC Journal,2014,65(6): 2042-2048. DOI: 10.3969/j.issn.0438-1157. 2014.06.013.

    [17] 董麗寧,全曉軍,鄭平. 微針肋陣列通道中水的層流摩擦系數(shù)的模擬[J]. 上海交通大學(xué)學(xué)報(bào),2010,44(11): 1561-1565. DONG L N,QUAN X J,ZHENG P. Numerical simulation for friction factor of laminar flow of water in microchannels with micro-pin-fin arrays [J]. Journal of Shanghai Jiao Tong University,2010,44(11): 1561-1565.

    [18] MOHAMMAD M M,KOK-CHEONG W,MANSOOR S. Numerical investigation of fluid flow and heat transfer under high heat flux using rectangular micro-channels [J]. Int. J. Comm. Heat Mass Transfer,2012,39: 291-297.

    [19] OU J,BLAIR P,ROTHSTEIN J P. Laminar drag reduction in microchannels using ultrahydrophobic surfaces [J]. Physics of Fluids,2004,16(12): 4635-4643.

    [20] SONG D,DANIELLO R J,ROTHSTEIN J P. Drag reduction using superhydrophobic sanded Teflon surfaces [J]. Exp. Fluids,2014,55: 1783.

    [21] 宋善鵬,于志家,劉興華,等. 超疏水表面微通道內(nèi)水的傳熱特征[J]. 化工學(xué)報(bào),2008,59(10): 2465-2469. SONG S P,YU Z J,LIU X H,et al. Heat transfer characteristics of water flowing in microchannels with super-hydrophobic inner surface [J]. Journal of Chemical Industry and Engineering (China),2008,59(10): 2465-2469.

    [22] LI B W,YAO Z H,HAO P F. Incompressible LBGK simulation of flow characteristics in a micro-channel with patterned superhydrophobic surfaces [J]. Applied Mathematical Modelling,2015,39: 300-308.

    [23] 霍素斌,于志家,李艷峰,等. 超疏水表面微通道內(nèi)水的流動(dòng)特征[J]. 化工學(xué)報(bào),2007,58(11): 2721-2726. HUO S B,YU Z J,LI Y F,et al. Flow characteristics of water in microchannel with super-hydrophobic surface [J]. Journal of Chemical Industry and Engineering (China),2007,58(11): 2721-2726.

    [24] QIAN B T,SHEN Z Q. Fabrication of superhydrophobic surface by dislocation-selective chemical on aluminum,copper and zinc substrate [J]. Langmuir,2005,21(3): 9007-9009.

    [25] 姜桂林,張承武,管寧,等. 水在不同接觸角微柱群內(nèi)的流動(dòng)特征[J]. 化工學(xué)報(bào),2015,66(5): 1704-1709. DOI: 10.11949/j.issn. 0438-1157.20141799 JIANG G L,ZHANG C W,GUAN N,et al. Flow characteristics of water on hydrophobic micro cylinders group with different contact angle [J]. CIESC Journal,2015,66(5): 1704-1709. DOI: 10.11949/j.issn. 0438-1157.20141799

    [26] MOFFAT R J. Describing the uncertainties in experimental results [J]. Exp. Thermal Fluid Sci.,1988,1(1): 3-17.

    [27] ZUKAUSKAS A A. Heat transfer from tubes in cross flow [J]. Advances in Heat Transfer,1972,8: 93-160.

    Characteristics of drag reduction in different cross sections of hydrophobic micro pin fins

    JIANG Guilin,GUAN Ning,ZHANG Chengwu,LIU Zhigang
    (Energy Research Institute of Shandong Academy of Sciences,Jinan 250014,Shandong,China)

    Abstract:The hydrophobic coatings with different properties are coated on the surfaces of micro pin fins with different cross sections of circular,diamond and elliptical to obtain hydrophobic micro pin fins with contact angles of 99.5°,119.5° and 151.5°,respectively. The reduction coefficient of the pressure drop and the friction resistance are measured experimentally when water flows through the hydrophobic micro pin fins at different flow rates. Experimental results show that the change rate of pressure drop in ellipse and diamond micro pin fins increases gradually,however the change rate of pressure drop in circular micro pin fins increases at first and then decreases when the contact angle changes from 83° (the red copper smooth surface) to 151.5°. With the increase of the flow rate,the change rate of pressure drop in ellipse micro pin fins gradually decreases,while it decreases at first and then almost keeps constant in diamond and circular micro pin fins. The investigation also illustrates that the resistance reduction rates in all of test sections become large with the increase of the contact angle at the same Re. Under the same contact angle conditions,the resistance reduction rate gradually decreases in elliptical micro pin fins with the increase of Re,while it decreases at first and then keeps almost constant in diamond shaped and circular test sections. The minimal resistance reduction rates in diamond shaped and circular micro pin fins are50.81% and 58.68% at contact angle of 151.5°,respectively.

    Key words:hydrophobic micro pin fins; contact angel; friction resistance coefficient; pressure drop; resistance reduction rate

    DOI:10.11949/j.issn.0438-1157.20150890

    中圖分類(lèi)號(hào):TK 124

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):0438—1157(2016)04—1258—11

    基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(51176105);山東省自然科學(xué)基金項(xiàng)目(ZR2012EEQ015);山東省科技發(fā)展計(jì)劃項(xiàng)目(2014GGX104008)。

    Corresponding author:Prof .LIU Zhigang,zgliu9322@163.com

    猜你喜歡
    壓力降
    陽(yáng)極壓力降在PEMFC故障診斷中的應(yīng)用
    非牛頓流體腈綸原液管道壓力降計(jì)算
    壓縮空氣泡沫流動(dòng)壓力降研究
    天津科技(2022年5期)2022-05-31 02:18:06
    上下分體式全年制冷風(fēng)冷冷水機(jī)組系統(tǒng)設(shè)計(jì)
    環(huán)境溫度對(duì)車(chē)用空調(diào)濾清器的影響研究
    管柱式氣-液分離器溢流壓力降計(jì)算模型
    熱交換器壓降改進(jìn)設(shè)計(jì)
    機(jī)載火焰抑制器流通性能研究
    SK型靜態(tài)混合器的液液分散數(shù)值模擬
    乳膠基質(zhì)管路的壓力降研究與安全性分析?
    爆破器材(2016年2期)2016-04-26 01:45:19
    欧美一级a爱片免费观看看| 国产精品 国内视频| 性欧美人与动物交配| 国产97色在线日韩免费| 女人被狂操c到高潮| 精品午夜福利视频在线观看一区| xxxwww97欧美| 在线观看美女被高潮喷水网站 | 黑人欧美特级aaaaaa片| 亚洲精品影视一区二区三区av| 深爱激情五月婷婷| 极品教师在线免费播放| 18禁在线播放成人免费| 久久久久久国产a免费观看| 精品久久久久久久毛片微露脸| 国产精品免费一区二区三区在线| 日韩欧美三级三区| 18禁国产床啪视频网站| 一进一出好大好爽视频| 国产91精品成人一区二区三区| 五月玫瑰六月丁香| 亚洲av成人av| 亚洲国产欧美人成| 免费大片18禁| 天堂√8在线中文| 十八禁网站免费在线| 99在线人妻在线中文字幕| 亚洲av一区综合| 亚洲精华国产精华精| 一区二区三区免费毛片| 久久精品91无色码中文字幕| 欧美在线一区亚洲| 黄色视频,在线免费观看| 免费观看人在逋| 免费看十八禁软件| 啪啪无遮挡十八禁网站| 男女视频在线观看网站免费| 成年人黄色毛片网站| 国产欧美日韩一区二区精品| 热99re8久久精品国产| 无人区码免费观看不卡| 两人在一起打扑克的视频| 热99re8久久精品国产| 亚洲av不卡在线观看| 精品日产1卡2卡| 欧美中文综合在线视频| 国产精品 欧美亚洲| 日本在线视频免费播放| 99久国产av精品| 2021天堂中文幕一二区在线观| 国产aⅴ精品一区二区三区波| 国产精品亚洲美女久久久| 一个人看的www免费观看视频| 国产老妇女一区| 日韩欧美国产一区二区入口| av国产免费在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产精品免费一区二区三区在线| 亚洲内射少妇av| 女生性感内裤真人,穿戴方法视频| 天堂动漫精品| 久久久国产精品麻豆| 少妇裸体淫交视频免费看高清| 免费大片18禁| 成年人黄色毛片网站| 热99re8久久精品国产| 精品久久久久久成人av| 国产久久久一区二区三区| a级一级毛片免费在线观看| 成人国产综合亚洲| 亚洲国产高清在线一区二区三| 九九在线视频观看精品| 欧美高清成人免费视频www| 少妇熟女aⅴ在线视频| 亚洲av五月六月丁香网| 黄色日韩在线| 全区人妻精品视频| 午夜福利高清视频| 最好的美女福利视频网| 国产精品香港三级国产av潘金莲| 亚洲激情在线av| 综合色av麻豆| 天天躁日日操中文字幕| 成年人黄色毛片网站| 少妇的丰满在线观看| 国产在视频线在精品| 一区二区三区国产精品乱码| 在线观看舔阴道视频| 欧美日韩瑟瑟在线播放| 全区人妻精品视频| 亚洲自拍偷在线| 免费av毛片视频| 亚洲av日韩精品久久久久久密| 女人十人毛片免费观看3o分钟| av视频在线观看入口| 在线观看免费午夜福利视频| 精品午夜福利视频在线观看一区| 日韩欧美国产在线观看| 亚洲一区二区三区不卡视频| 国产精品一区二区三区四区免费观看 | 搡女人真爽免费视频火全软件 | 国产精品亚洲一级av第二区| 久久久久亚洲av毛片大全| 午夜福利欧美成人| 99在线人妻在线中文字幕| 叶爱在线成人免费视频播放| 美女cb高潮喷水在线观看| 成人午夜高清在线视频| a在线观看视频网站| 欧美色欧美亚洲另类二区| 欧美大码av| 又黄又粗又硬又大视频| 午夜福利在线观看免费完整高清在 | 亚洲av免费在线观看| 天美传媒精品一区二区| 久久久久九九精品影院| 18禁黄网站禁片午夜丰满| 亚洲第一电影网av| 成熟少妇高潮喷水视频| 国产av麻豆久久久久久久| 又黄又爽又免费观看的视频| 日本成人三级电影网站| 日本与韩国留学比较| 亚洲av免费在线观看| 啪啪无遮挡十八禁网站| 在线观看免费午夜福利视频| 欧美极品一区二区三区四区| 淫妇啪啪啪对白视频| 欧美成人一区二区免费高清观看| 日韩欧美一区视频在线观看 | 国产av码专区亚洲av| 国产成人精品一,二区| 日产精品乱码卡一卡2卡三| 国产午夜精品论理片| 午夜免费激情av| 精品人妻熟女av久视频| 男人和女人高潮做爰伦理| 国产国拍精品亚洲av在线观看| 国产精品一及| 人妻少妇偷人精品九色| 熟女电影av网| 亚洲精品一区蜜桃| 青青草视频在线视频观看| 国产 一区精品| 日韩欧美国产在线观看| 丝袜美腿在线中文| 国产黄a三级三级三级人| 真实男女啪啪啪动态图| 一级毛片久久久久久久久女| 麻豆久久精品国产亚洲av| 好男人在线观看高清免费视频| 国产一区二区三区av在线| 亚洲人成网站高清观看| 免费人成在线观看视频色| 少妇人妻一区二区三区视频| 亚洲精品成人久久久久久| 成人毛片60女人毛片免费| 日韩av在线大香蕉| 九草在线视频观看| 成年人午夜在线观看视频 | 国产精品蜜桃在线观看| 国内精品美女久久久久久| 高清在线视频一区二区三区| 国产淫片久久久久久久久| 中文字幕亚洲精品专区| 免费不卡的大黄色大毛片视频在线观看 | 一级av片app| 插阴视频在线观看视频| av一本久久久久| 91久久精品国产一区二区三区| 国产亚洲精品av在线| 日日啪夜夜爽| 91久久精品电影网| 国产不卡一卡二| 亚洲怡红院男人天堂| 欧美精品国产亚洲| 岛国毛片在线播放| 国产国拍精品亚洲av在线观看| 国产精品麻豆人妻色哟哟久久 | 日韩av在线大香蕉| 网址你懂的国产日韩在线| 一个人观看的视频www高清免费观看| 国内少妇人妻偷人精品xxx网站| 亚洲三级黄色毛片| freevideosex欧美| 麻豆久久精品国产亚洲av| 久久久久久久久久成人| 成年人午夜在线观看视频 | 69av精品久久久久久| 亚洲欧美日韩东京热| 国产美女午夜福利| 啦啦啦中文免费视频观看日本| 91久久精品国产一区二区三区| 亚洲成人中文字幕在线播放| 老司机影院成人| 精品久久久久久久末码| 精品国产露脸久久av麻豆 | 91精品伊人久久大香线蕉| 中国美白少妇内射xxxbb| 国模一区二区三区四区视频| 97在线视频观看| 欧美人与善性xxx| 中文字幕免费在线视频6| 亚洲精品久久久久久婷婷小说| 少妇的逼水好多| 日韩欧美精品免费久久| 搡女人真爽免费视频火全软件| 岛国毛片在线播放| 可以在线观看毛片的网站| 国产在视频线在精品| 卡戴珊不雅视频在线播放| 亚洲欧美一区二区三区黑人 | 久久久成人免费电影| 国内少妇人妻偷人精品xxx网站| 国产爱豆传媒在线观看| 22中文网久久字幕| 欧美成人午夜免费资源| 国产激情偷乱视频一区二区| 国产成人a∨麻豆精品| 国产精品国产三级国产av玫瑰| 亚洲精品日本国产第一区| 深爱激情五月婷婷| 最新中文字幕久久久久| 久久久久久久大尺度免费视频| 婷婷色av中文字幕| 婷婷色av中文字幕| 成人毛片60女人毛片免费| 免费大片黄手机在线观看| 三级男女做爰猛烈吃奶摸视频| 久久精品久久久久久久性| 亚洲高清免费不卡视频| av国产久精品久网站免费入址| 日韩av在线免费看完整版不卡| 99久久精品一区二区三区| 99re6热这里在线精品视频| 成人午夜精彩视频在线观看| 一级毛片我不卡| 色吧在线观看| 久久久久久久午夜电影| 日本猛色少妇xxxxx猛交久久| 免费电影在线观看免费观看| 久久久久久久久久黄片| 久久久a久久爽久久v久久| 国产乱来视频区| 久久精品夜色国产| av在线亚洲专区| 伦精品一区二区三区| 内射极品少妇av片p| 欧美日本视频| 亚洲精品成人av观看孕妇| 国产亚洲精品久久久com| av免费在线看不卡| 国产免费一级a男人的天堂| 九九久久精品国产亚洲av麻豆| 中文字幕亚洲精品专区| 亚洲欧美中文字幕日韩二区| 大香蕉97超碰在线| 51国产日韩欧美| 亚洲久久久久久中文字幕| 亚洲av免费高清在线观看| 免费在线观看成人毛片| 97人妻精品一区二区三区麻豆| 国产精品综合久久久久久久免费| 亚洲国产精品专区欧美| 有码 亚洲区| 国产 一区 欧美 日韩| 国产三级在线视频| 18禁动态无遮挡网站| 午夜精品一区二区三区免费看| 亚洲欧美成人精品一区二区| 国产 一区 欧美 日韩| 国产激情偷乱视频一区二区| 欧美+日韩+精品| 久久久久久久久久久免费av| 亚洲av国产av综合av卡| 美女大奶头视频| 成人特级av手机在线观看| 国产在视频线精品| 欧美成人午夜免费资源| 国产黄a三级三级三级人| 美女xxoo啪啪120秒动态图| 美女黄网站色视频| 国产av国产精品国产| 亚洲熟妇中文字幕五十中出| 一级毛片久久久久久久久女| 两个人视频免费观看高清| 免费观看的影片在线观看| 亚洲精品中文字幕在线视频 | 国产av在哪里看| 国产精品一及| 成人无遮挡网站| 国产男人的电影天堂91| 国产精品美女特级片免费视频播放器| 久久精品熟女亚洲av麻豆精品 | 日韩一区二区视频免费看| 国产色爽女视频免费观看| 国产精品久久久久久久电影| 91午夜精品亚洲一区二区三区| 国产精品蜜桃在线观看| 亚洲精品一区蜜桃| 22中文网久久字幕| 高清日韩中文字幕在线| 国产精品1区2区在线观看.| 99热这里只有是精品50| 人妻一区二区av| 极品少妇高潮喷水抽搐| 久久亚洲国产成人精品v| 久久久久久久久中文| 成人午夜精彩视频在线观看| 身体一侧抽搐| 中文精品一卡2卡3卡4更新| ponron亚洲| 久久这里只有精品中国| 免费大片18禁| www.色视频.com| 欧美激情久久久久久爽电影| 色网站视频免费| 97人妻精品一区二区三区麻豆| 久久国产乱子免费精品| 国国产精品蜜臀av免费| 国产成人a∨麻豆精品| 美女内射精品一级片tv| 亚洲国产av新网站| 观看美女的网站| 干丝袜人妻中文字幕| 日本熟妇午夜| 搡女人真爽免费视频火全软件| 天美传媒精品一区二区| 中国国产av一级| 欧美丝袜亚洲另类| 国产精品人妻久久久影院| 日韩欧美精品免费久久| 又爽又黄a免费视频| 日韩电影二区| 高清视频免费观看一区二区 | 国产精品三级大全| 乱系列少妇在线播放| 精品99又大又爽又粗少妇毛片| 91久久精品国产一区二区成人| 国产v大片淫在线免费观看| 亚洲欧美成人精品一区二区| 亚洲成人精品中文字幕电影| 久久久久久国产a免费观看| 丰满人妻一区二区三区视频av| 一个人看的www免费观看视频| 国产白丝娇喘喷水9色精品| 国产av不卡久久| 久久久久久国产a免费观看| 精品一区二区三区视频在线| 亚洲图色成人| 国产有黄有色有爽视频| 日日撸夜夜添| 国产高清国产精品国产三级 | 免费大片黄手机在线观看| 久久久成人免费电影| 69av精品久久久久久| 亚洲最大成人手机在线| 久久久久久久国产电影| 天堂中文最新版在线下载 | 少妇人妻精品综合一区二区| 99热6这里只有精品| 特级一级黄色大片| 国产精品一区二区性色av| 日韩强制内射视频| 好男人在线观看高清免费视频| 久久草成人影院| 免费看光身美女| 你懂的网址亚洲精品在线观看| 99久久中文字幕三级久久日本| 99久国产av精品国产电影| 国产成人精品一,二区| 国产高潮美女av| 国产精品国产三级国产av玫瑰| 春色校园在线视频观看| 看黄色毛片网站| 日本黄大片高清| a级毛片免费高清观看在线播放| 亚洲精品国产av蜜桃| 免费av不卡在线播放| 久久亚洲国产成人精品v| 亚洲精品视频女| 午夜福利在线观看吧| 秋霞伦理黄片| 黑人高潮一二区| 乱系列少妇在线播放| 久久久久久久久中文| 日本wwww免费看| 国产片特级美女逼逼视频| 美女高潮的动态| 天天一区二区日本电影三级| 日产精品乱码卡一卡2卡三| 国模一区二区三区四区视频| 亚洲精华国产精华液的使用体验| 一级片'在线观看视频| 超碰97精品在线观看| 18+在线观看网站| 亚洲国产av新网站| 久99久视频精品免费| 成人亚洲精品一区在线观看 | 久久国产乱子免费精品| 欧美成人a在线观看| av专区在线播放| 国产日韩欧美在线精品| 国产日韩欧美在线精品| 久久久久久九九精品二区国产| 色综合站精品国产| 亚洲图色成人| 亚洲精品,欧美精品| 午夜精品在线福利| 免费av观看视频| 久久99蜜桃精品久久| 又粗又硬又长又爽又黄的视频| 最近2019中文字幕mv第一页| 中文字幕av在线有码专区| 国产av在哪里看| 日本猛色少妇xxxxx猛交久久| 男插女下体视频免费在线播放| 别揉我奶头 嗯啊视频| 男的添女的下面高潮视频| 国产成人a∨麻豆精品| 午夜精品国产一区二区电影 | 国产精品久久久久久久久免| 亚洲国产成人一精品久久久| 午夜视频国产福利| 亚洲成人久久爱视频| 伊人久久精品亚洲午夜| 寂寞人妻少妇视频99o| 午夜福利高清视频| 国产欧美另类精品又又久久亚洲欧美| 国产精品美女特级片免费视频播放器| 免费av不卡在线播放| 亚洲国产精品成人综合色| 十八禁国产超污无遮挡网站| 别揉我奶头 嗯啊视频| kizo精华| 国产不卡一卡二| 亚洲四区av| 久久久久久久久久久免费av| 日韩成人av中文字幕在线观看| 免费观看精品视频网站| 自拍偷自拍亚洲精品老妇| 久久久久九九精品影院| 色综合色国产| or卡值多少钱| 国产激情偷乱视频一区二区| xxx大片免费视频| 尤物成人国产欧美一区二区三区| 最近的中文字幕免费完整| 精品国产露脸久久av麻豆 | 简卡轻食公司| 在线天堂最新版资源| 91在线精品国自产拍蜜月| 中文资源天堂在线| 国产av不卡久久| 国产精品麻豆人妻色哟哟久久 | 亚洲精品乱码久久久v下载方式| 久久99蜜桃精品久久| 久99久视频精品免费| 国产免费福利视频在线观看| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区在线观看99 | 成人鲁丝片一二三区免费| 亚洲av中文av极速乱| 好男人在线观看高清免费视频| 看十八女毛片水多多多| 色哟哟·www| 天堂影院成人在线观看| 婷婷色综合大香蕉| 午夜爱爱视频在线播放| 亚洲第一区二区三区不卡| 两个人视频免费观看高清| 秋霞伦理黄片| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品国产av蜜桃| 日韩精品青青久久久久久| 欧美日韩国产mv在线观看视频 | 国产一区二区亚洲精品在线观看| 久久久久久久亚洲中文字幕| 国产精品国产三级国产av玫瑰| 久久精品国产亚洲av天美| 久久精品久久久久久噜噜老黄| 99热6这里只有精品| 中文字幕av成人在线电影| 久久久色成人| 大陆偷拍与自拍| 性插视频无遮挡在线免费观看| 黄色一级大片看看| 成人漫画全彩无遮挡| 国产激情偷乱视频一区二区| 九色成人免费人妻av| 亚洲在线观看片| 亚洲伊人久久精品综合| 亚洲精品日本国产第一区| 精品国产露脸久久av麻豆 | 国产精品久久久久久久久免| 国产精品综合久久久久久久免费| 久久韩国三级中文字幕| 日韩一区二区视频免费看| 亚洲av免费在线观看| 天堂影院成人在线观看| 草草在线视频免费看| 蜜桃亚洲精品一区二区三区| 国产精品福利在线免费观看| 在线免费观看不下载黄p国产| 国产乱人视频| 国产美女午夜福利| 天天躁夜夜躁狠狠久久av| 国产人妻一区二区三区在| 免费观看av网站的网址| 99久久精品国产国产毛片| 免费av不卡在线播放| 亚洲精品影视一区二区三区av| 国产精品一区二区三区四区免费观看| 乱人视频在线观看| 蜜桃亚洲精品一区二区三区| 国产黄频视频在线观看| 成人鲁丝片一二三区免费| 久久久亚洲精品成人影院| 久久国内精品自在自线图片| 欧美最新免费一区二区三区| 亚洲婷婷狠狠爱综合网| 亚洲av二区三区四区| 18禁裸乳无遮挡免费网站照片| 久久亚洲国产成人精品v| 天堂√8在线中文| 亚洲在线自拍视频| 国产精品人妻久久久影院| 久久久久久久亚洲中文字幕| 少妇的逼好多水| 久久久久久久久中文| 黑人高潮一二区| 久久久a久久爽久久v久久| 五月玫瑰六月丁香| 日本一本二区三区精品| 床上黄色一级片| 欧美一级a爱片免费观看看| 久久精品夜夜夜夜夜久久蜜豆| 免费看日本二区| 亚洲欧美成人综合另类久久久| 老女人水多毛片| 亚洲综合色惰| 欧美高清性xxxxhd video| 亚洲精品成人久久久久久| 一级av片app| 亚洲精品,欧美精品| 最近2019中文字幕mv第一页| 91久久精品电影网| 亚洲精品成人av观看孕妇| 黄片无遮挡物在线观看| 91精品国产九色| 久久99热这里只有精品18| 亚洲美女搞黄在线观看| 黄片wwwwww| 国产成人aa在线观看| 国产女主播在线喷水免费视频网站 | 亚洲自偷自拍三级| 又爽又黄a免费视频| 国产一区二区在线观看日韩| 久久久久久久国产电影| 亚洲欧洲国产日韩| 久久精品久久久久久久性| av在线观看视频网站免费| 国产麻豆成人av免费视频| 我的老师免费观看完整版| 老司机影院成人| 欧美不卡视频在线免费观看| 欧美成人一区二区免费高清观看| 在线天堂最新版资源| 男人和女人高潮做爰伦理| 91久久精品电影网| 久久精品国产亚洲网站| 久久国产乱子免费精品| 别揉我奶头 嗯啊视频| 精品午夜福利在线看| 国产探花在线观看一区二区| 国产精品一区二区三区四区久久| 十八禁国产超污无遮挡网站| 亚洲熟妇中文字幕五十中出| 国产淫语在线视频| 亚洲av电影在线观看一区二区三区 | 午夜福利视频1000在线观看| 日韩av不卡免费在线播放| 三级国产精品片| 亚洲欧美成人综合另类久久久| 国产色婷婷99| 床上黄色一级片| 国产在线男女| 黑人高潮一二区| 日本一本二区三区精品| 国产精品国产三级国产av玫瑰| 啦啦啦韩国在线观看视频| 日韩强制内射视频| 国产毛片a区久久久久| 亚洲最大成人手机在线| 麻豆久久精品国产亚洲av| 国产人妻一区二区三区在| 亚洲精品aⅴ在线观看| 精品一区二区三区视频在线| 纵有疾风起免费观看全集完整版 | 久久久久久久亚洲中文字幕| 搡老乐熟女国产| 国产亚洲精品av在线| 六月丁香七月| 免费观看在线日韩| 久久久久久久久大av| 国产亚洲av嫩草精品影院| 亚洲av中文av极速乱| 亚洲不卡免费看| 免费电影在线观看免费观看| 亚洲在线观看片| 国产黄a三级三级三级人| 国产激情偷乱视频一区二区| 亚洲av一区综合| 人人妻人人澡人人爽人人夜夜 | 免费看日本二区|