• <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
    欧美成人一区二区免费高清观看 | 欧美中文综合在线视频| 国产主播在线观看一区二区| 12—13女人毛片做爰片一| 日韩精品中文字幕看吧| 一区二区三区国产精品乱码| 老司机福利观看| 久99久视频精品免费| 久久久久九九精品影院| 黄色视频,在线免费观看| 国产一区二区三区综合在线观看| 无遮挡黄片免费观看| 每晚都被弄得嗷嗷叫到高潮| 变态另类成人亚洲欧美熟女 | 国产精品,欧美在线| 午夜影院日韩av| 欧美+亚洲+日韩+国产| 一区二区三区高清视频在线| 精品一区二区三区视频在线观看免费| 好看av亚洲va欧美ⅴa在| 国产在线精品亚洲第一网站| 欧美大码av| 欧美日韩亚洲国产一区二区在线观看| 久久久国产成人精品二区| 性少妇av在线| 三级毛片av免费| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲五月婷婷丁香| 久久精品人人爽人人爽视色| 婷婷六月久久综合丁香| 丁香欧美五月| 亚洲成人精品中文字幕电影| 亚洲全国av大片| 国产精品综合久久久久久久免费 | 亚洲成人国产一区在线观看| 亚洲精品久久成人aⅴ小说| 国产不卡一卡二| 久久久久久免费高清国产稀缺| 久久香蕉精品热| 久热爱精品视频在线9| 少妇熟女aⅴ在线视频| 麻豆一二三区av精品| 看片在线看免费视频| 日韩大尺度精品在线看网址 | 午夜a级毛片| 免费少妇av软件| 精品人妻1区二区| 亚洲成人久久性| 精品电影一区二区在线| 纯流量卡能插随身wifi吗| 久久婷婷成人综合色麻豆| 麻豆成人av在线观看| 成人国产综合亚洲| 欧美日韩瑟瑟在线播放| 亚洲自偷自拍图片 自拍| 韩国av一区二区三区四区| 热re99久久国产66热| 热re99久久国产66热| 国产亚洲精品综合一区在线观看 | 色综合亚洲欧美另类图片| 久久精品国产清高在天天线| 丰满的人妻完整版| 99热只有精品国产| 国产欧美日韩精品亚洲av| 国产午夜精品久久久久久| 亚洲av日韩精品久久久久久密| 久久久国产欧美日韩av| 久久精品91蜜桃| 亚洲精品在线美女| 亚洲精品中文字幕在线视频| 亚洲中文字幕一区二区三区有码在线看 | 99re在线观看精品视频| 乱人伦中国视频| 欧美激情久久久久久爽电影 | 欧美一区二区精品小视频在线| 久热爱精品视频在线9| 色播在线永久视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲色图av天堂| 夜夜看夜夜爽夜夜摸| 亚洲成人免费电影在线观看| 国产亚洲精品久久久久久毛片| 国产一区二区三区视频了| 国产一级毛片七仙女欲春2 | 久久久久久免费高清国产稀缺| 每晚都被弄得嗷嗷叫到高潮| 欧美+亚洲+日韩+国产| 亚洲国产毛片av蜜桃av| avwww免费| 国产免费av片在线观看野外av| 亚洲,欧美精品.| 国产一区二区在线av高清观看| 国产一区二区在线av高清观看| 91精品三级在线观看| 99在线人妻在线中文字幕| 亚洲精品美女久久av网站| 女性生殖器流出的白浆| 亚洲精品美女久久久久99蜜臀| 国产私拍福利视频在线观看| 十八禁网站免费在线| 久久久久久大精品| 欧美日韩中文字幕国产精品一区二区三区 | 99re在线观看精品视频| 国产蜜桃级精品一区二区三区| www.999成人在线观看| 欧美日本中文国产一区发布| 国产午夜精品久久久久久| 精品免费久久久久久久清纯| 高清毛片免费观看视频网站| 香蕉国产在线看| 日韩一卡2卡3卡4卡2021年| 首页视频小说图片口味搜索| 色av中文字幕| 人妻久久中文字幕网| 我的亚洲天堂| 久久国产亚洲av麻豆专区| 看黄色毛片网站| 熟妇人妻久久中文字幕3abv| 岛国在线观看网站| 精品福利观看| 777久久人妻少妇嫩草av网站| 一级黄色大片毛片| 丝袜在线中文字幕| 丝袜在线中文字幕| 美女大奶头视频| 国产欧美日韩一区二区三| 一二三四社区在线视频社区8| 日本欧美视频一区| 国产精品1区2区在线观看.| 午夜免费鲁丝| 久久久国产成人精品二区| 男人的好看免费观看在线视频 | 久久久久久大精品| 日本五十路高清| av电影中文网址| 熟女少妇亚洲综合色aaa.| 国产精品久久久久久人妻精品电影| 一a级毛片在线观看| 777久久人妻少妇嫩草av网站| 国产成人av激情在线播放| x7x7x7水蜜桃| 69av精品久久久久久| 天堂√8在线中文| 久久人人精品亚洲av| 国产精品一区二区精品视频观看| 两个人看的免费小视频| 50天的宝宝边吃奶边哭怎么回事| 男男h啪啪无遮挡| 欧美乱色亚洲激情| 亚洲人成网站在线播放欧美日韩| 女人被躁到高潮嗷嗷叫费观| 欧美日本亚洲视频在线播放| 国产不卡一卡二| 人人妻人人爽人人添夜夜欢视频| 女生性感内裤真人,穿戴方法视频| av在线天堂中文字幕| 日韩有码中文字幕| 狂野欧美激情性xxxx| 热re99久久国产66热| 老司机靠b影院| 一区二区日韩欧美中文字幕| 国产精品精品国产色婷婷| 18禁黄网站禁片午夜丰满| 成人亚洲精品av一区二区| 午夜免费鲁丝| 中文字幕精品免费在线观看视频| 国产视频一区二区在线看| 日韩欧美国产一区二区入口| 久久人人精品亚洲av| 国产av一区二区精品久久| 国产精品一区二区三区四区久久 | 麻豆久久精品国产亚洲av| 亚洲成人精品中文字幕电影| 国产精品,欧美在线| 黄色视频,在线免费观看| 在线观看午夜福利视频| 操美女的视频在线观看| 国产高清videossex| 国产1区2区3区精品| 国产精品98久久久久久宅男小说| 妹子高潮喷水视频| 国产成人系列免费观看| 99精品欧美一区二区三区四区| 精品免费久久久久久久清纯| 午夜精品久久久久久毛片777| 久久久国产成人免费| 亚洲性夜色夜夜综合| 一区二区三区精品91| 亚洲人成伊人成综合网2020| 成人18禁在线播放| 91精品国产国语对白视频| 老司机在亚洲福利影院| 国产精品久久久久久亚洲av鲁大| 日本精品一区二区三区蜜桃| 亚洲一码二码三码区别大吗| 可以在线观看的亚洲视频| 99riav亚洲国产免费| 免费女性裸体啪啪无遮挡网站| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲欧洲精品一区二区精品久久久| 久久这里只有精品19| 母亲3免费完整高清在线观看| 91成年电影在线观看| 宅男免费午夜| 两个人免费观看高清视频| 在线观看一区二区三区| 美女大奶头视频| 韩国精品一区二区三区| 亚洲全国av大片| 麻豆久久精品国产亚洲av| 亚洲中文字幕一区二区三区有码在线看 | 在线观看日韩欧美| 在线播放国产精品三级| 岛国在线观看网站| 男女午夜视频在线观看| 欧美日本中文国产一区发布| 日韩大码丰满熟妇| 亚洲男人的天堂狠狠| 中文字幕精品免费在线观看视频| 国产成人欧美| 欧美激情 高清一区二区三区| 一个人观看的视频www高清免费观看 | 长腿黑丝高跟| 后天国语完整版免费观看| 级片在线观看| 欧美精品亚洲一区二区| 一a级毛片在线观看| 亚洲七黄色美女视频| 免费在线观看亚洲国产| 国产1区2区3区精品| 12—13女人毛片做爰片一| 久久精品aⅴ一区二区三区四区| 亚洲狠狠婷婷综合久久图片| 日韩一卡2卡3卡4卡2021年| 黑人操中国人逼视频| 国产精品影院久久| 亚洲av电影在线进入| 老司机午夜十八禁免费视频| 亚洲精品国产色婷婷电影| 岛国在线观看网站| 最好的美女福利视频网| 国产精品一区二区三区四区久久 | 免费在线观看完整版高清| 免费人成视频x8x8入口观看| 亚洲男人的天堂狠狠| 中文字幕精品免费在线观看视频| 人人妻,人人澡人人爽秒播| 久久久久久大精品| 精品无人区乱码1区二区| 午夜免费鲁丝| 妹子高潮喷水视频| 免费在线观看完整版高清| 日日干狠狠操夜夜爽| 日本vs欧美在线观看视频| 看片在线看免费视频| 三级毛片av免费| 一卡2卡三卡四卡精品乱码亚洲| 国产成+人综合+亚洲专区| 午夜免费成人在线视频| 91成年电影在线观看| 中文字幕人妻熟女乱码| 夜夜躁狠狠躁天天躁| 亚洲av五月六月丁香网| 97人妻精品一区二区三区麻豆 | 欧美+亚洲+日韩+国产| 叶爱在线成人免费视频播放| 两性夫妻黄色片| 国产高清videossex| 亚洲熟女毛片儿| 精品乱码久久久久久99久播| 99riav亚洲国产免费| 级片在线观看| 国产精品亚洲一级av第二区| 久久精品国产亚洲av香蕉五月| 亚洲一码二码三码区别大吗| 自拍欧美九色日韩亚洲蝌蚪91| 岛国在线观看网站| 在线观看一区二区三区| 亚洲黑人精品在线| 欧美日韩亚洲国产一区二区在线观看| 99国产精品一区二区三区| 午夜福利18| 久久精品国产亚洲av高清一级| 欧美日韩亚洲国产一区二区在线观看| 色婷婷久久久亚洲欧美| 最好的美女福利视频网| 俄罗斯特黄特色一大片| 日本一区二区免费在线视频| 久久久久久久久中文| 悠悠久久av| 十八禁网站免费在线| 青草久久国产| 久久人妻av系列| 亚洲自偷自拍图片 自拍| 久久人人精品亚洲av| 久久中文字幕人妻熟女| 免费在线观看黄色视频的| 国产激情欧美一区二区| 老司机靠b影院| 欧美黑人欧美精品刺激| 一本久久中文字幕| 黑人欧美特级aaaaaa片| 久久香蕉精品热| 国产精品久久久人人做人人爽| 在线观看午夜福利视频| 91大片在线观看| 国产欧美日韩一区二区三区在线| 亚洲精品粉嫩美女一区| 99精品欧美一区二区三区四区| 国产伦人伦偷精品视频| 久久久国产成人精品二区| 午夜日韩欧美国产| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品亚洲av一区麻豆| 天天躁夜夜躁狠狠躁躁| 日韩欧美国产在线观看| 法律面前人人平等表现在哪些方面| 美女扒开内裤让男人捅视频| 欧美日韩瑟瑟在线播放| 在线观看一区二区三区| 久久天躁狠狠躁夜夜2o2o| 欧洲精品卡2卡3卡4卡5卡区| 不卡av一区二区三区| avwww免费| 丝袜人妻中文字幕| 怎么达到女性高潮| 久99久视频精品免费| 黄片大片在线免费观看| 亚洲欧美精品综合久久99| 成人精品一区二区免费| 美女大奶头视频| 18美女黄网站色大片免费观看| 色av中文字幕| 一级a爱视频在线免费观看| 国产精品免费一区二区三区在线| 两性夫妻黄色片| 99精品在免费线老司机午夜| 波多野结衣高清无吗| 一夜夜www| 伦理电影免费视频| 日韩免费av在线播放| 制服丝袜大香蕉在线| 欧美日韩亚洲国产一区二区在线观看| 国语自产精品视频在线第100页| 欧美黑人欧美精品刺激| 黄色片一级片一级黄色片| 欧美 亚洲 国产 日韩一| 久久婷婷人人爽人人干人人爱 | 夜夜夜夜夜久久久久| 久久久久国产精品人妻aⅴ院| 美女高潮喷水抽搐中文字幕| 黄色片一级片一级黄色片| 99国产精品一区二区三区| 1024香蕉在线观看| 亚洲人成电影免费在线| 成人永久免费在线观看视频| av视频免费观看在线观看| 久久九九热精品免费| 看片在线看免费视频| 欧美在线黄色| 中亚洲国语对白在线视频| 丁香六月欧美| 国产xxxxx性猛交| 欧美色欧美亚洲另类二区 | 亚洲第一av免费看| 日韩大码丰满熟妇| 丝袜美足系列| 乱人伦中国视频| 人成视频在线观看免费观看| 窝窝影院91人妻| 免费少妇av软件| 好男人在线观看高清免费视频 | 变态另类成人亚洲欧美熟女 | 久久久久精品国产欧美久久久| 日韩中文字幕欧美一区二区| 国产午夜精品久久久久久| 欧美色欧美亚洲另类二区 | 亚洲国产欧美一区二区综合| 久久国产精品人妻蜜桃| 涩涩av久久男人的天堂| 国产精品久久视频播放| 日日夜夜操网爽| avwww免费| 日韩有码中文字幕| 国产精品乱码一区二三区的特点 | 日韩欧美在线二视频| 久久久久久免费高清国产稀缺| 俄罗斯特黄特色一大片| 成人18禁高潮啪啪吃奶动态图| 人人妻,人人澡人人爽秒播| 欧美成人一区二区免费高清观看 | 女同久久另类99精品国产91| 视频区欧美日本亚洲| 成人18禁在线播放| 黑人巨大精品欧美一区二区蜜桃| av超薄肉色丝袜交足视频| 脱女人内裤的视频| 国产又爽黄色视频| 后天国语完整版免费观看| 在线十欧美十亚洲十日本专区| 男女床上黄色一级片免费看| 亚洲成国产人片在线观看| 亚洲国产高清在线一区二区三 | 婷婷六月久久综合丁香| 国产麻豆成人av免费视频| 久久香蕉精品热| 两个人视频免费观看高清| 日日摸夜夜添夜夜添小说| 午夜a级毛片| 非洲黑人性xxxx精品又粗又长| 精品少妇一区二区三区视频日本电影| 91成年电影在线观看| 久久久久精品国产欧美久久久| 欧美+亚洲+日韩+国产| 中亚洲国语对白在线视频| 国产亚洲精品综合一区在线观看 | 亚洲精品av麻豆狂野| 男女下面进入的视频免费午夜 | 日本a在线网址| 免费在线观看日本一区| avwww免费| 美女 人体艺术 gogo| 老司机在亚洲福利影院| 国产精品爽爽va在线观看网站 | 母亲3免费完整高清在线观看| 色综合欧美亚洲国产小说| 亚洲专区国产一区二区| 午夜a级毛片| bbb黄色大片| 国产精品亚洲一级av第二区| 国产亚洲精品久久久久5区| 欧美最黄视频在线播放免费| avwww免费| 免费在线观看亚洲国产| 欧美亚洲日本最大视频资源| 亚洲aⅴ乱码一区二区在线播放 | 两个人视频免费观看高清| 亚洲av片天天在线观看| 一级a爱片免费观看的视频| 精品人妻1区二区| 一边摸一边抽搐一进一小说| 色播在线永久视频| 日韩欧美三级三区| 成人三级做爰电影| 国产精品影院久久| 18禁裸乳无遮挡免费网站照片 | 成人亚洲精品av一区二区| 色哟哟哟哟哟哟| 精品免费久久久久久久清纯| 真人做人爱边吃奶动态| 黄色女人牲交| 操美女的视频在线观看| 亚洲五月色婷婷综合| а√天堂www在线а√下载| 中文字幕高清在线视频| 18禁国产床啪视频网站| 国产亚洲欧美98| 纯流量卡能插随身wifi吗| 在线观看午夜福利视频| 国产精品乱码一区二三区的特点 | 黄色毛片三级朝国网站| 在线观看免费视频网站a站| 女警被强在线播放| 久久精品91蜜桃| 一夜夜www| 可以在线观看毛片的网站| 午夜免费观看网址| 女人精品久久久久毛片| 亚洲国产欧美日韩在线播放| 久久亚洲真实| 免费人成视频x8x8入口观看| 韩国av一区二区三区四区| 色综合站精品国产| 午夜免费鲁丝| 美女高潮喷水抽搐中文字幕| 国产区一区二久久| 自线自在国产av| 青草久久国产| 精品久久久精品久久久| 琪琪午夜伦伦电影理论片6080| 午夜久久久久精精品| 国产xxxxx性猛交| 97超级碰碰碰精品色视频在线观看| 午夜精品国产一区二区电影| 久久精品成人免费网站| 久久精品亚洲熟妇少妇任你| 叶爱在线成人免费视频播放| 啦啦啦韩国在线观看视频| 免费在线观看完整版高清| 法律面前人人平等表现在哪些方面| 欧美亚洲日本最大视频资源| 一级,二级,三级黄色视频| 韩国av一区二区三区四区| 国产成年人精品一区二区| 国产亚洲欧美精品永久| 午夜影院日韩av| 欧美人与性动交α欧美精品济南到| 午夜激情av网站| 亚洲av第一区精品v没综合| 韩国av一区二区三区四区| 国产三级黄色录像| 欧美久久黑人一区二区| 日本 av在线| 亚洲第一欧美日韩一区二区三区| 神马国产精品三级电影在线观看 | 亚洲欧美日韩无卡精品| 很黄的视频免费| 91字幕亚洲| 亚洲色图 男人天堂 中文字幕| 欧美成人午夜精品| 老鸭窝网址在线观看| 欧美精品啪啪一区二区三区| 制服丝袜大香蕉在线| 91在线观看av| 黑丝袜美女国产一区| 亚洲 欧美一区二区三区| 1024香蕉在线观看| 亚洲精品国产一区二区精华液| 国产又爽黄色视频| 欧美丝袜亚洲另类 | av天堂久久9| 黄色女人牲交| 午夜成年电影在线免费观看| 美女大奶头视频| 两人在一起打扑克的视频| www.自偷自拍.com| 一区在线观看完整版| 亚洲一区二区三区色噜噜| 嫁个100分男人电影在线观看| 少妇 在线观看| 人成视频在线观看免费观看| 久久精品亚洲精品国产色婷小说| 午夜福利免费观看在线| 91精品国产国语对白视频| 久久久国产精品麻豆| 国产精品野战在线观看| 国产精品日韩av在线免费观看 | 欧美黑人欧美精品刺激| 高潮久久久久久久久久久不卡| www.999成人在线观看| 最近最新中文字幕大全电影3 | 91字幕亚洲| 日本vs欧美在线观看视频| www.自偷自拍.com| 桃色一区二区三区在线观看| 神马国产精品三级电影在线观看 | 亚洲中文字幕日韩| 久久性视频一级片| 少妇 在线观看| 国产精华一区二区三区| 亚洲 欧美 日韩 在线 免费| 欧美乱妇无乱码| 色哟哟哟哟哟哟| 久久香蕉激情| 精品电影一区二区在线| 精品熟女少妇八av免费久了| 男人舔女人的私密视频| 美女 人体艺术 gogo| 国产野战对白在线观看| 老司机午夜福利在线观看视频| 一区福利在线观看| 午夜久久久久精精品| 亚洲久久久国产精品| 啦啦啦韩国在线观看视频| 成人免费观看视频高清| 亚洲午夜理论影院| 日本免费a在线| 国产精品国产高清国产av| av超薄肉色丝袜交足视频| 欧美绝顶高潮抽搐喷水| 男女下面插进去视频免费观看| 亚洲色图av天堂| 黄片大片在线免费观看| 少妇粗大呻吟视频| 桃色一区二区三区在线观看| 在线视频色国产色| 成人18禁在线播放| 国产蜜桃级精品一区二区三区| 国产高清videossex| 欧美老熟妇乱子伦牲交| 亚洲国产精品久久男人天堂| 国产精品免费一区二区三区在线| 亚洲欧美日韩高清在线视频| 亚洲中文av在线| 成人永久免费在线观看视频| 免费人成视频x8x8入口观看| 亚洲中文字幕日韩| 欧美人与性动交α欧美精品济南到| 亚洲av成人不卡在线观看播放网| 老司机午夜十八禁免费视频| 亚洲少妇的诱惑av| 色综合欧美亚洲国产小说| 色综合亚洲欧美另类图片| 我的亚洲天堂| 亚洲人成77777在线视频| 国产伦一二天堂av在线观看| 国产精品久久久久久人妻精品电影| 亚洲人成网站在线播放欧美日韩| 亚洲色图 男人天堂 中文字幕| 乱人伦中国视频| 国产亚洲欧美精品永久| 美女高潮喷水抽搐中文字幕| 亚洲精品av麻豆狂野| 免费看美女性在线毛片视频| 曰老女人黄片| 亚洲五月婷婷丁香| ponron亚洲| 性欧美人与动物交配| 国产蜜桃级精品一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 日韩视频一区二区在线观看| 一边摸一边抽搐一进一出视频| 在线观看免费日韩欧美大片| 乱人伦中国视频| 国产精品久久电影中文字幕|