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

    基于WIMS格式多群核數(shù)據(jù)制作過程的優(yōu)化分析

    2015-12-02 03:39:26謝明亮陳玉清
    核技術(shù) 2015年11期
    關(guān)鍵詞:共振誤差加工

    謝明亮 陳玉清 于 雷 時 浩

    (海軍工程大學(xué) 核能科學(xué)與工程系 武漢 430033)

    基于WIMS格式多群核數(shù)據(jù)制作過程的優(yōu)化分析

    謝明亮 陳玉清 于 雷 時 浩

    (海軍工程大學(xué) 核能科學(xué)與工程系 武漢 430033)

    基于最新釋放的ENDF/B-VII.1核評價庫,采用核數(shù)據(jù)加工處理程序NJOY-99制作基于WIMS格式的多

    ENDF/B-VII.1,WIMS,NJOY-99,核數(shù)據(jù)庫

    NJOY[1]是國際上廣泛應(yīng)用的核數(shù)據(jù)加工處理程序,可以將ENDF/B格式核數(shù)據(jù)庫轉(zhuǎn)換為WIMS格式的多群截面庫,從而為燃料組件均勻化計算提供可直接應(yīng)用的核截面數(shù)據(jù)。ENDF/B庫包含有核反應(yīng)堆核設(shè)計所需各種材料的核素數(shù)據(jù),庫更新速度較快,新的核數(shù)據(jù)庫在評價核、多溫度點下群截面與共振處理等方面不斷改進(jìn),數(shù)據(jù)精度不斷提高,目前已發(fā)布至ENDF/B-VII.1[2]。為從源頭上減少核設(shè)計分析的誤差,本文基于最新釋放的程序輸入模塊參數(shù)變量開展了廣泛的優(yōu)化分析,制作了基于WIMS格式的多群截面數(shù)據(jù)庫。

    1 計算模型建立

    1.1 NJOY程序數(shù)據(jù)處理模塊簡介

    NJOY程序包含一系列數(shù)據(jù)處理模塊[3],每個模塊都有它特定的功能,每個模塊都是分開的計算程序,常用數(shù)據(jù)處理模塊的基本功能如表1所示。ENDF/B-VII.1評價庫,以輕水堆(Light Water Reactor, LWR)基本燃料柵元基準(zhǔn)題為分析對象,對NJOY-99

    表1 NJOY程序常用模塊Table 1 Commonly used modules in NJOY procedures.

    (續(xù)表1)

    1.2 WIMS格式多群數(shù)據(jù)庫

    WIMS格式多群數(shù)據(jù)庫廣泛應(yīng)用于輕水堆各類燃料組件的均勻化計算。數(shù)據(jù)庫一般包含吸收截面、裂變截面、ν·fission截面、輸運(yùn)及散射截面等。對共振吸收體,數(shù)據(jù)庫給出屏蔽共振積分以及本底截面和溫度的關(guān)系,同時還包含裂變產(chǎn)物產(chǎn)額及衰變常數(shù)等。采用NJOY程序制作核素WIMS格式核截面數(shù)據(jù)需要的處理模塊主要有:MODER、RECONR、BROADR、UNRESR、PURR、THERMR、GROUPR、WIMSR[4],制作的基本流程如圖1所示。

    圖1 WIMS格式數(shù)據(jù)生成簡圖Fig.1 WIMS format data processing sequence.

    1.3 模型建立及基準(zhǔn)題驗證分析

    使用NJOY程序生成核素WIMS格式數(shù)據(jù)時,需要經(jīng)過多個處理模塊,處理模塊參數(shù)選擇合理與否,將影響數(shù)據(jù)庫制作的精度與效率。本文以LWR基本燃料柵元均勻化基準(zhǔn)題[5]為應(yīng)用研究對象,該基準(zhǔn)題研究的是一個無限長燃料柵元零燃耗的兩種狀態(tài),即熱態(tài)零功率(Hot zero power, HZP)和熱態(tài)滿功率(Hot full power, HFP),兩種狀態(tài)下慢化劑溫度均為600 K,HZP狀態(tài)的包殼和燃料溫度為600 K,HFP狀態(tài)下的包殼及燃料溫度為900 K。該基準(zhǔn)題UO2燃料包含5類富集度,分別對應(yīng)為0.711%、1.6%、2.4%、3.1%和3.9%。該基準(zhǔn)題最初是由Rahnema等[5]提出,主要是驗證反應(yīng)性多普勒系數(shù)計算能力的柵元基準(zhǔn)例題,其幾何描述如圖2所示。

    針對上述基準(zhǔn)例題,重點以235U、238U核數(shù)據(jù)制作過程為優(yōu)化分析對象,通過更改多群數(shù)據(jù)庫制作過程不同處理模塊的特性輸入?yún)?shù),得到不同輸入變量下的核素多群截面數(shù)據(jù),將制作完成的多群數(shù)據(jù)庫導(dǎo)入組件均勻化程序計算。通過廣泛的對比分析,實現(xiàn)數(shù)據(jù)庫制作過程的優(yōu)化,其計算流程如圖3所示。

    圖2 柵元基準(zhǔn)題幾何描述Fig.2 Geometrical description of cell benchmark.

    圖3 計算流程示意圖Fig.3 Flow diagram of calculation.

    2 數(shù)據(jù)處理模塊輸入?yún)?shù)的影響分析

    對于一個核裝置(包括堆芯或?qū)嶒炑b置)可以定義一個宏觀特征參數(shù)R,它隨著變量q(可以是共振重造誤差、勢散射截面、參考邦達(dá)連柯本底截面、裂變譜、權(quán)重譜等)的函數(shù),通過改變變量q從而引起變量P的變化[6],可以表示為:

    式中,ε表示輸入?yún)?shù)引起變量q改變的矢量或微元,因此多群核數(shù)據(jù)庫制作過程靈敏度分析主要根據(jù)輸入模塊的變量控制和對組件有效增殖因子Keff的影響。為避免單一算例對輸出結(jié)果的偶然影響,設(shè)定ΔKeff?Ave=(ΔKeff1+ΔKeff2+…+ΔKeffN)/N作為結(jié)果

    優(yōu)劣判據(jù),其中ΔKeffN= (KeffN(Cal.)?KeffN(Exp.))/KeffN(Exp.),因該燃料柵元基準(zhǔn)題包含兩種工況(HZP、HFP)和5類不同富集度,因此需要計算每個對應(yīng)算例的ΔKeffN(即選取N=10)與ΔKeff?Ave,改變輸入模塊對應(yīng)的參數(shù),重復(fù)上述過程,得到相應(yīng)的ΔKeff?Ave(下文簡寫ΔKeff)值。

    2.1 RECONR共振重造誤差限影響分析

    共振重造模塊RECONR通過ENDF中的共振參數(shù)和插值方法重構(gòu)能量相關(guān)的點截面,在重造過程中改變誤差限,從而得到關(guān)于不同誤差限下的輸出文件,并重復(fù)上述模型進(jìn)行計算,以238U截面庫制作加工時間為例,加工過程中選取6個溫度點(293 K、600 K、900 K、1200 K、1600 K、2000 K),表2給出加工時間及ΔKeff隨RECONR誤差限的影響情況。

    表2 NJOY制作截面加工時間及ΔKeff隨RECONR誤差限的影響Table 2 Results of processing time and ΔKeffusing different errors in RECONR.

    從表2可以看出,隨RECONR誤差限的增加,MODER、PURR、THERMR、WIMSR模塊加工時間基本保持不變,而RECONR、BROADR、GROUPR模塊的加工時間逐漸減少,總的加工時間呈現(xiàn)下降的趨勢,且誤差限較小時影響較大,但在0.05附近以后加工時間下降不明顯;當(dāng)RECONR誤差限在0.005附近以下時,對于有效增殖系數(shù)計算結(jié)果沒有明顯影響,但隨著誤差限的增加,ΔKeff計算結(jié)果逐漸發(fā)散,因此在數(shù)據(jù)處理過程中,RECONR誤差限選取控制在0.005內(nèi),得到的計算結(jié)果更加接近真實值。

    2.2 BROADR誤差限影響分析

    BROADR模塊用于處理多普勒(Doppler)展寬[7],讀取共振重造PENDF點截面文件,固定其余模塊參數(shù),改變BROADR模塊輸入的誤差限,重復(fù)上述計算過程,表3給出加工時間及ΔKeff隨BROADR誤差限的影響情況。

    表3 NJOY制作截面加工時間及ΔKeff隨BROADR誤差限的影響Table 3 Results of processing time and ΔKeffusing different errors in BROADR.

    可以看出,隨著BROADR誤差限的逐漸增大,除了BROADR、GROUPR模塊加工時間逐漸減少外,其他輸入模塊加工時間基本保持不變,截面總的加工時間呈現(xiàn)下降的趨勢,但影響相對不是很大;選取BROADR誤差限在0.01附近以下,對有效增殖系數(shù)無明顯影響,具有較好的計算精度,但隨BROADR誤差限的增加,計算結(jié)果ΔKeff偏差變大。

    2.3 更改PURR共振級數(shù)及截面σ0影響分析

    采用PURR模塊處理共振自屏數(shù)據(jù),雖然速度較UNRESR慢,但計算精度較前者高,表4給出隨PURR共振級數(shù)變化、數(shù)據(jù)加工時間及ΔKeff的影響情況,可以看出隨著PURR模塊處理不可分辨共振自屏效應(yīng)的共振級數(shù)逐級增加,PURR加工時間逐漸增加,而其他模塊加工時間基本不變,ΔKeff影響偏差逐漸減少,計算精度逐漸增加;當(dāng)共振級數(shù)選擇100?1000時,對計算結(jié)果ΔKeff沒有明顯影響,但需要很長的加工時間(92.2?915.6 s),因此選取PURR共振級數(shù)100即可達(dá)到精度要求,同時保證了較快的計算速度。

    表4 加工時間及ΔKeff隨PURR共振級數(shù)的影響Table 4 Results of processing time and ΔKeffusing different resonance ladders in PURR.

    不同的核素處理選取不同的本底截面σ0值,通常根據(jù)核素維數(shù)及組成來進(jìn)行估計,常采用下面公式計算柵元的σ0值[8]:

    式中,NA為核素i的核子密度;N(t)為同一區(qū)核素i的核密度;λg(i)為核素第g群對應(yīng)的參數(shù)λG值;σp(i)為核素i的勢散射截面;α為對應(yīng)的Bell因子;l為對應(yīng)典型柵元圓柱直徑。改變本底截面σ0值,選取4組不同的σ0參考值并保持其他模塊輸入?yún)?shù)不變,可以得到對加工時間及ΔKeff的影響情況(表5)。由表5可以看出,對于4組不同截面σ0值,加工時間基本不發(fā)生變化,但對ΔKeff值影響較大,因此截面σ0值的選取對積分量Keff有較大的靈敏度。

    表5 加工時間及ΔKeff隨截面σ0的影響Table 5 Results of processing time and ΔKeffusing different cross section σ0.

    2.4 THERMR模塊熱化能量上限影響分析

    在進(jìn)行WIMS格式數(shù)據(jù)的制作過程中,中子熱化處理THERMR根據(jù)不同的能群結(jié)構(gòu),選取對應(yīng)的熱化能量上限,以基于WIMS格式的CASMO-70[9]群為例,熱化上限能量選取4.0 eV,保持其他模塊輸入?yún)?shù)不變,表6給出不同的熱化能量上限對加工時間及ΔKeff影響??梢钥闯鯰HERMR模塊熱化能量上限的變化對加工時間影響不大,對于偏離熱能群的熱化能量上限對積分量Keff有少量的影響。

    表6 加工時間及ΔKeff隨熱化能量上限的影響Table 6 Results of processing time and ΔKeffusing different upper boundary of thermal energy range.

    2.5 權(quán)重譜選擇影響分析

    在采用GROUPR模塊產(chǎn)生多群截面過程中,權(quán)重函數(shù)在制作能群常數(shù)之前通常是未知的,因此需要選擇合適的權(quán)重譜適用特定的裝置。表7給出常見的4種權(quán)重譜下,截面制作加工時間及ΔKeff的變化情況。由表7可見,選擇不同的權(quán)重譜對積分量Keff有較大的影響,權(quán)重譜選擇l/E+FISSION譜+MAXWELLIAN譜時,較其他譜而言需要更長的加工時間,因此權(quán)重譜的選擇具有較大的靈敏度,制作數(shù)據(jù)截面需要合理選取。IWT1:CONSTANT;IWT2:l/E;IWT3:l/E + FISSION SPECTRU+ MAXWELLIAN SPECTRU;IWT4:EPRI-CELL LWR SPECTRU

    表7 加工時間及ΔKeff隨權(quán)重譜選取的影響Table 7 Results of processing time and ΔKeffusing different weight spectrum.

    2.6 Goldstein-Cohen參數(shù)λG選擇影響分析

    通過WIMSR模塊產(chǎn)生核素WIMS格式數(shù)據(jù)過

    程中,需要選取合適的Goldstein-Cohen參數(shù)λG,λG

    值通常由特定柵元結(jié)構(gòu)、溫度等參數(shù)計算而來,主

    要是影響共振能群下的共振數(shù)據(jù)。表8給出不同參

    數(shù)λG選擇對結(jié)果的影響情況,對于不同的輸入?yún)?shù)λG值(13×λG),數(shù)據(jù)處理加工時間基本保持不變,對計算結(jié)果ΔKeff有一定影響,但靈敏度不是很大,ΔKeff變化量基本都在5×10?5以內(nèi),同時當(dāng)輸入?yún)?shù)僅輸入一個0.2時,與相同輸入13×0.2所得ΔKeff有少量差異。

    表8 加工時間及ΔKeff隨參數(shù)λG選擇的影響Table 8 Results of processing time and ΔKeffusing different λG.

    3 中子截面庫制作及驗證

    根據(jù)NJOY程序控制模塊輸入?yún)?shù)選擇分析,可以初步確定輸入?yún)?shù)的最佳選擇。在進(jìn)行WIMS格式數(shù)據(jù)的制作過程中,選取RECONR共振重造誤差0.001,最大重構(gòu)誤差0.005,多普勒(Doppler)展

    寬BROADR誤差限0.001,選取PURR模塊處理共振自屏數(shù)據(jù),共振級數(shù)選取100,表長取20,中子熱化處理THERMR誤差限0.001,散射角個數(shù)選取12,熱化的能量上限選取4.0 eV,GROUPR中子能群結(jié)構(gòu)選擇EPRI-CPM-69群,權(quán)重函數(shù)選擇EPRI-CELL LWR,此處γ能群結(jié)構(gòu)選擇NONE,勒讓德級數(shù)選取1,最后通過WIMSR模塊生成核素WIMS格式數(shù)據(jù),Goldstein-Cohen參數(shù)λG=0.2,輸入文件為GROUPR模塊的輸出文件GENDF TAPE,輸出WIMS-D格式數(shù)據(jù)WIMS TAPE文件。為簡便起見,圖4給出基于WIMS格式數(shù)據(jù)制作的235U多群中子截面數(shù)據(jù),并與IAEA核數(shù)據(jù)公布的WIMS-69群截面參考值進(jìn)行比對,兩者吻合較好。

    圖4 基于WIMS格式69群235U制作截面對比圖Fig.4 Cross-section of235U based on WIMS-69.

    將制作的中子數(shù)據(jù)庫導(dǎo)入組件均勻化計算程序(CASMO)中,選取基本燃料柵元與AP1000燃料組件作為基準(zhǔn)題驗證,AP1000堆芯組件采用改進(jìn)型17×17 XL Robust燃料(UO2)組件,其基本結(jié)構(gòu)如圖5所示。很好,絕對誤差控制在3‰以內(nèi),驗證了模塊輸入?yún)?shù)選擇的合理性以及基于WIMS格式制作的中子數(shù)據(jù)庫的正確性。

    圖5 AP1000燃料組件示意圖Fig.5 AP1000 fuel assembly.

    計算燃料柵元基準(zhǔn)題在HZP與HFP兩類工況,AP1000燃料組件600 K HFP工況下,不同富集度對應(yīng)的柵元、組件有效增殖系數(shù)Keff,并與基準(zhǔn)參考值進(jìn)行比較,結(jié)果如表9所示。根據(jù)結(jié)果可以看出,制作的中子數(shù)據(jù)庫計算結(jié)果與基準(zhǔn)參考值吻合

    表9 基準(zhǔn)題有效增殖系數(shù)KeffTable 9 Keffresults of benchmark problem.

    4 結(jié)語

    基于最新釋放的ENDF/B-VII.1核評價庫,采用核數(shù)據(jù)加工處理程序NJOY-99,制作基于WIMS格式的多群中子數(shù)據(jù)庫。針對LWR基本燃料柵元與AP1000燃料組件均勻化計算基準(zhǔn)題,以235U、238U核數(shù)據(jù)制作過程為優(yōu)化分析對象,通過分析制作過程不同處理模塊輸入?yún)?shù)選擇對截面庫加工時間及ΔKeff變化量的影響可以看出:(1) 共振重造模塊RECONR誤差限的改變、PURR共振級數(shù)逐級增加對加工時間有明顯的影響;(2) RECONR誤差限選取0.005以下,BROADR誤差限在0.01附近以下具有較好的計算精度,共振級數(shù)PURR選取100可以達(dá)到精度要求,同時保證較快的計算速度;(3) σ0值、熱化能量上限、權(quán)重譜、Goldstein-Cohen參數(shù)λG的選取對組件均勻化計算積分量Keff有較大的靈敏度,因此在多群數(shù)據(jù)庫的制作過程,處理模塊輸入?yún)?shù)的合理選擇非常重要。

    1 MacFarlane R E, Muir D W. The NJOY nuclear data processing system[R]. Version 91. Los Alamos Laboratory report LA-12740-M, October 1994

    2 Chadwick M B, Herman M, Oblozinsky P, et al. ENDF/B-Vll.1 nuelear data for science and technology: cross sections, covariances, fission product yields and decay data[J]. Nuelear Data Sheets, 2011, 112: 2887?2996

    3 鄒旸. 基于不同ENDF格式數(shù)據(jù)庫的MCNP溫度相關(guān)中子截面庫研究[J]. 核動力工程, 2012, 33(3): 12?16 ZOU Yang. Study on MCNP related to temperaturedependent neutron cross section library based on different ENDF format catabases[J]. Nuclear Power Engineering, 2012, 33(3): 12?16

    4 李松陽, 王侃, 余綱林, 等. MCNP溫度相關(guān)中子截面庫制作方法[J]. 核動力工程, 2008, 29(4): 10?14 LI Songyang, WANG Kan, YU Ganglin, et al. Processing methods for temperature-dependent MCNP libraries[J]. Nuclear Power Engineering, 2008, 29(4): 10?14

    5 Rahnema F, Gheorghiu H N M. Benchmark calculations for the Doppler coefficient of reactivity[J]. Ann Nuclear Energy, 1996, 23(12): 1011?1019

    6 劉萍. 核數(shù)據(jù)不確定性對ADS系統(tǒng)keff等積分量的影響和WIMS82群庫的研制[D]. 北京: 中國原子能科學(xué)研究院, 2004 LIU Ping. The effects on subcritical system of ADS caused by the uncertainty from nuclear data and development of WIMS82-group library[D]. Beijing: China Institute of Atomic Energy, 2004

    7 劉萍, 王耀清. TP2008核數(shù)據(jù)庫研制[J]. 核科學(xué)與工程, 2010, 30(3): 195?198 LIU Ping, WANG Yaoqing. The development of multigroup constants library TP2008[J]. Chinese Journal of Nuclear Science and Engineering, 2010, 30(3): 195?198

    8 劉萍. WIMS多群截面庫更新時NJOY輸入?yún)?shù)研究[J].原子核物理評論, 2001, 18(3): 188?191 LIU Ping. Study of NJOY input parameters for update of WIMS library[J]. Nuclear Physics Review, 2001, 18(3): 188?191

    9 Forssen B H. CASMO-3 user's manual[Z]. Studsvikofamerica, INC, 1994

    Optimization analysis of the production process based on WIMS format multigroup nuclear data

    XIE Mingliang CHEN Yuqing YU Lei SHI Hao

    (Department of Nuclear Energy Science and Engineering, Naval University of Engineering, Wuhan 430033, China)

    Background: NJOY is used widely nuclear data processing program that can convert the data format of ENDF/B nuclear database into WIMS format of the multi-group cross section database. Purpose: This study aims to produce a multi-group database of WIMS-D format on the basis of the latest release of nuclear evaluation library ENDF/B-VII.1 by using nuclear data processing program NJOY-99. Methods: In allusion to basic fuel cell homogenization calculation benchmark problem of the light water reactor (LWR), taking the235U and238U nuclide as the main objects of analysis, the effect on processing time of making cross section library, as well as integral parameters ΔKeffand data sensitivity of selecting input parameters of the NJOY program were compared and analyzed, which optimized the option of input parameter. Results and Conclusion: The results of validating on benchmark problems showed that the production of the multi-group database was correct, which embodied the high accuracy of calculation and provided the basis data of fuel assembly homogenization calculation for pressurized water reactor.

    ENDF/B-VII.1, WIMS, NJOY-99, Nuclear data library

    TL99

    10.11889/j.0253-3219.2015.hjs.38.110502

    謝明亮,男,1990年出生,2013年畢業(yè)于南華大學(xué),現(xiàn)為碩士研究生,從事核反應(yīng)堆安全分析方向工作

    陳玉清,E-mail: chenyuqing301@163.com

    2015-08-28,

    2015-10-11

    群數(shù)據(jù)庫,針對輕水堆(Light Water Reactor, LWR)基本燃料柵元均勻化計算基準(zhǔn)題,以235U、238U核素為主要分析對象,對比研究了NJOY程序輸入模塊參數(shù)的選擇對截面庫制作加工時間、積分量ΔKeff及靈敏度的影響,得到優(yōu)化的輸入?yún)?shù)選擇方案?;鶞?zhǔn)例題驗證結(jié)果表明:所制作的多群數(shù)據(jù)庫是正確的,Keff計算精度較高,可為壓水堆燃料組件均勻化計算提供數(shù)據(jù)基礎(chǔ)。

    CLC TL99

    猜你喜歡
    共振誤差加工
    認(rèn)識“超加工食品”
    中老年保健(2022年5期)2022-08-24 02:36:10
    后期加工
    角接觸球軸承接觸角誤差控制
    哈爾濱軸承(2020年2期)2020-11-06 09:22:26
    Beidou, le système de navigation par satellite compatible et interopérable
    壓力容器制造誤差探究
    安然 與時代同頻共振
    選硬人打硬仗——紫陽縣黨建與脫貧同頻共振
    CTA 中紡院+ 化纖聯(lián)盟 強(qiáng)強(qiáng)聯(lián)合 科技共振
    菱的簡易加工幾法
    九十億分之一的“生死”誤差
    山東青年(2016年2期)2016-02-28 14:25:41
    久久精品人妻少妇| 在线十欧美十亚洲十日本专区| av黄色大香蕉| 757午夜福利合集在线观看| 国内精品美女久久久久久| 亚洲美女黄片视频| 99国产综合亚洲精品| 亚洲成av人片免费观看| 国产色爽女视频免费观看| 国产精品乱码一区二三区的特点| 亚洲国产精品sss在线观看| 国产日本99.免费观看| 国产真实乱freesex| 国产真实乱freesex| 国内精品久久久久久久电影| aaaaa片日本免费| 国产单亲对白刺激| 国产中年淑女户外野战色| 很黄的视频免费| 99热6这里只有精品| 97超级碰碰碰精品色视频在线观看| 国产中年淑女户外野战色| 国产精品亚洲av一区麻豆| 精品日产1卡2卡| 欧美一级a爱片免费观看看| 午夜视频国产福利| 欧美绝顶高潮抽搐喷水| 国产毛片a区久久久久| 欧美区成人在线视频| 日本五十路高清| 国产欧美日韩精品亚洲av| 欧美成人一区二区免费高清观看| 蜜桃久久精品国产亚洲av| 伊人久久精品亚洲午夜| 五月伊人婷婷丁香| 亚洲aⅴ乱码一区二区在线播放| 操出白浆在线播放| 午夜视频国产福利| 国产三级在线视频| 性色av乱码一区二区三区2| 99热6这里只有精品| 国产av麻豆久久久久久久| 久久精品91蜜桃| 亚洲成人中文字幕在线播放| 99热只有精品国产| 国产高清激情床上av| 国产亚洲精品av在线| 免费看光身美女| 国产午夜福利久久久久久| 日本与韩国留学比较| 97人妻精品一区二区三区麻豆| 日韩欧美三级三区| 狠狠狠狠99中文字幕| 日韩亚洲欧美综合| 亚洲一区高清亚洲精品| 天堂影院成人在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利在线观看免费完整高清在 | 久久久久久大精品| 成人无遮挡网站| 亚洲人与动物交配视频| 午夜久久久久精精品| 免费av观看视频| 午夜免费男女啪啪视频观看 | 国产精品,欧美在线| 亚洲熟妇中文字幕五十中出| 免费高清视频大片| 国产伦人伦偷精品视频| 99在线视频只有这里精品首页| 色综合欧美亚洲国产小说| 成人18禁在线播放| 国产 一区 欧美 日韩| 亚洲人成网站高清观看| 国产av麻豆久久久久久久| 国产探花极品一区二区| 三级男女做爰猛烈吃奶摸视频| 欧美日韩黄片免| 最好的美女福利视频网| www.色视频.com| 丰满乱子伦码专区| 欧美xxxx黑人xx丫x性爽| 久99久视频精品免费| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看日本一区| 成年版毛片免费区| 国产午夜精品久久久久久一区二区三区 | 99久久无色码亚洲精品果冻| av黄色大香蕉| 日韩欧美 国产精品| 少妇的逼好多水| 久久久国产成人免费| 亚洲在线自拍视频| 12—13女人毛片做爰片一| 国产一区二区在线av高清观看| 精品一区二区三区视频在线 | 亚洲av第一区精品v没综合| 国产精品亚洲av一区麻豆| 亚洲最大成人中文| 日本 av在线| 亚洲aⅴ乱码一区二区在线播放| 熟女少妇亚洲综合色aaa.| 网址你懂的国产日韩在线| 男女做爰动态图高潮gif福利片| 中国美女看黄片| 尤物成人国产欧美一区二区三区| 欧美日韩综合久久久久久 | 国产午夜福利久久久久久| 国产精品一区二区三区四区久久| 日本撒尿小便嘘嘘汇集6| 欧美日本视频| 久久精品夜夜夜夜夜久久蜜豆| 一个人观看的视频www高清免费观看| 精品99又大又爽又粗少妇毛片 | 国产精品99久久久久久久久| 日韩 欧美 亚洲 中文字幕| 久久6这里有精品| 亚洲精品乱码久久久v下载方式 | 女人被狂操c到高潮| 久久精品91蜜桃| 91麻豆精品激情在线观看国产| 男人的好看免费观看在线视频| 九色国产91popny在线| 少妇高潮的动态图| 亚洲精品色激情综合| 成年人黄色毛片网站| 久久久久久久亚洲中文字幕 | 宅男免费午夜| 一本一本综合久久| 久久精品国产综合久久久| 99国产极品粉嫩在线观看| a级毛片a级免费在线| 午夜激情福利司机影院| 老司机深夜福利视频在线观看| 亚洲欧美日韩无卡精品| 午夜老司机福利剧场| 真实男女啪啪啪动态图| 人人妻人人澡欧美一区二区| 亚洲欧美日韩卡通动漫| 黄色丝袜av网址大全| 全区人妻精品视频| 淫秽高清视频在线观看| 熟女少妇亚洲综合色aaa.| 女同久久另类99精品国产91| 成年女人永久免费观看视频| 婷婷丁香在线五月| 国产精品98久久久久久宅男小说| 99热只有精品国产| 欧美3d第一页| 亚洲 欧美 日韩 在线 免费| 舔av片在线| 免费看日本二区| 手机成人av网站| 亚洲激情在线av| 日韩亚洲欧美综合| 欧美精品啪啪一区二区三区| 少妇的逼好多水| 久久久国产成人精品二区| 亚洲专区国产一区二区| 午夜福利免费观看在线| 国产精品久久电影中文字幕| 大型黄色视频在线免费观看| 亚洲精品影视一区二区三区av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | av中文乱码字幕在线| 天天添夜夜摸| 精品人妻一区二区三区麻豆 | 高清毛片免费观看视频网站| 天堂av国产一区二区熟女人妻| 亚洲精品一区av在线观看| 18禁美女被吸乳视频| 亚洲男人的天堂狠狠| 国产精品久久久久久久电影 | 亚洲午夜理论影院| www.色视频.com| 女警被强在线播放| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美一区二区三区在线观看| 免费在线观看成人毛片| 日本黄色视频三级网站网址| 一个人免费在线观看电影| 在线播放国产精品三级| 最好的美女福利视频网| 久久久久九九精品影院| 又黄又粗又硬又大视频| 免费av毛片视频| 久99久视频精品免费| 又黄又爽又免费观看的视频| 国产亚洲欧美在线一区二区| 成人18禁在线播放| 亚洲欧美日韩无卡精品| 国产成人aa在线观看| 日韩国内少妇激情av| 欧美黑人欧美精品刺激| 欧美日韩亚洲国产一区二区在线观看| 亚洲av成人不卡在线观看播放网| 中国美女看黄片| 色综合婷婷激情| 九九久久精品国产亚洲av麻豆| 国产欧美日韩一区二区精品| 老熟妇乱子伦视频在线观看| 国产精品爽爽va在线观看网站| 亚洲avbb在线观看| ponron亚洲| 最近最新中文字幕大全免费视频| 天堂av国产一区二区熟女人妻| 一区福利在线观看| 老司机午夜十八禁免费视频| 亚洲专区国产一区二区| 18禁美女被吸乳视频| 亚洲在线自拍视频| 日本 av在线| 少妇人妻精品综合一区二区 | 婷婷精品国产亚洲av在线| 母亲3免费完整高清在线观看| 脱女人内裤的视频| 午夜a级毛片| 国产真实乱freesex| 国产精品国产高清国产av| 女生性感内裤真人,穿戴方法视频| 夜夜夜夜夜久久久久| 亚洲国产高清在线一区二区三| 国产精品女同一区二区软件 | 高清在线国产一区| 亚洲成人中文字幕在线播放| 免费在线观看亚洲国产| 国产99白浆流出| 欧美日韩乱码在线| 色综合婷婷激情| 国产精品免费一区二区三区在线| 欧美另类亚洲清纯唯美| 亚洲真实伦在线观看| 十八禁人妻一区二区| 小蜜桃在线观看免费完整版高清| 亚洲国产欧洲综合997久久,| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品在线观看二区| 禁无遮挡网站| 麻豆成人av在线观看| 十八禁网站免费在线| 欧洲精品卡2卡3卡4卡5卡区| aaaaa片日本免费| 亚洲午夜理论影院| 亚洲内射少妇av| 中出人妻视频一区二区| 日韩中文字幕欧美一区二区| 91在线精品国自产拍蜜月 | 亚洲成人中文字幕在线播放| 好男人在线观看高清免费视频| 国产精品野战在线观看| 他把我摸到了高潮在线观看| 国产美女午夜福利| 国产一级毛片七仙女欲春2| 天天添夜夜摸| 国产精品一区二区免费欧美| 丰满的人妻完整版| 精品一区二区三区视频在线 | 久久精品91无色码中文字幕| 日韩欧美国产在线观看| 村上凉子中文字幕在线| 色精品久久人妻99蜜桃| 国产欧美日韩一区二区三| 禁无遮挡网站| 一级毛片高清免费大全| 九九在线视频观看精品| 美女高潮喷水抽搐中文字幕| 最近最新中文字幕大全免费视频| 波多野结衣高清作品| 国产精华一区二区三区| av片东京热男人的天堂| 日韩欧美在线乱码| 99视频精品全部免费 在线| 嫩草影视91久久| 两个人的视频大全免费| 亚洲国产精品sss在线观看| 久久精品91无色码中文字幕| 狂野欧美白嫩少妇大欣赏| 国产高清videossex| 激情在线观看视频在线高清| 亚洲片人在线观看| 听说在线观看完整版免费高清| 国产v大片淫在线免费观看| 亚洲无线观看免费| 国产蜜桃级精品一区二区三区| a在线观看视频网站| 中文字幕久久专区| 精品国产美女av久久久久小说| 给我免费播放毛片高清在线观看| 精品福利观看| 国产一区在线观看成人免费| 亚洲一区高清亚洲精品| 日韩精品青青久久久久久| 最近最新中文字幕大全免费视频| 国产午夜精品论理片| 国产主播在线观看一区二区| 色av中文字幕| 麻豆久久精品国产亚洲av| 日韩欧美精品免费久久 | 一级毛片女人18水好多| 人人妻,人人澡人人爽秒播| 欧美成狂野欧美在线观看| 可以在线观看毛片的网站| 精品国产美女av久久久久小说| 男人的好看免费观看在线视频| 看黄色毛片网站| 久久精品影院6| 国产精品亚洲av一区麻豆| 亚洲av成人不卡在线观看播放网| 中国美女看黄片| 国产精品国产高清国产av| 久久99热这里只有精品18| 欧美成人免费av一区二区三区| 国产免费男女视频| 免费高清视频大片| 一本精品99久久精品77| 村上凉子中文字幕在线| 女人被狂操c到高潮| 亚洲最大成人中文| 亚洲欧美日韩卡通动漫| 内射极品少妇av片p| 久久婷婷人人爽人人干人人爱| 国产单亲对白刺激| 免费在线观看影片大全网站| 99riav亚洲国产免费| av专区在线播放| 十八禁网站免费在线| 午夜久久久久精精品| 人人妻人人澡欧美一区二区| 亚洲人成伊人成综合网2020| 在线观看66精品国产| 亚洲黑人精品在线| 免费观看精品视频网站| 五月玫瑰六月丁香| 美女 人体艺术 gogo| 色老头精品视频在线观看| 国产精品女同一区二区软件 | 亚洲成人精品中文字幕电影| 亚洲精品一区av在线观看| 久久久久久久久大av| 99国产综合亚洲精品| 亚洲av不卡在线观看| 好看av亚洲va欧美ⅴa在| 一区二区三区高清视频在线| 色吧在线观看| 久久天躁狠狠躁夜夜2o2o| 99热这里只有精品一区| 国产毛片a区久久久久| 高清日韩中文字幕在线| 久久亚洲真实| 国产亚洲av嫩草精品影院| 午夜免费男女啪啪视频观看 | 亚洲电影在线观看av| 午夜免费观看网址| 成人无遮挡网站| 欧美色欧美亚洲另类二区| 老司机午夜十八禁免费视频| 99精品久久久久人妻精品| 亚洲av中文字字幕乱码综合| 欧美日韩瑟瑟在线播放| ponron亚洲| 深爱激情五月婷婷| 国产精品久久久久久久久免 | 99热只有精品国产| 日本 欧美在线| 美女高潮的动态| 日韩欧美在线乱码| 在线播放无遮挡| 日韩 欧美 亚洲 中文字幕| 男插女下体视频免费在线播放| 国产成人av激情在线播放| 最好的美女福利视频网| 琪琪午夜伦伦电影理论片6080| 女人被狂操c到高潮| 窝窝影院91人妻| 伊人久久大香线蕉亚洲五| 免费搜索国产男女视频| 亚洲avbb在线观看| 亚洲成人久久性| 久久精品综合一区二区三区| 久久草成人影院| 亚洲专区国产一区二区| 露出奶头的视频| 日本撒尿小便嘘嘘汇集6| 国产伦一二天堂av在线观看| 深夜精品福利| 女警被强在线播放| h日本视频在线播放| 丁香欧美五月| 成熟少妇高潮喷水视频| 此物有八面人人有两片| 国产亚洲欧美98| 村上凉子中文字幕在线| 国产成人aa在线观看| 男人舔奶头视频| 在线观看免费午夜福利视频| 黄色日韩在线| 午夜免费激情av| 国产极品精品免费视频能看的| 51午夜福利影视在线观看| 亚洲国产精品合色在线| 18禁在线播放成人免费| 天天一区二区日本电影三级| 亚洲内射少妇av| 国产av一区在线观看免费| 丰满的人妻完整版| 免费看日本二区| 最新中文字幕久久久久| 亚洲国产色片| 神马国产精品三级电影在线观看| 老司机午夜十八禁免费视频| 久久久久久人人人人人| 成年免费大片在线观看| 在线天堂最新版资源| 精品欧美国产一区二区三| ponron亚洲| 国产精品久久视频播放| 久久精品亚洲精品国产色婷小说| 亚洲精品国产精品久久久不卡| 九色国产91popny在线| 色噜噜av男人的天堂激情| 国产精品永久免费网站| 一进一出好大好爽视频| 亚洲人成网站在线播| 婷婷六月久久综合丁香| 99久国产av精品| 男女下面进入的视频免费午夜| 日韩欧美在线二视频| 无遮挡黄片免费观看| 一卡2卡三卡四卡精品乱码亚洲| 婷婷精品国产亚洲av在线| 免费看日本二区| 丁香欧美五月| 久久精品国产99精品国产亚洲性色| 一级a爱片免费观看的视频| 日韩有码中文字幕| 少妇人妻一区二区三区视频| 欧美色视频一区免费| 国语自产精品视频在线第100页| 最后的刺客免费高清国语| 日韩国内少妇激情av| 亚洲精华国产精华精| 国产午夜精品论理片| 女生性感内裤真人,穿戴方法视频| 蜜桃亚洲精品一区二区三区| 中文字幕av成人在线电影| 一区二区三区国产精品乱码| 国产精品久久久人人做人人爽| 亚洲在线观看片| 亚洲国产欧美人成| 日本三级黄在线观看| 亚洲,欧美精品.| 成人亚洲精品av一区二区| 免费电影在线观看免费观看| 看免费av毛片| 中文字幕人妻丝袜一区二区| 韩国av一区二区三区四区| 国产成人系列免费观看| 99久久成人亚洲精品观看| 国产精品一区二区三区四区久久| 久久国产精品人妻蜜桃| 国产精品爽爽va在线观看网站| 一个人看视频在线观看www免费 | 久久久久免费精品人妻一区二区| 国产黄a三级三级三级人| 国产亚洲精品一区二区www| 一个人观看的视频www高清免费观看| 韩国av一区二区三区四区| 日韩av在线大香蕉| 熟女电影av网| 欧美黑人巨大hd| 国产午夜福利久久久久久| 99久久无色码亚洲精品果冻| 在线天堂最新版资源| 18禁黄网站禁片午夜丰满| 1000部很黄的大片| 少妇的丰满在线观看| 国产免费av片在线观看野外av| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久久黄片| 欧美中文综合在线视频| 国产黄a三级三级三级人| 又黄又粗又硬又大视频| 欧美黄色片欧美黄色片| 在线国产一区二区在线| 岛国在线免费视频观看| 亚洲国产日韩欧美精品在线观看 | 欧美日韩综合久久久久久 | 在线观看av片永久免费下载| 欧美中文日本在线观看视频| 国产精品影院久久| www.熟女人妻精品国产| 成年女人永久免费观看视频| a级一级毛片免费在线观看| 国产高清有码在线观看视频| 精品不卡国产一区二区三区| 午夜免费激情av| 亚洲精品在线观看二区| 国产亚洲欧美98| 精品久久久久久,| 国产aⅴ精品一区二区三区波| 日本与韩国留学比较| 精品久久久久久久久久久久久| 国产亚洲欧美在线一区二区| 脱女人内裤的视频| 亚洲av熟女| 三级国产精品欧美在线观看| 99在线视频只有这里精品首页| 久久国产乱子伦精品免费另类| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久香蕉精品热| 欧美日韩福利视频一区二区| 亚洲国产精品999在线| av片东京热男人的天堂| 国产精品,欧美在线| 午夜免费男女啪啪视频观看 | 女人被狂操c到高潮| 国产精品久久久久久久电影 | 99热精品在线国产| 日韩成人在线观看一区二区三区| 欧美3d第一页| 好男人电影高清在线观看| 99在线视频只有这里精品首页| 色老头精品视频在线观看| 女人十人毛片免费观看3o分钟| netflix在线观看网站| 综合色av麻豆| 日韩欧美国产在线观看| 色综合亚洲欧美另类图片| 国产欧美日韩一区二区精品| 我的老师免费观看完整版| 18禁美女被吸乳视频| 欧美激情久久久久久爽电影| 亚洲欧美一区二区三区黑人| 国产毛片a区久久久久| 99精品久久久久人妻精品| 成人av一区二区三区在线看| 婷婷丁香在线五月| 长腿黑丝高跟| 看片在线看免费视频| 69av精品久久久久久| 国产爱豆传媒在线观看| 精品无人区乱码1区二区| 日韩欧美一区二区三区在线观看| 国产精品野战在线观看| 一个人免费在线观看的高清视频| 国产老妇女一区| 亚洲欧美激情综合另类| 久久人人精品亚洲av| 国产综合懂色| 俄罗斯特黄特色一大片| 国产毛片a区久久久久| 熟女人妻精品中文字幕| avwww免费| 日本免费a在线| 国产精品亚洲美女久久久| 一夜夜www| 成人国产综合亚洲| 毛片女人毛片| 99国产精品一区二区三区| 久9热在线精品视频| 18禁黄网站禁片免费观看直播| 欧美激情在线99| 大型黄色视频在线免费观看| 国产高清videossex| 午夜福利在线在线| 中文字幕人妻熟人妻熟丝袜美 | 久久久久久久亚洲中文字幕 | 天天躁日日操中文字幕| 国产乱人视频| 亚洲国产精品久久男人天堂| 男女做爰动态图高潮gif福利片| netflix在线观看网站| 禁无遮挡网站| 老汉色∧v一级毛片| 精品久久久久久成人av| 国产一区二区三区在线臀色熟女| 国产免费男女视频| 欧美日韩一级在线毛片| 草草在线视频免费看| www.熟女人妻精品国产| 成人特级av手机在线观看| 99久久九九国产精品国产免费| 亚洲国产色片| 国产野战对白在线观看| 国产综合懂色| 成人午夜高清在线视频| 久久精品影院6| 欧洲精品卡2卡3卡4卡5卡区| 国产精品精品国产色婷婷| 亚洲无线观看免费| 国产欧美日韩一区二区精品| 99国产极品粉嫩在线观看| 麻豆成人午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡| 日本a在线网址| 午夜视频国产福利| 国产高清有码在线观看视频| 97碰自拍视频| 国产成人a区在线观看| 久久久久久久午夜电影| 国产老妇女一区| 欧美一级毛片孕妇| av欧美777| 好看av亚洲va欧美ⅴa在| 国产淫片久久久久久久久 | 搡老熟女国产l中国老女人| 男人和女人高潮做爰伦理| 精品国产三级普通话版| 丰满人妻熟妇乱又伦精品不卡| 我要搜黄色片| 日本免费一区二区三区高清不卡| 老司机福利观看| 18禁国产床啪视频网站| 国产精品国产高清国产av| 日本撒尿小便嘘嘘汇集6| 一a级毛片在线观看|