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

    非對(duì)稱(chēng)通道內(nèi)親疏水結(jié)構(gòu)影響下的納米氣泡滑移效應(yīng)

    2022-03-03 01:10:28戴雙武盧艷高友明李揚(yáng)帆
    表面技術(shù) 2022年2期
    關(guān)鍵詞:壁面氣泡流體

    戴雙武,盧艷,高友明,李揚(yáng)帆

    非對(duì)稱(chēng)通道內(nèi)親疏水結(jié)構(gòu)影響下的納米氣泡滑移效應(yīng)

    戴雙武,盧艷,高友明,李揚(yáng)帆

    (武漢科技大學(xué) a.冶金裝備及其控制教育部重點(diǎn)實(shí)驗(yàn)室 b.武漢科技大學(xué)機(jī)械傳動(dòng)與制造工程湖北省重點(diǎn)實(shí)驗(yàn)室,武漢 430081)

    研究非對(duì)稱(chēng)性通道中親疏水表面結(jié)構(gòu)影響下納米氣泡特征與邊界滑移之間的關(guān)系,以實(shí)現(xiàn)良好的流體減阻效果。采用二元體系分子動(dòng)力學(xué)方法,研究納米氣泡在通道流動(dòng)中產(chǎn)生的滑移減阻效應(yīng)。首先建立上下壁面非對(duì)稱(chēng)微通道模型,通過(guò)考慮微通道流動(dòng)傳熱過(guò)程,探究納米氣泡影響下的微通道界面速度滑移現(xiàn)象。保持親水下壁面高度以及上下壁面溫差不變的情況下,增加上壁面納米結(jié)構(gòu)高度,對(duì)通道中納米氣泡體積產(chǎn)生促進(jìn)作用。另外,當(dāng)上壁面為疏水壁面時(shí),氣泡呈現(xiàn)為壁面形式,并且隨著體積增大,相對(duì)應(yīng)通道中上壁面滑移長(zhǎng)度增大;當(dāng)上壁面為親水壁面時(shí),納米氣泡呈現(xiàn)為體相形式,并且隨著體積增大,對(duì)應(yīng)上壁面滑移長(zhǎng)度減小。非對(duì)稱(chēng)性通道內(nèi),在上壁面結(jié)構(gòu)高度影響下,壁面形式的氣泡體積增大對(duì)通道內(nèi)減阻具有促進(jìn)作用,而體相形式的氣泡體積增大對(duì)通道內(nèi)減阻具有抑制作用。

    微通道;溫度階躍;速度滑移;納米氣泡;滑移減阻;潤(rùn)濕性

    近年來(lái),隨著科技的發(fā)展,微納科學(xué)已經(jīng)應(yīng)用到了各個(gè)工程領(lǐng)域[1-4]。當(dāng)熱流體流過(guò)納米通道時(shí),通道內(nèi)的能量傳遞和物質(zhì)運(yùn)輸都會(huì)在受限制的微小尺度空間內(nèi)進(jìn)行,因此納米通道的固液界面會(huì)產(chǎn)生與宏觀尺度不同的現(xiàn)象,宏觀尺度固液界面的理論不再適用。因?yàn)樵诤暧^尺度下,表面效應(yīng)對(duì)通道固液界面流體流動(dòng)性的影響較小,并且宏觀通道流體流動(dòng)性的研究都是基于無(wú)邊界滑移情況下進(jìn)行的,而微納尺度通道中,表面效應(yīng)則不可忽視,而且會(huì)影響到整個(gè)通道中流體的流動(dòng)性[5-7]。研究發(fā)現(xiàn),由于表面效應(yīng)作用,流體流過(guò)微通道固體壁面時(shí),會(huì)存在速度滑移[8],且具有一定的界面速度滑移長(zhǎng)度[9]。Zhang等人[10-11]對(duì)納米通道中流體的研究發(fā)現(xiàn),固體壁面的潤(rùn)濕性會(huì)影響流體運(yùn)動(dòng),并發(fā)現(xiàn)固液相互作用力較差時(shí),界面會(huì)出現(xiàn)比較明顯的速度滑移。Voronov等[12]利用分子動(dòng)力學(xué)方法研究納米通道中的流體行為時(shí),發(fā)現(xiàn)固液界面的滑移長(zhǎng)度與固體表面的潤(rùn)濕性相關(guān),對(duì)于潤(rùn)濕性較差的固體表面,滑移長(zhǎng)度會(huì)隨潤(rùn)濕性的降低而增加。曹炳陽(yáng)等人[13]通過(guò)研究通道中流體滑移現(xiàn)象發(fā)現(xiàn),固體壁面的表面勢(shì)能越強(qiáng),界面速度滑移長(zhǎng)度越小。表面潤(rùn)濕性不同,固液界面可能發(fā)生負(fù)滑移、無(wú)滑移以及正滑移等不同現(xiàn)象。另外,固體表面粗糙度大小也能影響固液之間的潤(rùn)濕性變化[14-15]。Pit等[16]研究發(fā)現(xiàn),固體表面粗糙度能夠影響固體表面的速度滑移。Choi等[17]在研究牛頓流體的減阻機(jī)理時(shí)發(fā)現(xiàn),固體表面粗糙度可以使界面滑移長(zhǎng)度增大,結(jié)構(gòu)變化會(huì)影響粗糙度的變化,進(jìn)而影響滑移長(zhǎng)度。張程賓等人[18]利用仿真方法對(duì)比光滑以及粗糙固體表面流體的流動(dòng)發(fā)現(xiàn),粗糙結(jié)構(gòu)的存在限制了流體的運(yùn)動(dòng),導(dǎo)致固液界面速度滑移長(zhǎng)度減小。Ou等人[19]通過(guò)合理設(shè)計(jì)疏水固體表面的結(jié)構(gòu)尺寸來(lái)探討固體表面的減阻性能時(shí)發(fā)現(xiàn),對(duì)比光滑表面,疏水表面微結(jié)構(gòu)的存在,使得界面速度滑移長(zhǎng)度較大。Rahmatipour等[20]利用分子動(dòng)力學(xué)方法對(duì)粗糙通道進(jìn)行研究時(shí)發(fā)現(xiàn),納米結(jié)構(gòu)高度越高,壁面滑移長(zhǎng)度越長(zhǎng)。其次,研究表明,當(dāng)流體流過(guò)粗糙固體表面時(shí),將會(huì)有氣體滯留在微織構(gòu)內(nèi)部,并在其表面形成穩(wěn)定的氣泡,從而使固液之間產(chǎn)生滑移減阻效應(yīng)[21-22]。Elias等人[23-24]在研究通道中流體的運(yùn)動(dòng)時(shí)發(fā)現(xiàn),氣泡的產(chǎn)生是實(shí)現(xiàn)通道固體表面滑移減阻的關(guān)鍵因素。

    溫度也是影響納米通道中氣泡成核的因素之一[25]。Lin等人[26-27]通過(guò)微通道中加熱絲加熱流體時(shí)發(fā)現(xiàn),只有當(dāng)微通道中液體溫度接近臨界溫度時(shí),才會(huì)有熱氣泡產(chǎn)生。Yang等[28]研究結(jié)果表明,升高疏水表面的溫度,在表面粗糙的地方更易形成納米氣泡。張雪花等[29]在研究固體表面上已經(jīng)形成的納米氣泡時(shí)發(fā)現(xiàn),當(dāng)液體溫度升高到某個(gè)定值時(shí),納米氣泡體積達(dá)到最大,然后隨溫度的升高而減小。劉樺等[30]發(fā)現(xiàn),固體表面溫度影響了通道中流體滑移減阻。

    綜上所述,關(guān)于納米通道中流體的滑移減阻已有大量研究,大部分是考慮納米氣泡影響下對(duì)稱(chēng)性管道中的流體運(yùn)動(dòng),且關(guān)于溫度對(duì)納米氣泡的研究也都處于實(shí)驗(yàn)階段。然而,非對(duì)稱(chēng)性管道的應(yīng)用卻十分廣泛[31-32],因此研究非對(duì)稱(chēng)通道內(nèi)納米氣泡的滑移減阻效應(yīng)具有重要意義。目前,通過(guò)調(diào)整非對(duì)稱(chēng)通道的微結(jié)構(gòu)來(lái)改變通道表面粗糙度和流體溫度,從而影響納米氣泡特性的通道減阻理論研究并不完善。本文首先采用分子動(dòng)力學(xué)模型,主要研究了非對(duì)稱(chēng)性通道內(nèi)上壁面結(jié)構(gòu)高度影響下納米氣泡與流體流動(dòng)之間的關(guān)系,通過(guò)控制上下壁面溫差、潤(rùn)濕性,并保持下壁面結(jié)構(gòu)高度不變,探討了改變疏水或者親水上壁面結(jié)構(gòu)高度,通道中納米氣泡體積的變化情況。然后,進(jìn)一步研究了納米氣泡體積與通道中界面滑移長(zhǎng)度之間的關(guān)系。

    1 微通道流體傳熱仿真模型

    本文采用分子動(dòng)力學(xué)(MD)模擬方法,利用開(kāi)源程序LAMMPS仿真軟件來(lái)研究納米通道內(nèi)納米氣泡成核體積與壁面能量傳遞之間的關(guān)系。非對(duì)稱(chēng)模型通道如圖1所示。

    圖1 非對(duì)稱(chēng)通道模型

    圖1a為非對(duì)稱(chēng)納米通道三視圖,模擬微通道在、和方向上的尺寸分別為9、2.6、6 nm。二維模型(圖1b)中,上壁面藍(lán)色粒子為疏水性鉑原子,下壁面紅色粒子為親水性鉑原子,中間黃色粒子為流體氬原子。流體氬原子被限制在兩個(gè)平行的鉑固體壁面之間,固體和液體原子的初始狀態(tài)都是面心立方體FCC晶格排列(固體粒子的晶格常數(shù)為0.392 nm,流體粒子的晶格常數(shù)為0.543 nm)。為了提高模型的真實(shí)性和保持系統(tǒng)的穩(wěn)定性,和方向設(shè)為周期性邊界,方向設(shè)為固定邊界。

    圖1b為非對(duì)稱(chēng)通道模型的正視圖,固體壁面分為內(nèi)外兩層,外層被固定,以保持系統(tǒng)的穩(wěn)定性(壁面外兩層原子標(biāo)記為固定層),內(nèi)層為浴熱層(下壁面內(nèi)兩層原子為高溫浴熱層,上壁面內(nèi)兩層原子為低溫浴熱層)。上壁面包括上壁面固定層、低溫浴熱層、結(jié)構(gòu),定義為top;下壁面包括下壁面固定層、高溫浴熱層、結(jié)構(gòu),定義為down。下壁面結(jié)構(gòu)高度定義為down,上壁面結(jié)構(gòu)高度定義為top。本文通過(guò)改變親疏水壁面結(jié)構(gòu)高度參數(shù)(top)來(lái)展開(kāi)通道減阻研究。另外,液體由2100個(gè)氬原子組成,上下固體壁面大約2800個(gè)鉑原子。Lennard-Jones(LJ)勢(shì)能描述了粒子之間的相互作用力,見(jiàn)式(1)。

    設(shè)置勢(shì)函數(shù)的截止距離=0.85 nm,仿真系統(tǒng)的時(shí)間步長(zhǎng)為0.01 ps。首先采用速度校定法將整個(gè)系統(tǒng)NVE馳豫,經(jīng)過(guò)2 ps達(dá)到平衡狀態(tài)。接著對(duì)系統(tǒng)浴熱,用NVT(粒子數(shù)恒定,體積、溫度保持不變)對(duì)上、下固體壁面進(jìn)行控溫,將上壁面溫度top控制在86 K,下壁面溫度down控制在100 K,同時(shí)對(duì)液體施加一個(gè)方向的恒定驅(qū)動(dòng)力(驅(qū)=0.001 eV/nm),并將系統(tǒng)加熱時(shí)間由2 ps延長(zhǎng)到9 ps,以達(dá)到平衡狀態(tài)。最后對(duì)整個(gè)系統(tǒng)進(jìn)行切片,并統(tǒng)計(jì)數(shù)據(jù)。

    表1 原子之間的相互作用力勢(shì)能參數(shù)及其作用距離[34]

    Tab.1 Potential parameters and interaction distancesof interaction force between atoms[34]

    2 疏水壁面結(jié)構(gòu)高度影響下的納米氣泡

    2.1 納米氣泡成核

    保持down=1 nm、top=86 K、down=100 K不變,改變疏水上壁面的結(jié)構(gòu)高度(0~0.8 nm)時(shí),上壁面納米氣泡的成核規(guī)律如圖2所示??梢钥闯?,上壁面的高度由0 nm增大到0.8 nm時(shí),近下壁面附近的流體被吸附,氣泡在近上壁面區(qū)域出現(xiàn),并且上壁面附近流體中納米氣泡的體積越來(lái)越大。這是因?yàn)槭杷员砻娼Y(jié)構(gòu)高度增加,壁面疏水性增強(qiáng),使得上壁面的氣泡成核體積增大。

    2.2 通道流體密度分布

    改變結(jié)構(gòu)高度,對(duì)流體中氣泡的形成有明顯影響,致使流體密度的分布有差異。為了研究通道中僅疏水上壁面結(jié)構(gòu)高度變化時(shí),通道中納米氣泡與流體密度之間的關(guān)系,將模型進(jìn)行切層,統(tǒng)計(jì)密度。待系統(tǒng)粒子分布均勻后,再對(duì)系統(tǒng)進(jìn)行切層處理,并將模型分為70層。

    圖2 上壁面結(jié)構(gòu)高度影響下的氣泡成核

    通道中流體密度沿方向的分布如圖3所示。由于體系中液體與固體粒子之間的勢(shì)能作用域?yàn)榧{米級(jí)別,因此固體壁面只影響近壁面流體區(qū)域的密度分布特性。由圖3可知,由于固體表面效應(yīng),近壁面流體區(qū)域密度分布出現(xiàn)振蕩現(xiàn)象。在通道中心主流區(qū)域,流體粒子受固體粒子的影響較小,密度幾乎不發(fā)生變化。由于上下壁面的潤(rùn)濕性不同,以及上下壁面結(jié)構(gòu)不對(duì)稱(chēng),所以通道內(nèi)流體密度呈不對(duì)稱(chēng)分布。當(dāng)down=1 nm、top=86 K、down=100 K時(shí),通過(guò)改變上壁面的結(jié)構(gòu)高度,來(lái)分析流體密度的分布規(guī)律。結(jié)果發(fā)現(xiàn),上壁面結(jié)構(gòu)高度增加,近下壁面區(qū)域流體分布基本不變,近上壁面區(qū)域流體密度呈近乎直線的下降趨勢(shì)。這是由于增加上壁面結(jié)構(gòu)高度,近上壁面區(qū)域流體中納米氣泡成核體積變大。

    圖3 疏水上壁面結(jié)構(gòu)高度影響下流體沿z方向的密度分布

    2.3 疏水界面氣泡影響下微通道流體的流動(dòng)性

    通道內(nèi)壁結(jié)構(gòu)高度變化會(huì)影響流體中納米氣泡的生成,導(dǎo)致流體密度分布發(fā)生變化,從而影響通道中流體的流動(dòng)傳熱。為了研究?jī)H通道上壁面結(jié)構(gòu)高度變化時(shí),納米氣泡體積與流體速度分布之間的關(guān)系,取3組數(shù)據(jù)進(jìn)行分析。當(dāng)top=86 K、down=100 K、down=1 nm時(shí),改變上壁面結(jié)構(gòu)高度,通道中流體的速度分布如圖4所示。結(jié)果顯示,在通道的中心區(qū)域,流體速度整體呈不對(duì)稱(chēng)分布。在上壁面結(jié)構(gòu)高度變化影響下,近上下壁面區(qū)域的流體速度、流體平均速度和納米氣泡體積見(jiàn)表2??梢苑治龅贸觯S著上壁面結(jié)構(gòu)高度的增加,上壁面納米氣泡體積變大,近上下壁面的流體速度正向增大,通道流體平均速度變大。這是由于上壁面流體中氣泡成核體積增大,減小了下壁面固液接觸面積,使得流體與固體之間的運(yùn)動(dòng)阻力降低。

    圖4 納米氣泡影響下流體沿z方向的速度分布

    表2 不同疏水上壁面高度影響下近上下壁面區(qū)域流體速度與板面的速度差、流體平均速度大小和納米氣泡體積

    Tab.2 Difference between the velocity of the fluid near the upper and lower walls and the velocity of the platesurface, the average velocity of the fluid and the volume of the nanobubbles on different height of drain upper wall

    2.4 氣泡影響下流體溫度分布

    在恒定外力驅(qū)動(dòng)下,流體在納米通道中流動(dòng)時(shí),流體粒子之間、流體粒子與粗糙固體表面之間都會(huì)存在摩擦,致使流體溫度變化,進(jìn)而影響納米氣泡的變化。為了討論改變疏水上壁面結(jié)構(gòu)高度變化時(shí),納米氣泡體積與流體溫度分布之間的關(guān)系,采用類(lèi)似的模型處理方式,選取具有代表性的幾組數(shù)據(jù)進(jìn)行分析。當(dāng)top=86 K、down=100 K、down=1 nm時(shí),改變上壁面結(jié)構(gòu)高度,流體溫度的分布狀況如圖5所示。可以看到,增加上壁面結(jié)構(gòu)高度,通道中流體溫度的平均水平下降。在上下壁面溫度以及下壁面結(jié)構(gòu)高度不變時(shí),改變上壁面結(jié)構(gòu)高度,通道中近上下壁面的流體溫差、平均流體溫度和納米氣泡體積見(jiàn)表3。結(jié)合表3與圖5分析可以得到,隨著上壁面結(jié)構(gòu)高度的增加,上壁面氣泡體積增大,對(duì)應(yīng)近上下壁面流體溫差以及流體平均溫度都呈正向減小的趨勢(shì)。其原因是,疏水上壁面結(jié)構(gòu)高度增加,一方面促進(jìn)了通道中的傳熱,流體溫度降低;另一方面,增大了壁面疏水性,壁面氣泡體積增大。

    圖5 納米氣泡影響下流體沿z方向的溫度分布

    表3 不同疏水上壁面高度影響下近上下壁面流體溫度差、流體平均溫度、納米氣泡的體積

    Tab.3 Fluid temperature difference near the upper and lower walls, average fluid temperature, volume of nanobubbles on different height of drain upper wall

    3 親水上壁面結(jié)構(gòu)高度影響下的納米氣泡

    3.1 上壁面結(jié)構(gòu)高度對(duì)氣泡成核影響

    模擬中,模型處理與第2節(jié)一樣,僅將上壁面的疏水性鉑壁面換為親水性鉑壁面。將模型運(yùn)行至系統(tǒng)能量平衡后,輸出并記錄仿真結(jié)果。在top=86 K、down=100 K、down=1 nm下,改變親水上壁面結(jié)構(gòu)高度時(shí),微通道流體中納米氣泡的成核規(guī)律如圖6所示。從圖6中可以看出,top為0、0.2、0.4、0.6、0.8 nm時(shí),氣泡成核出現(xiàn)在微通道主流區(qū)域,并且隨著上壁面結(jié)構(gòu)高度的增加,通道流體中納米氣泡體積越來(lái)越大。這是由于上壁面的結(jié)構(gòu)高度增加,使得壁面對(duì)流體粒子的吸引面積增大,通道中央流體出現(xiàn)密度不均勻的現(xiàn)象越來(lái)越明顯,因而氣泡越來(lái)越大。

    3.2 氣泡影響下通道流體密度分布

    通過(guò)在模型方向施加恒定力,保持流體運(yùn)動(dòng),同時(shí)改變上壁面結(jié)構(gòu)高度,通道中氣泡的生成就會(huì)發(fā)生變化,那么通道中粒子分布也會(huì)發(fā)生改變。為了研究?jī)H親水上壁面結(jié)構(gòu)高度變化時(shí),通道中氣泡與流體密度之間的關(guān)系,將模型進(jìn)行切層,統(tǒng)計(jì)密度,待系統(tǒng)粒子分布均勻后,再對(duì)系統(tǒng)進(jìn)行切層處理。當(dāng)top=86 K、down=100 K、down=1 nm時(shí),改變上壁面結(jié)構(gòu)高度,通道中流體密度的分布如圖7所示??梢钥闯?,當(dāng)上壁面為親水壁面時(shí),近壁面區(qū)域流體密度同樣存在震蕩現(xiàn)象,通道中央出現(xiàn)密度減小的現(xiàn)象。上下壁面的結(jié)構(gòu)不對(duì)稱(chēng),使得通道流體密度分布不對(duì)稱(chēng)。上壁面高度越高時(shí),相對(duì)應(yīng)近上壁面區(qū)域流體密度波動(dòng)越明顯。這是由于親水上壁面結(jié)構(gòu)高度增加,增大了近上壁面區(qū)域流體與親水鉑結(jié)構(gòu)之間的接觸面積。另外,上壁面結(jié)構(gòu)較高的通道內(nèi),通道中央流體密度大幅度減小。這是因?yàn)樯媳诿娼Y(jié)構(gòu)高度增加,通道中央氣泡的成核體積變大。

    圖6 親水上壁面結(jié)構(gòu)高度變化影響下的納米氣泡

    圖7 納米氣泡影響下的流體沿z方向的密度分布

    3.3 氣泡影響下通道流體的流動(dòng)性

    當(dāng)氣泡形成形式發(fā)生改變時(shí),流體密度分布發(fā)生突變,粒子間接觸形式發(fā)生改變,通道中的動(dòng)量傳遞也會(huì)發(fā)生變化。在top=86 K、down=100 K、down=1 nm時(shí),改變親水上壁面結(jié)構(gòu)高度,通道中的流體速度分布如圖8所示??梢钥吹?,速度整體呈不對(duì)稱(chēng)分布,并且上壁面高度增加,通道中流體平均速度降低。這是由于親水上壁面結(jié)構(gòu)高度增加,增大了固體壁面對(duì)流體粒子的吸附能力,干擾了流體運(yùn)動(dòng)。當(dāng)上下壁面溫差以及下壁面納米結(jié)構(gòu)高度不變時(shí),在上壁面結(jié)構(gòu)高度變化的影響下,近上下壁面區(qū)域的流體速度、通道流體平均速度和納米氣泡體積見(jiàn)表4??梢钥闯觯S著上壁面結(jié)構(gòu)高度的增加,通道中央納米氣泡體積不斷變大,近上下壁面區(qū)域流體速度正向減小,流體平均速度也不斷變小。這是由于通道流體中納米氣泡的成核體積增大,阻礙了流體流動(dòng)所致。

    圖8 納米結(jié)構(gòu)高度影響下流體沿z方向的速度分布

    3.4 氣泡影響下通道溫度分布

    流體的密度分布以及粒子間的動(dòng)量傳遞發(fā)生變化,流體與流體的黏性摩擦就會(huì)發(fā)生改變,通道中流體溫度也會(huì)發(fā)生變化。當(dāng)top=86 K、down=100 K、down=1 nm時(shí),改變上壁面結(jié)構(gòu)高度,流體溫度的分布狀況如圖9所示??梢钥吹剑黾由媳诿娼Y(jié)構(gòu)高度,通道中流體的平均溫度下降。這是因?yàn)楸诿娼Y(jié)構(gòu)高度增加,流體粒子與固體壁面的接觸面積增加,傳熱效率增大。在上下壁面溫度以及下壁面結(jié)構(gòu)高度不變時(shí),改變上壁面結(jié)構(gòu)高度,通道中近上下壁面的流體溫差、流體平均溫度和納米氣泡體積見(jiàn)表5。結(jié)合表5與圖9分析可得,隨著上壁面納米結(jié)構(gòu)高度的增加,通道流體中納米氣泡體積增大,對(duì)應(yīng)近上下壁面流體溫差及流體平均溫度都呈正向減小趨勢(shì)。這是因?yàn)橛H水上壁面結(jié)構(gòu)高增加,壁面的潤(rùn)濕性增大,通道中氣泡的成核體積增大。同時(shí),結(jié)構(gòu)高度增加,流體粒子的吸附面積增大,通道中界面的傳熱效率提高。

    表4 不同親水上壁面高度影響下近上下壁面流體速度與板面的速度差、流體平均速度大小和納米氣泡體積

    Tab.4 Difference between the velocity of the fluid near the upper and lower walls and the velocity of plate surface, average fluid velocity and the volume of the nanobubble on different hydrophilic upper wall height

    圖9 上壁面結(jié)構(gòu)高度影響下流體沿z方向的溫度分布

    表5 不同親水上壁面高度影響下近上下壁面流體溫度差、流體平均溫度、納米氣泡體積

    Tab.5 Relationship between fluid temperature difference near the upper and lower walls, the average fluid temperature, and the volume of the nanobubble on different hydrophilic upper wall height

    4 氣泡影響下流體中滑移

    改變通道中結(jié)構(gòu)高度,納米氣泡體積大小和流體密度、速度、溫度分布就會(huì)發(fā)生變化,界面能量傳導(dǎo)也會(huì)發(fā)生變化。對(duì)于疏水表面氣泡影響下滑移長(zhǎng)度與表面結(jié)構(gòu)粗糙度的關(guān)系,本文對(duì)比了文獻(xiàn)[35]的結(jié)果(如圖10所示),實(shí)驗(yàn)與仿真所得結(jié)果的變化規(guī)律一致,但存在差異。這是由于本文仿真通道中固液界面?zhèn)鳠岬拇嬖冢约拔墨I(xiàn)采用的表面材料等不同。

    為了探究親、疏水上壁面結(jié)構(gòu)高度變化影響下,通道中納米氣泡體積與界面速度滑移之間的關(guān)系,將所有數(shù)據(jù)進(jìn)行統(tǒng)計(jì),并按照式(2)計(jì)算速度滑移長(zhǎng)度。

    式中:Δv為近固體壁下的液體速度;為流體在固液界面下的速度梯度。

    當(dāng)top=86 K、down=100 K、down=1 nm時(shí),改變下壁面結(jié)構(gòu)高度,通道中納米氣泡體積與上下壁面速度滑移之間的關(guān)系如圖11所示。由圖11可知,保持下壁面結(jié)構(gòu)高度不變,當(dāng)上壁面為疏水壁面時(shí),增加上壁面結(jié)構(gòu)高度,上壁面附近流體中氣泡的體積增大,上壁面的速度滑移長(zhǎng)度增大。這是由于壁面納米氣泡的存在,降低了上壁面固液之間的接觸面積,動(dòng)量傳遞效率提高。當(dāng)上壁面為親水壁面時(shí),增加上壁面結(jié)構(gòu)高度,上壁面速度滑移的長(zhǎng)度減小。這是因?yàn)橥ǖ乐醒爰{米體積越大,對(duì)通道中流體流動(dòng)的阻礙越強(qiáng)。

    圖11 不同壁面結(jié)構(gòu)高度影響下氣泡體積與速度滑移的關(guān)系

    5 結(jié)論

    本文使用MD模擬研究了非對(duì)稱(chēng)納米通道內(nèi),上壁面為親、疏水壁面影響下納米氣泡成核體積的大小,然后進(jìn)一步探索了通道中納米氣泡與界面滑移行為之間的關(guān)系,得出以下結(jié)論:

    1)非對(duì)稱(chēng)通道中,當(dāng)保持下壁面結(jié)構(gòu)高度以及上下壁面溫差不變時(shí),增加上壁面結(jié)構(gòu)高度,通道中流體的平均溫度下降,納米氣泡的體積增大。

    2)非對(duì)稱(chēng)通道內(nèi),當(dāng)保持下壁面結(jié)構(gòu)高度以及上下壁面溫差不變,上壁面為疏水壁面時(shí),納米氣泡呈壁面形式,且疏水壁面附近納米氣泡體積的增大,對(duì)應(yīng)上壁面的速度滑移長(zhǎng)度正向增大,而下壁面的速度滑移長(zhǎng)度不變,通道中流體運(yùn)動(dòng)阻力變小。

    3)非對(duì)稱(chēng)通道內(nèi),當(dāng)保持下壁面結(jié)構(gòu)高度以及上下壁面溫差不變,上壁面為親水壁面時(shí),納米氣泡呈體相形式,通道中納米氣泡體積增大,上壁面的速度滑移長(zhǎng)度減小,而下壁面的速度滑移長(zhǎng)度不變,通道中流體運(yùn)動(dòng)阻力增大。

    [1] ZHANG Guang, XU Yan-tong, DUAN Zheng, et al. Enh-an-cement of Evaporative Heat Transfer on Carbon Nano-tube Sponges by Electric Field Reinforced Wettability[J]. Applied Surface Science, 2018, 454: 262-269.

    [2] CHOI S M, JANG S H, LEE Dong hun, et al. Spray Chara-cteristics of the Rotating Fuel Injection System of a Micro-Jet Engine[J]. Journal of Mechanical Science and Technology, 2010, 24(2): 551-558.

    [3] SOHRABI S, LIU Ya-ling. Modeling Thermal Inkjet and Cell Printing Process Using Modified Pseudopotential and Thermal Lattice Boltzmann Methods[J]. Physical Review E, 2018, 97(3-1): 033105.

    [4] ALTINTAS Z, AKGUN M, KOKTURK G, et al. A Fully Automated Microfluidic-Based Electrochemical Sensor for Real-Time Bacteria Detection[J]. Biosensors and Bioelec-tronics, 2018, 100: 541-548.

    [5] 黃橋高, 潘光, 宋保維. 疏水表面滑移流動(dòng)及減阻特性的格子Boltzmann方法模擬[J]. 物理學(xué)報(bào), 2014, 63(5): 236-242.

    HUANG Qiao-gao, PAN Guang, SONG Bao-wei. Lattice Boltzmann Simulation of Slip Flow and Drag Reduction Characteristics of Hydrophobic Surfaces[J]. Acta Physica Sinica, 2014, 63(5): 236-242.

    [6] 葛宋, 陳民. 接觸角與液固界面熱阻關(guān)系的分子動(dòng)力學(xué)模擬[J]. 物理學(xué)報(bào), 2013, 62(11): 110204.

    GE Song, CHEN Min. A Molecular Dynamics Simulation on the Relationship between Contact Angle and Solid- Liquid Interfacial Thermal Resistance[J]. Acta Physica Sinica, 2013, 62(11): 110204.

    [7] 顧驍坤, 陳民. 納米硅通道內(nèi)滑移現(xiàn)象的分子動(dòng)力學(xué)模擬[J]. 工程熱物理學(xué)報(bào), 2010, 31(10): 1724-1726.

    GU Xiao-kun, CHEN Min. Molecular Dynamics Simula-tions of the Slip Phenomena in Silicon Nano-channels[J]. Journal of Engineering Thermophysics, 2010, 31(10): 1724-1726.

    [8] JOSEPH P, TABELING P. Direct Measurement of the Apparent Slip Length[J]. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 2005, 71(3 Pt 2A): 035303.

    [9] HUANG Qiao-gao, PAN Guang, SONG Bao-wei, et al. Lattice Boltzmann Simulation of Slip Flow and Drag Re-duction Characteristics of Hydrophobic Surfaces[J]. Acta Physica Sinica, 2014, 63(5): 054701.

    [10] NAGAYAMA G, CHENG Ping. Effects of Interface Wet-tability on Microscale Flow by Molecular Dynamics Sim-ulation[J]. International Journal of Heat and Mass Tran-sfer, 2004, 47(3): 501-513.

    [11] ZHANG Hong-wu, ZHANG Zhong-qiang, YE Hong-fei. Molecular Dynamics-based Prediction of Boundary Slip of Fluids in Nanochannels[J]. Microfluidics and Nanoflui-dics, 2012, 12(1-4): 107-115.

    [12] VORONOV R S, PAPAVASSILIOU D V, LEE L L. Slip Length and Contact Angle over Hydrophobic Surfaces[J]. Chemical Physics Letters, 2007, 441(4-6): 273-276.

    [13] 曹炳陽(yáng), 陳民, 過(guò)增元. 納米通道內(nèi)液體流動(dòng)的滑移現(xiàn)象[J]. 物理學(xué)報(bào), 2006, 55(10): 5305-5310.

    CAO Bing-yang, CHEN Min, GUO Zeng-yuan. Velocity Slip of Liquid Flow in Nanochannels[J]. Acta Physica Sinica, 2006, 55(10): 5305-5310.

    [14] 王甜. 固著納米水滴潤(rùn)濕行為的分子動(dòng)力學(xué)模擬[D]. 大連: 大連理工大學(xué), 2015.

    WANG Tian. Wetting Behaviors of Sessile Water Nano- Droplets Using Molecular Dynamics Simulation[D]. Da-lian: Dalian University of Technology, 2015.

    [15] MARMUR A, BITTOUN E. When Wenzel and Cassie are Right: Reconciling Local and Global Considerations[J]. Langmuir: The ACS Journal of Surfaces and Colloids, 2009, 25(3): 1277-1281.

    [16] PIT R, HERVET H, LéGER L. Direct Experimental Evi-d-ence of Slip in Hexadecane: Solid Interfaces[J]. Physical Review Letters, 2000, 85(5): 980-983.

    [17] CHOI C H, ULMANELLA U, KIM J, et al. Effective Slip and Friction Reduction in Nanograted Superhydrophobic Microchannels[J]. Physics of Fluids, 2006, 18(8): 087105.

    [18] 張程賓, 許兆林, 陳永平. 粗糙納通道內(nèi)流體流動(dòng)與傳熱的分子動(dòng)力學(xué)模擬研究[J]. 物理學(xué)報(bào), 2014, 63(21): 263-270.

    ZHANG Cheng-bin, XU Zhao-lin, CHEN Yong-ping. Mo--lecular Dynamics Simulation on Fluid Flow and Heat Transfer in Rough Nano Channels[J]. Acta Physica Sinica, 2014, 63(21): 263-270.

    [19] OU Jia, PEROT B, ROTHSTEIN J P. Laminar Drag Red-u-ction in Microchannels Using Ultrahydrophobic Surfa-ces[J]. Physics of Fluids, 2004, 16(12): 4635-4643.

    [20] RAHMATIPOUR H, AZIMIAN A. Molecular Dynamics Simulation of Couette Flow Behavior in Smooth and Ro-ugh Nanochannels[J]. Physica A, 2015, 15(7): 119-130.

    [21] LUM K, CHANDLER D, WEEKS J D. Hydrophobicity at Small and Large Length Scales[J]. The Journal of Phy-sical Chemistry B, 1999, 103(22): 4570-4577.

    [22] OU Jia, PEROT B, ROTHSTEIN J P. Laminar Drag Re-duction in Microchannels Using Ultrahydrophobic Sur-faces[J]. Physics of Fluids, 2004, 16(12): 4635-4643.

    [23] ALJALLIS E, SARSHAR M A, DATLA R, et al. Exper-imental Study of Skin Friction Drag Reduction on Super-h-y-drophobic Flat Plates in High Reynolds Number Boun-dary Layer Flow[J]. Physics of Fluids, 2013, 25(2): 025103.

    [24] JAGDISH B N, XIAN BRANDON T Z, KWEE T J, et al. Experimental Study of Air Layer Sustainability for Frict-ional Drag Reduction[J]. Journal of Ship Research, 2014, 58(1): 30-42.

    [25] LOHSE D, ZHANG Xue-hua. Pinning and Gas Oversatu-ration Imply Stable Single Surface Nanobubbles[J]. Phy-sical Review E, Statistical, Nonlinear, and Soft Matter Ph-y-sics, 2015, 91(3): 031003.

    [26] LIN Li-wei. Microscale Thermal Bubble Formation: The-r-m-ophysical Phenomena and Applications[J]. Microscale Thermophysical Engineering, 1998, 2(2): 71-85.

    [27] LIN L, UDELL K S, PISANO A P. Liquid-V Apor Phase Transition and Bubble Formation in Micro Structures[J]. Thermal science engineering, 1994, 2(1): 52-59.

    [28] YANG Shang-jiong, DAMMER S M, BREMOND N, et al.Characterization of Nanobubbles on Hydrophobic Surfaces in Water[J]. Langmuir: The ACS Journal of Surfaces and Colloids, 2007, 23(13): 7072-7077.

    [29] ZHANG X H, LI G, WU Z H, et al. Effect of Temperature on the Morphology of Nanobubbles at Mica/Water Inter-face[J]. Chinese Physics B, 2005, 14(9): 1774-1778.

    [30] 劉樺, 孫仁. 控溫、系統(tǒng)溫度和外載荷對(duì)納米通道中水流性質(zhì)的影響[J]. 水動(dòng)力學(xué)研究與進(jìn)展, 2020, 35(6): 719-725.

    LIU Hua, SUN Ren. Effects of Temperature Control, Sys-tem Temperature and External Load on Properties of Water Flow Inside Nanochannel[J]. Chinese Journal of Hydro-d-y-na-mics, 2020, 35(6): 719-725.

    [31] CHEN Qin-wen, MENG Ling-yi, LI Qi-kai, et al. Water Transport and Purification in Nanochannels Controlled by Asymmetric Wettability[J]. Small, 2011, 7(15): 2225-2231.

    [32] DE GANS B J, KITA R, MüLLER B, et al. Negative Thermodiffusion of Polymers and Colloids in Solvent Mixtures[J]. The Journal of Chemical Physics, 2003, 118(17): 8073-8081.

    [33] 張龍艷. 微尺度下流體的流動(dòng)換熱及核化沸騰的分子動(dòng)力學(xué)研究[D]. 北京: 華北電力大學(xué)(北京), 2019.

    ZHANG Long-yan. Molecular Dynamics Simulation of Fluid Flow and Heat Transfer and Nucleate Boiling in Microscale[D]. Beijing: North China Electric Power Uni-versity, 2019.

    [34] 王艷紅. 潤(rùn)濕性對(duì)納米尺度薄液膜爆炸沸騰現(xiàn)象影響的分子動(dòng)力學(xué)模擬[D]. 北京: 華北電力大學(xué)(北京), 2018.

    WANG Yan-hong. Molecular Dynamics Simulations of the Effect of Wettability on Explosive Boiling of a Nan-oscale Thin Liquid Film[D]. Beijing: North China Electric Power University, 2018.

    [35] 李大勇. 固液界面納米氣泡及其對(duì)流體邊界滑移影響的研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2014.

    LI Da-yong. Study of Nanobubbles at Solid-Liquid Inter-face and the Influence of Nanobubbles on Boundary Slip of Fluids[D]. Harbin: Harbin Institute of Technology, 2014.

    Slip Effect of Nanobubbles under the Influence of Hydrophilic and Hydrophobic Structures in Asymmetric Channels

    ,,,

    (a. Key Laboratory of Metallurgical Equipment and Control Technology of Ministry of Education, b. Key Laboratory of Mechanical Transmission and Manufacturing Engineering, Wuhan University of Science and Technology, Wuhan 430081, China)

    The purpose of this paper is to achieve a good fluid drag reduction effect by studying the relationship between the characteristics of nanobubbles and boundary slip under the influence of hydrophilic and hydrophobic surface structures in asymmetric channels. In this paper, the binary system molecular dynamics method is used to study the slip drag reduction effect of nanobubbles flowing in the channel. Firstly, the asymmetric microchannel model on the upper and lower walls is established, and the microchannel interface velocity slip phenomenon under the influence of nanobubbles is explored by considering the flow and heat transfer process of the microchannel. The simulation results show that when keeping the height of the hydrophilic lower wall and the temperature difference between the upper and lower walls unchanged, the increase in the height of the nanostructures on the upper wall promotes the volume increase of nanobubbles in the channel; in addition, when the upper wall is a hydrophobic wall, the nanobubbles appear in the form of a wall surface. As the volume increases, the slip length of the upper wall in the corresponding channel increases; when the upper wall is a hydrophilic wall, the nanobubbles are in bulk form, and as the volume increases, the corresponding upper wall’s slip length is reduced. In an asymmetrical channel, under the influence of the height of the upper wall structure, the increase of the nanobubble volume in the form of a wall surface promotes the drag reduction in the channel, while the increase of the nanobubble volume in the bulk form inhibits the drag reduction in the channel. Therefore, the research results in this paper provide a theoretical basis for exploring the application of nanobubble engineering in drag reduction.

    microchannel; temperature step; velocity slip; nanobubbles; slip drag reduction; wettability

    TH117

    A

    1001-3660(2022)02-0202-09

    10.16490/j.cnki.issn.1001-3660.2022.02.019

    2021-03-16;

    2021-07-05

    2021-03-16;

    2021-07-05

    國(guó)家自然科學(xué)基金(51875417,51975425)

    The National Natural Science Foundation of China (51875417, 51975425)

    戴雙武(1995—),男,碩士研究生,主要研究方向?yàn)槟Σ翆W(xué)理論。

    DAI Shuang-wu (1995—), Male, Postgraduate, Research focus: tribological theories.

    盧艷(1984—),女,博士,副教授,主要研究方向?yàn)楸砻婺Σ翆W(xué)理論。

    LU Yan (1984—), Female, Doctor, Associate professor, Research focus: surface tribology theory.

    戴雙武, 盧艷, 高友明, 等.非對(duì)稱(chēng)通道內(nèi)親疏水結(jié)構(gòu)影響下的納米氣泡滑移效應(yīng)[J]. 表面技術(shù), 2022, 51(2): 202-210.

    DAI Shuang-wu, LU Yan, GAO You-ming, et al. Slip Effect of Nanobubbles under the Influence of Hydrophilic and Hydrophobic Structures in Asymmetric Channels[J]. Surface Technology, 2022, 51(2): 202-210.

    猜你喜歡
    壁面氣泡流體
    檸檬氣泡水
    欣漾(2024年2期)2024-04-27 15:19:49
    二維有限長(zhǎng)度柔性壁面上T-S波演化的數(shù)值研究
    流體壓強(qiáng)知多少
    SIAU詩(shī)杭便攜式氣泡水杯
    新潮電子(2021年7期)2021-08-14 15:53:12
    浮法玻璃氣泡的預(yù)防和控制對(duì)策
    山雨欲來(lái)風(fēng)滿樓之流體壓強(qiáng)與流速
    冰凍氣泡
    等效流體體積模量直接反演的流體識(shí)別方法
    壁面溫度對(duì)微型內(nèi)燃機(jī)燃燒特性的影響
    顆粒—壁面碰撞建模與數(shù)據(jù)處理
    一二三四社区在线视频社区8| 国产一级毛片七仙女欲春2| 少妇的丰满在线观看| 亚洲av第一区精品v没综合| 丰满人妻一区二区三区视频av | 在线免费观看不下载黄p国产 | 久久精品国产清高在天天线| 亚洲 国产 在线| 一个人看的www免费观看视频| 亚洲av片天天在线观看| 九色国产91popny在线| 午夜福利成人在线免费观看| 日日干狠狠操夜夜爽| 国产91精品成人一区二区三区| 成年女人看的毛片在线观看| 亚洲无线观看免费| 午夜久久久久精精品| 亚洲人成电影免费在线| 叶爱在线成人免费视频播放| 午夜免费观看网址| 嫩草影视91久久| 小蜜桃在线观看免费完整版高清| 精品久久久久久久久久久久久| 亚洲一区高清亚洲精品| 女警被强在线播放| 精品国产三级普通话版| 一区二区三区国产精品乱码| 99久久精品国产亚洲精品| or卡值多少钱| 12—13女人毛片做爰片一| 久久热在线av| 国产精品久久视频播放| 亚洲精品美女久久久久99蜜臀| www日本黄色视频网| 欧美激情久久久久久爽电影| 91在线观看av| 一级毛片精品| xxx96com| 级片在线观看| 在线播放国产精品三级| 人妻久久中文字幕网| 久9热在线精品视频| 最近最新中文字幕大全免费视频| 国产蜜桃级精品一区二区三区| 欧美乱码精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 亚洲精华国产精华精| 亚洲精品粉嫩美女一区| 午夜福利18| 久久香蕉精品热| 欧美日本视频| 观看免费一级毛片| 精品午夜福利视频在线观看一区| 欧美成狂野欧美在线观看| 观看美女的网站| 国产精品九九99| avwww免费| 男女午夜视频在线观看| 国产毛片a区久久久久| 夜夜躁狠狠躁天天躁| 人妻久久中文字幕网| 国产日本99.免费观看| 国产真实乱freesex| 国产亚洲精品久久久com| 久久久久久大精品| 久久天堂一区二区三区四区| 久久久久国产精品人妻aⅴ院| 国产激情偷乱视频一区二区| 高清在线国产一区| 亚洲精品中文字幕一二三四区| 久久99热这里只有精品18| 一区二区三区激情视频| 欧美日韩国产亚洲二区| 国产成人精品久久二区二区免费| 国产av麻豆久久久久久久| 欧美大码av| 精品久久久久久久人妻蜜臀av| 最新在线观看一区二区三区| 在线观看日韩欧美| 久久午夜亚洲精品久久| 嫁个100分男人电影在线观看| 一进一出好大好爽视频| 夜夜看夜夜爽夜夜摸| 国产亚洲精品综合一区在线观看| 精品久久久久久成人av| 国产精品 国内视频| 哪里可以看免费的av片| 久久精品aⅴ一区二区三区四区| 国产精品乱码一区二三区的特点| 欧美日本视频| 亚洲 欧美一区二区三区| 国产又黄又爽又无遮挡在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产中文字幕在线视频| 欧美一级毛片孕妇| 精品无人区乱码1区二区| 久9热在线精品视频| 在线观看午夜福利视频| 给我免费播放毛片高清在线观看| 亚洲av熟女| 综合色av麻豆| 亚洲精品美女久久av网站| 国产高清激情床上av| 亚洲精品一区av在线观看| 午夜两性在线视频| 岛国在线观看网站| 中文字幕高清在线视频| 黑人欧美特级aaaaaa片| 欧美极品一区二区三区四区| 亚洲欧美激情综合另类| 精品欧美国产一区二区三| 男人和女人高潮做爰伦理| 免费搜索国产男女视频| 日本免费a在线| 久久久精品欧美日韩精品| 亚洲精品美女久久av网站| 哪里可以看免费的av片| 九九久久精品国产亚洲av麻豆 | 在线观看66精品国产| 亚洲最大成人中文| 白带黄色成豆腐渣| 九九在线视频观看精品| 久久国产精品影院| 99精品欧美一区二区三区四区| 国产精品久久久久久亚洲av鲁大| 欧美日韩黄片免| 中文亚洲av片在线观看爽| 99久久99久久久精品蜜桃| 欧美高清成人免费视频www| 国产av麻豆久久久久久久| 小蜜桃在线观看免费完整版高清| 9191精品国产免费久久| 久久久精品大字幕| 国内精品久久久久精免费| 日本一本二区三区精品| 人妻丰满熟妇av一区二区三区| 99国产精品99久久久久| 中文字幕精品亚洲无线码一区| 国产亚洲欧美98| 99久久无色码亚洲精品果冻| 九色成人免费人妻av| 中文字幕熟女人妻在线| 美女午夜性视频免费| 人人妻,人人澡人人爽秒播| 欧美日韩黄片免| 韩国av一区二区三区四区| 99热精品在线国产| 亚洲人成网站高清观看| 男女那种视频在线观看| 啦啦啦观看免费观看视频高清| 久久国产精品人妻蜜桃| 国产精品野战在线观看| 麻豆一二三区av精品| 亚洲一区二区三区不卡视频| 天堂av国产一区二区熟女人妻| 国产1区2区3区精品| 一个人看视频在线观看www免费 | 国产成人精品久久二区二区91| 长腿黑丝高跟| 亚洲熟妇中文字幕五十中出| 国产91精品成人一区二区三区| 久久婷婷人人爽人人干人人爱| 国产精品香港三级国产av潘金莲| 午夜福利欧美成人| 在线观看一区二区三区| 亚洲片人在线观看| av欧美777| 久久亚洲精品不卡| 久久亚洲精品不卡| 色综合亚洲欧美另类图片| 国产午夜精品久久久久久| 淫秽高清视频在线观看| 国产精品自产拍在线观看55亚洲| 午夜亚洲福利在线播放| 亚洲九九香蕉| 在线观看美女被高潮喷水网站 | 国产黄色小视频在线观看| 天天添夜夜摸| 亚洲精品久久国产高清桃花| 久久午夜亚洲精品久久| 特大巨黑吊av在线直播| 欧美成人免费av一区二区三区| 又爽又黄无遮挡网站| 亚洲人成伊人成综合网2020| 天天躁日日操中文字幕| 18禁裸乳无遮挡免费网站照片| 香蕉国产在线看| 真人一进一出gif抽搐免费| 国产成人一区二区三区免费视频网站| 成人鲁丝片一二三区免费| 亚洲欧洲精品一区二区精品久久久| 观看免费一级毛片| 日韩欧美在线乱码| 变态另类成人亚洲欧美熟女| ponron亚洲| 亚洲片人在线观看| 亚洲av电影在线进入| 久久婷婷人人爽人人干人人爱| 久久久久久大精品| 99精品欧美一区二区三区四区| 一个人免费在线观看的高清视频| 黄片小视频在线播放| 中文字幕人成人乱码亚洲影| 国产一区二区在线av高清观看| 国产蜜桃级精品一区二区三区| 欧美黑人欧美精品刺激| 国产不卡一卡二| 19禁男女啪啪无遮挡网站| 午夜亚洲福利在线播放| 美女cb高潮喷水在线观看 | 一二三四在线观看免费中文在| 12—13女人毛片做爰片一| 一a级毛片在线观看| 国产成年人精品一区二区| 国产一区二区三区在线臀色熟女| 亚洲国产欧美一区二区综合| 亚洲精品456在线播放app | 亚洲人成伊人成综合网2020| 国产免费男女视频| 日韩国内少妇激情av| 国产精品1区2区在线观看.| 老司机在亚洲福利影院| 精品国产超薄肉色丝袜足j| 大型黄色视频在线免费观看| 欧美日本亚洲视频在线播放| av欧美777| 亚洲熟妇中文字幕五十中出| 一个人免费在线观看的高清视频| 狠狠狠狠99中文字幕| 一边摸一边抽搐一进一小说| 亚洲人与动物交配视频| 欧洲精品卡2卡3卡4卡5卡区| 99精品在免费线老司机午夜| 免费无遮挡裸体视频| 国产精品综合久久久久久久免费| 午夜影院日韩av| 精品乱码久久久久久99久播| АⅤ资源中文在线天堂| 久久中文字幕人妻熟女| 精品一区二区三区视频在线 | 麻豆成人午夜福利视频| 法律面前人人平等表现在哪些方面| 香蕉久久夜色| 国产精品久久久久久久电影 | 国产极品精品免费视频能看的| 人人妻,人人澡人人爽秒播| 99热只有精品国产| 最近最新中文字幕大全免费视频| 午夜免费激情av| 国产亚洲欧美98| 成人欧美大片| 亚洲人成网站在线播放欧美日韩| 久久久国产欧美日韩av| 天堂动漫精品| 精品久久久久久久末码| e午夜精品久久久久久久| 国产在线精品亚洲第一网站| 国产精品乱码一区二三区的特点| 日本五十路高清| 全区人妻精品视频| 婷婷精品国产亚洲av在线| 岛国在线观看网站| 曰老女人黄片| av黄色大香蕉| 亚洲国产精品久久男人天堂| 18美女黄网站色大片免费观看| 国产精品久久久久久人妻精品电影| 啦啦啦免费观看视频1| 美女高潮喷水抽搐中文字幕| 国产午夜精品久久久久久| 免费av毛片视频| 国产一级毛片七仙女欲春2| 手机成人av网站| 香蕉久久夜色| 一本综合久久免费| 亚洲av片天天在线观看| 精品久久久久久久人妻蜜臀av| 99热这里只有是精品50| 久久午夜综合久久蜜桃| 国产成人啪精品午夜网站| 亚洲av电影不卡..在线观看| 久久久久国产一级毛片高清牌| 美女 人体艺术 gogo| 欧美中文日本在线观看视频| 国产av一区在线观看免费| 精品久久久久久久人妻蜜臀av| 老司机在亚洲福利影院| 久久人人精品亚洲av| 成年女人永久免费观看视频| 啪啪无遮挡十八禁网站| 一进一出抽搐动态| 久久伊人香网站| 小蜜桃在线观看免费完整版高清| 国产乱人伦免费视频| 国产乱人视频| 丰满的人妻完整版| 精品国产乱码久久久久久男人| 两性夫妻黄色片| 一卡2卡三卡四卡精品乱码亚洲| 99久久精品一区二区三区| 亚洲中文字幕日韩| 亚洲国产看品久久| 天天添夜夜摸| 午夜精品久久久久久毛片777| 久久人妻av系列| 免费av不卡在线播放| 久久久成人免费电影| 狠狠狠狠99中文字幕| 久久这里只有精品19| 高潮久久久久久久久久久不卡| 亚洲人成网站在线播放欧美日韩| 久久这里只有精品中国| 国产黄色小视频在线观看| 美女高潮的动态| 国产亚洲精品av在线| 国产三级中文精品| 精品久久久久久久毛片微露脸| 日韩有码中文字幕| 啦啦啦观看免费观看视频高清| 特级一级黄色大片| 日韩欧美国产一区二区入口| 国产精品日韩av在线免费观看| 精品欧美国产一区二区三| 国产又黄又爽又无遮挡在线| 婷婷精品国产亚洲av| 五月伊人婷婷丁香| 日韩欧美免费精品| 搞女人的毛片| 一区二区三区高清视频在线| 日韩三级视频一区二区三区| 丝袜人妻中文字幕| 在线播放国产精品三级| 国产精品1区2区在线观看.| 亚洲片人在线观看| 免费人成视频x8x8入口观看| 国产久久久一区二区三区| 久久精品国产综合久久久| 一个人看的www免费观看视频| 欧美+亚洲+日韩+国产| avwww免费| 搡老妇女老女人老熟妇| 禁无遮挡网站| 特大巨黑吊av在线直播| 国内精品久久久久精免费| 国产精品99久久99久久久不卡| 身体一侧抽搐| 亚洲美女视频黄频| 成人特级av手机在线观看| 桃红色精品国产亚洲av| 国产高清videossex| 小说图片视频综合网站| 国产蜜桃级精品一区二区三区| 亚洲精华国产精华精| 午夜成年电影在线免费观看| 国产精品一区二区三区四区免费观看 | 国产伦在线观看视频一区| 亚洲成人中文字幕在线播放| 99riav亚洲国产免费| 精品国产乱子伦一区二区三区| 亚洲真实伦在线观看| 国产精品av久久久久免费| aaaaa片日本免费| 精品久久久久久久人妻蜜臀av| 三级男女做爰猛烈吃奶摸视频| 国产不卡一卡二| 18禁美女被吸乳视频| 亚洲一区高清亚洲精品| 在线国产一区二区在线| 999精品在线视频| 黄片小视频在线播放| 这个男人来自地球电影免费观看| 法律面前人人平等表现在哪些方面| 欧美日韩瑟瑟在线播放| 成人三级黄色视频| 日本a在线网址| 欧美日韩亚洲国产一区二区在线观看| 欧美绝顶高潮抽搐喷水| 五月伊人婷婷丁香| 亚洲欧美一区二区三区黑人| 美女扒开内裤让男人捅视频| 一个人看的www免费观看视频| 久久天堂一区二区三区四区| 给我免费播放毛片高清在线观看| 亚洲天堂国产精品一区在线| 亚洲欧美精品综合一区二区三区| 在线免费观看的www视频| 午夜激情欧美在线| 一二三四在线观看免费中文在| 亚洲精品粉嫩美女一区| 日本黄色片子视频| 免费在线观看视频国产中文字幕亚洲| 色噜噜av男人的天堂激情| 色精品久久人妻99蜜桃| 日韩欧美国产在线观看| 国产极品精品免费视频能看的| 国产成人福利小说| 亚洲一区高清亚洲精品| 看黄色毛片网站| 12—13女人毛片做爰片一| a级毛片a级免费在线| 亚洲欧美激情综合另类| 亚洲无线观看免费| 国产一区二区在线观看日韩 | 国产精品98久久久久久宅男小说| 国产免费av片在线观看野外av| 一本久久中文字幕| 天堂av国产一区二区熟女人妻| 欧美又色又爽又黄视频| 18禁美女被吸乳视频| av黄色大香蕉| 一进一出抽搐gif免费好疼| 久久久久久大精品| 97超级碰碰碰精品色视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 午夜a级毛片| 毛片女人毛片| 超碰成人久久| 男女午夜视频在线观看| 国产精品日韩av在线免费观看| 亚洲欧美日韩高清专用| 两个人看的免费小视频| 丁香欧美五月| 叶爱在线成人免费视频播放| 欧美黄色淫秽网站| 亚洲国产日韩欧美精品在线观看 | 成人国产综合亚洲| 男女做爰动态图高潮gif福利片| 亚洲aⅴ乱码一区二区在线播放| 国产69精品久久久久777片 | 亚洲狠狠婷婷综合久久图片| 亚洲av电影不卡..在线观看| 国产亚洲精品久久久久久毛片| 天堂√8在线中文| e午夜精品久久久久久久| 亚洲av中文字字幕乱码综合| 91av网一区二区| 男人和女人高潮做爰伦理| 91在线观看av| 午夜福利欧美成人| tocl精华| 免费一级毛片在线播放高清视频| 午夜免费成人在线视频| 欧美成人性av电影在线观看| 日本在线视频免费播放| 制服人妻中文乱码| 网址你懂的国产日韩在线| 国产三级在线视频| 亚洲国产精品成人综合色| 精品久久久久久久久久免费视频| 久久精品人妻少妇| av视频在线观看入口| 变态另类丝袜制服| 欧美3d第一页| 亚洲专区国产一区二区| 亚洲av成人av| 久久久久国产一级毛片高清牌| 在线观看66精品国产| 噜噜噜噜噜久久久久久91| 亚洲avbb在线观看| 免费在线观看影片大全网站| а√天堂www在线а√下载| 亚洲欧美日韩高清在线视频| 午夜福利高清视频| 亚洲熟妇熟女久久| 精品久久久久久成人av| 一进一出抽搐gif免费好疼| 久久久久久久精品吃奶| 精品一区二区三区av网在线观看| 国产野战对白在线观看| 性色avwww在线观看| 波多野结衣巨乳人妻| 99久久精品热视频| 少妇的丰满在线观看| 我的老师免费观看完整版| 亚洲中文av在线| 精品久久蜜臀av无| 99视频精品全部免费 在线 | 精品日产1卡2卡| 免费在线观看影片大全网站| 国产高清有码在线观看视频| 99热精品在线国产| 日韩成人在线观看一区二区三区| 欧美成人一区二区免费高清观看 | 亚洲国产看品久久| 熟女人妻精品中文字幕| 欧美不卡视频在线免费观看| 国产黄a三级三级三级人| 一卡2卡三卡四卡精品乱码亚洲| 两个人的视频大全免费| 亚洲av熟女| 波多野结衣巨乳人妻| 国产免费av片在线观看野外av| 精品一区二区三区四区五区乱码| 变态另类成人亚洲欧美熟女| 露出奶头的视频| 国产伦精品一区二区三区四那| 色在线成人网| 丝袜人妻中文字幕| 亚洲自偷自拍图片 自拍| 色综合亚洲欧美另类图片| 精品一区二区三区视频在线 | 精品一区二区三区av网在线观看| 久久欧美精品欧美久久欧美| 日本免费a在线| 亚洲成av人片在线播放无| 啦啦啦免费观看视频1| 黄片大片在线免费观看| 国产伦一二天堂av在线观看| 久久99热这里只有精品18| 女警被强在线播放| 国产精品一区二区三区四区免费观看 | 色精品久久人妻99蜜桃| 又黄又爽又免费观看的视频| 国内毛片毛片毛片毛片毛片| www国产在线视频色| 免费看光身美女| 国产免费av片在线观看野外av| 麻豆av在线久日| 国产在线精品亚洲第一网站| 不卡一级毛片| 久久久久久久午夜电影| 两人在一起打扑克的视频| 一本综合久久免费| 亚洲国产精品sss在线观看| 天天躁日日操中文字幕| 国产精品永久免费网站| 一区二区三区激情视频| 国产av不卡久久| 国产69精品久久久久777片 | 亚洲精品一区av在线观看| 久久久精品欧美日韩精品| 美女高潮的动态| 国产午夜精品论理片| 亚洲一区二区三区色噜噜| 国产av不卡久久| 99精品久久久久人妻精品| 啦啦啦观看免费观看视频高清| 不卡一级毛片| 九色成人免费人妻av| bbb黄色大片| 香蕉久久夜色| 亚洲精品国产精品久久久不卡| 国内久久婷婷六月综合欲色啪| 久久久国产欧美日韩av| 一级作爱视频免费观看| 一级黄色大片毛片| 性色av乱码一区二区三区2| 波多野结衣巨乳人妻| 日本精品一区二区三区蜜桃| 欧美中文综合在线视频| 日本与韩国留学比较| 亚洲国产色片| 国产精品亚洲一级av第二区| 夜夜爽天天搞| 特级一级黄色大片| 一个人免费在线观看的高清视频| 啦啦啦观看免费观看视频高清| 在线永久观看黄色视频| 国产单亲对白刺激| 日韩欧美 国产精品| 日本一二三区视频观看| 欧美日韩综合久久久久久 | 成人特级黄色片久久久久久久| 欧美激情在线99| 香蕉久久夜色| 亚洲av日韩精品久久久久久密| 丁香六月欧美| 天堂动漫精品| 国产精品 欧美亚洲| 五月伊人婷婷丁香| 免费在线观看成人毛片| 一区二区三区高清视频在线| 高清毛片免费观看视频网站| 国产精品98久久久久久宅男小说| 狠狠狠狠99中文字幕| 香蕉久久夜色| 舔av片在线| 小蜜桃在线观看免费完整版高清| 一二三四社区在线视频社区8| 99久久成人亚洲精品观看| 国产精品香港三级国产av潘金莲| 最新在线观看一区二区三区| 此物有八面人人有两片| 熟女人妻精品中文字幕| 成人av一区二区三区在线看| 黄片大片在线免费观看| 女同久久另类99精品国产91| 天天一区二区日本电影三级| 欧美午夜高清在线| 国产 一区 欧美 日韩| 国产不卡一卡二| 色视频www国产| 午夜福利高清视频| 精品一区二区三区视频在线 | 最近最新中文字幕大全电影3| 亚洲国产看品久久| 夜夜夜夜夜久久久久| 国产精品久久久久久亚洲av鲁大| tocl精华| 亚洲精品在线美女| 后天国语完整版免费观看| 亚洲成人免费电影在线观看| 欧美黑人巨大hd| 岛国在线观看网站| 日日干狠狠操夜夜爽| 亚洲专区字幕在线| 久久精品综合一区二区三区| 成人欧美大片| 欧美日韩中文字幕国产精品一区二区三区| 真人做人爱边吃奶动态| 亚洲中文av在线| 久久久久久久久中文| 十八禁网站免费在线| 人妻丰满熟妇av一区二区三区| www.999成人在线观看|