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

    基于參數(shù)的可視化裂變堆芯蒙特卡羅自動(dòng)建模方法

    2016-08-25 01:00:27程夢(mèng)云胡麗琴
    核技術(shù) 2016年6期
    關(guān)鍵詞:蒙特卡羅層次化堆芯

    甘 佺 吳 斌 宋 婧 程夢(mèng)云 胡麗琴

    ?

    基于參數(shù)的可視化裂變堆芯蒙特卡羅自動(dòng)建模方法

    甘 佺1,2吳 斌2宋 婧2程夢(mèng)云2胡麗琴1,2

    1(中國(guó)科學(xué)技術(shù)大學(xué) 核科學(xué)技術(shù)學(xué)院 合肥 230027)2(中國(guó)科學(xué)院核能安全技術(shù)研究所 中子輸運(yùn)理論與輻射安全重點(diǎn)實(shí)驗(yàn)室 合肥 230031)

    蒙特卡羅程序已經(jīng)廣泛應(yīng)用在裂變反應(yīng)堆設(shè)計(jì)和驗(yàn)證過(guò)程中,快速獲得高效的計(jì)算模型可以有效縮短反應(yīng)堆的設(shè)計(jì)周期。本研究提出并實(shí)現(xiàn)了一種裂變堆芯快速蒙特卡羅建模的方法,該方法基于參數(shù)可視化和層次化兩種建模思想快速構(gòu)建出精細(xì)裂變堆芯計(jì)算機(jī)輔助設(shè)計(jì)(Computer Aided Design, CAD)模型且將其快速轉(zhuǎn)換成蒙特卡羅計(jì)算模型,同時(shí)采用一種新的堆芯分段管理方法實(shí)現(xiàn)了大規(guī)模裂變堆模型流暢交互。基于此方法快速構(gòu)建了加速器驅(qū)動(dòng)次臨界反應(yīng)堆(Accelerator Driven Sub-critical System, ADS)的精細(xì)堆芯模型,通過(guò)與蒙特卡羅程序計(jì)算的參考結(jié)果進(jìn)行對(duì)比,證明了此建模方法的高效性和可靠性。

    蒙特卡羅,裂變堆,參數(shù)可視化,計(jì)算機(jī)輔助設(shè)計(jì),計(jì)算模型,堆芯分段管理

    隨著計(jì)算機(jī)性能的不斷提高,在核反應(yīng)堆計(jì)算中,為獲得更加精確的結(jié)果,得到更高的設(shè)計(jì)裕度,降低設(shè)計(jì)成本,高保真的計(jì)算能力已經(jīng)成為新一代計(jì)算程序的發(fā)展方向[1?2]。蒙特卡羅計(jì)算程序由于其幾何適應(yīng)性強(qiáng)等特點(diǎn),更易于滿足反應(yīng)堆高保真計(jì)算,并逐漸在反應(yīng)堆設(shè)計(jì)和驗(yàn)證領(lǐng)域獲得更廣泛的應(yīng)用。蒙特卡羅計(jì)算程序一般使用構(gòu)造實(shí)體幾何 (Constructive Solid Geometry, CSG)作為其幾何描述方式,并需要專業(yè)人員手工文本建模。不過(guò),如今高保真的裂變反應(yīng)堆可能需要模擬千萬(wàn)量級(jí)的反應(yīng)區(qū)域,例如精細(xì)壓水堆(Benchmark for Evaluation and Validation of Reactor Simulations, BEAVRS)模型[3],且隨著新型反應(yīng)堆的不斷問(wèn)世,裂變堆的幾何結(jié)構(gòu)也開(kāi)始變得更加復(fù)雜,再加上在實(shí)際工程應(yīng)用中反應(yīng)堆的模型需要反復(fù)的迭代修改,以致手工建模這樣不直觀的方式已經(jīng)無(wú)法滿足目前需求,并且其正確性也難以得到有效保證。為縮短基于蒙特卡羅方法的反應(yīng)堆設(shè)計(jì)與建模周期,建立一種能快速、高效且交互的裂變堆建模程序,是蒙特卡羅計(jì)算方法在裂變領(lǐng)域應(yīng)用的重大需求。

    鑒于以上需求,很多基于三維計(jì)算機(jī)輔助設(shè)計(jì)(Computer Aided Design, CAD)模型的自動(dòng)建模系統(tǒng)都得到了發(fā)展,此種系統(tǒng)能夠?qū)⑷SCAD模型采用邊界表示法(Boundary Representation, BREP)描述的幾何自動(dòng)轉(zhuǎn)換成蒙特卡羅計(jì)算模型[4]。不過(guò)傳統(tǒng)的轉(zhuǎn)換方法都需要將每一個(gè)BREP方法描述的實(shí)體模型分解成凸實(shí)體,再進(jìn)一步分解成半空間面。這種傳統(tǒng)的方法始終專注于復(fù)雜模型的分解,經(jīng)過(guò)多年的研究和發(fā)展,目前國(guó)內(nèi)外提出了多種先進(jìn)的方法,例如FDS團(tuán)隊(duì)自主開(kāi)發(fā)的多物理耦合分析自動(dòng)建模軟件SuperMC/MCAM[5?8]中的樣條面處理方法等,并已成功應(yīng)用于聚變反應(yīng)堆設(shè)計(jì)和核分析領(lǐng)域[9?16]。但針對(duì)裂變堆這種幾何規(guī)模大并具有明顯層次結(jié)構(gòu)和重復(fù)結(jié)構(gòu)的模型而言[17],此種方法面臨著內(nèi)存消耗大、計(jì)算硬件要求高、轉(zhuǎn)換效率低和生成模型幾何表達(dá)過(guò)度冗余等問(wèn)題。尤其是對(duì)大規(guī)模的裂變堆芯模型,例如多于千萬(wàn)量級(jí)柵元數(shù)目的BEAVRES反應(yīng)堆模型,傳統(tǒng)的轉(zhuǎn)換方法已經(jīng)無(wú)法適用。因此本研究針對(duì)裂變堆芯的特點(diǎn)發(fā)展了一種新的自動(dòng)建模方法,此方法可以通過(guò)少量的基本體按照層級(jí)關(guān)系和陣列關(guān)系構(gòu)建整個(gè)模型,使蒙特卡羅計(jì)算模型的生成過(guò)程僅需要解析少量基本體,從而極大提升了裂變堆蒙特卡羅模型的生成效率,并且生成的計(jì)算模型中將自動(dòng)保留層級(jí)關(guān)系和鄰居關(guān)系,使得新構(gòu)建的模型能獲得更高的計(jì)算效率。

    1 基于參數(shù)的可視化快速建模方法

    在本文中所謂基于參數(shù)的可視化建模方法就是指通過(guò)與圖形界面交互的方式輸入模型參數(shù),然后基于輸入?yún)?shù)構(gòu)建三維可視化裂變堆模型,或進(jìn)行迭代修改并最終生成完整蒙特卡羅計(jì)算模型的自動(dòng)化建模方法。此方法的具體流程如圖1所示。首先采用了層次化的建模思想,在參數(shù)化的基礎(chǔ)上分層構(gòu)造模型部件,然后通過(guò)層次化關(guān)系將所有參數(shù)組合成一種層級(jí)結(jié)構(gòu)樹(shù)。層級(jí)結(jié)構(gòu)樹(shù)則對(duì)整體模型進(jìn)行控制使得建模過(guò)程更加方便和精確?;趯蛹?jí)結(jié)構(gòu)樹(shù)能對(duì)大模型進(jìn)行分段渲染,獲取更高效的交互性。用戶能采用靈活的分段可視化方法觀察模型,如果判定模型沒(méi)有錯(cuò)誤則使用多蒙特卡羅語(yǔ)義適配器[18],將層級(jí)結(jié)構(gòu)樹(shù)中保存的模型和參數(shù)轉(zhuǎn)換成蒙特卡羅計(jì)算模型,并進(jìn)行蒙特卡羅計(jì)算。以下各節(jié)將分4個(gè)方面對(duì)該方法進(jìn)行更加詳細(xì)描述。

    圖1 參數(shù)化可視化迭代建模方法Fig.1 Parameter-based and visual modeling method.

    1.1 層次化建模方法

    裂變堆模型按層次可劃分為堆芯層(包括外屏蔽模型)、組件層、元件層等,所謂的層次化建模方法,就是通過(guò)模型中的層級(jí)關(guān)系分層構(gòu)建模型,優(yōu)點(diǎn)是建模的過(guò)程更加靈活,且層級(jí)之間相互獨(dú)立,并且得到的模型能夠保存層級(jí)和鄰居關(guān)系,提高蒙特卡羅程序計(jì)算時(shí)幾何處理的速度[19?22]。

    圖2為基于層次化建模方法的堆芯模型的建模過(guò)程示例。在圖2中U表示整體部件,C表示子部件,LAT表示一種陣列組合關(guān)系[23],其中U0代表為整體堆芯模型,包含有C1(堆芯燃料)、C2(外層冷卻劑)和C3(堆芯外屏蔽)。在建模過(guò)程中,用戶可以逐層地對(duì)模型進(jìn)行設(shè)計(jì)。例如C1引用一個(gè)半徑為20 cm、高為100 cm的圓柱,在該區(qū)域內(nèi)陣列填充燃料組件U1、反射組件U2、屏蔽組件U3、吸收棒U4、調(diào)節(jié)棒U5、補(bǔ)償棒U6,并按照LAT0的規(guī)則進(jìn)行陣列。其中燃料組件U1中包含了C4、C5、C6區(qū)域,C4為燃料區(qū)域,C5為包殼,C6為外層的冷卻劑。完成構(gòu)建如上的層級(jí)結(jié)構(gòu)關(guān)系后,可以生成總體的第一層和第二層堆芯模型,并得到對(duì)應(yīng)的層級(jí)結(jié)構(gòu)。如果堆芯需要進(jìn)一步精細(xì)建模,則在C4區(qū)域中填充燃料棒陣列進(jìn)行第三層建模,將燃料棒U7按照新的陣列規(guī)則進(jìn)行填充。在圖2中可以看到,其燃料棒由上下端塞C1、氣室C2、核燃料C3、C4和包殼C5組成,U7通過(guò)自身的不斷陣列填充進(jìn)C4中,得到元件級(jí)別的組件模型。

    如圖2所示,各節(jié)點(diǎn)間關(guān)系為模型的完成構(gòu)建后的層級(jí)關(guān)系,并明顯被劃分為三層結(jié)構(gòu),分別是堆芯層、組件層和最后構(gòu)建的元件層。其中,LAT節(jié)點(diǎn)中保留有各組件間鄰居關(guān)系,在GDML格式的計(jì)算模型中,C節(jié)點(diǎn)將轉(zhuǎn)換為邏輯體用來(lái)描述部件的幾何,并可以引用下層物理體構(gòu)建層級(jí)關(guān)系,U節(jié)點(diǎn)為物理體,物理體中引用唯一的邏輯體并將其放置在指定位置,同時(shí)賦予相應(yīng)的物理屬性,通過(guò)邏輯體和物理體的相互調(diào)用,從而實(shí)現(xiàn)層級(jí)關(guān)系的表達(dá),而鄰居關(guān)系將通過(guò)陣列矩陣來(lái)進(jìn)行表達(dá),并給出鄰居之間的間距和方向。這樣在粒子位置跟蹤過(guò)程中能直接判斷出穿出組件的粒子到達(dá)哪個(gè)新的組件。

    圖2 層次化堆芯模型構(gòu)建過(guò)程Fig.2 Procedure of hierarchical modeling method.

    1.2 分段模型管理與交互

    高保真裂變堆模型包含大量幾何體,在圖形交互時(shí),一次性渲染整體模型的效率十分低下,且對(duì)計(jì)算機(jī)硬件資源的要求較高,一般情況下,無(wú)法滿足實(shí)時(shí)交互操作的需求。為此,本研究提出了一種基于層次化模型的分段管理與渲染方法,包括:層級(jí)分段、部件分段、材料分段、混合分段。圖3描述了在完成層次化參數(shù)建模之后,該方法基于構(gòu)建完成的層級(jí)結(jié)構(gòu),創(chuàng)建各分段的過(guò)程。

    圖3 層次化模型分段管理方法Fig.3 Method and segments of hierarchical management.

    由圖3可以看出,層級(jí)結(jié)構(gòu)中將模型劃分為不同層次,例如堆芯整體模型在第一層,組件模型在第二層,元件級(jí)模型在第三層。在渲染單元樹(shù)的構(gòu)建過(guò)程中首先將根據(jù)層次構(gòu)建分支,并在分支內(nèi)部構(gòu)建相應(yīng)的部件節(jié)點(diǎn)及陣列關(guān)系。例如第一層的部件節(jié)點(diǎn)包含堆芯組件排布區(qū)域幾何,和堆芯外層冷卻劑與屏蔽結(jié)構(gòu)幾何;第二層為組件層包含所有類型組件幾何和組件在堆芯內(nèi)的陣列關(guān)系;第三層為元件層,包含各種元件幾何和元件在上層組件內(nèi)的陣列關(guān)系。其中元件類型可以在不同的陣列關(guān)系中進(jìn)行復(fù)用。

    從圖3可見(jiàn),在劃分為各層之后,層級(jí)之間保持互相連接關(guān)系,也就是說(shuō),同一個(gè)部件除了在本層有幾何外,在下層仍有屬于其的幾何。完成所有的渲染單元構(gòu)建后,用戶可以通過(guò)選擇不同的層次進(jìn)行管理,所選層次以下的幾何和陣列關(guān)系均被略去。也可以同時(shí)選取部件和層級(jí)屬性,決定提取出該部件在第幾層的幾何,此外還可以增加材料屬性的判斷,繼續(xù)篩選出更少量的幾何。于是,分段管理的實(shí)質(zhì)就是實(shí)現(xiàn)不同條件下的模型快速篩選。

    在圖3中,每個(gè)類型部件可以被不同關(guān)系所復(fù)用,每個(gè)幾何體可以被不同類型所復(fù)用,而幾何體可以和任意的材料相組合。不過(guò)每層之間幾何互相獨(dú)立,以防止改變層級(jí)獨(dú)立性,組件類型在調(diào)用一個(gè)幾何體的同時(shí)必須要明確與之匹配的材料并將材料編號(hào)作為自身屬性?;诖朔椒ㄍㄟ^(guò)層次、材料和部件的篩選,盡可能地減少了幾何體數(shù)量,實(shí)現(xiàn)按需提供給用戶最高效的模型,滿足了高效的交互式建模。

    1.3 CAD模型和蒙特卡羅計(jì)算模型的復(fù)用

    CAD模型的復(fù)用支持將通用CAD格式的裂變堆部件模型直接導(dǎo)入實(shí)現(xiàn)本研究方法的程序中,并指定一個(gè)柵元節(jié)點(diǎn)(C)引用該導(dǎo)入的模型。

    蒙特卡羅計(jì)算模型的復(fù)用支持直接導(dǎo)入蒙特卡羅計(jì)算模型的文本文件,自動(dòng)構(gòu)建相應(yīng)的層次結(jié)構(gòu)和部件節(jié)點(diǎn)?;诖?,用戶可以使用分段模型管理方法進(jìn)行可視化交互迭代設(shè)計(jì)。在復(fù)用CAD模型和蒙特卡羅計(jì)算模型生成新的蒙特卡羅計(jì)算文件時(shí),對(duì)于不規(guī)則的幾何體將采用基于分解的模型轉(zhuǎn)換算法,確保任意精細(xì)、復(fù)雜的裂變堆計(jì)算模型的創(chuàng)建。

    1.4 蒙特卡羅計(jì)算模型的生成

    基于以上的分段式模型管理方法,裂變反應(yīng)堆中的模型參數(shù)和模型之間的關(guān)系都保存在各層節(jié)點(diǎn)中,與分段可視化交互過(guò)程比較類似的是,實(shí)現(xiàn)本研究方法的程序同樣可以選取不同層級(jí)和部件生成對(duì)應(yīng)的蒙特卡羅計(jì)算模型,并在計(jì)算模型之中保留了部件之間的層級(jí)關(guān)系、陣列關(guān)系和鄰居關(guān)系,使蒙特卡羅計(jì)算幾何跟蹤更高效。

    2 測(cè)試?yán)}與測(cè)試方法

    2011年,中國(guó)科學(xué)院?jiǎn)?dòng)戰(zhàn)略性先導(dǎo)科技專項(xiàng)“先進(jìn)核裂變能-ADS (Accelerator Driven Sub-critical System)嬗變系統(tǒng)”,致力于自主發(fā)展ADS系統(tǒng)從實(shí)驗(yàn)裝置到示范裝置的全部核心技術(shù)和系統(tǒng)集成技術(shù)[24?25],本研究參考該反應(yīng)堆設(shè)計(jì)資料,構(gòu)建了ADS反應(yīng)堆蒙特卡羅計(jì)算模型。圖4?7為該模型的構(gòu)建過(guò)程,其中圖4為中子源部分幾何構(gòu)建;圖5為該堆芯組件排布設(shè)計(jì)界面;圖6為燃料組件內(nèi)部燃料棒排布設(shè)計(jì);圖7為該模型完成構(gòu)建后的可視化結(jié)果。建模過(guò)程總耗時(shí)不超過(guò)1h,如果復(fù)用已經(jīng)完成構(gòu)建的CAD部件模型或者蒙特卡羅計(jì)算模型,建模時(shí)間可縮短至幾分鐘。而傳統(tǒng)的手工建模需要將近一周時(shí)間并很難保證模型正確性。在測(cè)試中,分別使用了本研究方法構(gòu)建SuperMC[26]計(jì)算模型后和已經(jīng)得到充分驗(yàn)證的手工構(gòu)建的MCNP計(jì)算模型進(jìn)行了3000代、每代10000粒子的臨界源計(jì)算,并對(duì)比了兩次計(jì)算的eff和不同組件的通量的結(jié)果。

    圖4 ADS中子源的構(gòu)建Fig.4 Creating of neutron source in the ADS reactor.

    圖5 ADS鉛基反應(yīng)堆堆芯組件排布設(shè)計(jì)Fig.5 Arrangement of assemblies in ADS Lead-basedreactor core.

    圖6 ADS鉛基反應(yīng)堆燃料組件內(nèi)部燃料棒排布Fig.6 Arrangement of fuel pins in ADS Lead-basedreactor core.

    圖7 ADS鉛基快堆精細(xì)模型Fig.7 Detailed models of ADS Lead-based reactor core.

    3 結(jié)果與討論

    表1、2顯示了基于本文§2中方法和例題得到的測(cè)試結(jié)果,其中表1顯示了eff計(jì)算結(jié)果,表2為通量計(jì)算結(jié)果對(duì)比。本研究方法與MCNP參考計(jì)算通量偏差低于1%,eff偏差低于4 pcm,在MCNP中一般采用面幾何的描述方式,而通過(guò)本研究方法自動(dòng)構(gòu)建的SuperMC模型一般為面體混合描述方式(與面幾何采用半空間的形式相比,面體混合幾何還能使用各種基本體來(lái)描述其幾何,例如長(zhǎng)方體、圓柱、球體和六棱柱等),于是在蒙特卡羅計(jì)算過(guò)程中,其幾何處理的方法會(huì)產(chǎn)生差異,帶來(lái)一定的系統(tǒng)性誤差,并且所有偏差都在統(tǒng)計(jì)誤差范圍內(nèi)。此外,從建模過(guò)程描述可以看出,通過(guò)本文方法構(gòu)建蒙特卡羅模型,方便快捷,出錯(cuò)可能小,且支持快速迭代修改,與手工建模相比極大地提高了建模速度,與傳統(tǒng)轉(zhuǎn)換方便相比,提高了轉(zhuǎn)換效率和得到了更高效的模型(保留了層級(jí)結(jié)構(gòu)和重復(fù)結(jié)構(gòu),消除了冗余的幾何)。綜上所述,可以證明本研究自動(dòng)快速構(gòu)建的計(jì)算模型的方法是高效和正確的。


    表1eff計(jì)算結(jié)果對(duì)比
    Table 1 Comparison ofeffcalculation result.

    計(jì)算量Tally

    MCNP參考值Reference result

    本研究計(jì)算結(jié)果Calculation result

    1.00027 (0.00010)

    1.00022 (0.00011)

    表2 通量計(jì)算結(jié)果對(duì)比
    Table 2 Comparison of flux calculation result.

    柵元Cell

    MCNP參考值Reference result

    本方法Calculation result

    偏差Deviation

    1

    1.87759×10(0.0039)

    1.89434×10(0.0039)

    0.00892

    2

    6.52745×10(0.0007)

    6.49945×10(0.0007)

    0.00429

    4

    3.80315×10(0.0007)

    3.78466×10(0.0007)

    0.00486

    5

    6.82774×10(0.0005)

    6.79260×10(0.0005)

    0.00515

    6

    1.02052×10(0.0001)

    1.02063×10(0.0001)

    0.00011

    7

    3.37869×10(0.0007)

    3.39096×10(0.0007)

    0.00363

    8

    9.49932×10(0.0009)

    9.52201×10(0.0009)

    0.00239

    10

    3.06813×10(0.0010)

    3.07518×10(0.0010)

    0.00230

    150

    7.82561×10(0.0008)

    7.80693×10(0.0008)

    0.00239

    151

    4.36463×10(0.0008)

    4.35360×10(0.0008)

    0.00253

    152

    1.22123×10(0.0008)

    1.21857×10(0.0008)

    0.00218

    153

    3.94259×10(0.0008)

    3.93440×10(0.0008)

    0.00208

    4 結(jié)語(yǔ)

    蒙特卡羅計(jì)算程序在反應(yīng)堆設(shè)計(jì)和驗(yàn)證中應(yīng)用愈加廣泛,新型反應(yīng)堆的建模難度也將越來(lái)越大,而且在工程應(yīng)用中需要對(duì)模型進(jìn)行反復(fù)的迭代設(shè)計(jì),手工文本建模耗時(shí)易出錯(cuò),難以滿足其需求。本研究提出了一種參數(shù)可視化的層次化建模方法,實(shí)現(xiàn)了對(duì)精細(xì)裂變堆模型的快速構(gòu)建與高效交互。本方法能夠便捷地構(gòu)建出棒級(jí)的ADS反應(yīng)堆堆芯模型,與手工創(chuàng)建的MCNP模型的計(jì)算結(jié)果進(jìn)行對(duì)比,誤差均在統(tǒng)計(jì)誤差范圍內(nèi),證明了此方法的準(zhǔn)確性。在將來(lái),基于本文方法將繼續(xù)研究出針對(duì)不同反應(yīng)堆堆型定制的一鍵式迭代設(shè)計(jì)功能,從根本上解決蒙特卡羅方法建模耗費(fèi)過(guò)多人力成本問(wèn)題。

    1 Smith K, Forget B. Challenges in the development of high-fidelity LWR core neutronics tools[C]. International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, Idaho, USA, May 5?9, 2013

    2 Vaz P. Monte Carlo methods and techniques status and prospect for future evolution[J]. Applied Radiation and Isotopes, 2010, 68(4): 536?541. DOI: 10.1016/j.apradiso. 2009.10.020

    3 Horelik N, Herman B. Benchmark for evaluation and validation of reactor simulations (BEAVRS)[C]. International Conference on Mathematics and Computational Methods Applied to Nuclear Science & Engineering, Sun Valley, Idaho, USA, May 5?9, 2013

    4 吳宜燦, 李瑩, 盧磊, 等. 蒙特卡羅粒子輸運(yùn)計(jì)算自動(dòng)建模程序系統(tǒng)的研究與發(fā)展[J]. 核科學(xué)與工程, 2006, 26(1): 20?27. DOI: 10.3321/j.issn:0258-0918.2006.01. 004
    WU Yican, LI Ying, LU Lei,. Research and development of the automatic modeling system for Monte Carlo particle transport simulation[J]. Nuclear Science and Engineering, 2006, 26(1): 20?27. DOI: 10.3321/j.issn: 0258-0918.2006.01.004

    5 吳宜燦, 胡麗琴, 羅月童, 等. 先進(jìn)核能系統(tǒng)設(shè)計(jì)分析軟件與數(shù)據(jù)庫(kù)研發(fā)進(jìn)展[J]. 核科學(xué)與工程, 2010, 30(1): 60?69. DOI: 0258-0918(2010)01-0055-09
    WU Yican, HU Liqin, LUO Yuetong,Development of design and analysis software for advanced nuclear systems[J]. Nuclear Science and Engineering, 2010, 30(1): 60?69. DOI: 0258-0918 (2010)01-0055-09

    6 Wu Y C, FDS Team. CAD-based interface programs for fusion neutron transport simulation[J]. Fusion Engineering and Design, 2009, 84(7?11): 1987?1992. DOI: 10.1016/j.fusengdes.2008.12.041

    7 吳宜燦, 李靜驚, 李瑩. 大型集成多功能中子學(xué)計(jì)算與分析系統(tǒng)VisualBUS的研究與發(fā)展[J]. 核科學(xué)與工程, 2007, 27(4): 365?373. DOI: 10.3321/j.issn:0258-0918. 2007.04.014
    WU Yican, LI Jingjing, LI Ying. An integrated multi-functional neutronics calculation and analysis code system: VisualBUS[J]. Nuclear Science and Engineering, 2007, 27(4): 365?373. DOI: 10.3321/j.issn:0258-0918. 2007.04.014

    8 Li Y, Lu L, Ding A P,Benchmarking of MCAM4.0 with the ITER 3D model[J]. Fusion Engineering and Design, 2007, 82(15): 2861?2866. DOI: 10.1016/j. fusengdes.2007.02.022

    9 Wu Y C, FDS Team. Conceptual design of the China fusion power plant FDS-II[J]. Fusion Engineering and Design, 2008, 83(10?12): 1683?1689. DOI: 10.1016/j. fusengdes.2008.06.048

    10 Wu Y C, FDS Team. Conceptual design activities of FDS series fusion power plants in China[J]. Fusion Engineering and Design, 2006, 81(23?24): 2713?2718. DOI: 10.1016/j.fusengdes.2006.07.068

    11 Wu Y C, FDS Team. Conceptual design and testing strategy of a dual functional lithium-lead test blanket module in ITER and EAST[J]. Nuclear Fusion, 2007, 47(11): 1533?1539

    12 Wu Y C, Zheng S, Zhu X,Conceptual design of the fusion-driven subcritical system FDS-I[J]. Fusion Engineering and Design, 2006, 81(8?14): 1305?1311. DOI: 10.1016/j.fusengdes.2005.10.015

    13 Wu Y C, FDS Team. Design analysis of the China dual-functional lithium lead (DFLL) test blanket module in ITER[J]. Fusion Engineering and Design, 2007, 82(15): 1893?1903. DOI: 10.1016/j.fusengdes.2007.08.012

    14 Wu Y C, Jiang J Q, Wang M H,. A fusion-driven subcritical system concept based on viable technologies[J]. Nuclear Fusion, 2011, 51(10): 103036

    15 Wu Y C, FDS Team. Design status and development strategy of China liquid lithium-lead blankets and related material technology[J]. Journal of Nuclear Materials, 2007, 367?370: 1410?1415. DOI: 10.1016/j.jnucmat. 2007.04.031

    16 Song J, Sun G Y, Chen Z P,. Benchmarking of CAD-based SuperMC with ITER benchmark model[J]. Fusion Engineering and Design, 2014, 89(11): 2499?2503. DOI: 10.1016/j.fusengdes.2014.05.003

    17 Jain R, Tautges T J. Generating unstructured nuclear reactor core meshes in parallel[J]. Procedia Engineering, 2014, 82: 351?364. DOI: 10.1016/ j.proeng.2014.10.396

    18 吳宜燦, 俞盛朋, 程夢(mèng)云, 等. 多物理耦合分析自動(dòng)建模軟件SuperMC/MCAM5.2設(shè)計(jì)與實(shí)現(xiàn)[J]. 原子能科學(xué)技術(shù), 2015, 49(z1): 1?6. DOI: 10.7538/yzk.2015.49. S0.0001
    WU Yican, YU Shengpeng, CHENG Mengyun,Design and implementation of multi-physics coupling analysis modeling program SuperMC/MCAM5.2[J]. Atomic Energy Science and Technology, 2015, 49(z1): 1?6. DOI: 10.7538/yzk. 2015.49.S0.0001

    19 Donovan T, Tyburski L. Geometrics representations in the developmental Monte Carlo transport code MC21[R]. New York, USA: Knolls Atomic Power Laboratory, 2006. DOI: http://dx.doi.org/10.1051/snamc/201402505

    20 Geant4 Collaboration.[OL]. 2015. http://geant4.web. cern.ch

    21 Nease R N, Millman D L, Griesheimer D P,Geometric templates for improved tracking performance in Monte Carlo codes[C]. Joint International Conference on Supercomputing in Nuclear Applications and Monte Carlo, Paris, France, October 27?31, 2013

    22 Brown F B. Fundamentals of Monte Carlo transport[R]. USA: Los Alamos National Laboratory, 2004

    23 X-5 Monte Carlo Team. MCNP - a general Monte Carlo N-particle transport code[R]. Version 5. Los Alamos National Laboratory, 2003

    24 吳宜燦, 柏云清, 宋勇, 等. 中國(guó)鉛基反應(yīng)堆概念設(shè)計(jì)研究[J]. 核科學(xué)與工程, 2014, 34(2): 201?208
    WU Yican, BAI Yunqing, SONG Yong,Conceptual design of China Lead-based Research Reactor CLEAR-I[J]. Nuclear Science and Engineering, 2014, 34(2): 201?208

    25 Wu Y C, Bai Y Q, Song Y,Development strategy and conceptual design of China Lead-based Research Reactor[J]. Annals Nuclear Energy, 2016, 87(2): 511?516. DOI: 10.1016/j.anucene.2015.08.015

    26 Wu Y C, Song J, Zheng H Q,. CAD-based Monte Carlo program for integrated simulation of nuclear system SuperMC[J]. Annals of Nuclear Energy, 2015, 82: 161?168. DOI: 10.1016/j.anucene.2014.08.058

    中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(No.XDA03040000)、國(guó)家自然科學(xué)基金、No.11305205、No.11305203、No.11405204)、國(guó)家磁約束核聚變能發(fā)展研究專項(xiàng)(No.2014GB112001)資助

    Supported by Strategic Priority Research Program of Chinese Academy of Sciences (No.XDA03040000), National Natural Science Foundation of China (No.91026004, No.11305205, No.11305203, No.11405204), National Special Program for ITER (No.2014GB112001)

    parameter-based and visual Monte Carlo modeling method of fission reactor core

    GAN Quan1,2WU Bin2SONG Jing2CHENG Mengyun2HU Liqin1,2

    1(School of Nuclear Science and Technology, University of Science and Technology of China, Hefei 230027, China)2(Key Laboratory of Neutronics and Radiation Safety, Institute of Nuclear Energy Safety Technology, Chinese Academy of Sciences, Hefei 230031, China)

    Background: The Monte Carlo (MC) method is widely used in fission reactor design, because of its stronger geometry adaptability and the precise calculation result. The high-fidelity full core simulation demands for detailed fission reactor models, which is hard to build by manual and the conventional Computer Aided Design (CAD)-based modeling method. Purpose: In order to support the rapid design of fission reactor core with MC method, and generate the detailed MC calculation models, a parameter-based and visual MC modeling method was developed in this study. Method: The method can create the detailed CAD models and convert them into MC models in a high-efficiency way. Meanwhile, the huge amounts of models in fission reactor core can be managed by different segments for supporting smooth interactions. Result: Furthermore, the method was validated by the test of creating the Accelerator Driven Sub-critical System (ADS) reactor models, and the results were agree very well by comparing the reference models. Conclusion: Depending on the test, the detailed fission models were created more conveniently than conventional method and the numerical calculation results proved the accuracy of the new method.

    Monte Carlo, Fission reactor, Parameter-based and visual modeling, CAD, Calculation model, Manage by segments

    GAN Quan, male, born in 1988, graduated from Shenyang Normal University in 2009, doctor student, engaged in radiation transport automatic modeling method, major in nuclear science and technology

    HU Liqin, E-mail: liqin.hu@fds.org.cn

    TL329.2

    10.11889/j.0253-3219.2016.hjs.39.060501

    甘佺,男,1988年出生,2009年畢業(yè)于沈陽(yáng)師范大學(xué),現(xiàn)為博士研究生,從事輻射輸運(yùn)自動(dòng)建模方法,核科學(xué)與技術(shù)專業(yè)

    胡麗琴,E-mail: liqin.hu@fds.org.cn

    2016-02-22,

    2016-03-21

    猜你喜歡
    蒙特卡羅層次化堆芯
    面向量化分塊壓縮感知的區(qū)域?qū)哟位A(yù)測(cè)編碼
    利用蒙特卡羅方法求解二重積分
    應(yīng)用CDAG方法進(jìn)行EPR機(jī)組的嚴(yán)重事故堆芯損傷研究
    基于Hoogenboom基準(zhǔn)模型的SuperMC全堆芯計(jì)算能力校驗(yàn)
    核技術(shù)(2016年4期)2016-08-22 09:05:32
    鐵路傳送網(wǎng)OTN設(shè)備互聯(lián)互通開(kāi)銷層次化處理研究
    壓水堆堆芯中應(yīng)用可燃毒物的兩個(gè)重要實(shí)驗(yàn)
    探討蒙特卡羅方法在解微分方程邊值問(wèn)題中的應(yīng)用
    艦船系統(tǒng)間電磁兼容性的層次化優(yōu)化方法
    基于層次化分類器的遙感圖像飛機(jī)目標(biāo)檢測(cè)
    復(fù)合型種子源125I-103Pd劑量場(chǎng)分布的蒙特卡羅模擬與實(shí)驗(yàn)測(cè)定
    同位素(2014年2期)2014-04-16 04:57:20
    在线观看三级黄色| 国产在视频线精品| 中文字幕免费在线视频6| 人人妻人人看人人澡| 肉色欧美久久久久久久蜜桃 | 极品少妇高潮喷水抽搐| 成人亚洲精品一区在线观看 | 国产成人91sexporn| 日本爱情动作片www.在线观看| 水蜜桃什么品种好| 男人添女人高潮全过程视频| av专区在线播放| 免费人成在线观看视频色| av国产精品久久久久影院| 国产视频首页在线观看| 超碰av人人做人人爽久久| 深爱激情五月婷婷| 亚洲精品国产成人久久av| 一边亲一边摸免费视频| 国产成人福利小说| 狂野欧美激情性xxxx在线观看| 亚洲av免费在线观看| 2021天堂中文幕一二区在线观| 22中文网久久字幕| 日本色播在线视频| 国产真实伦视频高清在线观看| 亚洲欧美日韩无卡精品| 我的老师免费观看完整版| 一级毛片黄色毛片免费观看视频| 人人妻人人看人人澡| 联通29元200g的流量卡| av在线老鸭窝| 乱码一卡2卡4卡精品| 欧美高清性xxxxhd video| 国产亚洲一区二区精品| 秋霞在线观看毛片| 午夜免费男女啪啪视频观看| 亚洲自拍偷在线| 99九九线精品视频在线观看视频| 香蕉精品网在线| 久久精品国产亚洲av涩爱| 亚洲av国产av综合av卡| 我要看日韩黄色一级片| 欧美成人精品欧美一级黄| 欧美+日韩+精品| 中文欧美无线码| 欧美一区二区亚洲| 狠狠精品人妻久久久久久综合| 成人高潮视频无遮挡免费网站| 久久国产乱子免费精品| 一个人观看的视频www高清免费观看| 亚洲精品中文字幕在线视频 | 亚洲va在线va天堂va国产| 直男gayav资源| 国产一区二区三区综合在线观看 | 国产欧美另类精品又又久久亚洲欧美| 69av精品久久久久久| 天堂俺去俺来也www色官网| 最近中文字幕高清免费大全6| 午夜精品一区二区三区免费看| tube8黄色片| 大又大粗又爽又黄少妇毛片口| 看十八女毛片水多多多| 三级国产精品欧美在线观看| 成人黄色视频免费在线看| 欧美bdsm另类| 亚洲怡红院男人天堂| 亚洲av福利一区| 午夜亚洲福利在线播放| 国产精品久久久久久精品电影小说 | 在线看a的网站| 2022亚洲国产成人精品| 日本wwww免费看| 校园人妻丝袜中文字幕| 男人爽女人下面视频在线观看| 18禁动态无遮挡网站| 晚上一个人看的免费电影| 国产精品国产av在线观看| 色5月婷婷丁香| 尾随美女入室| 身体一侧抽搐| 婷婷色综合大香蕉| 91aial.com中文字幕在线观看| 自拍偷自拍亚洲精品老妇| 久久99精品国语久久久| 看非洲黑人一级黄片| 91在线精品国自产拍蜜月| 欧美性感艳星| 国产免费一区二区三区四区乱码| 97热精品久久久久久| 亚洲欧美清纯卡通| 午夜老司机福利剧场| 毛片女人毛片| 久久久久精品性色| 色婷婷久久久亚洲欧美| 热re99久久精品国产66热6| 伊人久久国产一区二区| 亚洲精品一二三| 国产毛片a区久久久久| 插逼视频在线观看| 久久国产乱子免费精品| 日本wwww免费看| 大码成人一级视频| 国精品久久久久久国模美| 婷婷色av中文字幕| 夫妻性生交免费视频一级片| 久久久国产一区二区| 在线免费观看不下载黄p国产| 日韩成人伦理影院| 搡女人真爽免费视频火全软件| 嫩草影院新地址| 亚洲美女搞黄在线观看| 纵有疾风起免费观看全集完整版| 嫩草影院新地址| 超碰97精品在线观看| 免费av观看视频| 久久99精品国语久久久| 嫩草影院入口| 精品一区二区三卡| 成人黄色视频免费在线看| 亚洲国产精品成人综合色| av黄色大香蕉| 午夜福利在线在线| 毛片一级片免费看久久久久| 人妻 亚洲 视频| 国产精品偷伦视频观看了| 国产av国产精品国产| 蜜桃久久精品国产亚洲av| 美女脱内裤让男人舔精品视频| 国产av不卡久久| 国产精品国产三级国产专区5o| 久久久久国产网址| 亚洲av成人精品一二三区| 婷婷色综合www| a级一级毛片免费在线观看| 久久99蜜桃精品久久| 国产成人a∨麻豆精品| 国产精品麻豆人妻色哟哟久久| 久久热精品热| 卡戴珊不雅视频在线播放| 国产亚洲av片在线观看秒播厂| 国产成人精品久久久久久| 国产一区二区亚洲精品在线观看| 人妻系列 视频| www.av在线官网国产| 久久久久久久久大av| 久久久久网色| 看免费成人av毛片| 日韩av免费高清视频| 韩国高清视频一区二区三区| 欧美日韩综合久久久久久| 国产国拍精品亚洲av在线观看| 成年免费大片在线观看| 熟女电影av网| 深爱激情五月婷婷| 夜夜爽夜夜爽视频| 免费播放大片免费观看视频在线观看| 毛片女人毛片| 特级一级黄色大片| 国产亚洲91精品色在线| 国产91av在线免费观看| 亚洲一区二区三区欧美精品 | 日韩国内少妇激情av| 视频区图区小说| 在线观看免费高清a一片| 97精品久久久久久久久久精品| 一区二区三区乱码不卡18| 中文字幕制服av| 国产精品不卡视频一区二区| 黄片无遮挡物在线观看| 久久ye,这里只有精品| 干丝袜人妻中文字幕| 真实男女啪啪啪动态图| 夜夜看夜夜爽夜夜摸| 亚洲一级一片aⅴ在线观看| 97热精品久久久久久| 少妇 在线观看| 熟女av电影| 日本午夜av视频| 午夜福利视频精品| 直男gayav资源| 亚洲色图av天堂| xxx大片免费视频| 亚洲人与动物交配视频| 亚洲三级黄色毛片| 两个人的视频大全免费| 亚洲av免费在线观看| 午夜老司机福利剧场| 国产精品一区www在线观看| 建设人人有责人人尽责人人享有的 | 国产乱人视频| 91午夜精品亚洲一区二区三区| 搡女人真爽免费视频火全软件| 啦啦啦在线观看免费高清www| 免费观看性生交大片5| 丝袜脚勾引网站| 如何舔出高潮| 国产一区有黄有色的免费视频| 99久久九九国产精品国产免费| 视频区图区小说| 乱系列少妇在线播放| 日韩,欧美,国产一区二区三区| 五月玫瑰六月丁香| 欧美成人午夜免费资源| 亚洲真实伦在线观看| 欧美日韩国产mv在线观看视频 | av在线观看视频网站免费| 亚洲婷婷狠狠爱综合网| a级一级毛片免费在线观看| 日韩不卡一区二区三区视频在线| 秋霞在线观看毛片| av在线蜜桃| 国产高清有码在线观看视频| 91精品一卡2卡3卡4卡| 久久久成人免费电影| 少妇熟女欧美另类| 国产精品一区二区性色av| 午夜免费鲁丝| 97超碰精品成人国产| 少妇人妻 视频| 水蜜桃什么品种好| 成年女人看的毛片在线观看| 国产高清有码在线观看视频| 赤兔流量卡办理| 少妇人妻久久综合中文| 欧美成人精品欧美一级黄| 亚洲精品aⅴ在线观看| 欧美zozozo另类| 日韩欧美一区视频在线观看 | 国产午夜精品久久久久久一区二区三区| 插阴视频在线观看视频| 欧美激情久久久久久爽电影| 国产精品女同一区二区软件| 小蜜桃在线观看免费完整版高清| 少妇人妻久久综合中文| 婷婷色av中文字幕| 亚洲精品久久久久久婷婷小说| 欧美成人一区二区免费高清观看| 久久久久国产网址| 超碰av人人做人人爽久久| 欧美97在线视频| 26uuu在线亚洲综合色| 51国产日韩欧美| 在线a可以看的网站| 听说在线观看完整版免费高清| videos熟女内射| 亚洲成人精品中文字幕电影| 少妇人妻精品综合一区二区| 欧美精品一区二区大全| 天天躁日日操中文字幕| 看免费成人av毛片| 激情五月婷婷亚洲| 亚洲丝袜综合中文字幕| 免费不卡的大黄色大毛片视频在线观看| 香蕉精品网在线| 日韩人妻高清精品专区| 好男人视频免费观看在线| 国产成人一区二区在线| 久久人人爽人人爽人人片va| 亚洲不卡免费看| 久久99蜜桃精品久久| 麻豆成人午夜福利视频| 在线天堂最新版资源| 黄色怎么调成土黄色| a级一级毛片免费在线观看| 美女国产视频在线观看| 亚洲人成网站在线播| 免费高清在线观看视频在线观看| 亚洲四区av| 看十八女毛片水多多多| 日日摸夜夜添夜夜爱| 国产在线男女| 有码 亚洲区| 国产大屁股一区二区在线视频| 五月伊人婷婷丁香| 国产69精品久久久久777片| 三级国产精品欧美在线观看| 高清欧美精品videossex| 女人久久www免费人成看片| 听说在线观看完整版免费高清| 国产探花在线观看一区二区| 国产精品秋霞免费鲁丝片| 免费看av在线观看网站| 亚洲av国产av综合av卡| 简卡轻食公司| 永久免费av网站大全| 天天躁日日操中文字幕| 国产毛片a区久久久久| 欧美性猛交╳xxx乱大交人| 国产探花在线观看一区二区| 国产亚洲91精品色在线| 夜夜爽夜夜爽视频| 日日撸夜夜添| 国产乱人偷精品视频| eeuss影院久久| 成年免费大片在线观看| 免费看a级黄色片| 国产精品久久久久久精品电影小说 | 国产高潮美女av| 欧美成人一区二区免费高清观看| 国产精品一区二区在线观看99| 日韩成人av中文字幕在线观看| av专区在线播放| 各种免费的搞黄视频| 国产精品伦人一区二区| 国产 一区精品| 少妇人妻精品综合一区二区| 亚洲一区二区三区欧美精品 | 久久国产乱子免费精品| 国产精品久久久久久精品古装| 在线观看av片永久免费下载| 人妻少妇偷人精品九色| 在线观看av片永久免费下载| 最近2019中文字幕mv第一页| 男女啪啪激烈高潮av片| 最近2019中文字幕mv第一页| 男女边摸边吃奶| 热99国产精品久久久久久7| 欧美最新免费一区二区三区| 大片免费播放器 马上看| 一个人看视频在线观看www免费| 狂野欧美激情性xxxx在线观看| 精品视频人人做人人爽| 国产精品三级大全| 80岁老熟妇乱子伦牲交| 啦啦啦在线观看免费高清www| 观看免费一级毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品乱码久久久久久按摩| 亚洲最大成人手机在线| 亚洲av日韩在线播放| 国产在视频线精品| 欧美xxⅹ黑人| av女优亚洲男人天堂| 美女脱内裤让男人舔精品视频| 毛片一级片免费看久久久久| 夫妻性生交免费视频一级片| 18禁在线无遮挡免费观看视频| 黄色配什么色好看| 亚洲,一卡二卡三卡| 一个人看视频在线观看www免费| av在线天堂中文字幕| 国产又色又爽无遮挡免| av在线天堂中文字幕| 免费av观看视频| 99久久精品国产国产毛片| 三级男女做爰猛烈吃奶摸视频| 色播亚洲综合网| 人妻制服诱惑在线中文字幕| 69人妻影院| 国产探花极品一区二区| 亚洲,一卡二卡三卡| 国产一区二区三区综合在线观看 | 国产v大片淫在线免费观看| 91aial.com中文字幕在线观看| 麻豆精品久久久久久蜜桃| av网站免费在线观看视频| 亚洲欧美中文字幕日韩二区| 国产爱豆传媒在线观看| 少妇熟女欧美另类| 日韩欧美精品v在线| 日本黄色片子视频| 婷婷色综合大香蕉| 丝瓜视频免费看黄片| 国产真实伦视频高清在线观看| 建设人人有责人人尽责人人享有的 | 高清视频免费观看一区二区| 精品一区在线观看国产| 亚洲无线观看免费| 久久97久久精品| 不卡视频在线观看欧美| 丝袜美腿在线中文| 色网站视频免费| 亚洲国产最新在线播放| 国产有黄有色有爽视频| 成年人午夜在线观看视频| 亚洲精品aⅴ在线观看| 国产淫片久久久久久久久| 大话2 男鬼变身卡| 狂野欧美激情性bbbbbb| 国产黄频视频在线观看| 18禁动态无遮挡网站| 久久精品国产自在天天线| 中文字幕免费在线视频6| 精品国产三级普通话版| 97在线人人人人妻| 日韩三级伦理在线观看| 日韩一区二区视频免费看| 国产男女内射视频| 国产精品不卡视频一区二区| 久久99热这里只频精品6学生| 亚洲国产日韩一区二区| 国产免费一区二区三区四区乱码| 欧美三级亚洲精品| 亚洲欧美一区二区三区国产| 国产亚洲5aaaaa淫片| 成人鲁丝片一二三区免费| 亚洲精品久久久久久婷婷小说| 亚洲精品成人av观看孕妇| 国产极品天堂在线| 又爽又黄无遮挡网站| 国产一区二区在线观看日韩| 久久精品国产鲁丝片午夜精品| 七月丁香在线播放| 成人一区二区视频在线观看| 综合色av麻豆| 高清视频免费观看一区二区| 91精品伊人久久大香线蕉| av天堂中文字幕网| 大香蕉97超碰在线| 午夜日本视频在线| 晚上一个人看的免费电影| 少妇人妻 视频| 涩涩av久久男人的天堂| 少妇的逼好多水| 久久久欧美国产精品| 欧美性猛交╳xxx乱大交人| 不卡视频在线观看欧美| videossex国产| 欧美日韩一区二区视频在线观看视频在线 | 久久精品国产亚洲av天美| 成人国产麻豆网| 新久久久久国产一级毛片| 亚洲欧洲国产日韩| 女的被弄到高潮叫床怎么办| 欧美区成人在线视频| 国产成人freesex在线| 一级毛片 在线播放| 国产精品国产三级国产专区5o| 国产有黄有色有爽视频| av在线天堂中文字幕| 久久久a久久爽久久v久久| 国产有黄有色有爽视频| 偷拍熟女少妇极品色| 亚洲三级黄色毛片| 亚洲国产精品专区欧美| 国产精品久久久久久久久免| 日韩av不卡免费在线播放| 亚洲精品国产av成人精品| 精品酒店卫生间| 美女xxoo啪啪120秒动态图| 一级毛片aaaaaa免费看小| 亚洲av成人精品一二三区| 国产精品伦人一区二区| 久久久午夜欧美精品| 简卡轻食公司| 国产成人一区二区在线| 亚洲国产精品999| 在线观看三级黄色| 亚洲国产精品国产精品| 日韩欧美精品v在线| 高清欧美精品videossex| 中文资源天堂在线| 国产欧美另类精品又又久久亚洲欧美| 国产精品久久久久久精品古装| 亚洲欧美精品自产自拍| 六月丁香七月| 国产成人freesex在线| 亚洲最大成人av| 亚洲丝袜综合中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩三级伦理在线观看| 亚洲av不卡在线观看| 久久久久久九九精品二区国产| 高清欧美精品videossex| 国产高清不卡午夜福利| 性插视频无遮挡在线免费观看| 免费大片黄手机在线观看| 亚洲婷婷狠狠爱综合网| av网站免费在线观看视频| 欧美丝袜亚洲另类| 视频区图区小说| 成人美女网站在线观看视频| 丰满乱子伦码专区| 大香蕉久久网| 国产片特级美女逼逼视频| 日日啪夜夜撸| 蜜桃久久精品国产亚洲av| 亚洲av福利一区| 日韩成人伦理影院| av国产久精品久网站免费入址| 亚洲欧美一区二区三区黑人 | 在线免费十八禁| 欧美亚洲 丝袜 人妻 在线| 国产欧美日韩一区二区三区在线 | 人人妻人人爽人人添夜夜欢视频 | 成人国产av品久久久| 国产精品爽爽va在线观看网站| 亚洲,欧美,日韩| 日韩欧美精品免费久久| 亚洲av中文av极速乱| 中文字幕免费在线视频6| 亚洲自偷自拍三级| 欧美人与善性xxx| 国模一区二区三区四区视频| 亚洲精品亚洲一区二区| 国产精品蜜桃在线观看| 亚洲自偷自拍三级| 国产精品99久久久久久久久| 春色校园在线视频观看| 国产黄频视频在线观看| 国产亚洲av片在线观看秒播厂| 激情五月婷婷亚洲| 亚洲经典国产精华液单| 最近2019中文字幕mv第一页| av在线亚洲专区| 国国产精品蜜臀av免费| 成年版毛片免费区| 卡戴珊不雅视频在线播放| 国产白丝娇喘喷水9色精品| 久久精品综合一区二区三区| 亚洲av电影在线观看一区二区三区 | 别揉我奶头 嗯啊视频| 少妇丰满av| 精品99又大又爽又粗少妇毛片| 欧美bdsm另类| 18禁动态无遮挡网站| 少妇丰满av| 别揉我奶头 嗯啊视频| 精品一区二区免费观看| 涩涩av久久男人的天堂| 黑人高潮一二区| 欧美97在线视频| 一边亲一边摸免费视频| 精品久久久久久久人妻蜜臀av| 激情五月婷婷亚洲| 欧美老熟妇乱子伦牲交| 黄色怎么调成土黄色| 黄色配什么色好看| 91精品伊人久久大香线蕉| 亚洲av二区三区四区| 18禁裸乳无遮挡免费网站照片| 永久网站在线| 熟女人妻精品中文字幕| 久久久午夜欧美精品| 成年免费大片在线观看| 亚洲精品成人久久久久久| 极品教师在线视频| 搡女人真爽免费视频火全软件| 新久久久久国产一级毛片| 爱豆传媒免费全集在线观看| 秋霞伦理黄片| 久久影院123| 亚洲无线观看免费| 午夜福利视频1000在线观看| 欧美日韩视频精品一区| 日韩视频在线欧美| 一二三四中文在线观看免费高清| 国产永久视频网站| 亚洲av在线观看美女高潮| 亚洲最大成人av| 熟女电影av网| 美女视频免费永久观看网站| 少妇的逼好多水| h日本视频在线播放| 亚洲精品成人av观看孕妇| 日日撸夜夜添| 内射极品少妇av片p| 中文乱码字字幕精品一区二区三区| 久久精品国产自在天天线| 国产伦精品一区二区三区四那| 大片免费播放器 马上看| 午夜福利在线在线| tube8黄色片| 国产黄片视频在线免费观看| av网站免费在线观看视频| av线在线观看网站| 国产精品女同一区二区软件| 亚洲av日韩在线播放| 亚洲人成网站在线播| 综合色av麻豆| 亚洲电影在线观看av| 女人被狂操c到高潮| 久久人人爽人人爽人人片va| 国产黄片美女视频| 青春草国产在线视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲自偷自拍三级| 大香蕉久久网| 人妻系列 视频| 亚洲国产欧美在线一区| 色综合色国产| 最近2019中文字幕mv第一页| 七月丁香在线播放| 伊人久久精品亚洲午夜| 国产黄色视频一区二区在线观看| 搡女人真爽免费视频火全软件| 国产精品不卡视频一区二区| 男人和女人高潮做爰伦理| 日韩伦理黄色片| 日韩在线高清观看一区二区三区| 亚洲四区av| 亚洲欧美成人综合另类久久久| 熟女av电影| 成年人午夜在线观看视频| 久久99热这里只有精品18| 高清视频免费观看一区二区| 亚洲欧美中文字幕日韩二区| 免费av毛片视频| 伦精品一区二区三区| 麻豆成人午夜福利视频| 精品一区二区免费观看| 十八禁网站网址无遮挡 | 日本爱情动作片www.在线观看| 成人亚洲欧美一区二区av| 久久久久久久午夜电影| 又爽又黄无遮挡网站| 视频区图区小说| 尤物成人国产欧美一区二区三区| 我的女老师完整版在线观看| 中文字幕久久专区| 噜噜噜噜噜久久久久久91| 看免费成人av毛片| 婷婷色综合www| 婷婷色麻豆天堂久久|