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

    破碎與團聚對冰漿存儲過程冰晶粒徑演化的影響

    2013-11-30 04:59:28徐愛祥劉志強王小倩趙騰磊
    中南大學學報(自然科學版) 2013年11期
    關(guān)鍵詞:概率密度冰晶粒徑

    徐愛祥,劉志強,王小倩,趙騰磊

    (中南大學 能源科學與工程學院,湖南 長沙,410083)

    冰漿(ice slurry)是指含有大量懸浮冰晶粒子的固液兩相溶液,其具有較高的儲能密度和良好的流動與換熱特性而受到人們廣泛關(guān)注[1?2],并且在眾多領(lǐng)域具有廣闊的應(yīng)用前景[3?4]。冰漿存儲期間冰晶的數(shù)量和粒徑存在著顯著的變化,這些變化會對冰漿的流動與傳熱特性產(chǎn)生重要影響[5]。Kozawa等[6?7]從宏觀角度研究了冰漿在存儲槽內(nèi)富集與分層現(xiàn)象,而這無法精確闡述冰漿存儲中冰晶演化過程與機理。Pronk等[8]從微觀角度出發(fā),認為存儲過程冰晶演化受磨損、團聚及Oswald熟化效應(yīng)共同作用,而熟化效應(yīng)的影響較為顯著,但計算結(jié)果與實驗存在較大差異。最近,Kumano[5]等分析發(fā)現(xiàn)團聚與破碎作用對冰晶粒徑演化也有重要影響。上述研究表明,存儲過程中冰晶演化的研究還不夠全面,特別是微尺度下冰晶演化過程與機理還有待深入研究。群體平衡模型(population balance model,PBM)是描述離散體系中分散相實體尺度分布的有效方法[9]。PBM在化工領(lǐng)域的廣泛應(yīng)用[10?11]為冰漿存儲過程離散相間相互作用的研究帶來很大啟示。冰漿系統(tǒng)是典型的離散系統(tǒng),其基本特點是載流溶液攜帶大量的冰晶顆粒群。冰漿存儲過程中冰晶生長與磨損、團聚與破碎等微觀現(xiàn)象都是群體平衡方程(population balance equation, PBE)所要表達的內(nèi)容。為此,本文作者針對冰漿存儲過程離散特點建立冰晶顆粒群體平衡模型,通過Matlab自編程序,采用分區(qū)法對群體平衡方程進行離散求解,著重討論破碎作用與團聚作用對冰晶粒徑演化的影響。

    1 冰晶群體平衡模型

    1.1 數(shù)學建模

    在建立冰漿存儲過程冰晶粒徑演化群體平衡數(shù)學模型時,進行如下假設(shè):(1) 冰晶粒子在空間上均勻分布,且具有各向同一特性;(2) 存儲過程冰晶的核化現(xiàn)象已不明顯,過冷度隨著冰晶的形成而逐漸變小,忽略冰晶的核化與生長;(3) 忽略外部空間對冰晶分布與演化的影響。因此,冰漿存儲過程中冰晶粒徑演化的零維空間群體平衡模型為

    式中:Bbr為由于破碎動力學行為造成冰晶數(shù)值密度的增加量,表示粒徑為x′的冰晶破碎成粒徑為x子顆粒的生產(chǎn)率;Dbr為由于破碎動力學行為造成冰晶數(shù)值密度的減小量,表示冰晶破碎成子顆粒后粒徑為x原顆粒的損失率;b(x)為破碎率(breakage frequency),表示單位時間粒徑為x的冰晶碰撞后發(fā)生破碎事件的速率;p(x|x′)為粒子概率密度函數(shù)(particle distribution function, PDF),表示粒徑為x′冰晶破碎后會形成粒徑為x冰晶的概率;v(x′)為粒徑為x′的冰晶破碎一次產(chǎn)生粒徑為x冰晶粒子的平均個數(shù),本文取v(x′)=2。

    冰晶的團聚項為

    式中:Bag為由于團聚動力學行為造成冰晶數(shù)值密度的增加量,表示粒徑為x?x′的冰晶與粒徑為x′冰晶團聚形成粒徑為x冰晶的生產(chǎn)率;Dag為由于團聚動力學行為造成冰晶數(shù)值密度的減小量,表示粒徑為x的冰晶在參與下一次的團聚事件后會產(chǎn)生更大粒徑的冰晶,從而造成粒徑為x原顆粒的損失率;a(x,x′)為團聚率(aggregation frequency),表示單位時間粒徑為x與x′冰晶碰撞后發(fā)生團聚事件的速率。需要注意的是:式(2)與式(3)積分項的上限∞指的是在計算時冰晶粒子粒徑的最大值。

    根據(jù)模型假設(shè),冰漿中冰晶粒子在慣性區(qū)域受流體擾動作用仍保持各向同性,故不同粒徑的冰晶碰撞時破碎率[12]的表達式如下:

    式中:ε為耗損率(dissipation rate),表示單位質(zhì)量的冰漿在外界擾動或內(nèi)部粒子間相互作用下所耗散的功率,W/kg;γ為冰晶與載流溶液間的表面張力,mN/m;φice為含冰率,%;ρ為載流溶液的密度,kg/m3;C1與C2為破碎率中常數(shù),其取值分別為0.004 81和0.08。

    粒徑為x′的冰晶破碎為粒徑x的冰晶概率分布函數(shù)[13]為

    式中:c為形狀因子,當冰晶為球形時,c可取為1。

    根據(jù)上述假設(shè),模型只考慮布朗碰撞作用下的團聚,因此,描述冰晶粒子在液相中的團聚率[12]為

    式中:KB為Boltzman常數(shù),J/K;T為冰漿存儲時的平衡溫度,K;μ為冰漿中載流溶液的動力黏度,mPa·s。

    模型中物性參數(shù)設(shè)置參照實驗條件,計算所需要的各物性參數(shù)如表1所示[12]。實驗中采用輕柔的攪拌方式以保持冰晶空間上的均勻分布,數(shù)值計算過程將其換算成耗損率。

    表1 載流溶液(9.2% NaCl)的物性參數(shù)Table 1 Physical parameters of ice slurry solutions(9.2% NaCl)

    1.2 求解方法

    群體平衡方程(PBE)是一個典型的混合積分與微分方程,其解呈非線性雙曲形式,求解過程需采取特殊的處理方式以封閉方程。PBE的求解方法主要有以下3種[14]:(1) 分區(qū)法;(2) 矩法;(3) Monte Carlo法。其中,分區(qū)法的離散方式思路直觀、簡潔, 且能夠得到顆粒尺度譜的動力學演變過程細節(jié)信息,計算精度與計算代價較適中。本文采用分區(qū)法對PBE進行離散求解:將冰晶尺度分布曲線進行離散,劃分為有限數(shù)目個區(qū)間,在每個區(qū)間內(nèi)針對某個冰晶粒徑的分布函數(shù),建立平衡方程,再聯(lián)立求解這些平衡方程。需要指出的是:求解過程中利用基于核函數(shù)的積分中值定理封閉PBE右邊的團聚與破碎項積分形式。上述求解過程通過Matlab自編程序來實現(xiàn)。

    1.3 邊界條件與初始條件

    1.3.1 邊界條件

    在 0~10?11m 范圍的冰晶粒徑概率密度在方程中無法表達與計算,而據(jù)成核理論可知[15],0~10?11m的冰晶粒子會有消融。因此,在任意時刻計算,可添加 1個與初始時刻相同的概率密度作為左邊界條件,即。

    1.3.2 初始條件

    考慮實驗中測定的初始狀態(tài)冰晶粒徑概率密度近似呈正態(tài)分布[16],在數(shù)值計算中,初始條件參照實驗結(jié)果中冰晶的初始分布形式設(shè)定,并通過實驗結(jié)果擬合而得到,擬合關(guān)聯(lián)式如下:

    其中:x的取值范圍為 10?11<x<vM+1。

    2 模型驗證

    為驗證本文建立的模型可靠性,將冰晶粒徑概率密度及冰晶平均粒徑的模擬結(jié)果與 Pronk的實驗結(jié)果[17]進行對比,如圖1和圖2所示。從圖1可以看出:冰晶粒徑實驗結(jié)果要比模擬結(jié)果稍偏高。從圖2可以看出:冰漿存儲22 h后,粒徑分布實驗結(jié)果與模擬結(jié)果在曲線谷峰有一定偏差,其余位置較吻合。實驗結(jié)果與模擬結(jié)果存在一定的偏差,其主要原因是:(1) 由于實驗過程中冰晶并非規(guī)則球狀,所獲取的冰晶圖片也非規(guī)則的圓形,因此,冰晶粒徑是經(jīng)過專門的粒徑圖像處理而得到,而模擬中冰晶形狀因子則按規(guī)則的球形來計算;(2) 冰漿實驗存儲過程可能涉及多種復雜微觀過程,而模擬計算忽略了冰晶生長與磨損,這也會給計算結(jié)果帶來偏差。通過定量比較,冰晶平均粒徑的實驗結(jié)果與模擬結(jié)果相對誤差在 5%以內(nèi),冰漿存儲22 h后冰晶粒徑概率密度基本吻合,說明本文模型準確、可靠。

    圖1 存儲期間冰晶平均粒徑的變化Fig.1 Variations of average ice crystal size during storage

    圖2 存儲22 h后粒徑概率密度的模擬值與實驗值Fig.2 Probability of ice crystal size distributions after 22 h storage obtained from model and experiment

    3 模擬結(jié)果與分析

    3.1 純破碎作用

    首先考察純破碎工況下冰漿存儲過程冰晶粒徑分布與演化。圖3所示為理想狀態(tài)下冰晶碰撞只發(fā)生破碎現(xiàn)象時的粒徑概率密度。由圖3可知:隨存儲時間的推移,冰晶分布曲線整體有向左移動的趨勢,也就是說,小粒徑冰晶出現(xiàn)的概率增大,大粒徑冰晶出現(xiàn)的概率減?。煌瑫r,峰值的左側(cè)更遠離初始分布,而峰值的右側(cè)更靠近初始分布,這表明破碎作用對小粒徑的冰晶影響更大,而對大粒徑的影響較小。這主要是由于冰漿在存儲過程中破碎作用使大粒徑冰晶減少及顆粒冰晶不斷增加。總體上,分布曲線移動幅度不大,且變化的幅度由快到慢,即在0~2 h內(nèi)冰晶粒徑有較小變化,2~22 h則無明顯變化。這表明在存儲過程中,破碎作用對冰晶的影響較小,且在短時間內(nèi),其作用有一定的影響;隨著存儲時間的增加,其作用逐漸削弱。

    圖3 冰晶在破碎作用下的粒徑概率密度Fig.3 Probability of ice crystal size distributions in time under breakage

    3.2 純團聚作用

    圖4 所示為理想狀態(tài)下冰晶碰撞只發(fā)生團聚現(xiàn)象時的粒徑概率密度。由圖4可知:與破碎作用下相反,隨著存儲時間的推移,團聚作用下冰晶分布曲線整體有向右移動,即小粒徑冰晶出現(xiàn)的概率減少,大粒徑冰晶出現(xiàn)概率增加;此外,峰值的左側(cè)與右側(cè)都遠離初始分布,這表明團聚作用對所有粒徑的冰晶都有影響。究其原因,主要是在存儲過程冰晶粒子群在團聚作用下,顆粒與顆粒之間不斷的黏附而導致冰晶整體呈現(xiàn)變大趨勢??傮w上來看,分布曲線移動明顯,且變化的幅度由慢到快,即在0~2 h內(nèi)冰晶粒徑有較小變化,2~22 h則變化明顯。這說明在存儲過程中,團聚作用對冰晶的影響是顯著的,且隨著存儲時間的增加,其作用也更加明顯。

    圖4 冰晶在團聚作用下的粒徑概率密度Fig.4 Probability of ice crystal size distributions in time under aggregation

    3.3 破碎與團聚共同作用

    在實際存儲過程中,破碎作用與團聚作用同時發(fā)生,共同作用,因此,很有必要分析破碎與團聚共同作用工況下冰漿存儲過程冰晶粒徑分布與演化。圖 5所示為冰晶在破碎與團聚共同作用下的粒徑概率密度,如圖5所示。從圖5可見:隨著冰漿存儲過程的推進,冰晶粒徑概率密度曲線向右移動且分布曲線坡度逐漸變緩,小粒徑冰晶出現(xiàn)的概率減小,大粒徑冰晶出現(xiàn)概率增大,隨時間推移冰晶粒徑分布范圍變廣。圖6所示為共同作用下冰晶粒徑分布函數(shù)。從圖6可以看出:初始時刻冰晶粒徑主要分布在100~500 μm,隨著存儲時間的增加,大粒徑冰晶比例明顯增加。對比純破碎作用、純團聚作用及破碎與團聚共同作用對粒徑分布與演化的影響可看出:冰漿存儲經(jīng)歷了破碎與團聚這2個不同作用的微觀過程,2個“作用方向”相反,破碎作用下大粒徑冰晶產(chǎn)生小粒徑冰晶,而團聚作用下小粒徑冰晶形成大粒徑冰晶;此外,在存儲初期,破碎與團聚單獨作用下(見圖3與圖4)粒徑概率密度曲線的變化比兩者共同作用下(見圖 5)粒徑概率密度曲線的變化大,表明在共同作用時破碎與團聚對冰晶粒徑的影響相互抵消;隨著存儲時間的增加,破碎作用不斷削弱,而團聚作用依然明顯,因此,二者的相互作用由最初的破碎與團聚“均衡作用”發(fā)展成團聚“主導”。這一現(xiàn)象的具體表現(xiàn)是冰晶在存儲期間平均粒徑會增大,圖1證明了這一結(jié)論,冰晶的平均粒徑從初始時刻的 300 μm 增長到存儲 22 h后的392 μm。

    圖5 冰晶在破碎與團聚共同作用下粒徑概率密度Fig.5 Probability of ice crystal size distributions in time under breakage and aggregation

    圖6 冰晶在破碎與團聚共同作用下粒徑分布函數(shù)Fig.6 Distribution function of ice crystal size distributions in time under breakage and aggregation

    4 結(jié)論

    (1) 引入群體平衡模型可較好地模擬冰漿存儲過程的粒徑分布及演化,獲得了冰晶粒徑概率密度隨時間的變化關(guān)系:純破碎作用下小粒徑冰晶出現(xiàn)的概率增大,大粒徑冰晶出現(xiàn)概率減小,純團聚及兩者共同作用下小粒徑冰晶出現(xiàn)的概率減少,大粒徑冰晶出現(xiàn)概率增加。

    (2) 破碎與團聚是冰晶粒徑演化的2個重要作用,在存儲過程中,兩者“作用方向”相反,且破碎與團聚經(jīng)歷從最初的破碎與團聚“均衡作用”發(fā)展成團聚“主導”的過程。

    (3) 在冰晶存儲過程中,冰晶平均粒徑不斷長大,即從初始時刻的300 μm增長到存儲22 h后的392 μm,且模擬值與實驗值相對誤差在5%以內(nèi)。

    [1]Egolf P W, Kauffeld M. From physical properties of ice slurries to industrial ice slurry applications[J]. International Journal of Refrigeration, 2005, 28(1): 4?12.

    [2]Ayel V, Lottin O, Peerhossaini H. Rheology, flow behaviour and heat transfer of ice slurries: a review of the state of the art[J].International Journal of Refrigeration, 2003, 26(1): 95?107.

    [3]Kauffeld M, Wang M J, Goldstein V, et al. Ice slurry applications[J]. International Journal of Refrigeration, 2010,33(8): 1491?1505.

    [4]Bellas I, Tassou S A. Present and future applications of ice slurries[J]. International Journal of Refrigeration, 2005, 28(1):115?121.

    [5]Kumano H, Hirata T, Hagiwara Y, et al. Effects of storage on flow and heat transfer characteristics of ice slurry[J].International Journal of Refrigeration, 2012, 35(1): 122?129.

    [6]Kozawa Y, Aizawa N, Tanino M. Study on ice storing characteristics in dynamic-type ice storage system by using supercooled water: Effects of the supplying conditions of ice-slurry at deployment to district heating and cooling system[J].International Journal of Refrigeration, 2005, 28(1): 73?82.

    [7]青春耀, 肖睿, 宋文吉, 等. 冰漿在蓄冰槽內(nèi)的蓄冰特性及其均勻度研究[J]. 低溫與超導, 2009, 37(5): 41?46.QING Chunyao, XIAO Rui, SONG Wenji, et al. Study on ice slurry storing characteristics and its uniformity in ice storage tank[J]. Cryogenics and Superconductivity, 2009, 37(5): 41?46.

    [8]Pronk P, Hansen T M, Ferreira C A I, et al. Time-dependent behavior of different ice slurries during storage[J]. International Journal of Refrigeration, 2005, 28(1): 27?36.

    [9]Ramkrishna D. Population balances: Theory and applications to particulate systems in engineering[M]. San Diego: Academic Press, 2000: 30?45.

    [10]顧兆林, 蘇軍偉, 李云, 等. 兩相及多相體系的離散相行為與群體平衡模型[J]. 化學反應(yīng)工程與工藝, 2007, 23(2): 162?167.GU Zhaolin, SU Junwei, LI Yun, et al. Behaviors of the dispersed phase in the multiphase system and population balance model[J]. Chemical Reaction Engineering and Technology, 2007,23(2): 162?167.

    [11]Gerstlauer A, Motz S, Mitrovic A, et al. Development, analysis and validation of population models for continuous and batch crystallizers[J]. Chemical Engineering Science, 2002, 57(20):4311?4327.

    [12]Marchisio D L, Vigil R D, Fox R O. Implementation of the quadrature method of moments in CFD codes for aggregation-breakage problems[J]. Chemical Engineering Science, 2003, 58(15): 3337?3351.

    [13]Hill P J, Ng K M. New discretization procedure for the breakage equation[J]. AIChE Journal, 1995, 41(5): 1204?1216.

    [14]趙海波, 鄭楚光. 離散系統(tǒng)動力學演變過程的顆粒群平衡模擬[M]. 北京: 科學出版社, 2008: 70?79.ZHAO Haibo, ZHENG Chuguang. Population balance modeling for the process of the dynamic evolution in dispersed systems[M]. Beijing: Science Press, 2008: 70?79.

    [15]Chégnimonhan V, Josset C, Peerhossaini H. Ice slurry crystallization based on kinetic phase-change modeling[J].International Journal of Refrigeration, 2010, 33(8): 1559?1568.

    [16]Hansen T M, Radosevic M, Kauffeld M. Behaviour of ice slurry in thermal storage systems[R]. Atlanta: ASHRAE, 2002:93?107.

    [17]Pronk P, Ferreira C A I, Witkamp G J. Effect of long-term ice slurry storage on crystal size distribution[C]// Melinder A.5th Workshop on Ice Slurries of the IIR. Stockholm: International Institute of Refrigeration, 2002: 151?160.

    猜你喜歡
    概率密度冰晶粒徑
    冰晶世界,凝凍微觀的美麗
    木屑粒徑對黑木耳栽培的影響試驗*
    雪為什么只在冬天下
    連續(xù)型隨機變量函數(shù)的概率密度公式
    小泥人冰晶畫 蘊藏大商機
    基于近場散射的顆粒粒徑分布測量
    冰晶奇域
    琴童(2016年5期)2016-05-14 15:43:58
    Hunt過程在Girsanov變換下的轉(zhuǎn)移概率密度的表示公式
    隨機變量線性組合的分布的一個算法
    Oslo結(jié)晶器晶體粒徑分布特征的CFD模擬
    国产爽快片一区二区三区| 最近中文字幕高清免费大全6| 成人手机av| 亚洲美女黄色视频免费看| 两个人免费观看高清视频| 十分钟在线观看高清视频www| 国产色爽女视频免费观看| 免费黄色在线免费观看| 国产视频内射| 国产色爽女视频免费观看| 免费看光身美女| 天天影视国产精品| 亚洲精品日韩av片在线观看| 日韩av在线免费看完整版不卡| freevideosex欧美| 国产精品 国内视频| 久久久久久久久久久免费av| 人体艺术视频欧美日本| 亚洲精品日韩av片在线观看| 26uuu在线亚洲综合色| 欧美bdsm另类| 一本色道久久久久久精品综合| 亚洲av电影在线观看一区二区三区| 久久午夜福利片| 亚洲成人一二三区av| 热99久久久久精品小说推荐| 午夜福利在线观看免费完整高清在| 国产精品成人在线| 人妻夜夜爽99麻豆av| 精品亚洲乱码少妇综合久久| 国产色婷婷99| 成年人免费黄色播放视频| 亚洲精品久久久久久婷婷小说| 最新中文字幕久久久久| 国产白丝娇喘喷水9色精品| 国产在线视频一区二区| 性高湖久久久久久久久免费观看| 97超碰精品成人国产| 国产av一区二区精品久久| 日韩在线高清观看一区二区三区| 亚洲性久久影院| 国产白丝娇喘喷水9色精品| 老司机影院毛片| 考比视频在线观看| 午夜福利,免费看| 人妻制服诱惑在线中文字幕| 国产成人精品福利久久| 欧美日韩综合久久久久久| 在线观看免费日韩欧美大片 | 成人毛片60女人毛片免费| 国产亚洲最大av| 国产精品人妻久久久影院| 国产精品三级大全| 各种免费的搞黄视频| 久久久久视频综合| 一边亲一边摸免费视频| 美女xxoo啪啪120秒动态图| 欧美性感艳星| 久久国产亚洲av麻豆专区| 能在线免费看毛片的网站| 2021少妇久久久久久久久久久| 国产老妇伦熟女老妇高清| 国产成人aa在线观看| 狠狠精品人妻久久久久久综合| 不卡视频在线观看欧美| 91aial.com中文字幕在线观看| 国产成人av激情在线播放 | 日本黄色日本黄色录像| 日日撸夜夜添| 亚洲精品日本国产第一区| 免费少妇av软件| 亚洲欧美精品自产自拍| 国产片内射在线| 97精品久久久久久久久久精品| 久久久精品免费免费高清| 久久热精品热| 日本爱情动作片www.在线观看| 亚洲av综合色区一区| 91精品三级在线观看| 国产精品一国产av| 王馨瑶露胸无遮挡在线观看| videos熟女内射| 人体艺术视频欧美日本| a级毛片免费高清观看在线播放| 少妇人妻精品综合一区二区| 国产成人一区二区在线| 成人黄色视频免费在线看| 久久青草综合色| 亚洲av不卡在线观看| 精品少妇黑人巨大在线播放| 极品人妻少妇av视频| 嫩草影院入口| 亚洲三级黄色毛片| 欧美日韩成人在线一区二区| 一二三四中文在线观看免费高清| 日本黄色片子视频| 少妇的逼水好多| a级毛片在线看网站| 婷婷色综合大香蕉| 成人无遮挡网站| 国产亚洲精品第一综合不卡 | 桃花免费在线播放| 水蜜桃什么品种好| 国产探花极品一区二区| 久久精品国产自在天天线| 美女福利国产在线| av在线老鸭窝| 国产精品 国内视频| 亚洲欧洲精品一区二区精品久久久 | 91精品国产九色| 青春草视频在线免费观看| 国产免费一区二区三区四区乱码| 夫妻午夜视频| a级毛片免费高清观看在线播放| 国产成人a∨麻豆精品| 97在线人人人人妻| 人妻夜夜爽99麻豆av| 丰满乱子伦码专区| 中文字幕亚洲精品专区| 成年美女黄网站色视频大全免费 | 丁香六月天网| 亚洲精品,欧美精品| 五月天丁香电影| 亚洲第一av免费看| 一区二区三区精品91| 成人漫画全彩无遮挡| 久热久热在线精品观看| 国产精品久久久久久久电影| 久久午夜综合久久蜜桃| 中文字幕人妻熟人妻熟丝袜美| 夜夜看夜夜爽夜夜摸| 大片电影免费在线观看免费| 嘟嘟电影网在线观看| 女人久久www免费人成看片| 美女主播在线视频| 高清欧美精品videossex| 人妻制服诱惑在线中文字幕| 视频区图区小说| 国语对白做爰xxxⅹ性视频网站| 日韩视频在线欧美| 国产熟女午夜一区二区三区 | 肉色欧美久久久久久久蜜桃| 国产欧美日韩综合在线一区二区| 亚洲精品,欧美精品| 在线精品无人区一区二区三| 99精国产麻豆久久婷婷| 91午夜精品亚洲一区二区三区| xxx大片免费视频| 99九九在线精品视频| 一本一本综合久久| 欧美成人午夜免费资源| 国产精品麻豆人妻色哟哟久久| 两个人免费观看高清视频| 亚洲精品一二三| 亚洲欧美精品自产自拍| 亚洲色图 男人天堂 中文字幕 | 久久国产精品男人的天堂亚洲 | 成人午夜精彩视频在线观看| 亚洲中文av在线| 国产黄频视频在线观看| 成人国产av品久久久| 日韩不卡一区二区三区视频在线| 亚洲av福利一区| 亚洲精品av麻豆狂野| 51国产日韩欧美| 十八禁高潮呻吟视频| 国产精品不卡视频一区二区| 免费观看在线日韩| 日韩中字成人| 欧美国产精品一级二级三级| 91精品伊人久久大香线蕉| 欧美日韩视频高清一区二区三区二| 91久久精品国产一区二区成人| 亚洲国产精品一区三区| 在线观看免费高清a一片| 丰满饥渴人妻一区二区三| 久热久热在线精品观看| 国产日韩欧美亚洲二区| 日日摸夜夜添夜夜添av毛片| 国产免费福利视频在线观看| 亚洲美女视频黄频| 观看美女的网站| 老司机影院成人| 高清毛片免费看| 尾随美女入室| 99视频精品全部免费 在线| 欧美人与性动交α欧美精品济南到 | 亚洲欧洲日产国产| 亚洲精品日本国产第一区| 涩涩av久久男人的天堂| 九草在线视频观看| 免费高清在线观看视频在线观看| 最新中文字幕久久久久| 精品一区在线观看国产| 久久国内精品自在自线图片| 亚洲精品国产av蜜桃| 亚洲第一av免费看| 精品亚洲成a人片在线观看| 赤兔流量卡办理| 国产精品.久久久| 只有这里有精品99| 国产精品一区二区在线不卡| 极品少妇高潮喷水抽搐| av播播在线观看一区| 成人影院久久| 亚洲精品第二区| 精品亚洲乱码少妇综合久久| 国产精品国产av在线观看| 18+在线观看网站| 黄色配什么色好看| 人妻夜夜爽99麻豆av| 亚洲精品日韩av片在线观看| 亚洲情色 制服丝袜| 久久久久久人妻| 少妇高潮的动态图| 这个男人来自地球电影免费观看 | 久久久精品免费免费高清| kizo精华| 成年美女黄网站色视频大全免费 | 免费黄网站久久成人精品| 男女免费视频国产| 国产黄片视频在线免费观看| 欧美亚洲日本最大视频资源| 80岁老熟妇乱子伦牲交| 亚洲国产日韩一区二区| 欧美三级亚洲精品| 最新的欧美精品一区二区| 美女主播在线视频| 五月玫瑰六月丁香| 欧美日韩国产mv在线观看视频| 久久国产精品男人的天堂亚洲 | 多毛熟女@视频| 欧美亚洲日本最大视频资源| 一个人看视频在线观看www免费| 高清黄色对白视频在线免费看| 亚洲国产精品一区二区三区在线| 美女主播在线视频| 有码 亚洲区| 久久人人爽av亚洲精品天堂| 亚洲综合精品二区| 国产欧美日韩一区二区三区在线 | 大码成人一级视频| 水蜜桃什么品种好| 中国美白少妇内射xxxbb| 国产免费福利视频在线观看| 一级毛片 在线播放| 丰满乱子伦码专区| av不卡在线播放| 999精品在线视频| 80岁老熟妇乱子伦牲交| 亚洲欧美精品自产自拍| 啦啦啦在线观看免费高清www| 超碰97精品在线观看| 性高湖久久久久久久久免费观看| 国产精品无大码| 午夜激情久久久久久久| 久久午夜福利片| 欧美精品亚洲一区二区| 成人二区视频| 国产又色又爽无遮挡免| 99国产综合亚洲精品| 又大又黄又爽视频免费| 亚洲五月色婷婷综合| 熟女电影av网| 九色成人免费人妻av| 女性被躁到高潮视频| 飞空精品影院首页| 中文字幕久久专区| 成人午夜精彩视频在线观看| 黄色一级大片看看| 亚洲精品,欧美精品| 国产精品一区二区在线不卡| 三级国产精品片| 高清午夜精品一区二区三区| 免费观看无遮挡的男女| 欧美人与善性xxx| √禁漫天堂资源中文www| 国产深夜福利视频在线观看| 哪个播放器可以免费观看大片| 亚洲欧美中文字幕日韩二区| 亚洲av日韩在线播放| 久久 成人 亚洲| 欧美激情 高清一区二区三区| 大片电影免费在线观看免费| 亚洲第一av免费看| 久久久午夜欧美精品| 视频中文字幕在线观看| 十八禁高潮呻吟视频| 亚洲人成77777在线视频| 汤姆久久久久久久影院中文字幕| 九色成人免费人妻av| 能在线免费看毛片的网站| 国产免费又黄又爽又色| 国产成人精品无人区| 最近最新中文字幕免费大全7| 亚洲成人手机| 美女中出高潮动态图| 久久精品久久久久久噜噜老黄| 少妇 在线观看| 亚洲第一区二区三区不卡| 久久 成人 亚洲| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久97久久精品| 欧美丝袜亚洲另类| 乱码一卡2卡4卡精品| 岛国毛片在线播放| 久久99蜜桃精品久久| 午夜免费鲁丝| 黑人巨大精品欧美一区二区蜜桃 | 97超碰精品成人国产| 伊人亚洲综合成人网| 国产精品一区二区在线不卡| 特大巨黑吊av在线直播| 人人妻人人澡人人爽人人夜夜| 精品久久久久久电影网| 卡戴珊不雅视频在线播放| 国产精品久久久久久久久免| 人妻一区二区av| 日本-黄色视频高清免费观看| 七月丁香在线播放| av电影中文网址| 免费少妇av软件| 熟女电影av网| 亚洲精品国产色婷婷电影| 国产精品蜜桃在线观看| 免费高清在线观看日韩| 欧美精品国产亚洲| 国产乱来视频区| 欧美人与善性xxx| 视频在线观看一区二区三区| 国产亚洲最大av| 国产乱来视频区| 伊人久久精品亚洲午夜| 人人妻人人澡人人看| 热99国产精品久久久久久7| .国产精品久久| 91精品国产国语对白视频| 热99久久久久精品小说推荐| 日韩不卡一区二区三区视频在线| 美女cb高潮喷水在线观看| 在线观看一区二区三区激情| 国产男女超爽视频在线观看| 一区在线观看完整版| 日本爱情动作片www.在线观看| 欧美日韩亚洲高清精品| 九草在线视频观看| 亚洲一级一片aⅴ在线观看| 欧美日韩av久久| 老女人水多毛片| 嘟嘟电影网在线观看| 亚洲中文av在线| 日韩中文字幕视频在线看片| 国产精品久久久久久精品电影小说| 亚洲国产精品999| 成人影院久久| 亚洲久久久国产精品| 在线观看国产h片| 欧美老熟妇乱子伦牲交| 免费大片18禁| 欧美激情极品国产一区二区三区 | 2021少妇久久久久久久久久久| 高清午夜精品一区二区三区| 精品久久久精品久久久| 国产又色又爽无遮挡免| 精品一区在线观看国产| 亚洲精华国产精华液的使用体验| 人人澡人人妻人| 国产伦理片在线播放av一区| 最近中文字幕高清免费大全6| 十八禁网站网址无遮挡| 少妇猛男粗大的猛烈进出视频| 国产精品一区二区三区四区免费观看| 国产成人免费无遮挡视频| 伦理电影大哥的女人| 免费观看性生交大片5| 日韩,欧美,国产一区二区三区| 天堂8中文在线网| 久久精品国产亚洲av涩爱| 久久99热这里只频精品6学生| 日韩,欧美,国产一区二区三区| 亚洲色图综合在线观看| 国产男女超爽视频在线观看| 2021少妇久久久久久久久久久| 亚洲国产精品一区二区三区在线| 国产黄色视频一区二区在线观看| 在线观看免费日韩欧美大片 | 我的女老师完整版在线观看| 久久国产精品男人的天堂亚洲 | 欧美精品亚洲一区二区| 大香蕉久久成人网| 美女视频免费永久观看网站| 草草在线视频免费看| 十八禁高潮呻吟视频| 日韩大片免费观看网站| 成年女人在线观看亚洲视频| 久久久精品94久久精品| 久久亚洲国产成人精品v| 久久综合国产亚洲精品| 国产成人精品福利久久| 插阴视频在线观看视频| 一级毛片黄色毛片免费观看视频| 亚洲av免费高清在线观看| 午夜福利网站1000一区二区三区| 亚洲精品第二区| 国产黄频视频在线观看| 亚洲人成网站在线播| 视频区图区小说| 极品少妇高潮喷水抽搐| 一边亲一边摸免费视频| 欧美变态另类bdsm刘玥| 亚洲精品久久午夜乱码| 久久精品人人爽人人爽视色| 欧美 亚洲 国产 日韩一| 制服丝袜香蕉在线| 国产精品一区www在线观看| 这个男人来自地球电影免费观看 | 免费不卡的大黄色大毛片视频在线观看| 天天影视国产精品| 狂野欧美激情性xxxx在线观看| 麻豆成人av视频| 黄色怎么调成土黄色| 看免费成人av毛片| 午夜福利影视在线免费观看| 国产永久视频网站| 妹子高潮喷水视频| 国产极品粉嫩免费观看在线 | 精品国产乱码久久久久久小说| 丝瓜视频免费看黄片| 日韩大片免费观看网站| 午夜视频国产福利| 观看美女的网站| 男人操女人黄网站| 欧美3d第一页| 亚洲伊人久久精品综合| 日韩视频在线欧美| 欧美性感艳星| 少妇高潮的动态图| 日本欧美国产在线视频| 国产精品国产三级专区第一集| 丰满饥渴人妻一区二区三| 最黄视频免费看| 亚洲图色成人| 免费久久久久久久精品成人欧美视频 | 纯流量卡能插随身wifi吗| 亚洲综合色惰| 成年女人在线观看亚洲视频| 在线观看美女被高潮喷水网站| 国产毛片在线视频| 少妇猛男粗大的猛烈进出视频| 熟女电影av网| 午夜激情福利司机影院| 国产欧美日韩综合在线一区二区| tube8黄色片| 有码 亚洲区| 伦精品一区二区三区| 热re99久久国产66热| 国产成人精品久久久久久| 水蜜桃什么品种好| 3wmmmm亚洲av在线观看| 精品亚洲成国产av| 亚洲欧洲精品一区二区精品久久久 | 国产 一区精品| 高清毛片免费看| 国产男人的电影天堂91| 国产成人一区二区在线| 日韩欧美精品免费久久| 精品国产一区二区三区久久久樱花| 亚洲精品国产av成人精品| a级毛片在线看网站| 欧美丝袜亚洲另类| 七月丁香在线播放| 国产又色又爽无遮挡免| 青春草亚洲视频在线观看| 久久这里有精品视频免费| 久久久午夜欧美精品| 一本—道久久a久久精品蜜桃钙片| 欧美亚洲日本最大视频资源| 成人手机av| 欧美3d第一页| 如日韩欧美国产精品一区二区三区 | 国产亚洲最大av| 精品一品国产午夜福利视频| 五月玫瑰六月丁香| 午夜福利视频在线观看免费| 日韩精品有码人妻一区| 欧美 日韩 精品 国产| 欧美日韩在线观看h| 妹子高潮喷水视频| 街头女战士在线观看网站| 国产亚洲一区二区精品| 亚洲av男天堂| 午夜视频国产福利| 毛片一级片免费看久久久久| 伊人亚洲综合成人网| 91精品国产九色| 亚洲精品色激情综合| 久久亚洲国产成人精品v| 18禁在线无遮挡免费观看视频| 纵有疾风起免费观看全集完整版| 欧美日韩av久久| 男的添女的下面高潮视频| 秋霞在线观看毛片| 亚洲欧美成人综合另类久久久| 老熟女久久久| 亚洲国产精品专区欧美| 免费黄频网站在线观看国产| 99九九线精品视频在线观看视频| 99视频精品全部免费 在线| 亚洲五月色婷婷综合| 亚洲欧洲国产日韩| 久久久精品免费免费高清| 国产av国产精品国产| 亚洲av日韩在线播放| 亚洲av综合色区一区| 中国三级夫妇交换| 人妻一区二区av| 午夜福利影视在线免费观看| 国产黄频视频在线观看| 免费黄频网站在线观看国产| 青青草视频在线视频观看| 熟女电影av网| 国产69精品久久久久777片| 日韩大片免费观看网站| 亚洲,欧美,日韩| 性色avwww在线观看| 超色免费av| 涩涩av久久男人的天堂| 亚洲av电影在线观看一区二区三区| 亚洲成人av在线免费| 最近2019中文字幕mv第一页| 亚洲国产av新网站| 97超碰精品成人国产| 成人无遮挡网站| 国产极品天堂在线| 久久久国产欧美日韩av| 视频中文字幕在线观看| 亚洲av日韩在线播放| 欧美成人精品欧美一级黄| 一级片'在线观看视频| 岛国毛片在线播放| 国产黄色视频一区二区在线观看| 九九久久精品国产亚洲av麻豆| 国产精品国产av在线观看| 18禁在线无遮挡免费观看视频| 日韩av在线免费看完整版不卡| 亚洲图色成人| 中国国产av一级| 你懂的网址亚洲精品在线观看| 韩国高清视频一区二区三区| 日韩三级伦理在线观看| 伊人久久精品亚洲午夜| 曰老女人黄片| 搡女人真爽免费视频火全软件| 久久毛片免费看一区二区三区| videossex国产| 大香蕉久久成人网| 天天影视国产精品| 日日摸夜夜添夜夜爱| 久久免费观看电影| 老熟女久久久| 成人漫画全彩无遮挡| av.在线天堂| 老熟女久久久| 中国国产av一级| 亚洲第一区二区三区不卡| 香蕉精品网在线| 午夜久久久在线观看| 精品酒店卫生间| 亚洲欧美日韩另类电影网站| 成人国产av品久久久| 在线观看美女被高潮喷水网站| 国产免费福利视频在线观看| 亚洲国产成人一精品久久久| 精品久久久久久久久av| 老司机影院毛片| 妹子高潮喷水视频| 日韩av免费高清视频| 亚洲中文av在线| 亚洲精品成人av观看孕妇| 免费不卡的大黄色大毛片视频在线观看| 国产成人freesex在线| 水蜜桃什么品种好| 七月丁香在线播放| 亚洲精品美女久久av网站| 欧美精品亚洲一区二区| 七月丁香在线播放| 国产一区二区三区综合在线观看 | 两个人免费观看高清视频| 亚洲国产最新在线播放| 大香蕉久久成人网| 男人添女人高潮全过程视频| 美女内射精品一级片tv| 一区二区三区精品91| 欧美97在线视频| 好男人视频免费观看在线| 亚洲美女黄色视频免费看| 成年女人在线观看亚洲视频| av国产久精品久网站免费入址| 亚洲美女黄色视频免费看| h视频一区二区三区| 亚洲婷婷狠狠爱综合网| 如日韩欧美国产精品一区二区三区 | 边亲边吃奶的免费视频| 十八禁高潮呻吟视频| 久久99热这里只频精品6学生| 久久久亚洲精品成人影院| 精品亚洲成国产av| 国精品久久久久久国模美| 2022亚洲国产成人精品| 又黄又爽又刺激的免费视频.| 91在线精品国自产拍蜜月| freevideosex欧美| 欧美97在线视频| 最后的刺客免费高清国语| 国产一区二区在线观看日韩|