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

    納米制冷劑在微通道中對流換熱特性及熵產(chǎn)分析

    2018-11-21 03:49:32郭文華彭浩
    制冷技術(shù) 2018年4期
    關(guān)鍵詞:雷諾數(shù)制冷劑對流

    郭文華,彭浩

    ?

    納米制冷劑在微通道中對流換熱特性及熵產(chǎn)分析

    郭文華,彭浩

    (上海海事大學(xué)商船學(xué)院,上海 201306)

    本文對R134a納米制冷劑在微通道中的對流換熱特性和熵產(chǎn)情況進(jìn)行了數(shù)值研究,其中納米顆粒包括氧化銅(CuO)和碳納米管(CNT)、體積分?jǐn)?shù)范圍為0~4%、雷諾數(shù)()范圍為400~2,000。結(jié)果表明,對流換熱系數(shù)和摩擦壓降隨著體積分?jǐn)?shù)和雷諾數(shù)的增加而增大;=1,200時,體積分?jǐn)?shù)為4%的R134a-CuO和R134a-CNT納米制冷劑的對流換熱系數(shù)比純R134a分別高出了45.7%、125.7%,摩擦壓降分別高出了192.4%、161.8%。熱熵產(chǎn)隨體積分?jǐn)?shù)和雷諾數(shù)的增加而減小,摩擦熵產(chǎn)隨體積分?jǐn)?shù)和雷諾數(shù)的增加而增大;R134a-CuO和R134a-CNT納米制冷劑的總熵產(chǎn)最小產(chǎn)生條件分別為體積分?jǐn)?shù)2%、=1,600和體積分?jǐn)?shù)4%、=800。

    納米制冷劑;對流換熱;熵產(chǎn);微通道;數(shù)值模擬

    0 引言

    微通道結(jié)構(gòu)是強(qiáng)化換熱的一種有效手段,廣泛應(yīng)用于制冷系統(tǒng)和電子芯片的冷卻裝置中[1-2]。制冷劑是其中的常用工質(zhì)[3]。在制冷劑中添加納米粒子形成納米制冷劑,研究表明其導(dǎo)熱和換熱性能顯著提高[4-5]。對于微通道中的制冷劑工質(zhì),對流換熱是基本的傳熱形式。納米粒子的種類、雷諾數(shù)和體積分?jǐn)?shù)會顯著影響納米制冷劑的對流換熱系數(shù)和摩擦壓降。為了選取合適的納米粒子以及確定使對流換熱過程的不可逆損失最小的工作條件,需要明確微通道中的納米制冷劑的對流換熱特性及進(jìn)行熵產(chǎn)分析。

    微通道中對流換熱及熵產(chǎn)特性的研究目前主要集中在水基納米流體[6-7],結(jié)果顯示,納米粒子的種類、體積分?jǐn)?shù)以及雷諾數(shù)等對熵產(chǎn)影響顯著。由于納米制冷劑的熱物性與水基納米流體差異明顯[8],需要對其在微通道中的對流換熱及熵產(chǎn)特性進(jìn)行專門的研究。

    本文以R134a-氧化銅(CuO)、R134a-碳納米管(CNT)兩種納米制冷劑為研究對象,擬建立納米制冷劑在矩形微通道中對流換熱的熵產(chǎn)計算模型;并明確納米粒子的種類、體積分?jǐn)?shù)和雷諾數(shù)等因素對微通道中對流換熱過程及熵產(chǎn)的影響規(guī)律。

    1 數(shù)值模擬

    1.1 模擬對象

    本文模擬對象是R134a-氧化銅(CuO)、R134a-碳納米管(CNT)兩種納米制冷劑三維矩形微通道中的對流換熱。選取的納米粒子的體積分?jǐn)?shù)為0、1%、2%、3%、4%,雷諾數(shù)范圍為400~2,000。微通道的長度為40 mm,寬度ch為382mm,深度ch為470mm,水力直徑h為421mm。微通道上表面絕熱,其余3個表面加恒熱流密度w= 3,530.01 W/m2。

    1.2 控制方程

    對物理模型進(jìn)行以下假設(shè)和簡化:1)納米制冷劑為單相不可壓縮流體;2)制冷劑和納米粒子處于平衡態(tài),且沒有速度滑移;3)考慮Boussinesq的假設(shè),即溫度變化會引起流體的密度變化。描述納米制冷劑對流換熱過程的連續(xù)性、動量和能量方程如下[9,12]。

    連續(xù)性方程:

    ?(nr) = 0 (1)

    動量方程:

    ?(nr) = -?+?(nr?) (2)

    能量方程:

    ?(nrVcT) =?(nr?) (3)

    式中:

    ——平均質(zhì)量流速,kg/(m2·s);

    ——壓力,Pa;

    nr——納米制冷劑的密度,kg/m3;

    nr——納米制冷劑的動力粘度,Pa·s;

    nr——納米制冷劑的熱導(dǎo)率,W/(m·K);

    c——納米制冷劑的比熱,J/(kg·K)。

    納米粒子CNT長為1.5mm,直徑為80 nm,CuO的直徑n為20 nm,R134a-CuO和R134a-CNT納米制冷劑的物性由公式(4)分別計算得出[10]。

    (5)

    (6)

    (8)

    式中:

    max——納米粒子影響形狀的最大填充因子,納米粒子CuO和CNT取值分別為0.632、0.268[11];

    ——固有粘度,CuO和CNT取值分別為2.5和9.25;

    ——CNT的熱導(dǎo)率參數(shù);

    ——球度參數(shù);

    ——納米制冷劑的體積分?jǐn)?shù);

    下角標(biāo)r——代表基液R134a;

    下角標(biāo)n——代表納米粒子。

    1.3 邊界條件

    模擬中,邊界條件使用速度入口,出口自由出流,底部及兩側(cè)加恒定熱流,上部絕熱。納米制冷劑進(jìn)入微通道的入口溫度為in為300 K,操作壓力為1,013,250 Pa。微通道內(nèi)的流動為層流,考慮重力場的影響,取x=0,y=-9.8 m/s2。

    1.4 數(shù)值方法

    模擬選用了穩(wěn)態(tài)條件,用有限體積法求解具有適當(dāng)邊界條件的控制方程。壓力場和速度場使用SIMPLE耦合,壓力修正方程采用PRESTO處理,能量和動量方程使用二階迎風(fēng)格式,壓力關(guān)聯(lián)方程、動量方程、能量方程的松弛因子分別為0.3、0.8、0.6;連續(xù)性方程、動量方程、能量方程的收斂殘差皆為10-6。

    1.5 網(wǎng)格劃分及無關(guān)性分析

    劃分4種不同的網(wǎng)格數(shù)(包括15,000、60,000、300,000和840,000)來進(jìn)行網(wǎng)格無關(guān)性分析。當(dāng)?shù)撞考雍愣崃鱳為3,530.01 W/m2,雷諾數(shù)為4,500,進(jìn)口溫度in為300 K,實(shí)驗(yàn)結(jié)果得R134a在矩形微通道內(nèi)流動的對流換熱系數(shù)為6,320 W/(m2·K)[12]。隨著網(wǎng)格數(shù)的增多,4種不同網(wǎng)格模擬得到的對流換熱系數(shù)與實(shí)驗(yàn)的誤差為8.90%、7.99%、6.69%和6.60%,由此可以看出,隨著網(wǎng)格數(shù)由300,000增加到840,000,誤差僅減小了0.09%。為了精確地反映對流換熱過程,又能減小計算量,模擬選取300,000網(wǎng)格數(shù)。

    1.6 數(shù)值方法的實(shí)驗(yàn)驗(yàn)證

    DALKILIC等[12]實(shí)驗(yàn)研究了R134a在矩形微通道對流換熱特性,其采用的矩形微通道的幾何尺寸和工況與本文一致,因此可以用來驗(yàn)證本文的數(shù)值方法準(zhǔn)確性。圖1為相同條件下,實(shí)驗(yàn)結(jié)果[12]與數(shù)值模擬結(jié)果的對比,可以看出最大誤差為12.9%,滿足精度要求。

    圖1 實(shí)驗(yàn)結(jié)果與數(shù)值模擬對比

    2 熵產(chǎn)計算模型

    在矩形微通道中,熵產(chǎn)包含由熱傳遞產(chǎn)生的熱熵產(chǎn)和壓降造成的摩擦熵產(chǎn)兩個不可逆因素。

    其中¢gen為總熵產(chǎn),¢gen, th為熱熵產(chǎn),¢gen, fr為摩擦熵產(chǎn)。

    由熱力學(xué)第一和第二定律分析得,穩(wěn)態(tài)系統(tǒng)中的熵產(chǎn)公式如下[13]:

    由上面式子化簡得:

    重新定義摩擦因子()和數(shù)():

    用上面的公式,熵產(chǎn)可以表示為:

    代入下面的特征數(shù):

    熵產(chǎn)公式變成:

    和的定義如下:

    式中:

    ——單位控制體的長度;

    ——流體的平均溫度,K;

    +D——壁面溫度,K;

    ——流體的速度,m/s;

    ¢——每個單位長度的熱流密度,W/(m·K);

    ——質(zhì)量流量,kg/s;

    ——橫截面積,m2;

    ——潤濕周長,m;

    ——運(yùn)動粘度,m2/s;

    ——熱擴(kuò)散率,m2/s。

    3 結(jié)果討論

    3.1 對流換熱及壓降特性

    圖2展示了R134a-CuO和R134a-CNT的對流換熱系數(shù)在不同體積分?jǐn)?shù)和不同雷諾數(shù)下的變化。對流換熱系數(shù)隨著體積分?jǐn)?shù)和雷諾數(shù)的增加而增加。在=1,200,納米顆粒體積分?jǐn)?shù)為1%和4%時,R134a-CuO對流換熱系數(shù)相比于純R134a分別增加了20.7%和45.7%,R134a-CNT比純R134a分別增加了31.1%和125.6%。從400增加到2,000時,體積分?jǐn)?shù)為4%的R134a-CuO和R134a-CNT的對流換熱系數(shù)分別增加了71.3%和65.8%。產(chǎn)生上述現(xiàn)象的原因與熱導(dǎo)率有關(guān)。隨著體積分?jǐn)?shù)的增加,納米制冷劑的熱導(dǎo)率增加;納米制冷劑熱導(dǎo)率與納米粒子周圍分子的吸附層傳熱有關(guān);液態(tài)里面的納米粒子會吸附一層液體分子層,液體層中的分子排列有序并且熱導(dǎo)率接近納米粒子的熱導(dǎo)率;納米粒子比表面積非常大,吸附的液體分子數(shù)量也比較大,因此納米粒子表面的液體吸附層對納米制冷劑導(dǎo)熱系數(shù)的強(qiáng)化有較大作用[8]。對流換熱系數(shù)還與納米粒子的熱遷移有關(guān)。熱遷移使得納米粒子更多的在微通道的中心處聚集,對流換熱系數(shù)增大[14];隨著雷諾數(shù)的增加,由于納米粒子對粘性底層的擾動隨著速度的增加而增加,從而邊界層變薄,對流換熱系數(shù)增加。

    圖3表示了納米制冷劑在不同體積分?jǐn)?shù)不同雷諾數(shù)下的壓降。由圖中可以看出,壓降隨體積分?jǐn)?shù)和雷諾數(shù)的增加而增加。在=1,200時,納米制冷劑4%的R134a-CuO和R134a-CNT比純R134a壓降分別增加了192.4%和161.8%。雷諾數(shù)從400增加到2,000時,體積分?jǐn)?shù)為2%的R134a-CuO和R134a-CNT納米制冷劑的壓降分別增加了546.2%和537.8%。隨著體積分?jǐn)?shù)的增大,納米制冷劑的粘度增加,摩擦阻力增大,從而壓降上升;由于納米粒子的移動,產(chǎn)生了局部渦流,增加了粘性耗散。

    圖2 對流換熱系數(shù)隨體積分?jǐn)?shù)和雷諾數(shù)的變化

    圖3 壓降隨體積分?jǐn)?shù)和雷諾數(shù)的變化

    3.2 熵產(chǎn)分析

    圖4和圖5顯示了摩擦熵產(chǎn)和熱熵產(chǎn)隨體積分?jǐn)?shù)和雷諾數(shù)的變化。由圖可以看出,摩擦熵產(chǎn)隨雷諾數(shù)的增加而增加,熱熵產(chǎn)隨雷諾數(shù)的增加而減小。雷諾數(shù)從400增加到2,000,體積分?jǐn)?shù)4%的R134a-CuO和R134a-CNT熱熵產(chǎn)分別減小了41.1%和39.2%;摩擦熵產(chǎn)隨體積分?jǐn)?shù)的增加而增加,熱熵產(chǎn)隨體積分?jǐn)?shù)的增加而減小。在=2,000時,4%的R134a-CuO和R134a-CNT的摩擦熵產(chǎn)相比純R134a分別增加了363.2%和316.1%,熱熵產(chǎn)分別減小了45.9%和55.2%。原因是隨著雷諾數(shù)的增加,速度增大導(dǎo)致了微通道內(nèi)壓降以及摩擦損失的增加,使得摩擦熵產(chǎn)增加。速度增大也使邊界層的厚度變薄,促進(jìn)了運(yùn)動的流體與固體壁面之間的對流換熱,從而熱熵產(chǎn)減小;隨著體積分?jǐn)?shù)的增加,雷諾數(shù)不變的情況下,納米制冷劑的粘度增加,從而摩擦熵產(chǎn)增加。體積分?jǐn)?shù)的增加提高了固體壁面與流體的對流換熱,減小了熱耗散,從而使熱熵產(chǎn)大大減小[15]。

    圖4 摩擦熵產(chǎn)隨體積分?jǐn)?shù)和雷諾數(shù)的變化

    圖5 熱熵產(chǎn)隨體積分?jǐn)?shù)和雷諾數(shù)的變化

    圖6表示總熵產(chǎn)隨著體積分?jǐn)?shù)和雷諾數(shù)的變化。純R134a的總熵產(chǎn)隨著雷諾數(shù)的增加而減小。體積分?jǐn)?shù)為1%的R134a-CuO和R134a-CNT隨著雷諾數(shù)的增加,總熵產(chǎn)減小且減小的幅度越來越小。對于體積分?jǐn)?shù)為2%~4%的R134a-CuO和R134a-CNT,隨雷諾數(shù)的增加,總熵產(chǎn)先減小后增加。產(chǎn)生上述現(xiàn)象的原因?yàn)椋弘S著雷諾數(shù)的增大,摩擦熵產(chǎn)的增加與熱熵產(chǎn)的減小出現(xiàn)了相互較量。低體積分?jǐn)?shù)時,摩擦熵產(chǎn)所占的份額遠(yuǎn)遠(yuǎn)小于熱熵產(chǎn);隨著雷諾數(shù)的增大,摩擦熵產(chǎn)的增加幅度小于熱熵產(chǎn)的減小幅度,從而總熵產(chǎn)減小。體積分?jǐn)?shù)較高時,存在一個最佳雷諾數(shù)值使得熵產(chǎn)達(dá)到最小且最佳雷諾數(shù)隨體積分?jǐn)?shù)的增加而減小,一旦超過對應(yīng)的最佳雷諾數(shù),摩擦熵產(chǎn)增加的幅度超過了熱熵產(chǎn)減小的幅度,熵產(chǎn)將會隨著雷諾數(shù)的增加而增加。為使微通道內(nèi)的熵產(chǎn)最小,需要同時考慮體積分?jǐn)?shù)和雷諾數(shù)作用[16]。本文模擬的層流條件下,R134a-CuO的熵產(chǎn)最小條件為體積分?jǐn)?shù)2%和=1,600;R134a-CNT的熵產(chǎn)最小條件為體積分?jǐn)?shù)4%和=800。

    圖6 總熵產(chǎn)隨體積分?jǐn)?shù)和雷諾數(shù)的變化

    4 結(jié)論

    1)對流換熱系數(shù)隨著體積分?jǐn)?shù)和雷諾數(shù)的增加而增加。壓降隨著體積分?jǐn)?shù)和雷諾數(shù)的增加而上升。在=1,200時,體積分?jǐn)?shù)為4%的R134a-CuO和R134a-CNT對流換熱系數(shù)相比于純R134a增加了45.7%和125.6%。雷諾數(shù)從400增加到2,000,體積分?jǐn)?shù)為2%的R134a-CuO和R134a-CNT納米制冷劑的壓降分別增加了546.2%和537.8%。

    2)隨著雷諾數(shù)的增加,熱熵產(chǎn)減小,摩擦熵產(chǎn)增加。低體積分?jǐn)?shù)時,總熵產(chǎn)隨著雷諾數(shù)的增加而減小。體積分?jǐn)?shù)較高時,存在一個最佳雷諾數(shù)值使得熵產(chǎn)達(dá)到最小且最佳雷諾數(shù)隨體積分?jǐn)?shù)的增加而減小。

    [1] 張蕾. 微通道換熱器空調(diào)系統(tǒng)性能試驗(yàn)分析[J]. 制冷技術(shù), 2010, 30(1): 33-36.

    [2] 汪年結(jié), 王穎, 李峰, 等. 微通道換熱器在家用空調(diào)上的應(yīng)用研究[J]. 制冷技術(shù), 2014, 34(2): 47-50.

    [3] 姚奕, 徐柏興, 邵翌旻, 等. R1234yf微通道平行流冷凝器仿真模型[J]. 制冷技術(shù), 2013, 33(3): 34-38.

    [4] SUN B, YANG D. Experimental study on the heat transfer characteristics of nanorefrigerants in an internal thread copper tube[J]. International Journal of Heat and Mass Transfer, 2013, 64(38): 559-566.

    [5] MAHBUBUL I M, SAIDUR R, AMALINA M A. Heat transfer and pressure drop characteristics of Al2O3-R141b nanorefrigerant in horizontal smooth circular tube[J]. Procedia Engineering, 2013, 56: 323-329.

    [6] MANAY E, AKYUREK E F, SAHIN B. Entropy generation of nanofluid flow in a microchannel heat sink[J]. Results in Physics, 2018, 9: 615-624.

    [7] XIE Y, ZHENG L, ZHANG D, et al. Entropy generation and heat transfer performances of Al2O3-water nanofluid transitional flow in rectangular channels with dimples and protrusions[J]. Entropy, 2016, 18(4): 148.

    [8] 姜未汀. 納米制冷劑的熱導(dǎo)率?穩(wěn)定性及納米流體電導(dǎo)率的實(shí)驗(yàn)與建模[D]. 上海: 上海交通大學(xué), 2009.

    [9] HELVACI H U, KHAN Z A. Heat transfer and entropy generation analysis of HFE 7000 based nanorefrigerants[J]. International Journal of Heat and Mass Transfer, 2017, 104: 318-327.

    [10] JIANG W T,DING G L,PENG H. Measurement and model on thermal conductivities of carbon nanotube nanorefrigerants[J]. International Journal of Thermal Sciences, 2009, 48(6): 1108-1115.

    [11] DAS N, TAKATA Y, KOHNO M, et al. Effect of carbon nano inclusion dimensionality on the melting of phase change nanocomposites in vertical shell-tube thermalenergy storage unit[J]. International Journal of Heat and Mass Transfer, 2017, 113: 423-431.

    [12] DALKILIC A S, MAHIAN O, YILMAZ S, et al. Experimental investigation of single-phase turbulent flow of R-134a in a multiport microchannel heat sink[J]. International Communications in Heat and Mass Transfer, 2017, 89: 47-56.

    [13] MOGHADDAMI M, MOHAMMADZADE A, ESFEHANI S A V. Second law analysis of nanofluid flow[J]. Energy Conversion and Management, 2011, 52(2): 1397-1405.

    [14] BAHIRAEI M, HOSSEINALIPOUR S M. Particle migration in nanofluids considering thermophoresis and its effect on convective heat transfer[J]. Thermochimica Acta, 2013, 574(2): 47-54.

    [15] SOHEL M R, SAIDUR R, HASSAN N H, et al. Analysis of entropy generation using nanofluid flow through the circular microchannel and minichannel heat sink[J]. International Communications in Heat and Mass Transfer, 2013, 46(46): 85-91.

    [16] MAHIAN O, KIANIFAR A, KLEINSTREUER C, et al. A review of entropy generation in nanofluid flow[J]. International Journal of Heat and Mass Transfer, 2013, 65(5): 514-532.

    Analysis on Convective Heat Transfer Characteristics and Entropy Generation of Nanorefrigerant Flow in Microchannel

    GUO Wenhua, PENG Hao*

    (Merchant Marine College, Shanghai Maritime University, Shanghai 201306, China)

    In the present study, the convective heat transfer and entropy generation of R134a nanorefrigerant flow in microchannel are numerically investigated. The nano-powders include cupric oxide (CuO) and carbon nanotube (CNT), the volume fraction range is from 0 to 4%, and Reynolds number () range is from 400 to 2,000. The results indicate that the convective heat transfer coefficient and frictional pressure drop increase with the increase of volume fraction and; whenis 1,200, compared to pure R134a, the convective heat transfer coefficients of R134a-CuO and R134a-CNT nanorefrigerants with volume fraction of 4% are increased by 45.7% and 125.7%, respectively, while the frictional pressure drops are increased by 192.4% and 161.8%, respectively. With the increase ofand volume fraction, the thermal entropy generation decreases while the frictional entropy generation increases; the minimum entropy generation for R134a-CuO and R134a-CNT nanorefrigerants occurs at conditions of 2% volume fraction combining=1,600, and 4% volume fraction combining=800, respectively.

    Nanorefrigerant; Convection heat transfer; Entropy generation; Microchannel; Computer simulation

    10.3969/j.issn.2095-4468.2018.04.102

    *彭浩(1982-),男,副教授,博士。研究方向:納米流體、相變傳熱。聯(lián)系地址:上海市浦東新區(qū)海港大道1550號上海海事大學(xué),郵編:201306。聯(lián)系電話:021-38282960。E-mail:hpeng@shmtu.edu.cn。

    上海高校青年東方學(xué)者人才計劃項(xiàng)目(No. QD2016045)。

    猜你喜歡
    雷諾數(shù)制冷劑對流
    揚(yáng)州市出臺“制冷劑六條”提升制冷劑安全管理水平
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    近共沸制冷劑R134a/R1234yf的PVTx性質(zhì)的實(shí)驗(yàn)研究
    基于ANSYS的自然對流換熱系數(shù)計算方法研究
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計研究
    民機(jī)高速風(fēng)洞試驗(yàn)的阻力雷諾數(shù)效應(yīng)修正
    二元驅(qū)油水界面Marangoni對流啟動殘余油機(jī)理
    霍尼韋爾新型制冷劑獲歐盟科研機(jī)構(gòu)認(rèn)可
    久久久久久九九精品二区国产| 黄色一级大片看看| 精品免费久久久久久久清纯| 波多野结衣高清作品| 麻豆久久精品国产亚洲av| 亚洲自偷自拍三级| 久久精品国产自在天天线| 动漫黄色视频在线观看| 在线播放国产精品三级| 最近最新免费中文字幕在线| netflix在线观看网站| 尤物成人国产欧美一区二区三区| 日本撒尿小便嘘嘘汇集6| 又爽又黄无遮挡网站| 成人亚洲精品av一区二区| 亚洲人成网站高清观看| 亚洲精品一卡2卡三卡4卡5卡| 伦精品一区二区三区| 亚洲国产精品久久男人天堂| 日韩人妻高清精品专区| 国产aⅴ精品一区二区三区波| 久久精品综合一区二区三区| av国产免费在线观看| 国产视频一区二区在线看| 国产欧美日韩精品一区二区| 天堂网av新在线| 精品一区二区三区av网在线观看| 亚洲av电影不卡..在线观看| 欧美中文日本在线观看视频| 麻豆精品久久久久久蜜桃| 内地一区二区视频在线| 深爱激情五月婷婷| 亚洲av成人av| 成人一区二区视频在线观看| 亚洲图色成人| 性欧美人与动物交配| 一级黄片播放器| 国产亚洲欧美98| 国产精品永久免费网站| 国产欧美日韩精品一区二区| 久久精品91蜜桃| 深夜a级毛片| 国产精品人妻久久久影院| 欧美成人性av电影在线观看| 国模一区二区三区四区视频| 成人亚洲精品av一区二区| 精品久久久久久久人妻蜜臀av| 少妇高潮的动态图| 男女做爰动态图高潮gif福利片| 免费大片18禁| 亚洲avbb在线观看| 亚洲av五月六月丁香网| 国产精品,欧美在线| 久久天躁狠狠躁夜夜2o2o| 91精品国产九色| 亚洲最大成人中文| 九九久久精品国产亚洲av麻豆| 免费看美女性在线毛片视频| 免费看av在线观看网站| 91在线观看av| 亚洲国产高清在线一区二区三| 久久欧美精品欧美久久欧美| 琪琪午夜伦伦电影理论片6080| 丝袜美腿在线中文| 久久久久久国产a免费观看| 俺也久久电影网| 九色成人免费人妻av| 久久精品久久久久久噜噜老黄 | 亚洲18禁久久av| 免费不卡的大黄色大毛片视频在线观看 | 国内毛片毛片毛片毛片毛片| 色精品久久人妻99蜜桃| 夜夜夜夜夜久久久久| 色视频www国产| 天堂网av新在线| 婷婷精品国产亚洲av| 国产三级在线视频| 免费高清视频大片| 亚洲四区av| 直男gayav资源| 国产精品99久久久久久久久| 日韩欧美精品v在线| 国国产精品蜜臀av免费| 人妻少妇偷人精品九色| 99精品在免费线老司机午夜| 成年女人看的毛片在线观看| 日韩欧美三级三区| 少妇猛男粗大的猛烈进出视频 | 国产高清三级在线| 波野结衣二区三区在线| 国产高清激情床上av| 国产真实伦视频高清在线观看 | 日韩高清综合在线| 女的被弄到高潮叫床怎么办 | 国产三级在线视频| 亚洲熟妇中文字幕五十中出| 日韩亚洲欧美综合| 又黄又爽又刺激的免费视频.| av专区在线播放| 中文亚洲av片在线观看爽| 色5月婷婷丁香| 在线观看av片永久免费下载| 亚洲国产欧美人成| 一进一出抽搐动态| 日本免费a在线| 一进一出抽搐动态| 啦啦啦观看免费观看视频高清| 亚洲av五月六月丁香网| 99久久精品热视频| 国模一区二区三区四区视频| 别揉我奶头 嗯啊视频| 欧美人与善性xxx| 日韩精品青青久久久久久| 波多野结衣高清作品| 99久久精品国产国产毛片| 精品无人区乱码1区二区| 别揉我奶头~嗯~啊~动态视频| 国产一区二区三区av在线 | 99在线视频只有这里精品首页| 色精品久久人妻99蜜桃| 免费人成在线观看视频色| 久久久久久九九精品二区国产| 88av欧美| 国产精品一区www在线观看 | 亚洲美女黄片视频| 久久久久久大精品| 久久久久性生活片| 国产成人一区二区在线| 少妇被粗大猛烈的视频| 色精品久久人妻99蜜桃| 亚洲av成人精品一区久久| 小说图片视频综合网站| 99久久精品热视频| 成年人黄色毛片网站| 成年人黄色毛片网站| 亚洲avbb在线观看| 少妇熟女aⅴ在线视频| 日韩欧美精品v在线| 淫秽高清视频在线观看| 美女高潮的动态| 一个人看视频在线观看www免费| 亚洲成人精品中文字幕电影| 性插视频无遮挡在线免费观看| 亚洲男人的天堂狠狠| 国产视频一区二区在线看| 日韩欧美三级三区| bbb黄色大片| 老司机深夜福利视频在线观看| 精华霜和精华液先用哪个| 69av精品久久久久久| 蜜桃久久精品国产亚洲av| 亚洲国产高清在线一区二区三| 丰满的人妻完整版| 一级毛片久久久久久久久女| 日本免费a在线| 好男人在线观看高清免费视频| 亚洲成人久久爱视频| 国产又黄又爽又无遮挡在线| 中文资源天堂在线| 国产精品亚洲一级av第二区| 91午夜精品亚洲一区二区三区 | 看黄色毛片网站| 欧美又色又爽又黄视频| 色播亚洲综合网| 18禁黄网站禁片午夜丰满| 黄色欧美视频在线观看| 精品一区二区三区视频在线观看免费| 国产精品98久久久久久宅男小说| 欧美另类亚洲清纯唯美| 中国美女看黄片| 嫁个100分男人电影在线观看| 热99re8久久精品国产| 日韩欧美在线二视频| 亚洲中文日韩欧美视频| 国产老妇女一区| 精品国产三级普通话版| 少妇裸体淫交视频免费看高清| 俺也久久电影网| 如何舔出高潮| 精品久久久久久久人妻蜜臀av| 日本黄大片高清| 男女下面进入的视频免费午夜| 麻豆国产97在线/欧美| 免费观看的影片在线观看| 午夜福利成人在线免费观看| 在线播放无遮挡| 97碰自拍视频| 天堂√8在线中文| 91久久精品国产一区二区成人| 久久精品国产亚洲网站| 十八禁国产超污无遮挡网站| 日韩中文字幕欧美一区二区| 一级a爱片免费观看的视频| 亚洲国产欧洲综合997久久,| 搡老妇女老女人老熟妇| 97超级碰碰碰精品色视频在线观看| 午夜久久久久精精品| 少妇熟女aⅴ在线视频| 国产三级在线视频| 舔av片在线| 精品无人区乱码1区二区| 国产又黄又爽又无遮挡在线| 欧美日韩国产亚洲二区| 久久热精品热| 色哟哟·www| a在线观看视频网站| 简卡轻食公司| 精品人妻一区二区三区麻豆 | a级毛片a级免费在线| 噜噜噜噜噜久久久久久91| 啪啪无遮挡十八禁网站| 精品一区二区三区视频在线| 日本在线视频免费播放| 成人精品一区二区免费| 又紧又爽又黄一区二区| 美女免费视频网站| 国产黄色小视频在线观看| 精品人妻一区二区三区麻豆 | 看免费成人av毛片| 亚洲精品日韩av片在线观看| 乱人视频在线观看| 久久99热6这里只有精品| 亚洲性久久影院| 亚洲av.av天堂| 两人在一起打扑克的视频| 亚洲av中文字字幕乱码综合| 人人妻人人澡欧美一区二区| 一级毛片久久久久久久久女| 99热这里只有是精品50| 九九久久精品国产亚洲av麻豆| 久9热在线精品视频| 欧美+日韩+精品| 88av欧美| 久久精品国产鲁丝片午夜精品 | 亚洲人成伊人成综合网2020| 国内精品一区二区在线观看| 我要看日韩黄色一级片| 大又大粗又爽又黄少妇毛片口| 日韩一本色道免费dvd| 欧美成人a在线观看| 天堂av国产一区二区熟女人妻| 少妇人妻一区二区三区视频| 亚洲人与动物交配视频| 欧美最新免费一区二区三区| 亚洲性夜色夜夜综合| 国产单亲对白刺激| 一区二区三区高清视频在线| www.www免费av| 欧美人与善性xxx| 深夜精品福利| 日韩精品中文字幕看吧| 亚洲第一区二区三区不卡| 俄罗斯特黄特色一大片| 国产精品久久久久久久久免| 1000部很黄的大片| 啪啪无遮挡十八禁网站| 中文字幕av在线有码专区| 搡女人真爽免费视频火全软件 | 免费观看的影片在线观看| 欧美绝顶高潮抽搐喷水| 久久久久久大精品| 国产欧美日韩精品一区二区| 国产高清激情床上av| 51国产日韩欧美| 国产亚洲精品久久久久久毛片| 少妇的逼好多水| 亚洲人成网站在线播| 日韩一本色道免费dvd| 最后的刺客免费高清国语| 色5月婷婷丁香| 哪里可以看免费的av片| 国产精品久久电影中文字幕| 亚洲美女黄片视频| 毛片女人毛片| 在线天堂最新版资源| 亚洲欧美日韩卡通动漫| 搞女人的毛片| 国产欧美日韩精品亚洲av| 久久草成人影院| 永久网站在线| 男女之事视频高清在线观看| 午夜久久久久精精品| www日本黄色视频网| bbb黄色大片| 欧美三级亚洲精品| 中国美女看黄片| 色5月婷婷丁香| 精品久久国产蜜桃| 啦啦啦韩国在线观看视频| 男女做爰动态图高潮gif福利片| 美女 人体艺术 gogo| 亚洲精品一卡2卡三卡4卡5卡| 日韩强制内射视频| 成人综合一区亚洲| 国产私拍福利视频在线观看| 又爽又黄a免费视频| 精华霜和精华液先用哪个| 一级黄片播放器| 1000部很黄的大片| 麻豆国产av国片精品| 少妇高潮的动态图| 99国产精品一区二区蜜桃av| 伊人久久精品亚洲午夜| 99热只有精品国产| 简卡轻食公司| 一区二区三区激情视频| 国模一区二区三区四区视频| 一个人看的www免费观看视频| 亚洲狠狠婷婷综合久久图片| 尤物成人国产欧美一区二区三区| 欧美成人a在线观看| 久久久国产成人精品二区| 欧美日韩黄片免| 欧美日韩乱码在线| 中文字幕av在线有码专区| 欧美精品国产亚洲| 国产色婷婷99| 国产精品,欧美在线| 欧美高清性xxxxhd video| 日韩欧美精品v在线| 在线免费观看不下载黄p国产 | 99精品久久久久人妻精品| 麻豆成人av在线观看| 亚洲无线在线观看| 亚洲av不卡在线观看| av在线蜜桃| 免费看美女性在线毛片视频| 久久久久久久久久久丰满 | 村上凉子中文字幕在线| 岛国在线免费视频观看| 又黄又爽又刺激的免费视频.| 在线观看av片永久免费下载| 一夜夜www| 亚洲不卡免费看| 日本黄大片高清| 免费观看人在逋| 精品人妻一区二区三区麻豆 | 赤兔流量卡办理| 97人妻精品一区二区三区麻豆| 亚洲黑人精品在线| 级片在线观看| 国产精品一及| 日本一本二区三区精品| 久久精品国产清高在天天线| 亚洲图色成人| 一区福利在线观看| 亚洲美女搞黄在线观看 | a级毛片a级免费在线| 99热精品在线国产| 男女边吃奶边做爰视频| 日韩一区二区视频免费看| 五月玫瑰六月丁香| 国产av一区在线观看免费| 999久久久精品免费观看国产| 淫秽高清视频在线观看| 高清在线国产一区| 男人舔奶头视频| 国产精品一区www在线观看 | 国内精品久久久久久久电影| 人人妻人人澡欧美一区二区| 久久久国产成人精品二区| 亚洲乱码一区二区免费版| 热99re8久久精品国产| 日韩 亚洲 欧美在线| 亚洲国产精品久久男人天堂| 成人国产麻豆网| 色av中文字幕| 国产 一区精品| 干丝袜人妻中文字幕| 国产精品福利在线免费观看| 成人鲁丝片一二三区免费| 国产av在哪里看| 麻豆精品久久久久久蜜桃| 神马国产精品三级电影在线观看| 精品久久久久久久末码| 女同久久另类99精品国产91| 欧美成人一区二区免费高清观看| 国产精品98久久久久久宅男小说| 夜夜爽天天搞| 亚洲av二区三区四区| 99久久无色码亚洲精品果冻| 婷婷精品国产亚洲av| 最近最新免费中文字幕在线| av天堂中文字幕网| 在线观看午夜福利视频| 国产精品98久久久久久宅男小说| 嫩草影院入口| 国产熟女欧美一区二区| 国产欧美日韩一区二区精品| 亚洲一区二区三区色噜噜| 在现免费观看毛片| 91久久精品电影网| 国产美女午夜福利| 自拍偷自拍亚洲精品老妇| 夜夜夜夜夜久久久久| 此物有八面人人有两片| 国产视频内射| 18禁黄网站禁片免费观看直播| xxxwww97欧美| 97人妻精品一区二区三区麻豆| 美女高潮的动态| 国产真实伦视频高清在线观看 | 精华霜和精华液先用哪个| 特大巨黑吊av在线直播| 国产蜜桃级精品一区二区三区| 国产精品久久电影中文字幕| 99久久中文字幕三级久久日本| 尤物成人国产欧美一区二区三区| 久久久久久伊人网av| 久久人妻av系列| 日韩欧美一区二区三区在线观看| 日本在线视频免费播放| 亚洲自偷自拍三级| 亚洲三级黄色毛片| 淫秽高清视频在线观看| 不卡视频在线观看欧美| 精品久久久久久久末码| 国产一区二区三区视频了| 午夜免费激情av| 午夜福利在线观看免费完整高清在 | 最好的美女福利视频网| 91麻豆av在线| 在线观看一区二区三区| 国产 一区 欧美 日韩| 免费不卡的大黄色大毛片视频在线观看 | 久久久久精品国产欧美久久久| 久久精品国产亚洲av香蕉五月| 校园春色视频在线观看| 日本 av在线| 桃红色精品国产亚洲av| 亚洲人成网站高清观看| 国国产精品蜜臀av免费| 国产伦一二天堂av在线观看| 99国产精品一区二区蜜桃av| 99视频精品全部免费 在线| 国内精品宾馆在线| 嫩草影视91久久| 色视频www国产| 3wmmmm亚洲av在线观看| 中文在线观看免费www的网站| 99热网站在线观看| 日日夜夜操网爽| 久久久久久国产a免费观看| 在线观看舔阴道视频| 免费看美女性在线毛片视频| h日本视频在线播放| 日本一二三区视频观看| 91久久精品电影网| 麻豆av噜噜一区二区三区| 97人妻精品一区二区三区麻豆| 亚洲精品456在线播放app | 亚洲不卡免费看| 亚洲一级一片aⅴ在线观看| 久久久久精品国产欧美久久久| 国产91精品成人一区二区三区| 真实男女啪啪啪动态图| 少妇人妻一区二区三区视频| 精品99又大又爽又粗少妇毛片 | 免费在线观看日本一区| 亚洲无线在线观看| 国产精品国产高清国产av| 在线观看午夜福利视频| 我的老师免费观看完整版| 欧美xxxx性猛交bbbb| 无遮挡黄片免费观看| 亚洲图色成人| 哪里可以看免费的av片| 熟妇人妻久久中文字幕3abv| 午夜福利在线观看吧| 亚洲人成网站在线播放欧美日韩| 国内精品宾馆在线| 最好的美女福利视频网| 直男gayav资源| АⅤ资源中文在线天堂| 欧美激情国产日韩精品一区| 最近最新免费中文字幕在线| 国产极品精品免费视频能看的| 久久久久九九精品影院| 91午夜精品亚洲一区二区三区 | 国产精品久久久久久精品电影| 级片在线观看| 丰满乱子伦码专区| 亚洲va日本ⅴa欧美va伊人久久| 99久久精品热视频| 日本一二三区视频观看| 99国产极品粉嫩在线观看| 免费在线观看成人毛片| 久久这里只有精品中国| 成人欧美大片| 最近中文字幕高清免费大全6 | 国产精品自产拍在线观看55亚洲| 人妻制服诱惑在线中文字幕| 男插女下体视频免费在线播放| 91久久精品电影网| 又爽又黄无遮挡网站| 99热网站在线观看| 在线国产一区二区在线| 日日干狠狠操夜夜爽| 亚洲电影在线观看av| 亚洲人成伊人成综合网2020| 精华霜和精华液先用哪个| av在线老鸭窝| 美女xxoo啪啪120秒动态图| 欧美最新免费一区二区三区| 一进一出抽搐动态| 精华霜和精华液先用哪个| 国产高清三级在线| 俄罗斯特黄特色一大片| h日本视频在线播放| 简卡轻食公司| 黄片wwwwww| 欧美一级a爱片免费观看看| 18禁在线播放成人免费| 熟女电影av网| 色视频www国产| 久久精品国产清高在天天线| a级一级毛片免费在线观看| 久久6这里有精品| 不卡一级毛片| 少妇人妻精品综合一区二区 | 中国美白少妇内射xxxbb| 午夜福利在线观看吧| 国产女主播在线喷水免费视频网站 | 动漫黄色视频在线观看| 精品无人区乱码1区二区| eeuss影院久久| 欧美一级a爱片免费观看看| 免费观看人在逋| 欧美三级亚洲精品| 男人狂女人下面高潮的视频| 免费电影在线观看免费观看| 亚洲精华国产精华液的使用体验 | 久久精品综合一区二区三区| 成年版毛片免费区| 真人一进一出gif抽搐免费| 国产精品无大码| 一本精品99久久精品77| av在线天堂中文字幕| 麻豆成人av在线观看| 国产精品美女特级片免费视频播放器| 一级a爱片免费观看的视频| 亚洲天堂国产精品一区在线| av在线天堂中文字幕| 在线天堂最新版资源| 桃色一区二区三区在线观看| 在线免费观看的www视频| 午夜久久久久精精品| 欧美性感艳星| 国产午夜福利久久久久久| 男人狂女人下面高潮的视频| 免费高清视频大片| 国产精品国产高清国产av| 国产精品一区www在线观看 | 国产在视频线在精品| 亚洲狠狠婷婷综合久久图片| 在线天堂最新版资源| 女人被狂操c到高潮| 亚洲美女视频黄频| 人人妻人人看人人澡| 欧美成人一区二区免费高清观看| 中文字幕av成人在线电影| 亚洲美女搞黄在线观看 | 久久久色成人| 午夜福利欧美成人| 久99久视频精品免费| 色综合婷婷激情| 91久久精品国产一区二区成人| 国产一区二区三区视频了| 亚洲一区高清亚洲精品| 久久中文看片网| 韩国av一区二区三区四区| 精品久久久久久久末码| 人人妻,人人澡人人爽秒播| 高清在线国产一区| a级毛片免费高清观看在线播放| 成年版毛片免费区| 能在线免费观看的黄片| 日韩欧美一区二区三区在线观看| 日韩大尺度精品在线看网址| 深爱激情五月婷婷| 麻豆久久精品国产亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久伊人网av| 伦精品一区二区三区| 99在线视频只有这里精品首页| 久久人人精品亚洲av| 成人美女网站在线观看视频| 女人十人毛片免费观看3o分钟| 三级国产精品欧美在线观看| 成熟少妇高潮喷水视频| 精品久久久久久久久av| 三级男女做爰猛烈吃奶摸视频| 18禁黄网站禁片午夜丰满| 国产高清激情床上av| 国产欧美日韩精品一区二区| 欧美xxxx黑人xx丫x性爽| 午夜福利18| 午夜精品一区二区三区免费看| 99久久中文字幕三级久久日本| 看片在线看免费视频| 欧美日韩乱码在线| 又紧又爽又黄一区二区| 午夜激情福利司机影院| 日日撸夜夜添| 成人高潮视频无遮挡免费网站| 99在线视频只有这里精品首页| 国产不卡一卡二| 啦啦啦韩国在线观看视频| 久久天躁狠狠躁夜夜2o2o| 人人妻,人人澡人人爽秒播| 精品久久久久久成人av| 国产一区二区亚洲精品在线观看| 色综合色国产| 午夜爱爱视频在线播放| 日韩精品有码人妻一区|