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

    基于多孔介質(zhì)模型的壓水堆堆芯溫場數(shù)值模擬

    2015-12-28 03:40:00田茂林田文喜秋穗正蘇光輝
    核技術(shù) 2015年9期
    關(guān)鍵詞:壓水堆冷卻劑熱工

    陳 森 劉 余 田茂林 田文喜 秋穗正 蘇光輝

    1(西安交通大學 核科學與技術(shù)學院 西安 710049)

    2(中國核動力研究設(shè)計院 核反應堆系統(tǒng)設(shè)計技術(shù)重點實驗室 成都 610041)

    基于多孔介質(zhì)模型的壓水堆堆芯溫場數(shù)值模擬

    陳 森1劉 余2田茂林1田文喜1秋穗正1蘇光輝1

    1(西安交通大學 核科學與技術(shù)學院 西安 710049)

    2(中國核動力研究設(shè)計院 核反應堆系統(tǒng)設(shè)計技術(shù)重點實驗室 成都 610041)

    針對壓水堆的復雜結(jié)構(gòu)特點,對堆芯采用多孔介質(zhì)模型,建立完整的壓力容器堆芯模型,使用商用軟件CFX對壓力容器堆芯的熱工水力特性進行數(shù)值模擬,得到偏環(huán)運行和典型事故工況下冷卻劑的熱工水力響應特性。計算結(jié)果表明:應用多孔介質(zhì)模型能有效正確直觀顯示堆芯的冷卻劑溫度分布情況,在偏環(huán)運行工況下堆芯會出現(xiàn)偏心現(xiàn)象,而通過瞬態(tài)事故工況計算結(jié)果表明堆芯中上部冷卻劑溫度最高,對壓水堆的熱工安全具有一定指導作用。

    壓水堆,堆芯,多孔介質(zhì)模型,瞬態(tài)工況

    隨著計算機性能的提高,計算流體動力學(Computational Fluid Dynamics, CFD)方法在反應堆研究中得到廣泛的應用。楊長江[1]采用CFD軟件對中國先進研究堆(China Advanced Research Reactor, CARR)堆芯的相關(guān)熱工水力研究,結(jié)果與實驗數(shù)據(jù)等比較符合較好,對堆芯的溫度場與流場進行了分析研究,得出的重要結(jié)論用于工程設(shè)計;桂學文等[2]利用CFD商用軟件CFX所提供的用戶接口編程技術(shù),對軟件CFX與三維中子動力學模型耦合這一方法進行了研究,并基于此方法對反應堆局部三維物理特性和流動進行了數(shù)值分析。壓水堆堆芯由很多組件(控制棒、燃料棒等)組成,而每個組件又由數(shù)量從十幾到上百不等的棒束組成,幾何結(jié)構(gòu)相當復雜,若直接模擬全堆芯的流動換熱特性,其網(wǎng)格的劃分極為復雜,此時采用多孔介質(zhì)模型簡化堆芯能有效簡化問題。目前,多孔介質(zhì)模型不僅能正確模擬堆芯冷卻劑的溫度、壓力和速度分布,還能分析冷卻劑的熱工水力特性[3–6]。

    本研究對壓水堆的壓力容器堆芯簡化建模,建立了進口段、下降段、壓力容器下腔室和上腔及堆芯的幾何模型,其中堆芯采用多孔介質(zhì)模型簡化堆芯,上腔室和下腔室設(shè)置阻力源項;然后基于商用軟件CFX多孔介質(zhì)模型計算堆芯冷卻劑的熱工水力特性。本文在穩(wěn)態(tài)的基礎(chǔ)之上進行偏環(huán)運行和蒸汽發(fā)生器傳熱管斷裂事故工況數(shù)值模擬計算,得到瞬態(tài)事故工況下堆芯冷卻劑的三維溫度響應特性。

    1 數(shù)學物理模型

    1.1 CFX多孔介質(zhì)模型介紹

    流動區(qū)域內(nèi)有固體物質(zhì)存在時,對流體的流動有兩方面的影響:一方面是幾何影響,即固體物質(zhì)的存在減少了流體流動的有效空間;另一方面是影響流體的動量和能量的轉(zhuǎn)換和傳遞,在ANSYS CFX多孔介質(zhì)模型[7]中,對標準的流體流動控制方程組中的動量方程增加了額外的源項來表示固體部分對流場的影響。該源項有兩部分組成粘性損失項(Darcy)和慣性損失項:

    式中,Si是某個坐標方向(x,y,z)動量方程中的源項;D和C是需要預先設(shè)定系數(shù)。由于多孔區(qū)域中存在壓力梯度,所形成的壓力降和控制體中的流體速度或平方根速度成比例,這就是多孔介質(zhì)模型動量方程中存在動量匯的原因。由于燃料組件呈有規(guī)律狀均勻排列,式(1)可簡化成:

    式中,α是滲透性參數(shù);C2是慣性阻力因子。與式(1)相比,分別取D為1/α和C為C2。

    1.2幾何模型的建立及網(wǎng)格劃分

    圖1為壓力容器模型示意圖,整個壓力容器堆芯十分復雜,本研究通過一定的簡化,建立了完整的壓力容器模型。由于堆芯采用多孔介質(zhì)模型,故只需畫出整個輪廓模型即可。在網(wǎng)格劃分中,對于壓力容器進口處、堆芯吊籃進行網(wǎng)格加密處理。而對于下腔室、堆芯和上腔室,由于該部分采用簡化方法計算,其網(wǎng)格相比其他部分較稀疏,同時壁面網(wǎng)格采用邊界層加密。本文中的網(wǎng)格在壓力容器上半分布采用四面體網(wǎng)格,其他分布均為六面體網(wǎng)格,圖2和3為幾何模型及模型的網(wǎng)格劃分。

    圖1壓力容器模型Fig.1 Pressure vessel model.

    圖2 壓力容器固體(a)及堆芯吊籃(b)網(wǎng)格Fig.2 Grid of pressure vessel solid (a) and core (b) barrel.

    圖3 下腔室(a)、堆芯(b)和上腔室(c)網(wǎng)格Fig.3 Grid of lower plenum (a), core (b) andupper plenum (c).

    1.3計算參數(shù)設(shè)置

    計算區(qū)域流體設(shè)為不可壓縮湍流,采用k-Epsilon湍流模型,壁面函數(shù)法選擇默認的Scalable。正確設(shè)置孔隙率、阻力系數(shù)等參數(shù)對多孔介質(zhì)模型的應用至關(guān)重要,下面為參數(shù)設(shè)置情況:

    (1) 孔隙率。由于堆芯內(nèi)固體空間和液體空間在軸向上的分布是均勻的,因此體空隙率等于面孔隙率。由堆芯的面積及有效流通面積,參考文獻[8]可計算得到孔隙率為0.492。

    (2) 阻力系數(shù)。在堆芯中由于文獻[8]給出堆芯壓降及堆芯高度,通過計算設(shè)置阻力損失系數(shù),在堆芯多孔介質(zhì)模型中設(shè)置阻力損失模型為定向損失,具體取值參考文獻[8]可得到。同時本研究中橫向阻力系數(shù)采用沿程方向阻力系數(shù)的倍數(shù)關(guān)系,其數(shù)值取為10。

    (3) 在本模型中考慮流固耦合傳熱[9],設(shè)置固體參數(shù)材料物性,同時設(shè)置進口溫度及質(zhì)量流量、堆芯功率、系統(tǒng)壓力邊界條件。

    2 數(shù)值計算結(jié)果分析

    2.1模型驗證

    首先對整個模型穩(wěn)態(tài)計算驗證,得到壓力容器內(nèi)冷卻劑的速度場、溫度場和壓力場。由圖4速度場可知冷卻劑從兩條進口接管流入壓力容器,沿著壓力容器與吊籃之間的環(huán)形空間向下流動,到壓力容器底部轉(zhuǎn)向向上流經(jīng)堆芯,帶出核反應堆放出的熱量,最后從兩條出口管道流出。計算的溫度場、壓強場如圖5和6所示。壓力容器冷卻劑出口溫度588.207K與設(shè)計值588.5K[8]相差很小,驗證了CFX多孔介質(zhì)模型在本研究的可行性和整個模型的可靠性與正確性。

    圖4 流線圖Fig.4 Streamline chart.

    圖5 壓力容器進出口截面溫度分布Fig.5 Temperature profile of pressure vessel inlet-section and outlet-section.

    圖6 壓力容器進出口截面壓力分布Fig.6 Pressure profile of pressure vessel inlet-section andoutlet-section.

    表1為壓力容器各部位壓降計算結(jié)果同設(shè)計值相比較,通過表1可以看出,整個模擬結(jié)果和水力實驗結(jié)果基本吻合,證明整個模型的可靠性。

    表1 壓降計算結(jié)果同實驗值比較Table1 Pressure drop calculation results compared with the experimental values.

    2.2偏環(huán)運行工況計算

    在兩個環(huán)路的核電站運行過程中,當其中的一個環(huán)路出現(xiàn)事故甚至停止運行時,會對整個反應堆的安全性有重要影響。本文針對有兩個環(huán)路的壓水堆核電站,在數(shù)值模擬時對兩個環(huán)路的進口設(shè)置不同的邊界條件,得到堆芯冷卻劑的溫度分布。本文采取不同的溫度邊界條件,即入口1和2設(shè)置不同的溫度進口,其中入口1的溫度為561.65K,而入口2的溫度為417.33K,兩入口質(zhì)量流量相同,計算得到溫度分布如圖7、8所示。

    由于壓力容器結(jié)構(gòu)的不對稱性,冷卻劑通過下降段進入下腔室后進入堆芯之前,冷卻劑在下腔室中沒有充分混合,導致不同溫度的冷卻劑進入堆芯后換熱量不一樣,堆芯溫度出現(xiàn)偏心現(xiàn)象,比較兩入口溫度一樣和入口溫度不一樣的工況下的壓力容器出口溫度平均值,對比結(jié)果如圖9所示。

    圖7 壓力容器進出口中心截面溫度分布Fig.7 Temperature distribution of pressure vessel inlet and outlet center section.

    圖8 堆芯溫度分布Fig.8 Core temperature distribution.

    圖9 不同邊界條件下出口溫度計算結(jié)果Fig.9 Outlet temperature calculation results under different boundary conditions.

    圖9 中曲線1和2代表入口溫度相同的條件下堆芯出口1和2的平均溫度,而曲線3和4代表入口溫度不同的條件下堆芯出口1和2的平均溫度。從圖9可以看出,在入口溫度相同的條件下,兩出口溫度最后基本相同,而不同的溫度進口條件下,堆芯出口平均溫度不同,且整體出口溫度相比入口溫度一樣工況時有明顯下降,說明后者運行時溫度不同的冷卻劑混合不充分,兩個環(huán)路偏環(huán)運行時堆芯會出現(xiàn)偏心現(xiàn)象。

    2.3蒸汽發(fā)生器傳熱管斷裂工況計算

    蒸汽發(fā)生器傳熱管斷裂事故(Steam Generator Tube Rupture, SGTR)是指有與蒸汽發(fā)生器(Steam Generator, SG)傳熱管破裂所造成的冷卻劑喪失率超過冷卻劑不給系統(tǒng)正常補水能力的冷卻劑裝量減少事故,在所有重大事故中,SGTR事故是最頻繁出現(xiàn)的。該瞬態(tài)過程事故邏輯為:100s蒸汽發(fā)生器一根傳熱管發(fā)生雙端剪切斷裂,448s低壓信號(13.0MPa)停堆,停堆信號發(fā)出主泵惰轉(zhuǎn),555s高壓安注在10.78MPa下投入運行,安注箱在4.91MPa下投入運行,低壓安注在0.98MPa下投入運行,停堆信號發(fā)出后給水由主給水切換到輔助給水,停堆信號發(fā)出后汽輪機脫扣。圖10為安注前后的壓力容器堆芯某些時刻的三維溫度場分布。

    圖10 SGTR工況不同時刻壓力容器進出口截面溫度分布Fig.10 Temperature distribution of pressure vessel inlet and outlet center section at different moments under SGTR condition.

    從圖10可以看出,停堆以后,冷卻劑流經(jīng)堆芯帶走堆芯余熱,冷卻劑在堆芯中部出口處溫度達到最高,且由中部到外圍冷卻劑溫度逐漸降低,隨著停堆深度的增加,堆芯余熱被冷卻劑帶出,冷卻劑的最高溫度在出口處所占份額也相應地降低;從出口溫度對比可以看出:發(fā)生SGTR事故后,由于冷卻劑出口溫度有較小溫升,停堆之后,隨著堆芯功率降低,冷卻劑的溫度隨之出現(xiàn)突降,當壓降到達10.78MPa,安全注入系統(tǒng)投入運行,此時堆芯余熱很小,冷卻劑溫度變化緩慢,RELAP5計算結(jié)果同CFX計算結(jié)果變化趨勢一致(圖11)。

    圖11 CFX計算與RELAP5計算結(jié)果對比Fig.11 Calculating results comparison between CFX and RELAP5.

    3 結(jié)語

    本文針對壓水堆堆芯采用多孔介質(zhì)模型,首先計算穩(wěn)態(tài)工況下冷卻劑的溫度分布、速度分布及壓力分布,將計算結(jié)果同實驗值進行對比,驗證了整個模型的可行性。然后針對溫度三維非對稱工況計算,計算結(jié)果表明在入口溫度不同的工況下,堆芯會出現(xiàn)偏心現(xiàn)象,對于堆芯的安全傳熱有重要的指導作用。最后模擬了蒸汽發(fā)生器傳熱管斷裂事故工況下在有安注條件下的堆芯冷卻劑的三維溫度分布情況,計算結(jié)果表明:冷卻劑最高溫度出現(xiàn)在堆芯中部出口處,冷卻劑的溫度分布與堆芯功率分布直接相關(guān)。本文研究結(jié)果對于堆芯的熱工水力安全性具有重要的參考價值。

    1 楊長江. CARR堆芯三維流場與溫場數(shù)值模擬研究[D].北京: 中國原子能科學研究院, 2006

    YANG Changjiang. Numerical simulation research of three-dimensional flow filed and temperature in CARR core[D]. Beijing: China Institute of Atomic Energy, 2006

    2 桂學文, 蔡琦, 陳玉清. 基于CFD的反應堆局部三維流動模型與時空中子動力學模型耦合研究[J]. 核科學與工程, 2010,30(3): 216–222

    GUI Xuewen, CAI Qi, CHEN Yuqing. Research of the coupling local three dimensional flow model with space-time reactor neutron kinetics model based on CFD[J]. Nuclear Science and Engineering, 2010,30(3): 216–222

    3 姚朝暉. 壓水堆核電站中若干熱工水力問題的研究[D].北京: 清華大學, 1995 YAO Chaohui. Research of some thermal hydraulic problems in PWR[D]. Beijing: Tsinghua University, 1995

    4 劉興民, 陸道綱, 劉天才, 等. 中國先進研究堆堆芯流量分配的數(shù)值模擬[J]. 核動力工程, 2003,24(增刊): 21–24 LIU Xingmin, LU Daogang, LIU Tiancai, et al. Numerical simulation of core flow distribution for China advanced reasearch reactor[J]. Nuclear Power Enginnering, 2003,24(Suppl): 21–24

    5 許義軍. 中國實驗快堆鈉池三維熱工水力分析[D]. 北京: 中國原子能科學研究院, 2003 XU Yijun. Sodium pool type three-dimensional thermal-hydraulic analysis of China experimental fast reactor[D]. Beijing: China Institute of Atomic Energy, 2003

    6 陳靜, 田文喜, 韋宏洋, 等. 基于多孔介質(zhì)模型的行波堆TP-1堆芯穩(wěn)態(tài)溫度場與流場數(shù)值模擬[J]. 原子能科學技術(shù), 2013,47(11): 1966–1970. DOI: 10.7538/ yzk.2013.47.11.1966 CHEN Jing, TIAN Wenxi, WEI Hongyang, et al. Numeircal simulation of temperature and flow fields under steady state condition in core for TP-1 sodium cooled TWR based on porous media model[J]. Atomic Energy Science and Technology, 2013,47(11): 1966–1970. DOI: 10.7538/yzk.2013.47.11.1966

    7 孫紀寧. ANSYS CFX 對流傳熱數(shù)值模擬基礎(chǔ)應用教程[M]. 北京: 國防工業(yè)出版社, 2010 SUN Jining. Numerical simulation of convection heat transfer basic application tutorial based on ANSYS CFX[M]. Beijing: National Defence Industry Press, 2010

    8 歐陽予. 秦山核電工程[M]. 北京: 原子能出版社, 2000 OUYANG Yu. Qinshan nuclear power engineering[M]. Beijing: Atomic Energy Press, 2000

    9 姚朝輝, 王學芳, 沈孟育. 堆芯冷卻劑流動和傳熱特性的數(shù)值模擬[J]. 核動力工程, 1997,18(4): 332–339 YAO Chaohui, WANG Xuefang, SHEN Mengyu. Numerical simulation on coolant flow and heat transfer in core[J]. Nuclear Power Enginnering, 1997,18(4): 332–339

    PWR core transient temperature numerical simulation based on porous media model

    CHEN Sen1LIU Yu2TIAN Maolin1TIAN Wenxi1QIU Suizheng1SU Guanghui1

    1(Institute of Nuclear Science and Technology,Xi’an Jiaotong University,Xi’an 710049,China)
    2(Science and Technology on Reactor System Design Technology Laboratory,Nuclear Power Institute of China,Chengdu 610041,China)

    Background:The coolant flow and heat transfer in the reactor are significant to the safety of pressure water reactor (PWR).Purpose:This study aims to obtain the coolant flow and heat transfer characteristics in the reactor core.Methods:The complete model of pressure vessel is built by Pro/E, and the thermal hydraulic characteristics of the core are simulated with commercial software CFX by using the porous medium model in the core.Results:The simulation results show that the application of porous media model can display the core coolant temperature distribution effectively and correctly, and the asymmetric operation conditions can cause asymmetric on reactor core. The transient accident calculation results show that the highest temperature of the coolant appears in the upper of the core.Conclusion:The porous media model can be used in the reactor core and it has a certain reference value for the thermal-hydraulic safety of pressurized water reactor.

    Pressurized water reactor, Reactor core, Porous media, Transient conditions

    TL33

    10.11889/j.0253-3219.2015.hjs.38.090601

    國家自然科學基金(No.11125522)資助

    陳森,男,1987年出生,2012年畢業(yè)于哈爾濱工程大學,現(xiàn)從事反應堆熱工水力學研究

    田文喜,E-mail: wxtian@mail.xjtu.edu.cn

    2014-07-29,

    2015-01-14

    CLCTL33

    猜你喜歡
    壓水堆冷卻劑熱工
    核電站主冷卻劑泵可取出部件一體化吊裝檢修工藝探索
    溶氫表在壓水堆核電站的應用
    壓水堆核電站α輻射的測量及防護
    熱工儀表自動化安裝探討的認識
    智能控制在電廠熱工自動化中的應用
    壓水堆核電站嚴重事故下移動泵的快速響應
    中國核電(2017年2期)2017-08-11 08:01:27
    反應堆冷卻劑pH對核電廠安全運行影響研究
    中國核電(2017年1期)2017-05-17 06:10:13
    冷卻劑泄漏監(jiān)測系統(tǒng)在核電廠的應用
    中國核電(2017年1期)2017-05-17 06:10:05
    智能控制在電廠熱工自動化中的應用
    小型壓水堆嚴重事故序列的篩選及模擬分析研究
    精品亚洲乱码少妇综合久久| 亚洲av男天堂| 一本色道久久久久久精品综合| 最黄视频免费看| 国产成人免费观看mmmm| 国产国拍精品亚洲av在线观看| 午夜激情久久久久久久| 免费不卡的大黄色大毛片视频在线观看| 午夜久久久在线观看| 一级毛片电影观看| 日日摸夜夜添夜夜爱| 久久午夜福利片| 午夜福利视频精品| 在线观看三级黄色| 欧美xxxx性猛交bbbb| 国产精品女同一区二区软件| 七月丁香在线播放| 色婷婷av一区二区三区视频| 插逼视频在线观看| 亚洲人成网站在线观看播放| 韩国高清视频一区二区三区| 国产亚洲最大av| 国产视频首页在线观看| 亚洲综合色网址| 国产极品天堂在线| 人成视频在线观看免费观看| 亚洲经典国产精华液单| 91精品三级在线观看| h视频一区二区三区| 97在线人人人人妻| 久久久久久久国产电影| 中国国产av一级| 成人毛片a级毛片在线播放| 熟女av电影| 久久精品国产亚洲av涩爱| 成年人免费黄色播放视频| 波野结衣二区三区在线| 18在线观看网站| 国产一区二区三区av在线| 91国产中文字幕| 亚洲国产av影院在线观看| 成人漫画全彩无遮挡| 下体分泌物呈黄色| 国国产精品蜜臀av免费| 免费av中文字幕在线| 18禁裸乳无遮挡动漫免费视频| 69精品国产乱码久久久| 中文欧美无线码| 日本免费在线观看一区| 色网站视频免费| 交换朋友夫妻互换小说| 五月玫瑰六月丁香| 最近手机中文字幕大全| 亚洲国产色片| 国产免费又黄又爽又色| 国产成人freesex在线| 狠狠婷婷综合久久久久久88av| 91午夜精品亚洲一区二区三区| 国产精品无大码| 日韩不卡一区二区三区视频在线| 黑人欧美特级aaaaaa片| 美女cb高潮喷水在线观看| 亚洲成人一二三区av| 国产亚洲av片在线观看秒播厂| 97超碰精品成人国产| 国产av精品麻豆| 亚洲天堂av无毛| a级毛片免费高清观看在线播放| 久久97久久精品| 午夜激情av网站| 老女人水多毛片| 色94色欧美一区二区| 亚洲精品乱码久久久久久按摩| 成人国产av品久久久| www.色视频.com| 久久青草综合色| 一边摸一边做爽爽视频免费| 欧美精品亚洲一区二区| 久久久久国产精品人妻一区二区| 中文天堂在线官网| 亚洲精品久久久久久婷婷小说| 性高湖久久久久久久久免费观看| 午夜精品国产一区二区电影| 99热这里只有是精品在线观看| 亚洲国产精品成人久久小说| 国模一区二区三区四区视频| 日韩欧美精品免费久久| 黄色毛片三级朝国网站| 在线观看免费高清a一片| 久久久久精品性色| 久久这里有精品视频免费| 国产男女超爽视频在线观看| 日韩精品有码人妻一区| 最近中文字幕高清免费大全6| 久久久久网色| 成人综合一区亚洲| 久久鲁丝午夜福利片| 国产免费福利视频在线观看| 这个男人来自地球电影免费观看 | 卡戴珊不雅视频在线播放| 久久久久精品性色| 免费少妇av软件| 久久女婷五月综合色啪小说| 亚洲精品久久久久久婷婷小说| 成人国语在线视频| 亚洲欧美清纯卡通| 久久久久久久亚洲中文字幕| 色94色欧美一区二区| 国产黄频视频在线观看| 亚洲精品国产色婷婷电影| 插阴视频在线观看视频| 久热久热在线精品观看| 男女无遮挡免费网站观看| 国产成人午夜福利电影在线观看| 天天躁夜夜躁狠狠久久av| 精品一区二区免费观看| 欧美日韩视频高清一区二区三区二| 国产精品熟女久久久久浪| 18禁动态无遮挡网站| 国产又色又爽无遮挡免| 国产免费一级a男人的天堂| 精品一区二区免费观看| 丝袜美足系列| 亚洲不卡免费看| 久久av网站| 国产黄频视频在线观看| 在线观看美女被高潮喷水网站| 美女国产高潮福利片在线看| 精品久久久噜噜| 国产成人精品久久久久久| 亚洲国产精品专区欧美| 婷婷色综合www| 免费黄频网站在线观看国产| 肉色欧美久久久久久久蜜桃| 桃花免费在线播放| 精品少妇内射三级| 99热这里只有精品一区| 伦理电影大哥的女人| 中文欧美无线码| 观看美女的网站| 久久精品国产亚洲网站| 亚洲精华国产精华液的使用体验| 18在线观看网站| 男女高潮啪啪啪动态图| 在线免费观看不下载黄p国产| 国产日韩欧美亚洲二区| 插逼视频在线观看| 国国产精品蜜臀av免费| 简卡轻食公司| 一级爰片在线观看| 日本av免费视频播放| 亚洲精品国产色婷婷电影| 99久久精品国产国产毛片| 久久精品国产亚洲av天美| 亚洲在久久综合| 91国产中文字幕| 青春草亚洲视频在线观看| 亚洲精品日本国产第一区| 国产在线一区二区三区精| 考比视频在线观看| 中文欧美无线码| 国产精品女同一区二区软件| 国产精品国产三级国产专区5o| tube8黄色片| 一个人免费看片子| 久久国产精品大桥未久av| 欧美最新免费一区二区三区| 欧美精品一区二区免费开放| 亚洲情色 制服丝袜| 曰老女人黄片| 赤兔流量卡办理| 18在线观看网站| 综合色丁香网| 九九爱精品视频在线观看| 国产一区二区三区综合在线观看 | 又大又黄又爽视频免费| 精品卡一卡二卡四卡免费| 999精品在线视频| 岛国毛片在线播放| 女人精品久久久久毛片| 成年人午夜在线观看视频| 国产高清三级在线| 国产黄色免费在线视频| 午夜久久久在线观看| 日韩亚洲欧美综合| 国产午夜精品一二区理论片| av网站免费在线观看视频| 亚洲怡红院男人天堂| 久久人人爽av亚洲精品天堂| 成人二区视频| 久久国产精品大桥未久av| 视频区图区小说| 日韩欧美精品免费久久| 男女无遮挡免费网站观看| 亚洲欧美一区二区三区黑人 | 国产亚洲av片在线观看秒播厂| 老司机影院毛片| 亚洲av男天堂| a 毛片基地| 免费看不卡的av| 精品国产乱码久久久久久小说| 天美传媒精品一区二区| 中文字幕免费在线视频6| 国产白丝娇喘喷水9色精品| 观看av在线不卡| 毛片一级片免费看久久久久| 高清毛片免费看| 精品国产一区二区三区久久久樱花| 欧美人与性动交α欧美精品济南到 | 中文天堂在线官网| 中文字幕av电影在线播放| 亚洲精品日本国产第一区| 少妇精品久久久久久久| 国产有黄有色有爽视频| 亚洲天堂av无毛| 插阴视频在线观看视频| 涩涩av久久男人的天堂| 91久久精品电影网| 久久青草综合色| 国产精品久久久久久精品古装| 欧美日本中文国产一区发布| 亚洲一区二区三区欧美精品| 丰满迷人的少妇在线观看| 最近2019中文字幕mv第一页| 中文字幕久久专区| 国产极品天堂在线| 性色avwww在线观看| 人妻少妇偷人精品九色| 少妇精品久久久久久久| 久久这里有精品视频免费| 狠狠婷婷综合久久久久久88av| 一级毛片 在线播放| 国产免费一区二区三区四区乱码| 天天影视国产精品| 国产精品偷伦视频观看了| 国产综合精华液| 蜜桃久久精品国产亚洲av| 午夜福利,免费看| 亚洲精品日韩av片在线观看| 久久韩国三级中文字幕| 国产精品麻豆人妻色哟哟久久| 国产精品99久久久久久久久| 肉色欧美久久久久久久蜜桃| 国产精品熟女久久久久浪| 欧美日韩视频高清一区二区三区二| 一区二区三区免费毛片| 熟女人妻精品中文字幕| 亚洲av成人精品一区久久| 精品少妇黑人巨大在线播放| 久久免费观看电影| 男女边吃奶边做爰视频| av电影中文网址| 老熟女久久久| 日韩,欧美,国产一区二区三区| 天堂中文最新版在线下载| 亚洲精品自拍成人| 午夜久久久在线观看| 春色校园在线视频观看| 春色校园在线视频观看| 久久精品国产亚洲网站| 久热久热在线精品观看| 一个人免费看片子| 日日摸夜夜添夜夜添av毛片| 国产一区二区三区综合在线观看 | 五月天丁香电影| 国产日韩欧美视频二区| av线在线观看网站| 蜜桃在线观看..| 久久久亚洲精品成人影院| 制服丝袜香蕉在线| 久久热精品热| 日本wwww免费看| 欧美性感艳星| 国模一区二区三区四区视频| 亚洲婷婷狠狠爱综合网| 日日啪夜夜爽| 热99久久久久精品小说推荐| 久久久a久久爽久久v久久| 亚洲成色77777| 日韩熟女老妇一区二区性免费视频| 国产一区二区三区综合在线观看 | 亚洲综合精品二区| 99热全是精品| 大又大粗又爽又黄少妇毛片口| 久久毛片免费看一区二区三区| 欧美日本中文国产一区发布| 亚洲精品日韩在线中文字幕| 一区在线观看完整版| 18禁裸乳无遮挡动漫免费视频| 日本猛色少妇xxxxx猛交久久| 最近最新中文字幕免费大全7| 午夜久久久在线观看| 在线观看美女被高潮喷水网站| av电影中文网址| 日本91视频免费播放| av黄色大香蕉| 久久av网站| av视频免费观看在线观看| 国产精品.久久久| 永久网站在线| 亚洲一区二区三区欧美精品| 最黄视频免费看| 色视频在线一区二区三区| 亚洲高清免费不卡视频| 亚洲欧美一区二区三区国产| 十八禁高潮呻吟视频| 另类亚洲欧美激情| 亚洲欧美一区二区三区黑人 | 午夜激情久久久久久久| av在线播放精品| 国产在视频线精品| 如何舔出高潮| 80岁老熟妇乱子伦牲交| 国产在线免费精品| 永久免费av网站大全| 久久久午夜欧美精品| 久久久久久伊人网av| 国产亚洲最大av| 国产精品女同一区二区软件| 久久久久视频综合| 丰满乱子伦码专区| 熟女人妻精品中文字幕| 2021少妇久久久久久久久久久| 免费观看的影片在线观看| 少妇人妻精品综合一区二区| 亚洲精品456在线播放app| 国产免费又黄又爽又色| 久久这里有精品视频免费| 亚洲欧美一区二区三区黑人 | 国产黄色免费在线视频| 69精品国产乱码久久久| 亚洲av福利一区| 精品一区二区三区视频在线| 免费看av在线观看网站| 少妇被粗大的猛进出69影院 | 亚洲av电影在线观看一区二区三区| 丰满少妇做爰视频| 国产亚洲av片在线观看秒播厂| 久久久欧美国产精品| 免费大片黄手机在线观看| 最黄视频免费看| 有码 亚洲区| 美女国产高潮福利片在线看| 国产高清三级在线| 永久网站在线| 亚洲精品自拍成人| 我的女老师完整版在线观看| 国产黄片视频在线免费观看| 简卡轻食公司| 日韩电影二区| 母亲3免费完整高清在线观看 | 亚洲精品乱码久久久v下载方式| 纵有疾风起免费观看全集完整版| 欧美精品国产亚洲| 啦啦啦中文免费视频观看日本| 国产 精品1| 飞空精品影院首页| 亚洲人成网站在线观看播放| 日本免费在线观看一区| 亚洲欧洲国产日韩| 国产老妇伦熟女老妇高清| 久久精品熟女亚洲av麻豆精品| 肉色欧美久久久久久久蜜桃| 少妇被粗大猛烈的视频| 9色porny在线观看| 狠狠精品人妻久久久久久综合| 熟女av电影| 美女主播在线视频| 一级毛片我不卡| 精品视频人人做人人爽| 精品酒店卫生间| 熟女电影av网| 九九久久精品国产亚洲av麻豆| 精品少妇久久久久久888优播| 美女内射精品一级片tv| 99热这里只有是精品在线观看| 精品人妻在线不人妻| 亚洲欧美精品自产自拍| 色婷婷av一区二区三区视频| 亚洲国产最新在线播放| 人妻系列 视频| 亚洲在久久综合| 99热国产这里只有精品6| 成人18禁高潮啪啪吃奶动态图 | 国产一区亚洲一区在线观看| 91在线精品国自产拍蜜月| av卡一久久| 久久久a久久爽久久v久久| 成人亚洲精品一区在线观看| 26uuu在线亚洲综合色| av在线播放精品| 久久久久久人妻| av.在线天堂| 激情五月婷婷亚洲| 简卡轻食公司| 国产精品99久久久久久久久| 街头女战士在线观看网站| 有码 亚洲区| 亚洲精品第二区| 欧美激情国产日韩精品一区| 成人黄色视频免费在线看| 国产片内射在线| 国产深夜福利视频在线观看| 国产男女超爽视频在线观看| 久久久久精品久久久久真实原创| 在线观看三级黄色| 国产免费一区二区三区四区乱码| 亚洲图色成人| 欧美精品人与动牲交sv欧美| 午夜久久久在线观看| 国产精品一区www在线观看| 国产乱人偷精品视频| 99视频精品全部免费 在线| 国产成人免费观看mmmm| 一区二区三区免费毛片| 亚洲国产精品国产精品| 熟女电影av网| 欧美人与性动交α欧美精品济南到 | 桃花免费在线播放| 2021少妇久久久久久久久久久| 国产精品麻豆人妻色哟哟久久| 日本黄大片高清| 男女高潮啪啪啪动态图| 亚洲高清免费不卡视频| 晚上一个人看的免费电影| 免费看不卡的av| 日本黄大片高清| 日韩av免费高清视频| 午夜激情av网站| 亚洲欧美成人精品一区二区| freevideosex欧美| 人人妻人人添人人爽欧美一区卜| 国产精品久久久久久久久免| 伦理电影大哥的女人| 欧美人与性动交α欧美精品济南到 | 亚洲精品国产av成人精品| 午夜福利视频在线观看免费| 午夜福利网站1000一区二区三区| 国产成人精品福利久久| 天天操日日干夜夜撸| 在线观看美女被高潮喷水网站| 中文字幕久久专区| 26uuu在线亚洲综合色| 国产精品熟女久久久久浪| 伦理电影大哥的女人| 中文字幕最新亚洲高清| 免费黄频网站在线观看国产| 3wmmmm亚洲av在线观看| 精品一区二区三区视频在线| 久久鲁丝午夜福利片| 免费久久久久久久精品成人欧美视频 | 精品酒店卫生间| 亚洲av电影在线观看一区二区三区| 中文字幕免费在线视频6| 在线看a的网站| 最后的刺客免费高清国语| 国产国语露脸激情在线看| 黑人欧美特级aaaaaa片| 少妇精品久久久久久久| 亚洲精品日韩av片在线观看| 少妇人妻精品综合一区二区| 亚洲国产精品一区二区三区在线| 免费观看无遮挡的男女| a级片在线免费高清观看视频| av免费在线看不卡| 丝袜在线中文字幕| 18在线观看网站| 欧美xxxx性猛交bbbb| 成年美女黄网站色视频大全免费 | 五月开心婷婷网| 美女xxoo啪啪120秒动态图| 乱码一卡2卡4卡精品| 欧美精品一区二区免费开放| 2021少妇久久久久久久久久久| 看免费成人av毛片| 国产片内射在线| 久久狼人影院| 99热网站在线观看| 新久久久久国产一级毛片| 久久久亚洲精品成人影院| 亚洲精品日本国产第一区| 国产成人精品在线电影| 欧美另类一区| 成人午夜精彩视频在线观看| 亚洲av.av天堂| 涩涩av久久男人的天堂| 亚洲天堂av无毛| 久久久精品区二区三区| av网站免费在线观看视频| 最近手机中文字幕大全| www.av在线官网国产| 丰满迷人的少妇在线观看| 久久国产精品男人的天堂亚洲 | 国产欧美日韩一区二区三区在线 | 亚洲精品第二区| 永久网站在线| 欧美 亚洲 国产 日韩一| 人妻制服诱惑在线中文字幕| 成人漫画全彩无遮挡| 三上悠亚av全集在线观看| 天天躁夜夜躁狠狠久久av| 日韩免费高清中文字幕av| www.av在线官网国产| 男女国产视频网站| 免费人成在线观看视频色| 成人国产av品久久久| 日韩大片免费观看网站| 国产欧美日韩一区二区三区在线 | 久久精品夜色国产| 男女边摸边吃奶| 久久人人爽人人片av| 亚洲国产精品一区二区三区在线| 国产高清有码在线观看视频| av专区在线播放| 久久影院123| 一个人看视频在线观看www免费| 男女边吃奶边做爰视频| 国产亚洲av片在线观看秒播厂| 少妇人妻久久综合中文| 丝瓜视频免费看黄片| 26uuu在线亚洲综合色| 99热这里只有是精品在线观看| 美女大奶头黄色视频| 插阴视频在线观看视频| av有码第一页| 91精品国产国语对白视频| 欧美激情 高清一区二区三区| 亚洲国产欧美在线一区| 午夜影院在线不卡| 亚洲欧洲日产国产| 精品久久久噜噜| 一级片'在线观看视频| 国产成人精品一,二区| 国产精品国产三级国产av玫瑰| 国产在线一区二区三区精| 最新中文字幕久久久久| 亚洲国产欧美日韩在线播放| 一区二区三区精品91| 在线观看www视频免费| 18在线观看网站| 国产精品国产三级国产专区5o| 交换朋友夫妻互换小说| 热99国产精品久久久久久7| 美女脱内裤让男人舔精品视频| 一区二区三区四区激情视频| 人妻 亚洲 视频| 国产精品一区二区在线不卡| 国产日韩欧美在线精品| av免费在线看不卡| 国产一区有黄有色的免费视频| 日韩制服骚丝袜av| 内地一区二区视频在线| av播播在线观看一区| 美女cb高潮喷水在线观看| 十八禁网站网址无遮挡| 日日爽夜夜爽网站| 波野结衣二区三区在线| 日日爽夜夜爽网站| 日韩熟女老妇一区二区性免费视频| 中国国产av一级| 热99久久久久精品小说推荐| 国产精品三级大全| 秋霞伦理黄片| 成人影院久久| a级片在线免费高清观看视频| √禁漫天堂资源中文www| 又黄又爽又刺激的免费视频.| 少妇人妻 视频| 丁香六月天网| 男女高潮啪啪啪动态图| 亚洲精品一二三| 亚洲色图综合在线观看| 欧美日本中文国产一区发布| 国产成人av激情在线播放 | 制服诱惑二区| 亚洲av不卡在线观看| 国产熟女午夜一区二区三区 | 久久精品国产自在天天线| 日本欧美视频一区| 久久国产亚洲av麻豆专区| 边亲边吃奶的免费视频| 18禁动态无遮挡网站| 久久人人爽人人爽人人片va| 免费看光身美女| 亚洲欧美一区二区三区国产| av播播在线观看一区| 91精品国产九色| 我的女老师完整版在线观看| 国产在线一区二区三区精| 国产国语露脸激情在线看| 国产欧美日韩一区二区三区在线 | 涩涩av久久男人的天堂| 精品久久国产蜜桃| 黄色配什么色好看| 男人爽女人下面视频在线观看| 蜜臀久久99精品久久宅男| 嘟嘟电影网在线观看| 草草在线视频免费看| 国产精品麻豆人妻色哟哟久久| 欧美国产精品一级二级三级| 亚洲人成网站在线观看播放| 高清av免费在线| 国产成人精品福利久久| 最近手机中文字幕大全| 天天操日日干夜夜撸| 一级黄片播放器| 草草在线视频免费看| 国产免费又黄又爽又色| 国产高清国产精品国产三级| 免费大片黄手机在线观看| 亚洲精品国产av成人精品| 男女免费视频国产| 午夜激情福利司机影院| 黄片无遮挡物在线观看| 久久精品国产亚洲av天美|