• <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
    久久精品91蜜桃| 全区人妻精品视频| 99热精品在线国产| 蜜桃久久精品国产亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 99热6这里只有精品| 久久久久九九精品影院| 在线观看美女被高潮喷水网站| 麻豆国产97在线/欧美| 久久久久久久久久黄片| 欧美bdsm另类| 久久国产乱子免费精品| 国产精品无大码| av在线观看视频网站免费| 亚洲精品aⅴ在线观看| 欧美三级亚洲精品| 99视频精品全部免费 在线| 国产老妇女一区| 尾随美女入室| 三级毛片av免费| 国产人妻一区二区三区在| 亚洲五月天丁香| 成人午夜高清在线视频| 亚洲人与动物交配视频| 女人十人毛片免费观看3o分钟| 日本免费a在线| 最近视频中文字幕2019在线8| 变态另类丝袜制服| 久久人人爽人人爽人人片va| 亚洲一级一片aⅴ在线观看| 亚洲怡红院男人天堂| 久久精品久久久久久噜噜老黄 | 性插视频无遮挡在线免费观看| 99国产精品一区二区蜜桃av| 一级毛片久久久久久久久女| 亚洲精品国产成人久久av| 男女边吃奶边做爰视频| 男女国产视频网站| 欧美性感艳星| 免费看美女性在线毛片视频| 偷拍熟女少妇极品色| 久久综合国产亚洲精品| 一个人免费在线观看电影| 青青草视频在线视频观看| 91在线精品国自产拍蜜月| 搡老妇女老女人老熟妇| 日本av手机在线免费观看| 日韩成人av中文字幕在线观看| 日韩成人av中文字幕在线观看| 成人午夜高清在线视频| 色播亚洲综合网| 嫩草影院新地址| 91精品伊人久久大香线蕉| videossex国产| 国产精品熟女久久久久浪| 亚洲欧美一区二区三区国产| 国产精品久久电影中文字幕| 婷婷色av中文字幕| 欧美3d第一页| 亚洲国产色片| 麻豆乱淫一区二区| 免费看av在线观看网站| 久久精品国产99精品国产亚洲性色| 国产单亲对白刺激| 免费黄网站久久成人精品| 超碰97精品在线观看| 国产av一区在线观看免费| 午夜久久久久精精品| 国产又黄又爽又无遮挡在线| 久久久精品大字幕| 精品国内亚洲2022精品成人| 亚洲丝袜综合中文字幕| 欧美高清成人免费视频www| 97超碰精品成人国产| 岛国毛片在线播放| 亚洲av二区三区四区| 夫妻性生交免费视频一级片| 午夜福利在线观看免费完整高清在| 国产免费福利视频在线观看| 亚洲av日韩在线播放| 成人毛片a级毛片在线播放| 日韩一本色道免费dvd| 午夜福利高清视频| 久久精品人妻少妇| 久久久午夜欧美精品| 日本午夜av视频| 99热网站在线观看| 国产一区亚洲一区在线观看| 久久精品综合一区二区三区| 偷拍熟女少妇极品色| 91在线精品国自产拍蜜月| 男插女下体视频免费在线播放| 人妻少妇偷人精品九色| 男人狂女人下面高潮的视频| 在线播放国产精品三级| 国产伦精品一区二区三区四那| 两性午夜刺激爽爽歪歪视频在线观看| 黄色一级大片看看| 亚洲经典国产精华液单| 日本免费在线观看一区| 成人毛片a级毛片在线播放| 久久久久久久久久黄片| 成人特级av手机在线观看| av视频在线观看入口| 毛片一级片免费看久久久久| 国内精品宾馆在线| 狂野欧美白嫩少妇大欣赏| 国产精品久久电影中文字幕| 黄色欧美视频在线观看| 亚洲性久久影院| 欧美最新免费一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 夜夜爽夜夜爽视频| 尾随美女入室| 欧美日韩综合久久久久久| 日韩人妻高清精品专区| 日本色播在线视频| 又粗又硬又长又爽又黄的视频| 午夜亚洲福利在线播放| 亚洲精品国产av成人精品| 少妇熟女欧美另类| videos熟女内射| 赤兔流量卡办理| 国产精品一区二区在线观看99 | 中文字幕久久专区| 看免费成人av毛片| 亚洲电影在线观看av| 日本与韩国留学比较| 国内揄拍国产精品人妻在线| 亚洲激情五月婷婷啪啪| 久久99热这里只有精品18| 色综合站精品国产| 69av精品久久久久久| 水蜜桃什么品种好| 美女高潮的动态| 亚洲人成网站高清观看| 精品久久久久久久人妻蜜臀av| www.av在线官网国产| 精品一区二区三区人妻视频| 村上凉子中文字幕在线| 七月丁香在线播放| 中文字幕熟女人妻在线| 男女视频在线观看网站免费| 老司机影院成人| 国产精品一二三区在线看| 乱码一卡2卡4卡精品| 国产伦理片在线播放av一区| 国产av码专区亚洲av| 久久久午夜欧美精品| 色吧在线观看| 在线观看美女被高潮喷水网站| 日本欧美国产在线视频| 亚洲精品456在线播放app| 亚洲天堂国产精品一区在线| 亚洲精华国产精华液的使用体验| 国产黄片美女视频| 男女啪啪激烈高潮av片| 丰满人妻一区二区三区视频av| 欧美性感艳星| 久久久国产成人免费| 一卡2卡三卡四卡精品乱码亚洲| 啦啦啦韩国在线观看视频| 亚洲美女视频黄频| or卡值多少钱| 午夜福利视频1000在线观看| 特大巨黑吊av在线直播| 夜夜爽夜夜爽视频| 看免费成人av毛片| 22中文网久久字幕| 国产伦精品一区二区三区视频9| 一级毛片电影观看 | 亚洲国产成人一精品久久久| 国产色婷婷99| 久久久久国产网址| 日韩一本色道免费dvd| 99热网站在线观看| 国产黄色小视频在线观看| 国产人妻一区二区三区在| 99久久成人亚洲精品观看| 国产亚洲精品av在线| 亚洲中文字幕日韩| 亚洲成人中文字幕在线播放| 国产一区二区在线av高清观看| 女人十人毛片免费观看3o分钟| 简卡轻食公司| 99久久成人亚洲精品观看| 两性午夜刺激爽爽歪歪视频在线观看| 91午夜精品亚洲一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 久久久亚洲精品成人影院| www.色视频.com| 成人漫画全彩无遮挡| videossex国产| 美女高潮的动态| 一级毛片aaaaaa免费看小| 一个人看的www免费观看视频| 日韩一本色道免费dvd| 色吧在线观看| 一夜夜www| 热99re8久久精品国产| 亚洲无线观看免费| 日本-黄色视频高清免费观看| 亚洲av成人精品一二三区| 欧美性猛交╳xxx乱大交人| 男人的好看免费观看在线视频| 国产一级毛片在线| 亚洲丝袜综合中文字幕| 人人妻人人澡人人爽人人夜夜 | 秋霞伦理黄片| 99视频精品全部免费 在线| 少妇被粗大猛烈的视频| 久久久精品欧美日韩精品| 99久国产av精品国产电影| 欧美成人午夜免费资源| 中国美白少妇内射xxxbb| 人妻夜夜爽99麻豆av| 99久久精品热视频| 久久精品久久精品一区二区三区| 国产精品美女特级片免费视频播放器| 韩国av在线不卡| 大香蕉97超碰在线| 国产精品无大码| 熟女人妻精品中文字幕| 欧美人与善性xxx| 男人的好看免费观看在线视频| 成人av在线播放网站| 国产一级毛片在线| 欧美变态另类bdsm刘玥| 精品一区二区三区人妻视频| 久久久久久久久久黄片| 丰满少妇做爰视频| 尤物成人国产欧美一区二区三区| 天天一区二区日本电影三级| 黄片无遮挡物在线观看| 免费看美女性在线毛片视频| 一级黄片播放器| 看黄色毛片网站| 91久久精品电影网| 亚洲av免费在线观看| 国产亚洲精品久久久com| 欧美高清性xxxxhd video| 特级一级黄色大片| 亚洲乱码一区二区免费版| 一级毛片电影观看 | 免费人成在线观看视频色| 九九在线视频观看精品| 久久99热这里只有精品18| 嫩草影院入口| 免费看av在线观看网站| 国产伦精品一区二区三区四那| 精品一区二区三区人妻视频| 亚洲国产欧美在线一区| 18禁裸乳无遮挡免费网站照片| 高清在线视频一区二区三区 | 国产伦精品一区二区三区视频9| 丝袜美腿在线中文| 午夜福利成人在线免费观看| 国产免费男女视频| 国产亚洲5aaaaa淫片| 有码 亚洲区| 国产麻豆成人av免费视频| 联通29元200g的流量卡| 成人毛片60女人毛片免费| 天堂√8在线中文| 精品一区二区三区视频在线| 亚洲av成人精品一区久久| 高清午夜精品一区二区三区| 国产伦精品一区二区三区视频9| 六月丁香七月| 欧美成人午夜免费资源| 一个人看的www免费观看视频| 97超视频在线观看视频| 国产淫片久久久久久久久| 久久久久久久亚洲中文字幕| 一级毛片久久久久久久久女| 麻豆乱淫一区二区| 免费av观看视频| 男的添女的下面高潮视频| 国产亚洲午夜精品一区二区久久 | 国内精品美女久久久久久| 国产av不卡久久| 久久久久网色| 大又大粗又爽又黄少妇毛片口| 午夜福利在线观看免费完整高清在| 日韩三级伦理在线观看| 99久久精品国产国产毛片| 十八禁国产超污无遮挡网站| 熟女电影av网| 美女黄网站色视频| 99久国产av精品国产电影| 免费人成在线观看视频色| av线在线观看网站| 夫妻性生交免费视频一级片| 99久久精品热视频| 亚洲国产最新在线播放| 高清在线视频一区二区三区 | 久久久成人免费电影| 一边摸一边抽搐一进一小说| 日本一本二区三区精品| 成年女人看的毛片在线观看| 长腿黑丝高跟| 最近的中文字幕免费完整| 日韩欧美在线乱码| 婷婷色麻豆天堂久久 | 赤兔流量卡办理| av播播在线观看一区| 亚洲av熟女| 亚洲精品一区蜜桃| 亚洲国产高清在线一区二区三| 国产私拍福利视频在线观看| 国产老妇伦熟女老妇高清| 国产伦精品一区二区三区四那| 午夜日本视频在线| 免费看美女性在线毛片视频| 欧美极品一区二区三区四区| 99久久九九国产精品国产免费| 婷婷六月久久综合丁香| 好男人视频免费观看在线| 久久精品国产鲁丝片午夜精品| 在线天堂最新版资源| 国产国拍精品亚洲av在线观看| 精品欧美国产一区二区三| 欧美97在线视频| 国内少妇人妻偷人精品xxx网站| 久久久精品欧美日韩精品| 国产v大片淫在线免费观看| 免费av不卡在线播放| 亚洲av男天堂| 久久久国产成人免费| 国产亚洲精品av在线| 亚洲精品成人久久久久久| 一区二区三区免费毛片| 国产av一区在线观看免费| 国内精品美女久久久久久| 精品人妻偷拍中文字幕| 最近视频中文字幕2019在线8| 成人鲁丝片一二三区免费| .国产精品久久| 国产私拍福利视频在线观看| 久久久久网色| 麻豆一二三区av精品| 91av网一区二区| 日韩高清综合在线| 久久精品熟女亚洲av麻豆精品 | 成人高潮视频无遮挡免费网站| 国产成人91sexporn| 精品一区二区三区视频在线| 国产伦在线观看视频一区| 精品久久久久久久久亚洲| 丰满少妇做爰视频| 久久精品熟女亚洲av麻豆精品 | 少妇熟女aⅴ在线视频| 亚洲怡红院男人天堂| 欧美色视频一区免费| 欧美日本视频| 男人舔女人下体高潮全视频| 美女cb高潮喷水在线观看| 99热精品在线国产| 亚洲乱码一区二区免费版| 精品久久久久久久末码| 国产亚洲精品av在线| 91aial.com中文字幕在线观看| 日本av手机在线免费观看| 国产久久久一区二区三区| 午夜福利高清视频| 精品免费久久久久久久清纯| 亚洲精品一区蜜桃| 午夜福利网站1000一区二区三区| 超碰av人人做人人爽久久| 中文在线观看免费www的网站| 人妻夜夜爽99麻豆av| 一个人免费在线观看电影| 国产综合懂色| 国产乱人视频| 国产精华一区二区三区| 成人午夜精彩视频在线观看| 91av网一区二区| 欧美色视频一区免费| 亚洲最大成人中文| av又黄又爽大尺度在线免费看 | 欧美变态另类bdsm刘玥| 国产亚洲最大av| 色尼玛亚洲综合影院| 两个人的视频大全免费| 最近中文字幕高清免费大全6| 啦啦啦观看免费观看视频高清| 久久国产乱子免费精品| 久久国内精品自在自线图片| 18+在线观看网站| 搡女人真爽免费视频火全软件| 好男人在线观看高清免费视频| 国产成人午夜福利电影在线观看| 成人国产麻豆网| 国产美女午夜福利| 黄色配什么色好看| 亚洲内射少妇av| 超碰97精品在线观看| 夫妻性生交免费视频一级片| 久久久久久久午夜电影| 18禁动态无遮挡网站| 日韩人妻高清精品专区| 亚洲高清免费不卡视频| 精品国内亚洲2022精品成人| 国产不卡一卡二| 六月丁香七月| 国产精品电影一区二区三区| 日本午夜av视频| 国产老妇伦熟女老妇高清| 尾随美女入室| 亚洲av免费在线观看| 久久精品国产亚洲网站| 国产精品久久久久久久电影| 综合色丁香网| 午夜精品一区二区三区免费看| 久久国产乱子免费精品| 亚洲av免费高清在线观看| 黄色欧美视频在线观看| 26uuu在线亚洲综合色| 51国产日韩欧美| 人体艺术视频欧美日本| 超碰av人人做人人爽久久| 亚洲人成网站高清观看| 亚洲乱码一区二区免费版| 久久久久九九精品影院| 亚洲性久久影院| 最近最新中文字幕免费大全7| 一边亲一边摸免费视频| 建设人人有责人人尽责人人享有的 | 18禁裸乳无遮挡免费网站照片| 亚洲精华国产精华液的使用体验| 草草在线视频免费看| 午夜激情福利司机影院| 一级黄色大片毛片| 午夜亚洲福利在线播放| 国产亚洲精品av在线| 免费av毛片视频| 青春草国产在线视频| 久久精品91蜜桃| 最近中文字幕高清免费大全6| 一级毛片aaaaaa免费看小| 男女下面进入的视频免费午夜| 三级毛片av免费| 永久网站在线| 亚洲国产精品成人综合色| 亚洲人成网站在线观看播放| 最新中文字幕久久久久| 国产美女午夜福利| 久久久久国产网址| 全区人妻精品视频| 日本与韩国留学比较| 午夜免费激情av| 91狼人影院| 天堂√8在线中文| 欧美成人精品欧美一级黄| 麻豆久久精品国产亚洲av| 午夜福利在线观看吧| 久久久久久久久久成人| 精品久久久久久久久av| 国产av一区在线观看免费| 精品人妻视频免费看| 亚洲国产色片| 国产精品一区二区三区四区免费观看| 丰满乱子伦码专区| 国产色爽女视频免费观看| 97人妻精品一区二区三区麻豆| 国产成人aa在线观看| 欧美三级亚洲精品| 一二三四中文在线观看免费高清| 久久久久久久久中文| 亚洲国产欧美在线一区| 亚洲国产精品久久男人天堂| 日本与韩国留学比较| 亚洲欧美一区二区三区国产| 亚洲欧美精品自产自拍| 一区二区三区免费毛片| 国产精品乱码一区二三区的特点| 精品久久久久久久久av| 精品99又大又爽又粗少妇毛片| 成人三级黄色视频| 狂野欧美白嫩少妇大欣赏| 欧美三级亚洲精品| 男女视频在线观看网站免费| 麻豆成人av视频| 欧美精品国产亚洲| 久久精品国产亚洲av天美| av播播在线观看一区| 免费黄色在线免费观看| 精品国产露脸久久av麻豆 | 国产 一区精品| 男插女下体视频免费在线播放| 国产综合懂色| 国产精品福利在线免费观看| 国产麻豆成人av免费视频| 国产精品一区二区性色av| 一本久久精品| 国语自产精品视频在线第100页| 日本av手机在线免费观看| 欧美日本视频| 成人特级av手机在线观看| 啦啦啦韩国在线观看视频| 国产精品人妻久久久影院| 最近最新中文字幕大全电影3| 亚洲欧美成人综合另类久久久 | 在线a可以看的网站| 一二三四中文在线观看免费高清| 欧美人与善性xxx| 欧美性猛交╳xxx乱大交人| 看黄色毛片网站| 有码 亚洲区| 免费观看精品视频网站| 最近最新中文字幕免费大全7| 99久久成人亚洲精品观看| 搡女人真爽免费视频火全软件| 如何舔出高潮| 国产 一区精品| 男人狂女人下面高潮的视频| 亚洲欧美精品专区久久| 国产精品国产三级国产专区5o | 在现免费观看毛片| 老司机影院毛片| 午夜老司机福利剧场| 亚洲成色77777| 69av精品久久久久久| 少妇熟女aⅴ在线视频| 观看美女的网站| 中国美白少妇内射xxxbb| 国产高清国产精品国产三级 | 国产精品一区二区三区四区免费观看| 久久99热这里只频精品6学生 | 精品一区二区三区人妻视频| 亚洲va在线va天堂va国产| av在线天堂中文字幕| 成人欧美大片| 精品久久国产蜜桃| 午夜精品一区二区三区免费看| 色吧在线观看| 国产一区有黄有色的免费视频 | 中文天堂在线官网| 午夜爱爱视频在线播放| 国产高潮美女av| 99在线视频只有这里精品首页| 三级国产精品欧美在线观看| 毛片女人毛片| 欧美日韩一区二区视频在线观看视频在线 | 国产精品嫩草影院av在线观看| 成年女人看的毛片在线观看| 禁无遮挡网站| 国产精品久久久久久久电影| 久久午夜福利片| 熟女电影av网| 久久久久九九精品影院| 亚洲中文字幕日韩| 最近手机中文字幕大全| 国产老妇女一区| 久久欧美精品欧美久久欧美| 在线免费观看不下载黄p国产| ponron亚洲| 国产成人aa在线观看| 久久欧美精品欧美久久欧美| 久久精品国产亚洲网站| 久久欧美精品欧美久久欧美| 国产中年淑女户外野战色| 亚洲精品亚洲一区二区| 建设人人有责人人尽责人人享有的 | 深爱激情五月婷婷| 午夜免费男女啪啪视频观看| 久久精品综合一区二区三区| 午夜视频国产福利| 丰满乱子伦码专区| 特级一级黄色大片| 69av精品久久久久久| 国语自产精品视频在线第100页| 美女内射精品一级片tv| 插逼视频在线观看| 性色avwww在线观看| 国产av码专区亚洲av| 男人和女人高潮做爰伦理| 国产免费一级a男人的天堂| 亚洲成人久久爱视频| 蜜臀久久99精品久久宅男| 成人二区视频| 亚洲欧美成人综合另类久久久 | 国内揄拍国产精品人妻在线| 高清在线视频一区二区三区 | 欧美成人一区二区免费高清观看| 99视频精品全部免费 在线| 1024手机看黄色片| 欧美一区二区精品小视频在线| 久99久视频精品免费| 在线观看av片永久免费下载| 秋霞在线观看毛片| 日韩一本色道免费dvd| 麻豆精品久久久久久蜜桃| 欧美日本亚洲视频在线播放| 欧美xxxx性猛交bbbb| 最近中文字幕2019免费版| 性插视频无遮挡在线免费观看| 在线天堂最新版资源| 美女内射精品一级片tv| 天天一区二区日本电影三级| 视频中文字幕在线观看| 黄色日韩在线| 2021少妇久久久久久久久久久| 国产免费一级a男人的天堂| 99热全是精品| 免费观看a级毛片全部| 直男gayav资源| 国产精品人妻久久久影院| 午夜免费男女啪啪视频观看| 观看美女的网站| 久久久久久久久久久免费av| 日韩欧美精品免费久久| 久久久久久九九精品二区国产| 久久精品国产亚洲网站| 3wmmmm亚洲av在线观看|