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

    柵徑比對環(huán)形燃料外周向溫度分布不均勻性影響

    2023-09-21 09:21:22曾付林趙鵬程向釗才
    核技術(shù) 2023年9期
    關(guān)鍵詞:周向六邊形正方形

    曾付林 唐 錨 趙鵬程 向釗才

    (南華大學(xué) 核科學(xué)技術(shù)學(xué)院 衡陽 421001)

    環(huán)形燃料采用內(nèi)、外冷卻流道增大了冷卻能力,相較于傳統(tǒng)棒狀燃料元件,降低了燃料芯塊的最大溫度,可以在保持或改進(jìn)現(xiàn)有反應(yīng)堆安全裕量的前提下,提升反應(yīng)堆功率密度20%~50%,具有優(yōu)良的安全性與經(jīng)濟(jì)性,已成為壓水堆先進(jìn)燃料元件的重要發(fā)展趨勢之一[1-2]。其自從MIT研究者提出以來,中國、美國、韓國以及伊朗等國家均已對其做了大量的熱工水力性能分析[2-5],且環(huán)形燃料棒束直徑明顯大于傳統(tǒng)棒狀燃料直徑,因此更加地穩(wěn)定,不容易產(chǎn)生振動,只要保持合理的壓降即可將棒束間間隙降低到1 mm[6]。中國原子能科學(xué)研究院初步明確了環(huán)形燃料組件概念設(shè)計(jì)以AFA-3G 組件為參考,在保持原有組件幾何尺寸和“水鈾比”不變的情況下,綜合考慮堆芯物理、熱工水力以及元件性能等因素,采用13×13的排列結(jié)構(gòu)[2]。因此,目前大多數(shù)已有試驗(yàn)將環(huán)形燃料柵徑比設(shè)置為1.07,且對環(huán)形燃料外徑尺寸的設(shè)置一般在15.3~15.6 mm,通過柵徑比、柵距、燃料外徑以及棒束間間隙寬度之間的關(guān)系可計(jì)算得棒束間間隙寬度僅有1.1~1.2 mm。這種緊密柵組件與柵徑比較大的燃料組件相比,其直觀的特征是燃料棒之間的間距較小,導(dǎo)致流動與換熱沿燃料棒周向存在明顯的不均勻性,使得包殼溫度沿燃料棒周向分布不均,進(jìn)而容易導(dǎo)致包殼在熱應(yīng)力下發(fā)生形變甚至破損,對反應(yīng)堆的安全性存在隱患。在燃料組件內(nèi),換熱沿周向存在不均勻性,且不均勻性隨著柵徑比的減小而急劇增大[7],通常認(rèn)為柵徑比大于1.3 時(shí)這種不均勻性可以忽略,因此,對于大多應(yīng)用于沸水堆及壓水堆的傳統(tǒng)棒狀燃料無須考慮周向溫度的分布,而對于緊密排列的環(huán)形燃料而言,其周向不同角度處的溫度可能差異很大,并且中心通道、邊通道、角通道的換熱規(guī)律并不相同,需要分別進(jìn)行考慮。

    Cheng等[8]對正方形、六邊形組件典型子通道熱工水力行為進(jìn)行模擬,發(fā)現(xiàn)緊密柵棒束中,尤其是正方形組件,存在強(qiáng)烈的周向傳熱不均勻性。Subbotin 等[9]對緊密排列的七棒束模型進(jìn)行實(shí)驗(yàn)研究,也觀測到了溫度沿周向分布的不均勻性,并提出了包殼表面溫度關(guān)于角度的半經(jīng)驗(yàn)公式。Krauss等[10-11]對緊密柵內(nèi)流動和傳熱的物理量進(jìn)行測算,發(fā)現(xiàn)了燃料棒間隙處流體流速和溫度存在周期性波動。楊婷[12]對緊密柵燃料組件進(jìn)行分析,提出影響溫度分布不均的主要無量綱數(shù)是柵徑比,并提出確定描述傳熱沿周向不均勻性的經(jīng)驗(yàn)關(guān)系式。魏君翰等[13]對環(huán)形燃料冷卻劑流動特性進(jìn)行數(shù)值模擬時(shí)發(fā)現(xiàn)兩環(huán)形燃料棒間隙處存在局部高溫區(qū)。以上多數(shù)研究均基于棒狀燃料的緊密柵組件,針對雙面冷卻的環(huán)形燃料組件其周向溫度分布的研究還較為有限。

    本文建立環(huán)形燃料組件計(jì)算流體力學(xué)分析模型,分析不同柵徑比下正方形與六邊形環(huán)形燃料組件流動換熱特性,研究組件內(nèi)環(huán)形燃料外周向溫度分布不均勻性。通過歸一化分析各類型子通道處的環(huán)形燃料外周向溫度,分別得出正方形、六邊形環(huán)形燃料組件的最佳柵徑比范圍。

    1 分析模型建立

    本文討論不同柵徑比下正方形與六邊形環(huán)形燃料組件的周向溫度分布,不考慮定位格架等的影響,建立3×3的正方形組件以及19棒束六邊形組件的光棒模型,其橫截面示意圖分別如圖1(a)、(b)所示。參考相關(guān)試驗(yàn)[14],對單根燃料棒的幾何尺寸參數(shù)的設(shè)置如表1 所示,并基于該尺寸參數(shù)建立柵徑比為1.03、1.06、1.09、1.12 的3×3 正方形和19 棒束六邊形燃料組件模型。

    表1 環(huán)形燃料幾何尺寸設(shè)計(jì)參數(shù)Table 1 Ring fuel geometry design parameters

    圖1 燃料棒束橫截面 (a) 正方形組件,(b) 六邊形組件Fig.1 Cross-section of the fuel rod bundle (a) Square component, (b) Hexagonal component

    為節(jié)約計(jì)算資源,本文基于燃料組件為徑向無限大棒束的假設(shè),根據(jù)幾何模型的對稱性,采用1/8正方形燃料組件和1/12六邊形燃料組件構(gòu)建環(huán)形燃料熱工水力分析模型。由于燃料組件中心通道與靠近壁面處的壁面通道冷卻劑流動特征和換熱規(guī)律均存在差異,按圖2所示,將環(huán)形燃料模型的外通道劃分為中心通道、邊通道以及角通道。以距離入口0.25 m、0.50 m 以及0.75 m 處最具代表性的三個(gè)截面來反映不同高度處周向溫度分布的差異,在各截面中的θ 范圍內(nèi)沿外包殼外表面每隔3°~4°取一個(gè)溫度測點(diǎn),正方形排布棒束各子通道均取了15個(gè)溫度測點(diǎn),六邊形排布棒束各子通道選取10個(gè)溫度測點(diǎn),并提取該處的溫度作為各子通道周向溫度的分布。

    圖2 燃料組件計(jì)算域 (a) 正方形組件,(b) 六邊形組件Fig.2 Component computing area of fuel assembly (a) Square component, (b) Hexagonal component

    2 網(wǎng)格劃分及網(wǎng)格敏感性分析

    利用ANSYS中的Mesh模塊對模型進(jìn)行全結(jié)構(gòu)化網(wǎng)格劃分,將內(nèi)、外流體區(qū)域、氣隙區(qū)域以及不同材料的固體區(qū)域通過交界面進(jìn)行連接,交界面網(wǎng)格尺寸相同,以保證不同區(qū)域的網(wǎng)格節(jié)點(diǎn)完全重合,提高計(jì)算的準(zhǔn)確性。關(guān)閉自動清除特征功能,以避免狹小氣隙的特征尺寸在網(wǎng)格劃分時(shí)被清除。由于燃料模型截面存在較多的圓弧線,不利于高質(zhì)量網(wǎng)格的生成,因而設(shè)置精細(xì)的捕獲曲率,此時(shí)在有曲率變化處網(wǎng)格會進(jìn)行自動加密,以控制曲面處網(wǎng)格的變化,使轉(zhuǎn)角處的曲邊網(wǎng)格細(xì)化。對邊緣尺寸進(jìn)行精細(xì)調(diào)整,根據(jù)不同材料的厚度調(diào)整網(wǎng)格數(shù)量,分別對內(nèi)包殼設(shè)置3 層網(wǎng)格、對外包殼設(shè)置4 層網(wǎng)格、對芯塊設(shè)置8 層網(wǎng)格、內(nèi)外氣隙分別設(shè)置一層網(wǎng)格。對內(nèi)通道采用全三角形的掃掠方法,包殼、芯塊、氣隙以及外通道采用全四邊形的掃掠方法,生成的網(wǎng)格如圖3所示。

    圖3 組件網(wǎng)格劃分 (a) 正方形組件,(b) 六邊形組件Fig.3 Mesh partitioning of the components(a) Square component, (b) Hexagonal component

    生成的網(wǎng)格質(zhì)量均大于0.6,不同柵徑比的正方形組件網(wǎng)格數(shù)量在1 588 000~1 708 000 之間,六邊形組件網(wǎng)格數(shù)量在3 032 000~3 263 000之間。

    將柵徑比為1.06的正方形和六邊形組件的網(wǎng)格數(shù)量分別增至5 232 000和8 245 854,通過計(jì)算外通道壓降的變化進(jìn)行網(wǎng)格敏感性分析。正方形組件在網(wǎng)格加密前的外通道壓降為74.109 kPa,進(jìn)行網(wǎng)格加密后為73.778 kPa,僅變化了0.45%;六邊形組件在網(wǎng)格加密前后的外通道壓降分別為150.802 kPa、151.872 kPa,僅變化了0.71%,可認(rèn)為數(shù)值解在未進(jìn)行網(wǎng)格加密前已逐漸收斂。為節(jié)省計(jì)算量與計(jì)算時(shí)間,選取初始網(wǎng)格數(shù)進(jìn)行計(jì)算。

    3 求解條件設(shè)置

    本文參考典型壓水堆的運(yùn)行參數(shù)[15],設(shè)置材料物性參數(shù)及計(jì)算工況分別如表2、3所示。

    表2 材料物性參數(shù)Table 2 Physical properties of the material

    表3 單棒模型運(yùn)行參數(shù)Table 3 Operating parameters for a single fuel rod model

    基于計(jì)算流體力學(xué)程序ANSYS FLUENT 開展模擬分析,參考國內(nèi)外相關(guān)學(xué)者研究結(jié)論,采用SSTk-ω湍流模型[16-18],使熱功率沿軸向呈余弦分布,設(shè)置冷卻劑、包殼、氣隙、芯塊之間的接觸面為耦合傳熱面,并設(shè)置為無滑移界面,組件外壁及頂、底部均與外界絕熱,組件對稱分割面設(shè)置為對稱面,各物理量在該面處法向梯度為零。內(nèi)、外通道入口設(shè)置為質(zhì)量流入口,質(zhì)量流量按橫截面積的大小比例進(jìn)行分配,出口均為壓力出口,所有方程的離散格式采用較高精度的二階迎風(fēng)格式,設(shè)置計(jì)算的迭代步數(shù)為200,所有計(jì)算的物理量都能在200 步以內(nèi)收斂,且殘差均小于10-5。

    4 結(jié)果分析及討論

    4.1 溫度、角度及不均勻系數(shù)的歸一化

    為便于對不同柵徑比、不同子通道類型,以及不同排布的燃料組件周向溫度分布進(jìn)行比較,將所取得的溫度數(shù)值及其對應(yīng)的角度進(jìn)行min-max 歸一化,令歸一化溫度為:

    角度的歸一化與此類似,將溫度與角度的值均固定在0~1 之間以得到不同角度處歸一化的溫度分布。

    并且為了更加準(zhǔn)確、直觀地描述溫度分布的不均勻性,引入顆粒分析方法中的篩分法[19]計(jì)算不均勻系數(shù)Cu,對不均勻系數(shù)Cu的說明如下:

    定義一個(gè)數(shù)列:X={x1,x2,x3,…,xn}以及兩個(gè)閾值C10、C80,將所有滿足xi≤C10(i= 1,2,3,…,n)的xi進(jìn)行求和,使其求得的和恰好不大于數(shù)列X之和的10%;將所有滿足xi≤C80(i= 1,2,3,…,n)的xi進(jìn)行求和,使其求得的和恰好不大于數(shù)列X之和的80%,得到數(shù)列的不均勻系數(shù):

    不均勻系數(shù)越大,則數(shù)據(jù)分布的不均勻性越大,不均勻系數(shù)越接近于1,則數(shù)據(jù)分布的不均勻性越小。

    對不同柵徑比下各子通道外包殼周向溫度分布及其不均勻系數(shù)進(jìn)行整理后觀察到如下現(xiàn)象。

    4.2 正方形組件中心通道周向溫度分布

    將不同柵徑比的各個(gè)正方形組件相同高度處的外包殼表面溫度以及對應(yīng)的角度作min-max 歸一化,得到如圖4 所示的正方形組件中心通道周向溫度分布圖。由圖4 發(fā)現(xiàn),外包殼周向溫度分布數(shù)據(jù)點(diǎn)存在些許震蕩與跳動,相關(guān)研究結(jié)果表明,外包殼周向溫度分布的不均勻性主要取決于流道局部水力直徑的變化[20],并且與冷卻劑的橫向流動有關(guān)[21],因此可認(rèn)為其是由緊密柵內(nèi)環(huán)形燃料棒束中復(fù)雜的流動特征以及強(qiáng)烈的交混所造成的。進(jìn)一步可看出正方形組件中心通道處周向溫度分布的均勻性受柵徑比的影響較大,兩燃料棒的間隙處溫度較高,而隨著θ角的增大,溫度逐漸降低,當(dāng)柵徑比為1.03 時(shí),中心通道處周向溫度分布最不均勻,而將柵徑比增大至1.06 后,中心通道處周向溫度分布的均勻性得到了明顯的改善,當(dāng)柵徑比增大至1.09,此時(shí)的周向溫度分布與柵徑比等于1.12 時(shí)的十分相近,此后若繼續(xù)增大柵徑比對周向溫度分布不均勻性的改善意義不大。

    圖4 正方形組件中心通道周向溫度分布 (a) 距入口0.25 m處,(b) 距入口0.50 m處,(c) 距入口0.75 m處Fig.4 The central channel temperature of a square component distributed along the circumference(a) 0.25 m from the entrance, (b) 0.50 m from the entrance, (c) 0.75 m from the entrance

    計(jì)算出所有組件三個(gè)子通道的三個(gè)典型高度處周向溫度分布的不均勻系數(shù),并對其作min-max 歸一化,將正方形組件中心通道的歸一化不均勻系數(shù)展示如圖5。可以看出,隨著高度的增大,歸一化不均勻系數(shù)也隨之增大,冷卻劑逐漸成為充分發(fā)展流,流道較窄處冷卻劑流量逐漸減小,較寬處冷卻劑流量逐漸增大。而隨著柵徑比的增大,歸一化不均勻系數(shù)隨之減小,柵徑比從1.03增至1.06時(shí),歸一化不均勻系數(shù)減小幅度最大,在0.25 m、0.50 m、0.75 m高度處分別由0.46 減小到0.14、0.60 減小到0.30、0.80減小到0.34,而將柵徑比從1.09增大到1.12時(shí),歸一化不均勻系數(shù)的減小量均低于0.1?;趯χ行耐ǖ赖姆治?,可以認(rèn)為,目前大多數(shù)環(huán)形燃料試驗(yàn)將柵徑比設(shè)置為1.07是合理的。

    圖5 正方形組件中心通道歸一化不均勻系數(shù)Fig.5 Normalized coefficient of nonuniformity of the center channel of the square component

    4.3 正方形組件邊通道周向溫度分布

    使用§4.2中相同的方法得出正方形組件邊通道歸一化周向溫度分布如圖6 所示,并將其歸一化不均勻系數(shù)展示如圖7所示。由邊通道溫度分布可以看出,燃料棒靠近壁面的區(qū)域與兩燃料棒間隙處溫度均較高,其中靠近壁面處溫度比間隙處溫度高出不少,邊通道周向溫度分布的不均勻性主要由燃料棒與壁面距離決定,因此其周向溫度分布的不均勻性較中心通道的大。通過觀察發(fā)現(xiàn),柵徑比越小的組件其棒束間隙處及近壁面處溫度越高,外通道流道較寬處溫度越低,由此可以推得柵徑比越小其外流道較窄處和較寬處冷卻劑流量分配差異越大。

    圖6 正方形組件邊通道周向溫度分布 (a) 距入口0.25 m處,(b) 距入口0.50 m處,(c) 距入口0.75 m處Fig.6 The side channel temperature of a square component distributed along the circumference(a) 0.25 m from the entrance, (b) 0.50 m from the entrance, (c) 0.75 m from the entrance

    圖7 正方形組件邊通道歸一化不均勻系數(shù)Fig.7 Normalized coefficient of nonuniformity of the side channel of the square component

    由圖6 可以看出,柵徑比從1.03 增至1.06 時(shí)周向溫度分布仍有較大的改善,其0.25 m、0.50 m、0.75 m高度處的歸一化不均勻系數(shù)分別從0.33減小到0.25、0.66 減小到0.48、1 減小到0.61,但降低程度較中心通道的小。

    4.4 正方形組件角通道周向溫度分布

    沿用上述方法得到歸一化周向溫度分布展示如圖8所示,歸一化不均勻系數(shù)展示如圖9所示。由角通道周向溫度分布可以看出,燃料棒靠近壁面處溫度較高,其溫度分布僅受壁面距離的影響,且分布趨勢與邊通道歸一化角度在0~0.5 范圍內(nèi)一致。由圖8可看出,柵徑比從1.03增至1.06時(shí),周向溫度分布不均勻性大幅改善,其0.25 m、0.50 m、0.75 m 高度處的歸一化不均勻系數(shù)分別從0.23 減小到0.20、0.46 減小到0.29、0.77 減小到0.51,但改善程度均小于中心通道和邊通道,而柵徑比從1.06增至1.09時(shí),周向溫度分布的不均勻性幾乎不變,這是由于此時(shí)增大柵徑比時(shí)僅增大了少許燃料棒與壁面間的距離,導(dǎo)致靠近壁面處的溫度僅有少許下降,對周向溫度的展平效果有限,與邊通道類似,可通過增大壁面距離對其進(jìn)行改善。

    圖8 正方形組件角通道周向溫度分布 (a) 距入口0.25 m處,(b) 距入口0.50 m處,(c) 距入口0.75 m處Fig.8 Angular channel temperature of a square component distributed along the circumference(a) 0.25 m from the entrance, (b) 0.50 m from the entrance, (c) 0.75 m from the entrance

    圖9 正方形組件邊通道歸一化不均勻系數(shù)Fig.9 Normalized coefficient of nonuniformity of the angular channel of the square component

    4.5 六邊形組件各子通道周向溫度分布

    按照上述方法得到六邊形組件中心通道、邊通道以及角通道的三個(gè)典型高度處歸一化周向溫度分布分別如圖10~12 所示,并將計(jì)算得的不均勻系數(shù)展示如圖13 所示。六邊形組件各子通道外包殼的周向溫度的分布規(guī)律與正方形組件的類似,且將柵徑比從1.06增大至1.09后對周向溫度的展平效果均不明顯。但與正方形組件不同的是,其柵徑比從1.03增大至1.06時(shí),中心通道、邊通道以及角通道的歸一化不均勻系數(shù)分別下降了0.21~0.37、0.16~0.24以及0.09~0.42,各子通道對柵徑比的敏感性從大到小依次為角通道、中心通道、邊通道,這是由于六邊形組件角通道面積太小造成的。在同一柵徑比下,六邊形組件各子通道的周向溫度分布均比正方形組件的均勻,各柵徑比下中心通道、邊通道和角通道的歸一化不均勻系數(shù)分別比正方形組件各子通道歸一化不均勻系數(shù)平均低了0.12、0.13 和0.1,這主要是由于六邊形組件各子通道水力直徑較小,各子通道的間隙區(qū)與中心區(qū)流道寬度差異不如正方形組件的大。綜合上述考慮可將六邊形組件柵徑比設(shè)在1.06~1.09之間。

    圖10 六邊形組件中心通道周向溫度分布 (a) 距入口0.25 m處,(b) 距入口0.50 m處,(c) 距入口0.75 m處Fig.10 Central channel temperature of a hexagonal component distributed along the circumference(a) 0.25 m from the entrance, (b) 0.50 m from the entrance, (c) 0.75 m from the entrance

    圖11 六邊形組件邊通道周向溫度分布 (a) 距入口0.25 m處,(b) 距入口0.50 m處,(c) 距入口0.75 m處Fig.11 Side channel temperature of a hexagonal component distributed along the circumference(a) 0.25 m from the entrance, (b) 0.50 m from the entrance, (c) 0.75 m from the entrance

    圖12 六邊形組件角通道周向溫度分布 (a) 距入口0.25 m處,(b) 距入口0.50 m處,(c) 距入口0.75 m處Fig.12 Angular channel temperature of a hexagonal component distributed along the circumference(a) 0.25 m from the entrance, (b) 0.50 m from the entrance, (c) 0.75 m from the entrance

    圖13 六邊形組件溫度分布不均勻系數(shù) (a) 中心通道,(b) 邊通道,(c) 角通道Fig.13 Normalized coefficient of nonuniformity of temperature distribution of hexagonal components(a) Central channel, (b) Side channel, (c) Angular channel

    5 結(jié)語

    本文基于ANSYS FLUENT 建立環(huán)形燃料組件計(jì)算流體力學(xué)分析模型,分析不同柵徑比下正方形與六邊形環(huán)形燃料組件流動換熱特性,研究組件內(nèi)環(huán)形燃料外周向溫度分布不均勻性,并分別得出正方形、六邊形環(huán)形燃料組件的最佳柵徑比范圍?;诓煌β拭芏?、流量大小會對包殼溫度大小進(jìn)行影響,本文在僅有柵徑比不同的條件下,對溫度結(jié)果分別進(jìn)行了min-max 歸一化的處理,最終得出無量綱溫度分布趨勢,因此本文結(jié)論在不同功率密度以及流量大小情況下仍然成立,主要結(jié)論如下:

    1)增大柵徑比有利于對環(huán)形燃料外周向溫度進(jìn)行展平,對正方形組件而言,柵徑比在1.07~1.09 之間較為合適,對于六邊形組件而言,其外周向溫度分布不均勻性略低于正方形組件,柵徑比可選擇在1.06~1.09之間。

    2)增大柵徑比對燃料棒間隙處溫度的改善最為明顯,對近壁面處溫度的改善次之,對正方形組件而言,中心通道的環(huán)形燃料棒外周向溫度分布對柵徑比最為敏感,其次是邊通道,最后是角通道;對六邊組件而言,角通道周向溫度分布對柵徑比最為敏感,其次是中心通道,最后是邊通道。

    作者貢獻(xiàn)聲明曾付林負(fù)責(zé)起草文章,分析/解釋數(shù)據(jù);唐錨負(fù)責(zé)實(shí)施研究,采集數(shù)據(jù);趙鵬程負(fù)責(zé)獲取研究經(jīng)費(fèi),對文章的知識性內(nèi)容作批評性審閱,指導(dǎo);向釗才負(fù)責(zé)支持性貢獻(xiàn)。

    猜你喜歡
    周向六邊形正方形
    周向拉桿轉(zhuǎn)子瞬態(tài)應(yīng)力分析與啟動曲線優(yōu)化
    知識快餐店 到處都是六邊形
    剪正方形
    剪拼正方形
    拼正方形
    拼正方形
    創(chuàng)意六邊形無限翻
    童話世界(2018年32期)2018-12-03 05:14:56
    周向定位旋轉(zhuǎn)分度鉆模設(shè)計(jì)
    怎樣剪拼
    怎樣剪拼
    欧美日本中文国产一区发布| 色婷婷av一区二区三区视频| 91aial.com中文字幕在线观看| 日韩伦理黄色片| 99热网站在线观看| 亚洲欧美日韩卡通动漫| 乱系列少妇在线播放| 人妻系列 视频| 亚洲美女搞黄在线观看| 免费看不卡的av| 成人亚洲精品一区在线观看| 亚洲精品日韩在线中文字幕| 欧美三级亚洲精品| 岛国毛片在线播放| 成人午夜精彩视频在线观看| 国产一级毛片在线| 日韩 亚洲 欧美在线| 欧美日韩av久久| 麻豆成人午夜福利视频| 男人爽女人下面视频在线观看| 五月天丁香电影| 免费观看的影片在线观看| 日日摸夜夜添夜夜爱| 青春草国产在线视频| 久久人人爽人人片av| 亚洲精品成人av观看孕妇| 最近的中文字幕免费完整| 自线自在国产av| 国产色婷婷99| 不卡视频在线观看欧美| av网站免费在线观看视频| 少妇被粗大猛烈的视频| 日韩成人av中文字幕在线观看| 男人爽女人下面视频在线观看| 伊人久久国产一区二区| 日韩熟女老妇一区二区性免费视频| 男女边摸边吃奶| 极品人妻少妇av视频| 青青草视频在线视频观看| 国产乱来视频区| 少妇被粗大的猛进出69影院 | 欧美日韩国产mv在线观看视频| 内地一区二区视频在线| av天堂中文字幕网| 亚洲自偷自拍三级| 日韩 亚洲 欧美在线| 一级毛片aaaaaa免费看小| 亚洲经典国产精华液单| 视频区图区小说| 久久综合国产亚洲精品| 国产日韩欧美视频二区| a级一级毛片免费在线观看| 黄色日韩在线| 哪个播放器可以免费观看大片| 高清av免费在线| 亚洲婷婷狠狠爱综合网| 97精品久久久久久久久久精品| 丁香六月天网| 黑丝袜美女国产一区| 亚洲美女黄色视频免费看| 色吧在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲精品456在线播放app| 中文在线观看免费www的网站| 99国产精品免费福利视频| 老司机影院毛片| 精华霜和精华液先用哪个| 激情五月婷婷亚洲| 日本欧美国产在线视频| 久久国产精品男人的天堂亚洲 | 熟妇人妻不卡中文字幕| www.av在线官网国产| 国产极品天堂在线| 久久影院123| 插阴视频在线观看视频| 免费观看无遮挡的男女| 久久久久久久久久久免费av| 老司机影院毛片| 黄色日韩在线| 老司机影院毛片| 亚洲成色77777| av.在线天堂| 交换朋友夫妻互换小说| 亚洲精品,欧美精品| 男人舔奶头视频| 精华霜和精华液先用哪个| 亚洲综合色惰| 一级片'在线观看视频| 99九九线精品视频在线观看视频| 久久国产亚洲av麻豆专区| 久久 成人 亚洲| 欧美成人精品欧美一级黄| 久久亚洲国产成人精品v| 亚洲av欧美aⅴ国产| 中国三级夫妇交换| 亚洲av免费高清在线观看| 又爽又黄a免费视频| 国产白丝娇喘喷水9色精品| 亚洲天堂av无毛| 9色porny在线观看| 久久精品夜色国产| 欧美+日韩+精品| 国产极品天堂在线| 午夜免费鲁丝| 国产精品久久久久久久久免| 国产白丝娇喘喷水9色精品| 高清黄色对白视频在线免费看 | 天堂8中文在线网| 亚洲精品第二区| 国产欧美日韩精品一区二区| 久久久亚洲精品成人影院| 国产国拍精品亚洲av在线观看| 国语对白做爰xxxⅹ性视频网站| 高清毛片免费看| 黄色配什么色好看| 99久久精品热视频| 国产精品久久久久成人av| 老熟女久久久| 黄色日韩在线| 国产成人免费观看mmmm| 午夜福利视频精品| 国产乱人偷精品视频| 国产白丝娇喘喷水9色精品| 婷婷色综合大香蕉| 日日摸夜夜添夜夜爱| av不卡在线播放| 国产中年淑女户外野战色| 美女xxoo啪啪120秒动态图| 亚洲欧美一区二区三区黑人 | 狂野欧美激情性xxxx在线观看| 精品一区在线观看国产| 久久6这里有精品| 日日爽夜夜爽网站| 国语对白做爰xxxⅹ性视频网站| 欧美精品亚洲一区二区| 在线观看www视频免费| 中文字幕免费在线视频6| 成年女人在线观看亚洲视频| 国产伦理片在线播放av一区| 亚洲欧美精品自产自拍| 美女xxoo啪啪120秒动态图| 久久久久久久国产电影| 边亲边吃奶的免费视频| 中文字幕av电影在线播放| 国产日韩欧美亚洲二区| av在线观看视频网站免费| 一区二区三区四区激情视频| 亚洲三级黄色毛片| 五月伊人婷婷丁香| 国产精品一区www在线观看| 久久av网站| 亚洲欧洲国产日韩| 亚洲欧美一区二区三区国产| 久久久久久久久久人人人人人人| 国产国拍精品亚洲av在线观看| 国产亚洲一区二区精品| 高清午夜精品一区二区三区| 18禁在线播放成人免费| tube8黄色片| 一级,二级,三级黄色视频| 观看美女的网站| 久久精品国产亚洲av天美| 国产免费福利视频在线观看| 91在线精品国自产拍蜜月| 精品人妻一区二区三区麻豆| 国产黄色免费在线视频| 六月丁香七月| 国产精品不卡视频一区二区| 日韩三级伦理在线观看| 如何舔出高潮| 亚洲av日韩在线播放| 成人18禁高潮啪啪吃奶动态图 | 观看免费一级毛片| 国产真实伦视频高清在线观看| 99久久综合免费| 免费人成在线观看视频色| 国产真实伦视频高清在线观看| 国模一区二区三区四区视频| 99re6热这里在线精品视频| 七月丁香在线播放| 超碰97精品在线观看| 美女脱内裤让男人舔精品视频| 色哟哟·www| 国产精品伦人一区二区| 一级毛片黄色毛片免费观看视频| av在线播放精品| 亚洲国产毛片av蜜桃av| 99热国产这里只有精品6| 一本一本综合久久| 五月开心婷婷网| 啦啦啦在线观看免费高清www| 亚洲成人一二三区av| 大陆偷拍与自拍| 亚洲欧洲精品一区二区精品久久久 | 精品卡一卡二卡四卡免费| 一级毛片aaaaaa免费看小| 精华霜和精华液先用哪个| 男人和女人高潮做爰伦理| 伦精品一区二区三区| av免费在线看不卡| 热re99久久精品国产66热6| 99精国产麻豆久久婷婷| 一级毛片久久久久久久久女| 一级毛片aaaaaa免费看小| 精品国产乱码久久久久久小说| 日韩亚洲欧美综合| 看十八女毛片水多多多| 18禁在线播放成人免费| 亚洲丝袜综合中文字幕| 少妇丰满av| av女优亚洲男人天堂| 国产成人freesex在线| 国产亚洲91精品色在线| av在线app专区| 欧美亚洲 丝袜 人妻 在线| 一区二区三区免费毛片| 欧美激情国产日韩精品一区| 久久久精品94久久精品| 少妇精品久久久久久久| 最新中文字幕久久久久| 天堂俺去俺来也www色官网| 99精国产麻豆久久婷婷| 99九九线精品视频在线观看视频| 成年女人在线观看亚洲视频| 最近最新中文字幕免费大全7| 国产精品一二三区在线看| 久久女婷五月综合色啪小说| av天堂久久9| 精品亚洲成国产av| 观看av在线不卡| 热99国产精品久久久久久7| 国产探花极品一区二区| 国产精品一区二区在线不卡| 男女无遮挡免费网站观看| 最黄视频免费看| 欧美日韩av久久| 亚洲久久久国产精品| 亚州av有码| 嘟嘟电影网在线观看| 国产av码专区亚洲av| 国产综合精华液| 精品久久久噜噜| 国产一区二区在线观看日韩| 久久精品久久久久久久性| 永久网站在线| 一区二区三区免费毛片| 亚洲中文av在线| 91久久精品国产一区二区三区| 十八禁高潮呻吟视频 | 精品人妻一区二区三区麻豆| 日韩伦理黄色片| 在线免费观看不下载黄p国产| 九九在线视频观看精品| 嫩草影院入口| 看免费成人av毛片| 亚洲情色 制服丝袜| 日韩电影二区| 国产精品人妻久久久影院| 久久久久久久国产电影| 亚洲天堂av无毛| 熟女人妻精品中文字幕| 日本黄大片高清| 卡戴珊不雅视频在线播放| kizo精华| 亚洲精品国产成人久久av| 又黄又爽又刺激的免费视频.| 哪个播放器可以免费观看大片| 丝袜脚勾引网站| 久久国产精品大桥未久av | 国产有黄有色有爽视频| 一区二区三区四区激情视频| 国产乱人偷精品视频| 精品熟女少妇av免费看| 人人妻人人澡人人看| 视频中文字幕在线观看| 国产乱来视频区| 欧美人与善性xxx| 国产极品粉嫩免费观看在线 | 成人影院久久| 国产精品久久久久久久久免| 欧美最新免费一区二区三区| av又黄又爽大尺度在线免费看| 狂野欧美白嫩少妇大欣赏| 欧美老熟妇乱子伦牲交| 91aial.com中文字幕在线观看| 久久精品国产亚洲av天美| 十八禁网站网址无遮挡 | 男人爽女人下面视频在线观看| 黄色日韩在线| 嫩草影院入口| 亚洲国产日韩一区二区| 久久人人爽av亚洲精品天堂| 日韩精品免费视频一区二区三区 | 免费观看在线日韩| 日日摸夜夜添夜夜爱| 尾随美女入室| 在线观看国产h片| 久久精品熟女亚洲av麻豆精品| 99久久中文字幕三级久久日本| 国产欧美日韩一区二区三区在线 | 国产精品女同一区二区软件| 国产又色又爽无遮挡免| 一级毛片电影观看| av视频免费观看在线观看| 久久久久视频综合| 久久久久久久亚洲中文字幕| 久久久久久久久大av| 中文字幕亚洲精品专区| 婷婷色综合大香蕉| 亚洲av不卡在线观看| 亚洲精品视频女| 在现免费观看毛片| 日本-黄色视频高清免费观看| 日日爽夜夜爽网站| 国产一区二区三区综合在线观看 | 97超视频在线观看视频| 91成人精品电影| 亚洲美女黄色视频免费看| 99热6这里只有精品| 边亲边吃奶的免费视频| 亚洲精品日本国产第一区| av福利片在线观看| 日日爽夜夜爽网站| 日韩熟女老妇一区二区性免费视频| 色视频在线一区二区三区| 免费看不卡的av| 亚洲欧洲日产国产| av在线老鸭窝| 大码成人一级视频| 日韩视频在线欧美| 国产成人精品婷婷| 久久精品熟女亚洲av麻豆精品| 一本—道久久a久久精品蜜桃钙片| 熟妇人妻不卡中文字幕| 亚洲怡红院男人天堂| 亚洲人成网站在线播| 亚洲精品日韩av片在线观看| 丝袜在线中文字幕| 精品国产乱码久久久久久小说| 国产精品伦人一区二区| 最黄视频免费看| 91精品国产九色| 国产一区二区三区av在线| 韩国av在线不卡| 一级,二级,三级黄色视频| a 毛片基地| 成人毛片60女人毛片免费| 日本欧美国产在线视频| 纵有疾风起免费观看全集完整版| 久久久久精品性色| 亚洲一区二区三区欧美精品| 青春草国产在线视频| 汤姆久久久久久久影院中文字幕| 韩国高清视频一区二区三区| 国产 一区精品| 亚洲av综合色区一区| 街头女战士在线观看网站| 亚洲欧美日韩另类电影网站| 一区二区三区精品91| 久久人人爽人人片av| 国产黄色免费在线视频| 性色av一级| 国产亚洲5aaaaa淫片| 亚洲欧美日韩另类电影网站| 亚洲欧美精品专区久久| 18禁动态无遮挡网站| 3wmmmm亚洲av在线观看| 亚洲人成网站在线播| 又大又黄又爽视频免费| 又爽又黄a免费视频| 国产成人91sexporn| 99久久人妻综合| 精品国产一区二区久久| 精品亚洲乱码少妇综合久久| 成人国产av品久久久| 尾随美女入室| 嫩草影院新地址| 国产探花极品一区二区| 国产成人a∨麻豆精品| 女人久久www免费人成看片| 日日摸夜夜添夜夜爱| 女人久久www免费人成看片| 免费观看无遮挡的男女| 99精国产麻豆久久婷婷| 色婷婷久久久亚洲欧美| 伊人亚洲综合成人网| 人妻少妇偷人精品九色| 免费高清在线观看视频在线观看| 婷婷色av中文字幕| 26uuu在线亚洲综合色| 国产爽快片一区二区三区| 又大又黄又爽视频免费| 乱码一卡2卡4卡精品| 色吧在线观看| 精品一品国产午夜福利视频| 亚洲欧洲精品一区二区精品久久久 | 在线观看人妻少妇| 大片电影免费在线观看免费| 亚洲人成网站在线播| 高清午夜精品一区二区三区| 成人亚洲欧美一区二区av| 国产亚洲91精品色在线| 亚洲国产成人一精品久久久| 我的老师免费观看完整版| 国产av国产精品国产| 久久婷婷青草| 亚洲一级一片aⅴ在线观看| 亚洲av不卡在线观看| 亚洲精品久久久久久婷婷小说| 插阴视频在线观看视频| 久久精品国产亚洲av涩爱| 成人漫画全彩无遮挡| 欧美一级a爱片免费观看看| 免费播放大片免费观看视频在线观看| 国产日韩欧美视频二区| 寂寞人妻少妇视频99o| 黄色毛片三级朝国网站 | 精品国产一区二区久久| 色视频www国产| 中文天堂在线官网| 人妻人人澡人人爽人人| 成年人免费黄色播放视频 | 成年av动漫网址| 美女国产视频在线观看| 男人添女人高潮全过程视频| 青春草亚洲视频在线观看| 国产成人freesex在线| 毛片一级片免费看久久久久| 人体艺术视频欧美日本| av女优亚洲男人天堂| 最新的欧美精品一区二区| 美女视频免费永久观看网站| 美女国产视频在线观看| 在线观看免费视频网站a站| 在线观看国产h片| 国产欧美亚洲国产| 免费观看的影片在线观看| 免费观看在线日韩| av国产久精品久网站免费入址| 久久国产精品男人的天堂亚洲 | 欧美另类一区| av在线app专区| 少妇人妻一区二区三区视频| 最近中文字幕高清免费大全6| 久久综合国产亚洲精品| 日韩一区二区视频免费看| 爱豆传媒免费全集在线观看| 91精品国产九色| 在线观看免费日韩欧美大片 | 午夜精品国产一区二区电影| 黄色日韩在线| 最黄视频免费看| 久久久久精品性色| av有码第一页| 成人影院久久| 乱码一卡2卡4卡精品| 日本色播在线视频| 狂野欧美激情性bbbbbb| 亚洲国产精品一区二区三区在线| 王馨瑶露胸无遮挡在线观看| 国产伦理片在线播放av一区| 亚洲欧洲日产国产| 精品久久久噜噜| 日本vs欧美在线观看视频 | 精品人妻一区二区三区麻豆| 三级国产精品欧美在线观看| 亚洲综合精品二区| 伦理电影大哥的女人| 亚洲国产精品一区三区| av视频免费观看在线观看| 亚洲精品色激情综合| av在线老鸭窝| 一级片'在线观看视频| 极品少妇高潮喷水抽搐| 午夜日本视频在线| 男人爽女人下面视频在线观看| 男女免费视频国产| 亚洲高清免费不卡视频| 国产男女超爽视频在线观看| 久久午夜福利片| 熟女电影av网| 我要看黄色一级片免费的| 国产成人精品久久久久久| 亚洲精品日本国产第一区| 成人毛片a级毛片在线播放| 婷婷色综合www| 欧美日韩一区二区视频在线观看视频在线| 国产精品女同一区二区软件| 亚洲av成人精品一区久久| 久久久久视频综合| 乱码一卡2卡4卡精品| 色婷婷av一区二区三区视频| 内地一区二区视频在线| 少妇人妻 视频| 99热这里只有精品一区| 亚洲国产av新网站| 乱人伦中国视频| 五月伊人婷婷丁香| 韩国高清视频一区二区三区| 欧美日韩视频精品一区| 成人二区视频| 五月伊人婷婷丁香| 一级黄片播放器| 国产成人精品久久久久久| 美女福利国产在线| 日韩熟女老妇一区二区性免费视频| 亚洲国产精品一区三区| 麻豆成人av视频| 日韩大片免费观看网站| 色婷婷久久久亚洲欧美| 亚洲国产日韩一区二区| √禁漫天堂资源中文www| 国产91av在线免费观看| 美女大奶头黄色视频| 午夜91福利影院| 国产成人freesex在线| 丰满乱子伦码专区| 精品一区二区三卡| 九色成人免费人妻av| 久久婷婷青草| av天堂久久9| 久久毛片免费看一区二区三区| 久久ye,这里只有精品| 中国美白少妇内射xxxbb| 亚洲欧洲精品一区二区精品久久久 | 人人妻人人看人人澡| 美女国产视频在线观看| 久久精品国产a三级三级三级| 国产熟女午夜一区二区三区 | 嘟嘟电影网在线观看| 国产熟女午夜一区二区三区 | 国产成人aa在线观看| 久久人妻熟女aⅴ| av国产久精品久网站免费入址| 日本黄色日本黄色录像| 亚洲av电影在线观看一区二区三区| 日韩精品免费视频一区二区三区 | 女性被躁到高潮视频| 精品国产露脸久久av麻豆| 久久97久久精品| 麻豆成人av视频| 黄色怎么调成土黄色| 日本黄色日本黄色录像| 国产日韩一区二区三区精品不卡 | 亚洲人与动物交配视频| 22中文网久久字幕| 狂野欧美激情性bbbbbb| 一级,二级,三级黄色视频| 最近2019中文字幕mv第一页| 亚洲真实伦在线观看| 免费在线观看成人毛片| 国产精品偷伦视频观看了| 国产精品99久久99久久久不卡 | 国产视频内射| 亚洲一区二区三区欧美精品| 2021少妇久久久久久久久久久| 肉色欧美久久久久久久蜜桃| 免费观看无遮挡的男女| 黑丝袜美女国产一区| 亚洲四区av| 日本黄大片高清| 成年人午夜在线观看视频| 我的老师免费观看完整版| 久久精品国产亚洲网站| 哪个播放器可以免费观看大片| 又爽又黄a免费视频| 最近最新中文字幕免费大全7| 97在线人人人人妻| 日日撸夜夜添| 久久狼人影院| 国产精品人妻久久久影院| 涩涩av久久男人的天堂| 日本wwww免费看| 女的被弄到高潮叫床怎么办| 免费高清在线观看视频在线观看| 久久久久精品性色| 高清黄色对白视频在线免费看 | 久久鲁丝午夜福利片| 天堂8中文在线网| 中国美白少妇内射xxxbb| 亚洲精品国产色婷婷电影| av视频免费观看在线观看| av天堂久久9| 国产成人午夜福利电影在线观看| 一边亲一边摸免费视频| 日韩精品免费视频一区二区三区 | 黄色毛片三级朝国网站 | 亚洲国产精品一区三区| 日本黄色日本黄色录像| 午夜91福利影院| 99视频精品全部免费 在线| 国产成人a∨麻豆精品| 女人久久www免费人成看片| 成人特级av手机在线观看| 久热这里只有精品99| 欧美另类一区| 蜜桃在线观看..| 国产精品欧美亚洲77777| 两个人免费观看高清视频 | 久久av网站| 亚洲精品久久午夜乱码| 久久99热6这里只有精品| 只有这里有精品99| 热re99久久精品国产66热6| 国产 一区精品| 欧美日韩亚洲高清精品| 国内揄拍国产精品人妻在线| av线在线观看网站| 亚洲av成人精品一二三区| 国产深夜福利视频在线观看| 日本wwww免费看| 精品亚洲成国产av| 精品久久久久久久久亚洲| 欧美日韩亚洲高清精品| 久久国产乱子免费精品|