• <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ù)處理
    成年人午夜在线观看视频| 男的添女的下面高潮视频| 狂野欧美激情性bbbbbb| 日韩熟女老妇一区二区性免费视频| 日韩一区二区三区影片| 国内精品宾馆在线| 国内精品宾馆在线| 在现免费观看毛片| 日产精品乱码卡一卡2卡三| 狂野欧美激情性bbbbbb| 日韩中文字幕视频在线看片| 免费高清在线观看视频在线观看| 国产免费又黄又爽又色| 狂野欧美白嫩少妇大欣赏| 啦啦啦中文免费视频观看日本| 26uuu在线亚洲综合色| 成人毛片a级毛片在线播放| 日本免费在线观看一区| 久久女婷五月综合色啪小说| 老熟女久久久| av又黄又爽大尺度在线免费看| 欧美人与善性xxx| 中文字幕免费在线视频6| 韩国av在线不卡| 国产精品一区www在线观看| 黑人高潮一二区| 日韩欧美一区视频在线观看| 国产片特级美女逼逼视频| 欧美精品国产亚洲| 制服诱惑二区| 国产黄频视频在线观看| 在线看a的网站| 亚洲欧美色中文字幕在线| 亚洲欧美日韩另类电影网站| 国产成人免费无遮挡视频| 日日撸夜夜添| 日韩不卡一区二区三区视频在线| 久久99精品国语久久久| 日韩av在线免费看完整版不卡| 成年av动漫网址| 十分钟在线观看高清视频www| 成人毛片60女人毛片免费| 国产av精品麻豆| 麻豆成人av视频| 国产日韩欧美在线精品| 下体分泌物呈黄色| 亚洲色图综合在线观看| 亚洲精品久久久久久婷婷小说| 午夜福利视频在线观看免费| 久久国产亚洲av麻豆专区| av一本久久久久| 免费观看性生交大片5| 少妇人妻 视频| 黄色配什么色好看| 亚洲av成人精品一区久久| 午夜老司机福利剧场| 亚洲av福利一区| 国产黄频视频在线观看| 精品人妻熟女毛片av久久网站| 丝袜在线中文字幕| 老司机影院成人| 日韩成人av中文字幕在线观看| 日韩大片免费观看网站| 欧美日韩成人在线一区二区| 精品国产一区二区三区久久久樱花| 久久久精品免费免费高清| 自线自在国产av| 午夜激情av网站| 中国国产av一级| 3wmmmm亚洲av在线观看| 色婷婷av一区二区三区视频| 在线观看三级黄色| 国产成人精品婷婷| 午夜激情久久久久久久| 欧美日韩综合久久久久久| 国产一级毛片在线| 精品久久久精品久久久| 一级毛片我不卡| 春色校园在线视频观看| 久久青草综合色| 亚洲欧洲国产日韩| 亚洲欧美成人综合另类久久久| 国产一区二区三区av在线| 黄片无遮挡物在线观看| 成人免费观看视频高清| 久久久久国产精品人妻一区二区| 一级片'在线观看视频| 26uuu在线亚洲综合色| 2022亚洲国产成人精品| 久久99精品国语久久久| 久久av网站| 亚洲五月色婷婷综合| 国语对白做爰xxxⅹ性视频网站| 黄片播放在线免费| 精品一区在线观看国产| 久久人人爽人人爽人人片va| 十分钟在线观看高清视频www| 91午夜精品亚洲一区二区三区| 三级国产精品欧美在线观看| 夜夜看夜夜爽夜夜摸| 美女中出高潮动态图| 久久久久久人妻| 欧美亚洲日本最大视频资源| 校园人妻丝袜中文字幕| 亚洲精品国产色婷婷电影| 国产亚洲欧美精品永久| 色吧在线观看| 日本与韩国留学比较| 亚洲精华国产精华液的使用体验| 一区二区三区乱码不卡18| 久久久精品94久久精品| 国产亚洲最大av| 精品久久久久久久久av| 亚洲第一区二区三区不卡| 日韩精品免费视频一区二区三区 | 在线观看一区二区三区激情| 蜜桃国产av成人99| 嫩草影院入口| 精品一区二区免费观看| 亚洲第一区二区三区不卡| 久久亚洲国产成人精品v| 18禁观看日本| 欧美xxxx性猛交bbbb| 婷婷色综合www| 在线观看免费高清a一片| 久久综合国产亚洲精品| 少妇的逼好多水| 中国三级夫妇交换| 国产一区二区在线观看日韩| 成人手机av| 亚洲综合色惰| 一级,二级,三级黄色视频| 狠狠婷婷综合久久久久久88av| 亚洲中文av在线| 亚洲性久久影院| 一级毛片aaaaaa免费看小| 丰满少妇做爰视频| 久久99热这里只频精品6学生| 国产高清不卡午夜福利| 午夜福利视频在线观看免费| 成年人午夜在线观看视频| 曰老女人黄片| 久久国产亚洲av麻豆专区| 国产成人免费观看mmmm| 亚洲精品乱码久久久久久按摩| 街头女战士在线观看网站| 日本vs欧美在线观看视频| 美女国产视频在线观看| 99久久中文字幕三级久久日本| av福利片在线| 亚洲欧美成人综合另类久久久| 18禁动态无遮挡网站| 亚洲欧洲精品一区二区精品久久久 | 国产又色又爽无遮挡免| 亚洲精华国产精华液的使用体验| 熟女电影av网| 日韩av不卡免费在线播放| 久久久久精品久久久久真实原创| 欧美另类一区| 男的添女的下面高潮视频| 国产成人一区二区在线| 精品少妇黑人巨大在线播放| 国产午夜精品一二区理论片| 久久精品国产亚洲av天美| 国产色婷婷99| 国产精品不卡视频一区二区| 久久久国产精品麻豆| 亚洲成人av在线免费| 国产毛片在线视频| 国产免费现黄频在线看| 一级,二级,三级黄色视频| 满18在线观看网站| 99精国产麻豆久久婷婷| 久久免费观看电影| 成人国产av品久久久| 国产国语露脸激情在线看| 亚洲av日韩在线播放| 日韩一本色道免费dvd| 亚洲怡红院男人天堂| 国产亚洲午夜精品一区二区久久| 亚州av有码| √禁漫天堂资源中文www| 超色免费av| 天天操日日干夜夜撸| 天天操日日干夜夜撸| 精品人妻偷拍中文字幕| 这个男人来自地球电影免费观看 | 中文字幕制服av| 麻豆乱淫一区二区| 乱人伦中国视频| 国产毛片在线视频| 久久狼人影院| 日本91视频免费播放| 国产乱人偷精品视频| 欧美成人午夜免费资源| 内地一区二区视频在线| 久久久久久久国产电影| 99热这里只有是精品在线观看| 日本黄色日本黄色录像| 亚洲怡红院男人天堂| 母亲3免费完整高清在线观看 | 大陆偷拍与自拍| 免费观看在线日韩| 高清黄色对白视频在线免费看| 2022亚洲国产成人精品| 国产熟女午夜一区二区三区 | 亚洲欧美中文字幕日韩二区| 国产精品一二三区在线看| 在线观看www视频免费| 国产亚洲一区二区精品| 一级黄片播放器| 三级国产精品片| 人妻系列 视频| 99热全是精品| 欧美丝袜亚洲另类| 久热久热在线精品观看| 中文字幕制服av| 亚洲精品乱码久久久v下载方式| 亚洲av在线观看美女高潮| 国产亚洲精品第一综合不卡 | 香蕉精品网在线| 国产精品久久久久久久久免| 午夜精品国产一区二区电影| 老熟女久久久| 精品亚洲成国产av| 少妇精品久久久久久久| 国产免费又黄又爽又色| 久久婷婷青草| 免费高清在线观看日韩| 亚洲色图 男人天堂 中文字幕 | 女性生殖器流出的白浆| 中文字幕人妻熟人妻熟丝袜美| 日韩视频在线欧美| 在线播放无遮挡| 国产一区二区三区综合在线观看 | 欧美日本中文国产一区发布| a级毛片在线看网站| 2018国产大陆天天弄谢| 能在线免费看毛片的网站| 亚洲少妇的诱惑av| 国产成人精品在线电影| 亚洲欧洲国产日韩| 国产av国产精品国产| 纯流量卡能插随身wifi吗| 赤兔流量卡办理| 午夜福利网站1000一区二区三区| 日日撸夜夜添| 啦啦啦中文免费视频观看日本| 麻豆乱淫一区二区| 国产精品不卡视频一区二区| 国产精品 国内视频| 九色成人免费人妻av| 99九九在线精品视频| 免费黄网站久久成人精品| 中国美白少妇内射xxxbb| 永久网站在线| 亚洲美女视频黄频| 国产亚洲欧美精品永久| 亚洲国产精品一区三区| 好男人视频免费观看在线| 中文乱码字字幕精品一区二区三区| 麻豆乱淫一区二区| 纯流量卡能插随身wifi吗| 又黄又爽又刺激的免费视频.| 大又大粗又爽又黄少妇毛片口| 日日摸夜夜添夜夜添av毛片| 亚洲经典国产精华液单| 永久网站在线| 国产熟女午夜一区二区三区 | 亚洲成人一二三区av| 久久精品国产亚洲av天美| 久久久久国产精品人妻一区二区| 纯流量卡能插随身wifi吗| 亚洲国产精品999| 国产精品无大码| 天美传媒精品一区二区| 久久99热6这里只有精品| 国产精品不卡视频一区二区| 精品国产一区二区三区久久久樱花| 一区二区三区精品91| 制服人妻中文乱码| 麻豆精品久久久久久蜜桃| 人妻一区二区av| 91国产中文字幕| 久久精品久久精品一区二区三区| 亚洲国产成人一精品久久久| 伦理电影大哥的女人| 国产极品粉嫩免费观看在线 | 国产爽快片一区二区三区| 国产高清不卡午夜福利| av又黄又爽大尺度在线免费看| 国产片特级美女逼逼视频| 日韩一区二区视频免费看| 亚洲av国产av综合av卡| 精品午夜福利在线看| 久久久久精品久久久久真实原创| av一本久久久久| 婷婷成人精品国产| 黑丝袜美女国产一区| 777米奇影视久久| 天堂俺去俺来也www色官网| 亚洲成色77777| 一级片'在线观看视频| 精品人妻熟女av久视频| 在线精品无人区一区二区三| 又大又黄又爽视频免费| 99久久精品国产国产毛片| av视频免费观看在线观看| 精品久久久久久久久亚洲| 国语对白做爰xxxⅹ性视频网站| 少妇的逼水好多| 欧美日韩精品成人综合77777| 大香蕉97超碰在线| 高清视频免费观看一区二区| 视频区图区小说| 国产精品熟女久久久久浪| 91精品一卡2卡3卡4卡| 国产黄色免费在线视频| 久久午夜综合久久蜜桃| 飞空精品影院首页| 日韩熟女老妇一区二区性免费视频| 亚洲av男天堂| 丝袜脚勾引网站| 亚洲国产精品999| 又大又黄又爽视频免费| 一本大道久久a久久精品| 成年女人在线观看亚洲视频| 国产黄频视频在线观看| 九九在线视频观看精品| 国产一区二区三区av在线| 国产精品久久久久久av不卡| 国产免费又黄又爽又色| av在线观看视频网站免费| 天美传媒精品一区二区| 亚洲欧美日韩另类电影网站| 国产在线视频一区二区| 91aial.com中文字幕在线观看| 久久久精品免费免费高清| 街头女战士在线观看网站| 日本爱情动作片www.在线观看| 免费黄频网站在线观看国产| 汤姆久久久久久久影院中文字幕| 91久久精品国产一区二区成人| 老熟女久久久| 亚洲第一区二区三区不卡| 国产亚洲最大av| 晚上一个人看的免费电影| 亚洲国产日韩一区二区| 日本vs欧美在线观看视频| 水蜜桃什么品种好| 国产探花极品一区二区| 日韩三级伦理在线观看| 99热6这里只有精品| 性高湖久久久久久久久免费观看| 亚洲精品乱久久久久久| 一级黄片播放器| 性高湖久久久久久久久免费观看| 亚洲国产欧美在线一区| 性色av一级| 99精国产麻豆久久婷婷| av一本久久久久| 蜜桃久久精品国产亚洲av| 伊人亚洲综合成人网| 中文字幕av电影在线播放| 午夜福利视频在线观看免费| 亚洲av国产av综合av卡| 日韩免费高清中文字幕av| 日韩一区二区视频免费看| 成人手机av| 青春草国产在线视频| a级毛片在线看网站| 日韩亚洲欧美综合| 成人毛片60女人毛片免费| 免费日韩欧美在线观看| 满18在线观看网站| 精品一区二区三区视频在线| 国产女主播在线喷水免费视频网站| 免费日韩欧美在线观看| 国产成人精品无人区| 2021少妇久久久久久久久久久| 九九在线视频观看精品| 热99国产精品久久久久久7| 欧美日韩综合久久久久久| 国模一区二区三区四区视频| 熟女电影av网| 亚洲av免费高清在线观看| 夜夜看夜夜爽夜夜摸| 免费观看的影片在线观看| 国产av码专区亚洲av| 国模一区二区三区四区视频| 国产精品久久久久久久电影| 春色校园在线视频观看| 国产精品久久久久久久电影| 热99国产精品久久久久久7| 国产男女超爽视频在线观看| 成人国产av品久久久| 少妇人妻久久综合中文| 亚洲,一卡二卡三卡| 精品视频人人做人人爽| 观看av在线不卡| 亚洲综合精品二区| 久久婷婷青草| 国产伦精品一区二区三区视频9| 中文字幕精品免费在线观看视频 | 人妻 亚洲 视频| 日韩不卡一区二区三区视频在线| 国产精品久久久久久精品电影小说| 大陆偷拍与自拍| 久久97久久精品| a级片在线免费高清观看视频| 国产成人精品婷婷| 亚洲精品乱码久久久v下载方式| 亚洲精品久久久久久婷婷小说| 能在线免费看毛片的网站| 日日摸夜夜添夜夜添av毛片| 亚洲av福利一区| 欧美 日韩 精品 国产| av网站免费在线观看视频| 亚洲国产精品一区三区| 国产爽快片一区二区三区| 久久免费观看电影| 国产免费一区二区三区四区乱码| 日韩成人伦理影院| 精品久久蜜臀av无| 亚洲国产成人一精品久久久| 日韩人妻高清精品专区| 性色avwww在线观看| 亚洲,欧美,日韩| 只有这里有精品99| 99久国产av精品国产电影| 久久久久久久久久久免费av| 搡女人真爽免费视频火全软件| 成人无遮挡网站| 午夜福利,免费看| 又粗又硬又长又爽又黄的视频| 内地一区二区视频在线| 亚洲精品日韩在线中文字幕| 亚洲欧美清纯卡通| 人妻制服诱惑在线中文字幕| 国产亚洲午夜精品一区二区久久| 中文字幕亚洲精品专区| av网站免费在线观看视频| 美女福利国产在线| 另类亚洲欧美激情| 日韩三级伦理在线观看| 嘟嘟电影网在线观看| 纵有疾风起免费观看全集完整版| 久久国产亚洲av麻豆专区| 91精品一卡2卡3卡4卡| 五月玫瑰六月丁香| 飞空精品影院首页| 日韩av在线免费看完整版不卡| 午夜91福利影院| 大又大粗又爽又黄少妇毛片口| 亚洲av成人精品一二三区| 人妻系列 视频| 三级国产精品片| 久久鲁丝午夜福利片| 日本vs欧美在线观看视频| 999精品在线视频| 日韩人妻高清精品专区| 嫩草影院入口| 国产成人一区二区在线| 国产成人精品婷婷| 国产亚洲一区二区精品| 最近中文字幕2019免费版| 日韩av免费高清视频| 久久久久久人妻| 国产极品天堂在线| videossex国产| 日本欧美视频一区| 国产极品粉嫩免费观看在线 | 在线观看一区二区三区激情| 观看av在线不卡| 免费日韩欧美在线观看| 九草在线视频观看| 男女边吃奶边做爰视频| 中文字幕免费在线视频6| 亚洲av国产av综合av卡| av电影中文网址| 亚洲无线观看免费| 丰满乱子伦码专区| 午夜视频国产福利| 亚洲精品色激情综合| 美女中出高潮动态图| 精品一区在线观看国产| 国产精品99久久久久久久久| 91久久精品国产一区二区成人| 色网站视频免费| 老司机亚洲免费影院| 亚洲精品乱码久久久久久按摩| 天堂俺去俺来也www色官网| 五月伊人婷婷丁香| 日日啪夜夜爽| 国产成人精品婷婷| 日韩成人伦理影院| 91精品伊人久久大香线蕉| 插逼视频在线观看| 最近的中文字幕免费完整| 久久99蜜桃精品久久| 九九爱精品视频在线观看| 国产一区亚洲一区在线观看| 在线播放无遮挡| 国产成人免费无遮挡视频| av.在线天堂| 亚洲精品第二区| 国产一区二区在线观看日韩| 亚洲欧美色中文字幕在线| 黑人猛操日本美女一级片| 欧美激情国产日韩精品一区| 午夜老司机福利剧场| av网站免费在线观看视频| 亚洲成人一二三区av| 一级毛片 在线播放| 秋霞伦理黄片| 久久97久久精品| 肉色欧美久久久久久久蜜桃| 丰满饥渴人妻一区二区三| 视频区图区小说| 精品人妻一区二区三区麻豆| av天堂久久9| videosex国产| 视频在线观看一区二区三区| 国产一级毛片在线| 各种免费的搞黄视频| 亚洲无线观看免费| 午夜91福利影院| 天天躁夜夜躁狠狠久久av| 久热久热在线精品观看| 中文精品一卡2卡3卡4更新| 亚洲无线观看免费| 最后的刺客免费高清国语| 久久久亚洲精品成人影院| 亚洲国产成人一精品久久久| 国产欧美另类精品又又久久亚洲欧美| 国产精品人妻久久久久久| 午夜福利视频精品| 美女脱内裤让男人舔精品视频| 一区二区三区乱码不卡18| 日本-黄色视频高清免费观看| 国产日韩欧美视频二区| 亚洲国产av影院在线观看| 99久久精品一区二区三区| 卡戴珊不雅视频在线播放| 国产女主播在线喷水免费视频网站| 婷婷成人精品国产| 精品久久蜜臀av无| 欧美日本中文国产一区发布| 99久国产av精品国产电影| 中文欧美无线码| 久久久a久久爽久久v久久| 国产综合精华液| 黄色一级大片看看| 丝袜脚勾引网站| 精品国产乱码久久久久久小说| 精品一品国产午夜福利视频| 亚洲欧美日韩卡通动漫| 永久免费av网站大全| 亚洲精品中文字幕在线视频| 午夜激情久久久久久久| 成年人免费黄色播放视频| 中文乱码字字幕精品一区二区三区| 亚洲天堂av无毛| 国产色爽女视频免费观看| 黄片播放在线免费| 亚洲第一区二区三区不卡| 亚洲成人一二三区av| 又粗又硬又长又爽又黄的视频| 国产在线视频一区二区| 中文乱码字字幕精品一区二区三区| 日本色播在线视频| 成人综合一区亚洲| 在线观看免费高清a一片| 日韩人妻高清精品专区| 高清毛片免费看| 久久影院123| 全区人妻精品视频| 韩国高清视频一区二区三区| 18禁裸乳无遮挡动漫免费视频| 不卡视频在线观看欧美| a级毛片黄视频| 国产一区二区在线观看日韩| 少妇的逼水好多| 欧美日韩国产mv在线观看视频| 欧美日韩成人在线一区二区| 久久久a久久爽久久v久久| 91在线精品国自产拍蜜月| 人妻夜夜爽99麻豆av| 狠狠精品人妻久久久久久综合| 男女高潮啪啪啪动态图| 色网站视频免费| 久久精品久久久久久久性| 亚洲精品色激情综合| 黄色一级大片看看| 狂野欧美白嫩少妇大欣赏| 精品一区二区免费观看| 欧美人与善性xxx| 女性被躁到高潮视频| 99热国产这里只有精品6| 考比视频在线观看| 最近手机中文字幕大全| 男人操女人黄网站| 一级毛片aaaaaa免费看小| av国产久精品久网站免费入址| 日产精品乱码卡一卡2卡三| tube8黄色片| 日本爱情动作片www.在线观看| 日韩av不卡免费在线播放| 免费黄频网站在线观看国产| 最近中文字幕高清免费大全6| 国产永久视频网站| 十八禁高潮呻吟视频| 亚洲情色 制服丝袜| 亚洲av欧美aⅴ国产| 精品亚洲乱码少妇综合久久|