• <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视频| 免费av不卡在线播放| 亚洲精品一二三| 成人二区视频| 国产精品熟女久久久久浪| 看免费av毛片| 免费大片18禁| 伦理电影免费视频| 国产欧美日韩综合在线一区二区| av播播在线观看一区| 亚洲一区二区三区欧美精品| 日本-黄色视频高清免费观看| 精品人妻熟女毛片av久久网站| 久久国产精品男人的天堂亚洲 | 午夜福利,免费看| 两个人免费观看高清视频| 日韩电影二区| 久久久久久久亚洲中文字幕| 男女午夜视频在线观看 | 久久人妻熟女aⅴ| 99视频精品全部免费 在线| 日韩熟女老妇一区二区性免费视频| 十八禁网站网址无遮挡| 亚洲久久久国产精品| 高清视频免费观看一区二区| 国产又色又爽无遮挡免| 久久婷婷青草| 国产淫语在线视频| 久久婷婷青草| 精品视频人人做人人爽| 巨乳人妻的诱惑在线观看| 久久久国产欧美日韩av| 国产日韩欧美在线精品| 成年av动漫网址| 中文字幕亚洲精品专区| 国产69精品久久久久777片| 韩国av在线不卡| 天堂俺去俺来也www色官网| 飞空精品影院首页| 婷婷成人精品国产| 久久久久精品性色| 久久精品夜色国产| 99热这里只有是精品在线观看| 午夜福利乱码中文字幕| 最近2019中文字幕mv第一页| 日本欧美国产在线视频| 高清欧美精品videossex| 又黄又粗又硬又大视频| 老司机亚洲免费影院| 人妻少妇偷人精品九色| 亚洲精品成人av观看孕妇| 激情视频va一区二区三区| 一级黄片播放器| 妹子高潮喷水视频| 狠狠精品人妻久久久久久综合| 中文字幕亚洲精品专区| 韩国精品一区二区三区 | 99热国产这里只有精品6| 国产精品国产三级专区第一集| 国产精品成人在线| 69精品国产乱码久久久| 色婷婷av一区二区三区视频| 亚洲精品色激情综合| 国产欧美另类精品又又久久亚洲欧美| www.色视频.com| 又粗又硬又长又爽又黄的视频| 欧美少妇被猛烈插入视频| 视频中文字幕在线观看| 久久韩国三级中文字幕| 亚洲精品国产av成人精品| 女性被躁到高潮视频| 一级黄片播放器| 狠狠精品人妻久久久久久综合| 亚洲av综合色区一区| 亚洲图色成人| 丁香六月天网| 两个人免费观看高清视频| 亚洲久久久国产精品| videosex国产| 你懂的网址亚洲精品在线观看| 欧美精品人与动牲交sv欧美| 一级,二级,三级黄色视频| 亚洲精品,欧美精品| 久久久久精品久久久久真实原创| 99香蕉大伊视频| 中国三级夫妇交换| 日本猛色少妇xxxxx猛交久久| 久久久精品94久久精品| 在线观看美女被高潮喷水网站| 99久久综合免费| a 毛片基地| a级毛色黄片| 婷婷色av中文字幕| 亚洲国产成人一精品久久久| 亚洲欧美色中文字幕在线| 免费看不卡的av| 不卡视频在线观看欧美| 人人妻人人澡人人爽人人夜夜| 欧美日韩一区二区视频在线观看视频在线| 久久99一区二区三区| freevideosex欧美| 26uuu在线亚洲综合色| 日日撸夜夜添| 精品久久国产蜜桃| 高清不卡的av网站| www日本在线高清视频| 99热这里只有是精品在线观看| 丝袜在线中文字幕| 亚洲伊人色综图| 夜夜爽夜夜爽视频| 99re6热这里在线精品视频| 亚洲性久久影院| 国产在线免费精品| 久久国产精品大桥未久av| 午夜福利网站1000一区二区三区| 汤姆久久久久久久影院中文字幕| 五月伊人婷婷丁香| 日韩精品有码人妻一区| 我的女老师完整版在线观看| 日韩av免费高清视频| 黄色怎么调成土黄色| 最近中文字幕高清免费大全6| 高清欧美精品videossex| 建设人人有责人人尽责人人享有的| 美女视频免费永久观看网站| 久久精品国产a三级三级三级| 亚洲伊人色综图| 国产高清三级在线| 日韩在线高清观看一区二区三区| 国产精品.久久久| 美女大奶头黄色视频| 久久99蜜桃精品久久| 日韩大片免费观看网站| 久久久精品94久久精品| 观看美女的网站| 亚洲综合精品二区| 亚洲国产精品成人久久小说| 亚洲精品456在线播放app| 丝袜喷水一区| 伦精品一区二区三区| 精品国产一区二区三区久久久樱花| 午夜福利网站1000一区二区三区| 国产精品国产三级专区第一集| 人人妻人人澡人人看| 日韩中文字幕视频在线看片| 午夜免费鲁丝| 婷婷色综合www| 蜜桃国产av成人99| 丁香六月天网| 寂寞人妻少妇视频99o| 超色免费av| 亚洲综合色网址| 久久精品国产综合久久久 | 国产成人一区二区在线| 免费高清在线观看日韩| 永久网站在线| 纵有疾风起免费观看全集完整版| 黄片无遮挡物在线观看| 最近最新中文字幕免费大全7| 亚洲av电影在线观看一区二区三区| 欧美日韩国产mv在线观看视频| 99热全是精品| 97在线人人人人妻| 亚洲精品一二三| 亚洲人成77777在线视频| 蜜臀久久99精品久久宅男| 26uuu在线亚洲综合色| 久久久国产一区二区| 日本爱情动作片www.在线观看| 99久久精品国产国产毛片| 黄网站色视频无遮挡免费观看| 97在线人人人人妻| 午夜免费男女啪啪视频观看| 亚洲国产日韩一区二区| 精品福利永久在线观看| 婷婷色综合大香蕉| 亚洲综合色惰| 亚洲欧洲精品一区二区精品久久久 | 国产女主播在线喷水免费视频网站| 欧美精品一区二区免费开放| 五月玫瑰六月丁香| 天美传媒精品一区二区| 热99国产精品久久久久久7| 中文字幕亚洲精品专区| 亚洲精华国产精华液的使用体验| 国产成人精品久久久久久| 中文字幕人妻熟女乱码| 2018国产大陆天天弄谢| 啦啦啦视频在线资源免费观看| 一本大道久久a久久精品| 蜜桃国产av成人99| 一级爰片在线观看| 看十八女毛片水多多多| 日产精品乱码卡一卡2卡三| 日韩制服骚丝袜av| 久久97久久精品| 一二三四在线观看免费中文在 | 久久精品人人爽人人爽视色| 国产永久视频网站| 一区二区av电影网| 9色porny在线观看| 国产片特级美女逼逼视频| 亚洲少妇的诱惑av| 亚洲欧美成人精品一区二区| 成人亚洲精品一区在线观看| 久久久久久久久久久久大奶| 精品一区二区三区四区五区乱码 | 亚洲av成人精品一二三区| 97在线视频观看| 精品久久久精品久久久| 欧美丝袜亚洲另类| 赤兔流量卡办理| 最近中文字幕2019免费版| 男人操女人黄网站| 久久久国产一区二区| 一本久久精品| 高清在线视频一区二区三区| 欧美激情 高清一区二区三区| 国产免费现黄频在线看| 久久久久久久久久久久大奶| 国产精品国产三级国产av玫瑰| 伦理电影免费视频| 亚洲色图综合在线观看| a级毛色黄片| 亚洲人成网站在线观看播放| 精品熟女少妇av免费看| 精品国产一区二区久久| 久久综合国产亚洲精品| 日韩一区二区三区影片| 午夜福利乱码中文字幕| 国产成人精品无人区| 七月丁香在线播放| videossex国产| 亚洲av在线观看美女高潮| 精品人妻熟女毛片av久久网站| 久久99精品国语久久久| 大码成人一级视频| 精品卡一卡二卡四卡免费| 久久久久久久大尺度免费视频| 国产成人一区二区在线| a级毛片在线看网站| 国产国拍精品亚洲av在线观看| 热re99久久国产66热| 国产欧美日韩一区二区三区在线| 国产探花极品一区二区| 视频在线观看一区二区三区| 亚洲一级一片aⅴ在线观看| 久久99热这里只频精品6学生| 久久精品夜色国产| 日韩一区二区视频免费看| 全区人妻精品视频| 色婷婷av一区二区三区视频| 日日摸夜夜添夜夜爱| 久久精品久久久久久久性| 亚洲精品一二三| 国产极品天堂在线| 色哟哟·www| 99久国产av精品国产电影| 亚洲少妇的诱惑av| 国产一区二区激情短视频 | 亚洲av免费高清在线观看| 香蕉精品网在线| 99re6热这里在线精品视频| 国产精品偷伦视频观看了| 成人影院久久| 九色亚洲精品在线播放| 日韩欧美一区视频在线观看| 纵有疾风起免费观看全集完整版| 亚洲成国产人片在线观看| 亚洲婷婷狠狠爱综合网| 亚洲国产欧美日韩在线播放| 免费在线观看完整版高清| 大香蕉久久成人网| 在线观看美女被高潮喷水网站| 久久久精品区二区三区| 少妇被粗大的猛进出69影院 | 日韩大片免费观看网站| 2021少妇久久久久久久久久久| 纯流量卡能插随身wifi吗| 高清视频免费观看一区二区| 亚洲av国产av综合av卡| 激情五月婷婷亚洲| 侵犯人妻中文字幕一二三四区| 久久精品aⅴ一区二区三区四区 | 超色免费av| 亚洲欧洲日产国产| 男人爽女人下面视频在线观看| 巨乳人妻的诱惑在线观看| 免费观看av网站的网址| 国产精品国产三级国产av玫瑰| 亚洲av免费高清在线观看| 国产成人精品一,二区| 国产免费福利视频在线观看| 国产国拍精品亚洲av在线观看| 春色校园在线视频观看| 丝袜喷水一区| 亚洲国产精品一区二区三区在线| 成人国产av品久久久| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产最新在线播放| 亚洲第一av免费看| 天天操日日干夜夜撸| 欧美少妇被猛烈插入视频| 少妇 在线观看| 宅男免费午夜| 国产一区二区三区av在线| xxxhd国产人妻xxx| 免费黄网站久久成人精品| 国产白丝娇喘喷水9色精品| 成人亚洲欧美一区二区av| 亚洲精品日韩在线中文字幕| 日韩伦理黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品女同一区二区软件| 好男人视频免费观看在线| 国语对白做爰xxxⅹ性视频网站| 在线 av 中文字幕| 天天躁夜夜躁狠狠躁躁| 精品亚洲成国产av| 人妻 亚洲 视频| 99国产精品免费福利视频| 中文字幕另类日韩欧美亚洲嫩草| 1024视频免费在线观看| 国产 一区精品| 高清av免费在线| 国产1区2区3区精品| 欧美日韩av久久| 久久久久久久国产电影| 精品一品国产午夜福利视频| freevideosex欧美| 久久国产亚洲av麻豆专区| 亚洲欧美色中文字幕在线| a级毛片在线看网站| 国产精品女同一区二区软件| 久久久久久久国产电影| 熟女电影av网| 亚洲欧洲日产国产| 建设人人有责人人尽责人人享有的| 深夜精品福利| 欧美人与性动交α欧美精品济南到 | 亚洲丝袜综合中文字幕| 男女边摸边吃奶| 天美传媒精品一区二区| 亚洲少妇的诱惑av| 亚洲精品乱码久久久久久按摩| 中文字幕最新亚洲高清| 中文字幕制服av| 制服丝袜香蕉在线| 色哟哟·www| 18禁动态无遮挡网站| 欧美日韩一区二区视频在线观看视频在线| 有码 亚洲区| 婷婷成人精品国产| 亚洲av电影在线进入| 少妇的丰满在线观看| 美女国产高潮福利片在线看| 亚洲精华国产精华液的使用体验| 看非洲黑人一级黄片| 大香蕉97超碰在线| 精品一品国产午夜福利视频| 国产精品无大码| 热re99久久国产66热| 18在线观看网站| 久久影院123| 欧美日韩亚洲高清精品| 中文字幕人妻丝袜制服| 校园人妻丝袜中文字幕| 新久久久久国产一级毛片| 免费观看在线日韩| 高清av免费在线| 欧美亚洲日本最大视频资源| 波野结衣二区三区在线| av又黄又爽大尺度在线免费看| 久久人人爽人人片av| 国产精品.久久久| 国产一区亚洲一区在线观看| 韩国高清视频一区二区三区| 成人亚洲精品一区在线观看| 国产精品人妻久久久影院| 一级黄片播放器| 夫妻性生交免费视频一级片| 国产精品免费大片| a级毛色黄片| 99国产综合亚洲精品| 有码 亚洲区| 熟妇人妻不卡中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 欧美少妇被猛烈插入视频| 亚洲精品国产av蜜桃| 边亲边吃奶的免费视频| 久久99热这里只频精品6学生| 乱码一卡2卡4卡精品| 午夜日本视频在线| av卡一久久| 少妇的逼水好多| 久久这里只有精品19| 国产黄色视频一区二区在线观看| 国产探花极品一区二区| 精品午夜福利在线看| 亚洲精品456在线播放app| 狠狠婷婷综合久久久久久88av| 一级毛片电影观看| 熟妇人妻不卡中文字幕| 国产成人精品久久久久久| 久久午夜综合久久蜜桃| 欧美成人午夜免费资源| 伊人亚洲综合成人网| 一区二区三区四区激情视频| 日本欧美视频一区| 欧美精品一区二区免费开放| 一级毛片我不卡| 亚洲,欧美精品.| 色吧在线观看| 一二三四在线观看免费中文在 | 亚洲美女视频黄频| 九色亚洲精品在线播放| 久久精品人人爽人人爽视色| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av综合色区一区| 满18在线观看网站| 黄色毛片三级朝国网站| 久久久欧美国产精品| 一区二区av电影网| 99热全是精品| 日本黄大片高清| 嫩草影院入口| kizo精华| 亚洲一级一片aⅴ在线观看| 超碰97精品在线观看| 天堂中文最新版在线下载| 日本av手机在线免费观看| 男女无遮挡免费网站观看| 极品人妻少妇av视频| 成人二区视频| 久久99一区二区三区| 久久97久久精品| 热99久久久久精品小说推荐| av免费在线看不卡| 久久99热6这里只有精品| 欧美xxⅹ黑人| 久久 成人 亚洲| 热99久久久久精品小说推荐| 日日撸夜夜添| 婷婷成人精品国产| 国产女主播在线喷水免费视频网站| 91精品伊人久久大香线蕉| 亚洲欧美成人精品一区二区| 18禁裸乳无遮挡动漫免费视频| 国产精品久久久久久久久免| 日韩中字成人| 久久久国产一区二区| 欧美成人午夜免费资源| 免费高清在线观看视频在线观看| 一区二区三区乱码不卡18| 久久久久精品人妻al黑| 大码成人一级视频| 18禁动态无遮挡网站| 男女啪啪激烈高潮av片| 国产亚洲精品久久久com| 欧美激情极品国产一区二区三区 | 亚洲欧美色中文字幕在线| 又大又黄又爽视频免费| 亚洲精品日本国产第一区| 女人久久www免费人成看片| 制服人妻中文乱码| 国产xxxxx性猛交| 极品人妻少妇av视频| 久久久久久人人人人人| 蜜臀久久99精品久久宅男| a级毛片黄视频| 天堂中文最新版在线下载| 在线观看免费高清a一片| 精品少妇黑人巨大在线播放| 国国产精品蜜臀av免费| 精品亚洲成国产av| 欧美bdsm另类| a级毛色黄片| 男的添女的下面高潮视频| 午夜福利影视在线免费观看| 亚洲久久久国产精品| 成人影院久久| av在线老鸭窝| 日日撸夜夜添| 成年av动漫网址| 91aial.com中文字幕在线观看| 国产精品99久久99久久久不卡 | 亚洲五月色婷婷综合| 大话2 男鬼变身卡| 一级毛片我不卡| 免费大片黄手机在线观看| 黄片播放在线免费| 日韩成人伦理影院| 免费久久久久久久精品成人欧美视频 | 国产成人aa在线观看| 日本黄大片高清| 最后的刺客免费高清国语| 亚洲成色77777| 中文字幕av电影在线播放| 久久亚洲国产成人精品v| 99精国产麻豆久久婷婷| 看免费av毛片| 人人妻人人爽人人添夜夜欢视频| av天堂久久9| 午夜福利乱码中文字幕| 色哟哟·www| 亚洲三级黄色毛片| 高清av免费在线| 夜夜爽夜夜爽视频| 一本—道久久a久久精品蜜桃钙片| av线在线观看网站| 男人操女人黄网站| 999精品在线视频| 一区二区三区精品91| 激情视频va一区二区三区| 22中文网久久字幕| 日韩 亚洲 欧美在线| 纵有疾风起免费观看全集完整版| 春色校园在线视频观看| 欧美精品国产亚洲| 各种免费的搞黄视频| 亚洲av在线观看美女高潮| 精品国产乱码久久久久久小说| 香蕉精品网在线| 男人操女人黄网站| 国产日韩欧美视频二区| 这个男人来自地球电影免费观看 | 汤姆久久久久久久影院中文字幕| 男人爽女人下面视频在线观看| 日本-黄色视频高清免费观看| 一本久久精品| 国产av码专区亚洲av| 亚洲高清免费不卡视频| 黄色配什么色好看| 校园人妻丝袜中文字幕| 人人妻人人爽人人添夜夜欢视频| 制服人妻中文乱码| 国产精品国产av在线观看| av天堂久久9| 国产伦理片在线播放av一区| 久久午夜综合久久蜜桃| 亚洲,欧美,日韩| 欧美性感艳星| 色5月婷婷丁香| 综合色丁香网| 国产在线视频一区二区| 中国国产av一级| a级片在线免费高清观看视频| 激情五月婷婷亚洲| 亚洲精品一二三| 国产精品女同一区二区软件| 亚洲久久久国产精品| 日本av手机在线免费观看| 久久久久精品人妻al黑| 久久精品久久精品一区二区三区| 亚洲久久久国产精品| 亚洲av电影在线观看一区二区三区| 亚洲精品一二三| 日本爱情动作片www.在线观看| 亚洲久久久国产精品| 日本爱情动作片www.在线观看| 免费在线观看完整版高清| 久久人人爽av亚洲精品天堂| 蜜桃在线观看..| 青春草亚洲视频在线观看| 人人妻人人爽人人添夜夜欢视频| 9色porny在线观看| 在线精品无人区一区二区三| 精品少妇久久久久久888优播| 日韩,欧美,国产一区二区三区| 国产日韩一区二区三区精品不卡| 久久亚洲国产成人精品v| 黑人欧美特级aaaaaa片| 十八禁网站网址无遮挡| 高清不卡的av网站| 亚洲人成77777在线视频| 精品国产国语对白av| 亚洲第一区二区三区不卡| 老司机影院成人| 天天操日日干夜夜撸| 秋霞在线观看毛片| 在线观看免费日韩欧美大片| 国产69精品久久久久777片| 观看av在线不卡| av黄色大香蕉| 女人精品久久久久毛片| 夜夜爽夜夜爽视频| 多毛熟女@视频| 97超碰精品成人国产| 人人妻人人添人人爽欧美一区卜| 黄色一级大片看看| 欧美日韩国产mv在线观看视频| 日本免费在线观看一区| 中文乱码字字幕精品一区二区三区| 男的添女的下面高潮视频| 国产乱来视频区| 日韩中字成人| 十八禁高潮呻吟视频| 亚洲一级一片aⅴ在线观看| 两个人免费观看高清视频| 成人18禁高潮啪啪吃奶动态图| 亚洲精品美女久久久久99蜜臀 | 在线 av 中文字幕| 国产极品天堂在线| 国产欧美日韩一区二区三区在线| 亚洲av成人精品一二三区| 久久久a久久爽久久v久久| 99热网站在线观看| 久久久久久人人人人人| av播播在线观看一区| 免费黄网站久久成人精品| 26uuu在线亚洲综合色| 激情视频va一区二区三区| 1024视频免费在线观看| 亚洲av综合色区一区|