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

    高溫CO2熱泵套管式氣冷器的仿真與優(yōu)化設(shè)計?

    2017-10-13 12:38:07劉和成趙建峰倪秒華
    制冷與空調(diào) 2017年3期
    關(guān)鍵詞:管長管程排氣

    劉和成 趙建峰 倪秒華 趙 銳

    ?

    高溫CO2熱泵套管式氣冷器的仿真與優(yōu)化設(shè)計?

    劉和成 趙建峰 倪秒華 趙銳

    (杭州佳力斯韋姆新能源科技有限公司 杭州 311241)

    以Petrov換熱關(guān)聯(lián)式為基礎(chǔ),建立了CO2氣冷器的數(shù)值仿真模型,并根據(jù)實驗數(shù)據(jù)對模型進行了修正,并通過聯(lián)立擬合的壓縮機模型,分析了CO2系統(tǒng)在不同工況下,排氣壓力、氣冷器單管長度和并聯(lián)管程數(shù)對其性能的影響:給定氣冷器布局,最優(yōu)排氣壓力隨供水溫度上升而升高;給定工況下,氣冷器CO2出口溫度隨著排氣壓力升高而降低,而熱水流量和制熱量都會增加;不同工況下,當排氣壓力低于最優(yōu)壓力時,管長的增加對系統(tǒng)COP值增大的影響非常顯著;在供水溫度低于70℃時,并聯(lián)管程數(shù)的增加使系統(tǒng)的COP值增大,但高于70℃時,并聯(lián)管程數(shù)多的系統(tǒng)COP值反而更低。

    高溫熱泵;套管式氣冷器;仿真優(yōu)化;CO2

    0 引言

    由于CO2被認為是CHC和HCHC類制冷劑的長期替代制冷劑,因而,近些年來CO2技術(shù)在制冷、熱泵和空調(diào)等領(lǐng)域的應(yīng)用研究成為當前的一個熱門課題[1]。在CO2跨臨界循環(huán)系統(tǒng)中,氣冷器的放熱過程存在較大的溫度滑移,能使水加熱到高于60℃以上的溫度,這種高溫CO2熱泵技術(shù)的開發(fā)對于工業(yè)用鍋爐的替代具有重要意義[2,3]。套管式氣冷器已廣泛應(yīng)用于CO2系統(tǒng)中,并有不少文獻以數(shù)值仿真方法對CO2套管式氣冷器的換熱性能進行了研究。Yu[4]建立了以Dang[5]換熱關(guān)聯(lián)式為基礎(chǔ)的仿真模型,并驗證螺旋套管式氣冷器沿管長方向溫度和壓力的分布與實驗的吻合程度。鐘瑜等[6]通過實驗設(shè)計4.5kW紫銅管氣冷器的CO2熱泵系統(tǒng),驗證了不同換熱關(guān)聯(lián)式的計算精度。通過文獻調(diào)研,發(fā)現(xiàn)已有的相關(guān)文獻都是局限于驗證不同計算模型的精度,而未揭示或總結(jié)高溫CO2氣冷器的優(yōu)化設(shè)計方向。

    在跨臨界CO2系統(tǒng)中,對于給定的工況,存在一個使系統(tǒng)效率最高的最優(yōu)排氣壓力。Liao[7]給出了最優(yōu)排氣壓力P的表達式:

    (1)

    式中:t為蒸發(fā)溫度,℃;t為氣冷器出口溫度,℃;,為與壓縮機性能相關(guān)的常數(shù)。

    其中氣冷器出口溫度t對最優(yōu)排氣壓力有很大的影響,而t又與氣冷器的布局設(shè)計相關(guān),因此合理的氣冷器布局設(shè)計對系統(tǒng)最優(yōu)壓力的控制以及效率的提高有重要的意義。本文正是基于這點,借助數(shù)值仿真技術(shù),研究高溫CO2系統(tǒng)中套管式氣冷器的設(shè)計參數(shù)(單管管長和管程數(shù))對最優(yōu)排氣壓力和系統(tǒng)效率的影響,并總結(jié)氣冷器設(shè)計的優(yōu)化方向。

    1 仿真模型的選擇

    1.1 換熱和壓降關(guān)聯(lián)式

    陶于兵等[8]總結(jié)了常用的換熱關(guān)聯(lián)式,它們主要分為兩類:常物性的關(guān)聯(lián)式和考慮流體物性隨溫度顯著變化的變物性關(guān)聯(lián)式。對于超臨界CO2流體,選擇變物性關(guān)聯(lián)式有利于提高計算精度。本文前期通過比較變物性關(guān)聯(lián)式Petrov[9],Pitla[10]和Yoon[11]模型的計算精度,發(fā)現(xiàn)Petrov模型的精度最高,因而選擇其作為數(shù)值仿真的基礎(chǔ)數(shù)學(xué)模型。Petrov關(guān)聯(lián)式是在Sieder-Tater[12]關(guān)聯(lián)式的基礎(chǔ)上演變過來,其表達式為:

    (2)

    式中:Nu為Sieder-Tater關(guān)聯(lián)式的努賽爾數(shù);c為流體的平均比熱,J/(kg·K);c為以壁溫為特征溫度的比熱,J/(kg·K);q為熱流密度,W/m2;G為質(zhì)量通量,kg/(s·m2);為常數(shù)指數(shù),取值由cc的關(guān)系決定。

    對于水側(cè)的換熱關(guān)聯(lián)式,本文選用Gnielinski[13]關(guān)聯(lián)式,其表達式為:

    (3)

    式中:摩擦系數(shù),選用Filonenko[14]公式;Re為雷諾數(shù);Pr為普朗特數(shù)。

    由于氣冷器的CO2側(cè)壓降相對于超臨界流體壓力來說,占比很小,所以本文用Darcy- Weisbach[15]方程來估算管路的總壓降。

    1.2 算法流程

    為了對整個氣冷器進行求解,本文基于體積微元法將氣冷器分成足夠多的子段,然后對每段采用平均對數(shù)溫差法進行換熱量校核。由于已有許多文獻對這種計算方法有詳細的介紹,因此本文不再做進一步陳述,只介紹基于Matlab仿真模型的算法流程。本文算法有別于其它文獻的是:由于所選的Petrov換熱關(guān)聯(lián)式含有壁溫,因此只有先求得CO2和水側(cè)的壁溫,然后才能計算出微元體的出口參數(shù)。如圖1所示,在經(jīng)過壁溫和換熱平衡雙重迭代算法程序后,CO2和水側(cè)的出口溫度和壓力等參數(shù)由分段子程序計算而得。主程序則通過從excel表格讀取工況參數(shù),再調(diào)用子程序而完成對整個換熱器的求解。對于模型仿真所需的CO2和水的流體物性數(shù)據(jù),本文是通過Matlab程序調(diào)用美國國家標準與技術(shù)研究院NIST的Refprop[16]數(shù)據(jù)庫,這是目前學(xué)術(shù)界認可程度最高的物性數(shù)據(jù)庫之一。

    圖1 氣冷器仿真模型的算法流程

    2 氣冷器分段仿真模型的實驗驗證

    2.1 實驗測試方法

    為了驗證氣冷器分段仿真模型的精確度,本文選取了幾何參數(shù)如表1所示的紫銅管套管式冷卻器作為實驗對象。

    表1 套管式氣冷器的幾何尺寸

    如圖2所示,把實驗測試的氣冷器分成兩段,分別測量每段CO2側(cè)和水則的進出口溫度和壓力參數(shù)。對于氣冷器分段仿真模型,分別輸入每段實驗測試所得的CO2側(cè)入口溫度tCO2和壓力PCO2以及質(zhì)量流量m,水側(cè)的出口溫度twater和壓力Pwater以及流量m,然后計算輸出CO2側(cè)的出口溫度tCO2和壓力PCO2、水側(cè)的入口溫度twater和壓力Pwater以及總換熱量。通過對比模型計算的輸出參數(shù)與實驗測試所得的參數(shù),就可以判定仿真模型的精確度。

    圖2 氣冷器分段實驗測試示意圖

    2.2 誤差分析與模型修正

    為了驗證仿真模型的適用性,本文選取了共92組實驗數(shù)據(jù)(包括供水溫度在45~80℃和熱源進水溫度在10~30℃之間變化的工況范圍,以及同一工況下由節(jié)流開度調(diào)節(jié)所引起的排氣壓力變化范圍為8~12MPa)作為參考,并根據(jù)模擬值與實驗值的總體偏差,對仿真模型進行修正。根據(jù)楊世銘等[17],對于螺旋式換熱管(由于直管過長,所以套管式氣冷器通常做成螺旋狀的換熱管),基于換熱關(guān)聯(lián)式的模型中,CO2和水側(cè)的對流換熱加強系數(shù)可以分別表示為

    為了讓模擬值與實驗值的總體偏差減少,本文在上式的基礎(chǔ)上,對CO2和水側(cè)的加強系數(shù)進行了進一步的修正。

    圖3 模擬計算溫度與實驗測試溫度值之間的誤差

    3 仿真計算結(jié)果與分析

    3.1 仿真工況與模型的輸入?yún)?shù)范圍

    由于不考慮蒸發(fā)器側(cè)的影響,本文假定蒸發(fā)溫度t和過熱度t不變,分別取值5℃。對于氣冷器側(cè),假定回水溫度t為15℃,供水溫度t在55~80℃之間每隔5℃取一個值。單管管長分別取值20m,25m和30m,并聯(lián)管程數(shù)分別取3和4作比較。取決于系統(tǒng)的節(jié)流方式和節(jié)流開度,本文假定壓縮機的排氣壓力在8~14MPa之間變化。下表2總結(jié)了氣冷器仿真模型輸入?yún)?shù)的取值范圍。

    表2 氣冷器仿真模型輸入?yún)?shù)范圍

    換熱量是評價氣冷器性能的一個重要指標,但不能反映氣冷器布局設(shè)計的改變對整個系統(tǒng)性能的影響。為了研究氣冷器的設(shè)計參數(shù)對系統(tǒng)性能的影響,本文根據(jù)壓縮機生產(chǎn)廠商的數(shù)據(jù)擬合了一個壓縮機模型,再聯(lián)立之前所建立的氣冷器仿真模型,從而就可以分析氣冷器設(shè)計布局的改變對系統(tǒng)性能(COP)的影響。

    式中:為氣冷器模型計算的制熱量,kW;為壓縮機模型計算的耗功,kW;其中制熱量是由氣冷器模型計算所得,而氣冷器模型的輸入?yún)?shù)包括排氣壓力P、排氣溫度t和CO2質(zhì)量流量m是由壓縮機模型所傳遞;壓縮機耗功完全由壓縮機模型計算所獲得。根據(jù)美國制冷工業(yè)協(xié)會壓縮機AHRI10[18]系數(shù)模型,對于本文選用的15匹活塞式壓縮機(型號為CD1500H的Dorin公司壓縮機),其性能參數(shù)可由含有4系數(shù)的二次函數(shù)計算所得:

    (7)

    式中:P為排氣壓力,MPa;為壓縮機的性能參數(shù)mt或。

    表3總結(jié)了壓縮機模型中不同性能參數(shù)計算式中的擬合系數(shù)。

    表3 壓縮機性能仿真模型系數(shù)

    3.2 仿真結(jié)果與優(yōu)化分析

    3.2.1 壓縮機排氣壓力對氣冷器性能的影響

    在3.1節(jié)中所定義的仿真輸入?yún)?shù)范圍內(nèi),下圖4中展示了管長=20m,并聯(lián)管程數(shù)=3時,不同供熱出水溫度工況下,系統(tǒng)性能參數(shù)(包括COP、制熱量、氣冷器CO2出口溫度和熱水流量)隨排氣壓力的變化。分析圖中曲線的變化趨勢,可以發(fā)現(xiàn):(1)不同的運行工況下,系統(tǒng)都存在一個最優(yōu)排氣壓力(圖4左上子圖所示),最優(yōu)排氣壓力隨供水溫度的提高而上升,供水溫度為60℃時約為10.5MPa,而80℃時則高達13.5MPa。在最優(yōu)排氣壓力兩側(cè),系統(tǒng)COP隨排氣壓力上升而增大的趨勢比之后的下降趨勢顯著得多;(2)系統(tǒng)制熱量(右上子圖)隨排氣壓力上升而增加,其變化趨勢是由陡變緩。不同運行工況下,排氣壓力由低到高變化時,換熱量由略高于45kW上升到接近65kW,增加了近40%;(3)隨著壓力的上升,氣冷器出口溫度(左下子圖)明顯降低,從工況的最低排氣壓力到最高壓力,出口溫度下降近1/2;(4)熱水流量(右下子圖)隨排氣壓力的升高而增大,從最低排氣壓力到最高壓力,熱水流量增加的幅度接近1/3。

    圖4 排氣壓力對氣冷器性能的影響

    3.2.2 管長對COP和最優(yōu)排氣壓力的影響

    圖5展示了當并聯(lián)管程數(shù)=3,而單管管長分別為20m,25m和30m時,不同供水溫度工況下,COP值隨排氣壓力的變化。通過比較系統(tǒng)COP值隨排氣壓力的變化曲線,可以發(fā)現(xiàn),管長增加(即增加換熱面積),系統(tǒng)COP增大,但增加的幅度會因排氣壓力和工況的不同而不同:當排氣壓力低于最優(yōu)排氣壓力時,增加幅度非常明顯,而當超過最優(yōu)壓力后,增加幅度變小。另外,增加管長可以降低最優(yōu)排氣壓力,從而會使系統(tǒng)有更多的節(jié)流調(diào)節(jié)余量。

    圖5 管長對COP和最優(yōu)排氣壓力的影響

    3.2.3管程數(shù)對COP和最優(yōu)排氣壓力的影響

    圖6則比較了單管管長為20m,而并聯(lián)管程數(shù)分別為3和4時,COP值隨排氣壓力的變化。在低供水溫度工況時(低于70℃),并聯(lián)管程數(shù)多的氣冷器布局,其系統(tǒng)COP值較高,最優(yōu)排氣壓力較低。但當供水溫度高于70℃時,情況恰好相反,并聯(lián)管程數(shù)少的氣冷器布局,其系統(tǒng)COP值反而更高。

    圖6 管程數(shù)對COP和最優(yōu)排氣壓力的影響

    4 結(jié)論

    氣冷器布局設(shè)計的優(yōu)化與設(shè)定的運行工況相關(guān),因而通過仿真程序分析不同工況下氣冷器布局對系統(tǒng)性能的影響,有利于前期系統(tǒng)的優(yōu)化設(shè)計。本文通過建立以Petrov換熱關(guān)聯(lián)式為基礎(chǔ)的微元分段數(shù)學(xué)模型,用Matlab軟件編寫了套管式氣冷器的仿真程序,并選取了不同工況下的92組實驗數(shù)據(jù)對模型的精確度進行了驗證。同時,本文根據(jù)Dorin公司提供的壓縮機性能數(shù)據(jù),擬合了基于AHRI10系數(shù)的壓縮機模型。通過聯(lián)立氣冷器和壓縮機模型,在假定CO2熱泵系統(tǒng)的蒸發(fā)溫度t和過熱度t分別為5℃(即不考慮蒸發(fā)器側(cè)的影響)以及回水溫度t為15℃時,分析并比較了供水溫度工況為55~80℃時,壓縮機排氣壓力對系統(tǒng)性能的影響。分析結(jié)果表明:最優(yōu)排氣壓力隨供水溫度升高而上升;在未達到最優(yōu)壓力之前,系統(tǒng)COP值隨排氣壓力升高而增加的趨勢比越過最優(yōu)壓力之后隨排氣壓力上升而下降的趨勢顯著得多;系統(tǒng)制熱量隨排氣壓力升高而增大,上升的趨勢在越過最優(yōu)壓力之后變得平緩;氣冷器出口溫度隨排氣壓力升高而下降,同樣這種趨勢也由陡峭變?yōu)槠骄?,從工況的最低排氣壓力到最高壓力,出口溫度下降近1/2;熱水流量隨排氣壓力升高而增大,從最低壓力到最高壓力,上升幅度達1/3。

    同時,本文也比較所設(shè)計氣冷器的單管管長和并聯(lián)管程數(shù)對系統(tǒng)性能和最優(yōu)排氣壓力的影響:不同供水溫度下,系統(tǒng)COP都會隨管長的增長(即換熱面積的增大)而增大,而最優(yōu)排氣壓力隨管長的增長而小幅度降低,但越過最優(yōu)壓力之后,管長的增長對系統(tǒng)COP值增大的影響變得不明顯;當供水溫度低于70℃時,并聯(lián)管程數(shù)的增加會使系統(tǒng)的COP值增加,最優(yōu)排氣壓力降低,但當溫度高于70℃時,并聯(lián)管程數(shù)少的系統(tǒng)COP值反而越高,即高供水溫度工況下,宜選擇并聯(lián)管程數(shù)少的氣冷器布局。

    [1] Cecchinato L, Corradi M, Fornasieri E, et al. Carbon dioxide as refrigerant for tap water heat pumps: a comparison with the traditional solution[J]. International Journal of Refrigeration, 2005,28(8):1250-1258.

    [2] 趙力.高溫熱泵在我國的應(yīng)用及研究進展[J].制冷學(xué)報,2005,26(2):8-13.

    [3] 胡斌,王文毅,王凱,等.高溫熱泵技術(shù)在工業(yè)制冷領(lǐng)域的應(yīng)用[J].制冷學(xué)報,2011,32(5):1-5.

    [4] Yu P Y, Lin W K, Wang C C. Performance evaluation of a tube-in-tube CO2gas cooler used in a heat pump water heater[J]. Experimental Thermal & Fluid Science, 2014,54(4):304-312.

    [5] Dang C, Hihara E. In-tube cooling heat transfer of supercritical carbon dioxide. Part 1. Experimental measurement[J]. International Journal of Refrigeration, 2004,27(7):736-747.

    [6] 鐘瑜,賈磊,劉期聶,等.CO2套管式氣冷器設(shè)計及試驗研究[J].制冷與空調(diào),2013,13(9):45-47.

    [7] Liao S M, Zhao T S, Jakobsen A. A correlation of optimal heat rejection pressures in transcritical carbon dioxide cycles[J]. Applied Thermal Engineering, 2000,20(9):831-841.

    [8] 陶于兵,吳志根,周俊杰,等.CO2超臨界管內(nèi)冷卻換熱和壓降關(guān)聯(lián)式的比較[J].暖通空調(diào),2006,36(3):25-28.

    [9] Petrov NE, Popov VN. Heat Transfer and Resistance of Carbon Dioxide Cooled in the Supercritical Region[J]. Thermal Engineering, 1985,32(3):131-134.

    [10] Pitla S S, Groll E A, Ramadhyani S. New correlation to predict the heat transfer coefficient during in-tube cooling of turbulent supercritical CO2[J]. International Journal of Refrigeration, 2002,25(7):887-895.

    [11] Seok Ho Yoon, Ju Hyok Kim, Yun Wook Hwang, et al. Heat transfer and pressure drop characteristics during the in-tube cooling process of carbon dioxide in the supercritical region[J]. International Journal of Refrigeration, 2003,26(8):857-864.

    [12] Sieder E N, Tate G E. Heat transfer and pressure drop of liquids in tubes[J]. Industrial & Engineering Chemistry, 1936,28(12):1429-1435.

    [13] Gnielinski V. new equations for heat and mass-transfer in turbulent pipe and channel flow[J]. International Chemical Engineering, 1976,16(2):359-368.

    [14] Filonenko G K. Hydraulic Resistance in Pipes[J]. Teploergetica, 1954,1(4):40-44.

    [15] Gilley J E, Kottwitz E R. Darcy-Weisbach Roughness Coefficients for Surfaces with Residue and Gravel Cover[J]. Transactions of the Asae, 1995,38(2):539-544.

    [16] McLinden M O, Klein S A, Lemmon E W, et al. Thermodynamic properties of refrigerants and refrigerant mixtures database (REFPROP)[J]. Gaithersburg, MD: NIST, 1998.

    [17] 楊世銘,陶文銓.傳熱學(xué)(第四版)[M].北京:高等教育出版社,2006:249-250.

    [18] ANSI/AHRI 540-2004: Performance Rating of Positive Displacement Refrigerant Compressors and Compressor Units[S]. 2004.

    Simulation and Optimal Design of Tube-in-tube Gas Cooler for High-temperature CO2Heat Pump

    Liu Hecheng Zhao Jianfeng Ni Miaohua Zhao Rui

    ( Zhejiang Jiali Technology Co., Ltd, Hangzhou, 311241 )

    Based on the Petrov model for supercritical CO2heat transfer, a numerical model for gas cooler was developed and then corrected according to the experimental data. By combining the fitted compressor model, the effect of discharge pressure, single tube length and tube passes on the system performance under different operating conditions was investigated: for a given gas cooler configuration, the optimal discharge pressure increased with increasing hot water supply temperature; Given the operating condition, the gas cooler CO2outlet temperature decreased, hot water flow rate and heat capacity increased as the discharge pressure increased; The trend of increasing COP with increasing tube length before approaching the optimal discharge pressure was much more significant than that exceeding the optimal discharge pressure under different operating conditions; A larger COP was achieved by increasing tube passes when the hot water supply temperature was set below 70℃, as opposed to that at a hot water temperature higher than 70℃.

    high-temperature heat pump; tube-in-tube gas cooler; simulation and optimization; CO2

    1671-6612(2017)03-249-06

    TK11+1

    A

    劉和成(1984.02-),男,博士,博士后研究人員,E-mail:lentiancai@163.com

    2016-05-24

    猜你喜歡
    管長管程排氣
    坡地滴灌順逆坡雙向布置毛管出水規(guī)律研究
    單壁碳納米管熱導(dǎo)率幾個問題的定量描述與分析
    噴射器不同管長設(shè)計
    多管程布置微通道分液冷凝器的熱力性能
    利用管程概念求解哲學(xué)家進餐問題
    堀場制作所的新型排氣流量計
    堀場制作所的新型排氣流量計
    為RTX51Tiny項目添加管程模塊※
    排氣歧管螺栓緊固力的測定
    縮缸強化汽油機用高效排氣消聲器
    91成人精品电影| 一进一出好大好爽视频| 国产精品永久免费网站| 少妇裸体淫交视频免费看高清 | 精品国产一区二区三区四区第35| 亚洲 欧美 日韩 在线 免费| 国产成人免费无遮挡视频| 狂野欧美激情性xxxx| 天天添夜夜摸| 亚洲片人在线观看| 两个人看的免费小视频| 日本 欧美在线| 色综合亚洲欧美另类图片| 男人舔女人的私密视频| 精品卡一卡二卡四卡免费| 国产精品 欧美亚洲| 电影成人av| 国产黄a三级三级三级人| 亚洲精品国产精品久久久不卡| 99国产极品粉嫩在线观看| 免费一级毛片在线播放高清视频 | 欧美日本视频| 国产激情久久老熟女| 亚洲色图 男人天堂 中文字幕| 嫩草影院精品99| 日本 av在线| 成人av一区二区三区在线看| 丝袜美足系列| 国产激情久久老熟女| 国产欧美日韩一区二区精品| 人成视频在线观看免费观看| 亚洲熟妇中文字幕五十中出| 91成人精品电影| 国产欧美日韩一区二区三区在线| 美女国产高潮福利片在线看| 久久国产亚洲av麻豆专区| 欧美成人免费av一区二区三区| 国产精华一区二区三区| 性色av乱码一区二区三区2| 日韩欧美一区视频在线观看| 亚洲美女黄片视频| 日韩免费av在线播放| 亚洲av五月六月丁香网| 天天一区二区日本电影三级 | 美女午夜性视频免费| 人人妻人人澡欧美一区二区 | a在线观看视频网站| 久久久久久久精品吃奶| 国产精品自产拍在线观看55亚洲| 亚洲最大成人中文| 亚洲av成人一区二区三| 日韩大尺度精品在线看网址 | 男女午夜视频在线观看| 人人妻,人人澡人人爽秒播| 老熟妇仑乱视频hdxx| 悠悠久久av| 激情在线观看视频在线高清| 韩国av一区二区三区四区| 久久精品国产综合久久久| 亚洲 欧美一区二区三区| 法律面前人人平等表现在哪些方面| 久久久久精品国产欧美久久久| 午夜福利成人在线免费观看| 精品午夜福利视频在线观看一区| 亚洲一区高清亚洲精品| 亚洲成国产人片在线观看| 美女 人体艺术 gogo| 免费在线观看亚洲国产| 精品国产一区二区久久| 黄色片一级片一级黄色片| 18禁黄网站禁片午夜丰满| 免费在线观看黄色视频的| 亚洲五月婷婷丁香| 久久性视频一级片| 国产精华一区二区三区| 一边摸一边抽搐一进一小说| a级毛片在线看网站| 叶爱在线成人免费视频播放| 国产精品自产拍在线观看55亚洲| 亚洲视频免费观看视频| 丰满的人妻完整版| 午夜a级毛片| 欧美在线黄色| 精品人妻1区二区| 国产欧美日韩一区二区三区在线| 亚洲欧美日韩另类电影网站| 手机成人av网站| 女人爽到高潮嗷嗷叫在线视频| 亚洲中文日韩欧美视频| 亚洲欧美日韩高清在线视频| 国内久久婷婷六月综合欲色啪| 亚洲成国产人片在线观看| 免费在线观看黄色视频的| 欧美人与性动交α欧美精品济南到| 久久久久九九精品影院| 18禁裸乳无遮挡免费网站照片 | 人人妻人人澡欧美一区二区 | 男人的好看免费观看在线视频 | 欧美性长视频在线观看| 色播在线永久视频| 亚洲国产欧美网| 麻豆一二三区av精品| 中文字幕高清在线视频| 国产激情欧美一区二区| 啦啦啦韩国在线观看视频| 91大片在线观看| 高清在线国产一区| 免费一级毛片在线播放高清视频 | 波多野结衣高清无吗| 久久久久久免费高清国产稀缺| 亚洲成人免费电影在线观看| 中文亚洲av片在线观看爽| 一级黄色大片毛片| 亚洲五月色婷婷综合| 中文亚洲av片在线观看爽| 97人妻天天添夜夜摸| 老熟妇乱子伦视频在线观看| 无限看片的www在线观看| 精品高清国产在线一区| 99久久99久久久精品蜜桃| 精品卡一卡二卡四卡免费| 老司机深夜福利视频在线观看| 亚洲午夜精品一区,二区,三区| 少妇 在线观看| 欧美激情久久久久久爽电影 | 免费观看人在逋| 制服丝袜大香蕉在线| а√天堂www在线а√下载| 久久精品91无色码中文字幕| 在线观看午夜福利视频| 亚洲精品久久国产高清桃花| 欧美乱妇无乱码| 99香蕉大伊视频| 日本vs欧美在线观看视频| 最新在线观看一区二区三区| 亚洲免费av在线视频| 午夜老司机福利片| 欧美中文日本在线观看视频| 老司机福利观看| 成人亚洲精品av一区二区| 校园春色视频在线观看| 午夜成年电影在线免费观看| 最新美女视频免费是黄的| 久久精品国产清高在天天线| 18禁国产床啪视频网站| 亚洲熟妇中文字幕五十中出| 国产精品日韩av在线免费观看 | 免费观看精品视频网站| 亚洲精品av麻豆狂野| 亚洲一区中文字幕在线| 麻豆成人av在线观看| 国产一级毛片七仙女欲春2 | 亚洲欧美日韩另类电影网站| 久久精品91蜜桃| 亚洲少妇的诱惑av| 日本 av在线| 免费高清视频大片| 丝袜美腿诱惑在线| cao死你这个sao货| 色精品久久人妻99蜜桃| av福利片在线| e午夜精品久久久久久久| 国产人伦9x9x在线观看| 在线国产一区二区在线| 久久精品亚洲熟妇少妇任你| 少妇 在线观看| 亚洲全国av大片| av视频免费观看在线观看| 久久精品人人爽人人爽视色| 51午夜福利影视在线观看| 久久人妻熟女aⅴ| 99在线人妻在线中文字幕| 国产精品精品国产色婷婷| av视频在线观看入口| 亚洲,欧美精品.| 在线观看一区二区三区| 女同久久另类99精品国产91| 久久婷婷成人综合色麻豆| 亚洲全国av大片| 女人高潮潮喷娇喘18禁视频| 国产又爽黄色视频| 久热这里只有精品99| 国产三级黄色录像| 最近最新中文字幕大全免费视频| 色精品久久人妻99蜜桃| 欧美色视频一区免费| 非洲黑人性xxxx精品又粗又长| 一级片免费观看大全| 又黄又粗又硬又大视频| 亚洲精品一卡2卡三卡4卡5卡| 人妻丰满熟妇av一区二区三区| 久久人妻av系列| 国产av又大| 欧美丝袜亚洲另类 | 此物有八面人人有两片| 国产一区在线观看成人免费| 亚洲成av片中文字幕在线观看| 免费观看人在逋| 18禁裸乳无遮挡免费网站照片 | 一个人免费在线观看的高清视频| 免费高清在线观看日韩| 母亲3免费完整高清在线观看| 免费看十八禁软件| 亚洲一区高清亚洲精品| av天堂久久9| 黑人操中国人逼视频| 亚洲专区国产一区二区| 伊人久久大香线蕉亚洲五| av在线播放免费不卡| 久久中文看片网| 久久人妻熟女aⅴ| 国产av又大| 精品久久蜜臀av无| 99精品在免费线老司机午夜| 午夜福利影视在线免费观看| 在线观看66精品国产| 91字幕亚洲| 久久亚洲精品不卡| 亚洲狠狠婷婷综合久久图片| 又大又爽又粗| av片东京热男人的天堂| 正在播放国产对白刺激| 精品久久久精品久久久| 精品不卡国产一区二区三区| www.自偷自拍.com| 欧美激情 高清一区二区三区| 久久久水蜜桃国产精品网| 嫩草影院精品99| 久久亚洲真实| 老司机午夜十八禁免费视频| √禁漫天堂资源中文www| 久久 成人 亚洲| 欧美激情久久久久久爽电影 | 满18在线观看网站| 久久久国产成人免费| 国产激情久久老熟女| 国产精品电影一区二区三区| 黑人操中国人逼视频| 一边摸一边抽搐一进一出视频| 日日爽夜夜爽网站| 两人在一起打扑克的视频| 美女扒开内裤让男人捅视频| 操美女的视频在线观看| 午夜成年电影在线免费观看| 中文字幕最新亚洲高清| 亚洲第一电影网av| 99国产精品99久久久久| 国语自产精品视频在线第100页| 亚洲三区欧美一区| 看片在线看免费视频| 亚洲自偷自拍图片 自拍| 999精品在线视频| 日韩精品青青久久久久久| 日本黄色视频三级网站网址| 亚洲欧美日韩无卡精品| 欧美激情极品国产一区二区三区| 久久久久久久久久久久大奶| 免费看十八禁软件| 国产亚洲欧美98| 成人三级做爰电影| 日日摸夜夜添夜夜添小说| 国产成人免费无遮挡视频| 日韩精品中文字幕看吧| 国产日韩一区二区三区精品不卡| 亚洲国产精品久久男人天堂| 国产精品九九99| 免费看美女性在线毛片视频| 精品国产乱码久久久久久男人| 国产精品爽爽va在线观看网站 | 欧美中文综合在线视频| www.精华液| 久久人人爽av亚洲精品天堂| 看黄色毛片网站| 欧美成人午夜精品| 悠悠久久av| 国产野战对白在线观看| www.999成人在线观看| 国产精品av久久久久免费| 麻豆国产av国片精品| 国产人伦9x9x在线观看| 美女国产高潮福利片在线看| 村上凉子中文字幕在线| 久久精品人人爽人人爽视色| 午夜福利视频1000在线观看 | 国产视频一区二区在线看| 91成人精品电影| 国产精品日韩av在线免费观看 | 久久久久九九精品影院| 亚洲欧美精品综合久久99| svipshipincom国产片| 91九色精品人成在线观看| 国产精品1区2区在线观看.| 99久久久亚洲精品蜜臀av| 精品国产美女av久久久久小说| www.999成人在线观看| 最近最新免费中文字幕在线| x7x7x7水蜜桃| 国产av一区在线观看免费| 午夜免费鲁丝| 一级a爱视频在线免费观看| 欧美亚洲日本最大视频资源| 亚洲国产毛片av蜜桃av| 久久久久久久午夜电影| 日日干狠狠操夜夜爽| 久久久久久久精品吃奶| 色综合欧美亚洲国产小说| av电影中文网址| 亚洲熟女毛片儿| 日韩欧美一区视频在线观看| 国产精品日韩av在线免费观看 | 精品午夜福利视频在线观看一区| 国产单亲对白刺激| 久久热在线av| 老熟妇乱子伦视频在线观看| 不卡一级毛片| 国产1区2区3区精品| 纯流量卡能插随身wifi吗| 成人精品一区二区免费| 国产精品亚洲一级av第二区| 啦啦啦免费观看视频1| 免费观看人在逋| 热99re8久久精品国产| 亚洲在线自拍视频| 色综合婷婷激情| 999精品在线视频| 国产私拍福利视频在线观看| 黑丝袜美女国产一区| 狠狠狠狠99中文字幕| 18禁黄网站禁片午夜丰满| 无人区码免费观看不卡| 成熟少妇高潮喷水视频| 亚洲一区中文字幕在线| 国产精品久久久久久亚洲av鲁大| 久久国产精品男人的天堂亚洲| 国内毛片毛片毛片毛片毛片| 久久这里只有精品19| 无限看片的www在线观看| 亚洲欧美日韩无卡精品| 老汉色av国产亚洲站长工具| 国产欧美日韩精品亚洲av| 精品国产乱码久久久久久男人| 99久久精品国产亚洲精品| 最新美女视频免费是黄的| 亚洲第一欧美日韩一区二区三区| 久久久久国产精品人妻aⅴ院| 国产高清视频在线播放一区| 亚洲九九香蕉| 日本免费a在线| 中出人妻视频一区二区| 咕卡用的链子| 男人舔女人的私密视频| 亚洲激情在线av| 日韩av在线大香蕉| 国产激情久久老熟女| 精品久久久精品久久久| 亚洲精品国产区一区二| 午夜福利,免费看| 亚洲一码二码三码区别大吗| 大香蕉久久成人网| √禁漫天堂资源中文www| 亚洲中文日韩欧美视频| 日韩欧美国产一区二区入口| 久久中文看片网| a级毛片在线看网站| 日本三级黄在线观看| 欧美乱色亚洲激情| 免费在线观看影片大全网站| 岛国视频午夜一区免费看| 丝袜美腿诱惑在线| 大码成人一级视频| 可以在线观看毛片的网站| 亚洲精品av麻豆狂野| 国产成+人综合+亚洲专区| 色老头精品视频在线观看| 国产精品99久久99久久久不卡| 波多野结衣av一区二区av| 亚洲精品中文字幕一二三四区| 免费不卡黄色视频| 亚洲专区字幕在线| 精品电影一区二区在线| 他把我摸到了高潮在线观看| 男人的好看免费观看在线视频 | 好看av亚洲va欧美ⅴa在| 九色亚洲精品在线播放| 亚洲欧美日韩无卡精品| 啪啪无遮挡十八禁网站| 97超级碰碰碰精品色视频在线观看| 欧美最黄视频在线播放免费| 国产高清视频在线播放一区| www.999成人在线观看| 一区在线观看完整版| 久久久国产成人精品二区| 亚洲在线自拍视频| 欧美性长视频在线观看| 日韩精品免费视频一区二区三区| 国产精品精品国产色婷婷| 国产亚洲av嫩草精品影院| 日韩大尺度精品在线看网址 | 90打野战视频偷拍视频| 国产成人精品久久二区二区91| 久久久久国产一级毛片高清牌| 亚洲午夜精品一区,二区,三区| 97人妻精品一区二区三区麻豆 | 别揉我奶头~嗯~啊~动态视频| 搡老熟女国产l中国老女人| 国产亚洲精品久久久久久毛片| 男女下面进入的视频免费午夜 | 国产精品久久久久久亚洲av鲁大| 日韩欧美一区二区三区在线观看| 日日干狠狠操夜夜爽| 色尼玛亚洲综合影院| 国产精品98久久久久久宅男小说| 亚洲视频免费观看视频| 69av精品久久久久久| 色综合欧美亚洲国产小说| 女同久久另类99精品国产91| 国产成人系列免费观看| 日日干狠狠操夜夜爽| 成人亚洲精品av一区二区| 不卡av一区二区三区| 人人妻人人爽人人添夜夜欢视频| 最近最新中文字幕大全免费视频| 丰满人妻熟妇乱又伦精品不卡| 欧美激情久久久久久爽电影 | 禁无遮挡网站| 精品久久久久久久久久免费视频| 亚洲av成人不卡在线观看播放网| 亚洲熟妇中文字幕五十中出| 一级,二级,三级黄色视频| 我的亚洲天堂| 男人舔女人的私密视频| 久久精品国产亚洲av高清一级| 亚洲片人在线观看| 色av中文字幕| av电影中文网址| 午夜福利高清视频| 啦啦啦 在线观看视频| av福利片在线| 99国产精品一区二区三区| 精品一品国产午夜福利视频| 亚洲五月婷婷丁香| 久久久久久久久免费视频了| 国产1区2区3区精品| 91字幕亚洲| 久热这里只有精品99| 免费看美女性在线毛片视频| 久久香蕉激情| 久久精品亚洲熟妇少妇任你| 国产在线观看jvid| 一区二区三区精品91| 一级毛片高清免费大全| 国产亚洲精品久久久久久毛片| 在线观看免费视频网站a站| 精品欧美国产一区二区三| 国产av一区在线观看免费| 亚洲中文字幕一区二区三区有码在线看 | 久久精品aⅴ一区二区三区四区| 黄色视频,在线免费观看| 国产亚洲欧美98| 国产精品1区2区在线观看.| 久久精品aⅴ一区二区三区四区| 人成视频在线观看免费观看| 欧美一级a爱片免费观看看 | 亚洲国产精品999在线| 国产成人精品在线电影| 后天国语完整版免费观看| 成人国产综合亚洲| 91麻豆精品激情在线观看国产| bbb黄色大片| 国产免费男女视频| 亚洲男人天堂网一区| 18禁裸乳无遮挡免费网站照片 | 精品人妻1区二区| 国产又爽黄色视频| 麻豆av在线久日| 最近最新中文字幕大全电影3 | 两人在一起打扑克的视频| 久久精品国产综合久久久| 国产一级毛片七仙女欲春2 | 女人高潮潮喷娇喘18禁视频| 日韩欧美三级三区| 日韩高清综合在线| av电影中文网址| 一级a爱片免费观看的视频| 午夜亚洲福利在线播放| 中文字幕久久专区| 男人舔女人的私密视频| av超薄肉色丝袜交足视频| 国产高清videossex| 国产成人欧美| 亚洲国产欧美一区二区综合| 成人18禁高潮啪啪吃奶动态图| 97人妻天天添夜夜摸| 老汉色av国产亚洲站长工具| 亚洲成av人片免费观看| 国产人伦9x9x在线观看| 日韩精品青青久久久久久| 制服丝袜大香蕉在线| 国产99久久九九免费精品| 人人妻,人人澡人人爽秒播| 欧美日韩亚洲综合一区二区三区_| 黑人欧美特级aaaaaa片| 日韩免费av在线播放| 亚洲性夜色夜夜综合| 国产成人av激情在线播放| 精品电影一区二区在线| 国产色视频综合| 国产亚洲精品综合一区在线观看 | 纯流量卡能插随身wifi吗| videosex国产| 国内精品久久久久久久电影| 欧美日韩黄片免| 成人免费观看视频高清| 欧美大码av| 乱人伦中国视频| 一级毛片精品| 69精品国产乱码久久久| 热re99久久国产66热| 国产精品影院久久| 欧美乱妇无乱码| 一区二区三区国产精品乱码| 欧美午夜高清在线| 欧美在线一区亚洲| 美国免费a级毛片| 一a级毛片在线观看| 国产精品秋霞免费鲁丝片| 久久国产亚洲av麻豆专区| ponron亚洲| 亚洲国产精品成人综合色| 亚洲专区字幕在线| 午夜影院日韩av| 免费在线观看黄色视频的| 一区二区三区国产精品乱码| 宅男免费午夜| 在线观看午夜福利视频| 欧美国产精品va在线观看不卡| 涩涩av久久男人的天堂| 国产麻豆成人av免费视频| www日本在线高清视频| 欧美激情极品国产一区二区三区| 精品一区二区三区av网在线观看| 婷婷精品国产亚洲av在线| 在线观看免费午夜福利视频| 免费在线观看视频国产中文字幕亚洲| 久久精品影院6| 91麻豆av在线| 亚洲成人精品中文字幕电影| 男人的好看免费观看在线视频 | 精品欧美一区二区三区在线| 亚洲精品在线观看二区| 免费看a级黄色片| 黄片大片在线免费观看| 欧美色视频一区免费| 日日夜夜操网爽| 丁香六月欧美| 香蕉久久夜色| 无限看片的www在线观看| 精品一区二区三区视频在线观看免费| 此物有八面人人有两片| 天堂动漫精品| 久久久久国产一级毛片高清牌| 欧美精品亚洲一区二区| 少妇熟女aⅴ在线视频| 久久人妻熟女aⅴ| 多毛熟女@视频| av中文乱码字幕在线| 亚洲精品国产精品久久久不卡| 人妻久久中文字幕网| 黄网站色视频无遮挡免费观看| av超薄肉色丝袜交足视频| 成人精品一区二区免费| 亚洲熟妇熟女久久| 人妻久久中文字幕网| 欧美在线黄色| 桃红色精品国产亚洲av| 少妇 在线观看| 成人亚洲精品一区在线观看| 亚洲色图综合在线观看| 久久精品aⅴ一区二区三区四区| 亚洲av电影不卡..在线观看| 亚洲 欧美 日韩 在线 免费| 国产亚洲欧美在线一区二区| 国产精品电影一区二区三区| 亚洲成人精品中文字幕电影| 国产区一区二久久| 老司机靠b影院| 久久九九热精品免费| 国产免费男女视频| 国产高清激情床上av| 久久久国产成人免费| 日日夜夜操网爽| 一二三四社区在线视频社区8| 国产成人啪精品午夜网站| 一a级毛片在线观看| 非洲黑人性xxxx精品又粗又长| 久久人人97超碰香蕉20202| 视频在线观看一区二区三区| 亚洲色图av天堂| 亚洲精品国产精品久久久不卡| 久久久久久亚洲精品国产蜜桃av| 国产精品综合久久久久久久免费 | 国产三级黄色录像| 免费在线观看日本一区| 久久天堂一区二区三区四区| 天堂影院成人在线观看| 嫩草影视91久久| 亚洲欧美激情综合另类| 亚洲中文字幕一区二区三区有码在线看 | 亚洲色图 男人天堂 中文字幕| 热re99久久国产66热| 免费高清在线观看日韩| www.精华液| 午夜精品国产一区二区电影| 午夜精品在线福利| 精品少妇一区二区三区视频日本电影|