• <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
    国产v大片淫在线免费观看| 日本一本二区三区精品| 美女高潮的动态| 51国产日韩欧美| 国产精品久久久久久精品电影小说 | 99久久精品热视频| 欧美日韩精品成人综合77777| 99国产极品粉嫩在线观看| 国内精品美女久久久久久| 在线观看免费视频日本深夜| 成人亚洲欧美一区二区av| 免费看光身美女| 成人高潮视频无遮挡免费网站| 亚洲国产精品合色在线| 老女人水多毛片| 亚洲国产精品久久男人天堂| 欧洲精品卡2卡3卡4卡5卡区| 国产一级毛片在线| 男人狂女人下面高潮的视频| 国产一级毛片在线| 99久国产av精品| 又爽又黄无遮挡网站| 中出人妻视频一区二区| 晚上一个人看的免费电影| 中文字幕制服av| 国产精品伦人一区二区| 99久国产av精品| 内地一区二区视频在线| 99久国产av精品| 亚洲久久久久久中文字幕| 亚洲图色成人| 老司机福利观看| 亚洲第一电影网av| 少妇熟女欧美另类| 麻豆成人午夜福利视频| 欧美xxxx性猛交bbbb| 亚洲最大成人手机在线| 亚洲欧美中文字幕日韩二区| 国产片特级美女逼逼视频| 免费观看的影片在线观看| 男插女下体视频免费在线播放| 午夜老司机福利剧场| 嫩草影院入口| 毛片女人毛片| 99在线视频只有这里精品首页| av女优亚洲男人天堂| 最近2019中文字幕mv第一页| 99热全是精品| 最近2019中文字幕mv第一页| 美女xxoo啪啪120秒动态图| 亚洲精品国产色婷婷电影| 激情五月婷婷亚洲| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产欧美在线一区| 国产男人的电影天堂91| 桃花免费在线播放| 最新中文字幕久久久久| 国产一区二区在线观看日韩| 精品国产一区二区久久| 精品久久久久久电影网| 亚州av有码| 成人国语在线视频| xxx大片免费视频| 国产精品欧美亚洲77777| 青春草亚洲视频在线观看| 亚洲精品久久午夜乱码| 插阴视频在线观看视频| 国产老妇伦熟女老妇高清| av女优亚洲男人天堂| 麻豆乱淫一区二区| 丝袜在线中文字幕| 18禁观看日本| 老女人水多毛片| 国产午夜精品一二区理论片| 久久亚洲国产成人精品v| 亚洲国产精品999| 高清视频免费观看一区二区| 欧美人与善性xxx| 男女免费视频国产| 亚洲熟女精品中文字幕| 国产精品秋霞免费鲁丝片| 亚洲欧美成人综合另类久久久| 亚洲av成人精品一二三区| 久久久久久伊人网av| 国产男女内射视频| 国产高清有码在线观看视频| 国产免费一区二区三区四区乱码| 99久国产av精品国产电影| 麻豆乱淫一区二区| 亚洲人成网站在线播| 日日啪夜夜爽| 男女无遮挡免费网站观看| 少妇被粗大的猛进出69影院 | a 毛片基地| 成人二区视频| 乱码一卡2卡4卡精品| 少妇被粗大的猛进出69影院 | 能在线免费看毛片的网站| 老司机亚洲免费影院| 久久人人爽人人片av| 美女视频免费永久观看网站| 国产在线一区二区三区精| 免费少妇av软件| 亚洲无线观看免费| 人人妻人人澡人人爽人人夜夜| 午夜福利在线观看免费完整高清在| 22中文网久久字幕| 成人国产麻豆网| 精品久久久久久久久av| 亚洲欧美日韩卡通动漫| 国语对白做爰xxxⅹ性视频网站| 简卡轻食公司| 九九爱精品视频在线观看| 成年女人在线观看亚洲视频| 亚洲国产日韩一区二区| 制服丝袜香蕉在线| 国产黄频视频在线观看| 一本色道久久久久久精品综合| 国产高清国产精品国产三级| 美女福利国产在线| 精品一区二区三区视频在线| 久热久热在线精品观看| 大又大粗又爽又黄少妇毛片口| 亚洲性久久影院| 久久人妻熟女aⅴ| 免费高清在线观看视频在线观看| 少妇人妻 视频| 国产视频首页在线观看| 一本大道久久a久久精品| 韩国高清视频一区二区三区| 内地一区二区视频在线| 国产色婷婷99| 久久久久国产网址| 欧美精品一区二区大全| 亚洲欧洲国产日韩| 欧美国产精品一级二级三级| videossex国产| 18禁在线无遮挡免费观看视频| 热99国产精品久久久久久7| 999精品在线视频| 青春草视频在线免费观看| 在线播放无遮挡| 欧美+日韩+精品| 久久精品国产亚洲av涩爱| 亚洲欧美一区二区三区黑人 | 黄色视频在线播放观看不卡| 最近的中文字幕免费完整| av线在线观看网站| 亚洲欧美成人精品一区二区| .国产精品久久| 蜜臀久久99精品久久宅男| 久久精品久久久久久噜噜老黄| 天堂中文最新版在线下载| 少妇高潮的动态图| 91久久精品电影网| 日韩制服骚丝袜av| 国产黄片视频在线免费观看| 国产黄色视频一区二区在线观看| 午夜精品国产一区二区电影| 亚洲丝袜综合中文字幕| 性高湖久久久久久久久免费观看| 亚洲激情五月婷婷啪啪| 最近最新中文字幕免费大全7| 国产亚洲最大av| 精品久久久精品久久久| 欧美亚洲日本最大视频资源| 国产精品麻豆人妻色哟哟久久| av视频免费观看在线观看| 亚洲av二区三区四区| 这个男人来自地球电影免费观看 | 视频区图区小说| 久久久久人妻精品一区果冻| 夜夜骑夜夜射夜夜干| 久热久热在线精品观看| 午夜免费男女啪啪视频观看| 免费观看在线日韩| 免费观看a级毛片全部| 夫妻午夜视频| 在现免费观看毛片| 亚洲不卡免费看| 久久久久久久久久成人| 国产免费福利视频在线观看| 婷婷色综合大香蕉| av一本久久久久| 亚洲国产精品成人久久小说| 精品亚洲成国产av| 伦理电影大哥的女人| 成人手机av| 有码 亚洲区| 色婷婷久久久亚洲欧美| 黑人巨大精品欧美一区二区蜜桃 | 91精品一卡2卡3卡4卡| 男女免费视频国产| 亚洲怡红院男人天堂| 人人妻人人添人人爽欧美一区卜| 国产老妇伦熟女老妇高清| 久久精品国产a三级三级三级| 久久精品国产亚洲av涩爱| 91午夜精品亚洲一区二区三区| 最黄视频免费看| 综合色丁香网| 国产亚洲一区二区精品| 少妇精品久久久久久久| 女人久久www免费人成看片| 2022亚洲国产成人精品| 高清av免费在线| 国产在视频线精品| 亚洲国产成人一精品久久久| 亚洲无线观看免费| 91久久精品电影网| 五月玫瑰六月丁香| 中文字幕精品免费在线观看视频 | 久久人人爽人人片av| 国产日韩欧美在线精品| 亚洲人成77777在线视频| 亚洲色图 男人天堂 中文字幕 | 精品人妻偷拍中文字幕| 国产精品国产三级国产av玫瑰| av在线播放精品| 午夜激情福利司机影院| 精品一区二区三卡| 精品国产露脸久久av麻豆| 国产国语露脸激情在线看| 视频中文字幕在线观看| 午夜激情福利司机影院| 久热这里只有精品99| 欧美精品一区二区大全| 国产一区亚洲一区在线观看| 国模一区二区三区四区视频| 精品99又大又爽又粗少妇毛片| 99re6热这里在线精品视频| 亚洲激情五月婷婷啪啪| 自拍欧美九色日韩亚洲蝌蚪91| 99久国产av精品国产电影| 亚洲人成网站在线播| 国产一区二区三区综合在线观看 | 中文欧美无线码| av一本久久久久| 亚洲国产精品成人久久小说| 夜夜爽夜夜爽视频| 老司机影院毛片| 日韩中字成人| 青春草亚洲视频在线观看| 下体分泌物呈黄色| 中文字幕久久专区| 老司机影院毛片| 国产成人免费观看mmmm| 制服人妻中文乱码| 妹子高潮喷水视频| 国产日韩欧美在线精品| 免费少妇av软件| 成人亚洲精品一区在线观看| 99国产综合亚洲精品| 精品国产乱码久久久久久小说| 日本爱情动作片www.在线观看| 天堂俺去俺来也www色官网| 亚洲精华国产精华液的使用体验| 美女大奶头黄色视频| av一本久久久久| 久久人人爽人人片av| 热99久久久久精品小说推荐| 母亲3免费完整高清在线观看 | 婷婷色综合www| 老女人水多毛片| 三级国产精品片| 18+在线观看网站| 亚洲av中文av极速乱| 97在线视频观看| 亚洲精品第二区| 少妇被粗大的猛进出69影院 | 亚洲欧美成人精品一区二区| 蜜臀久久99精品久久宅男| 亚洲欧洲精品一区二区精品久久久 | 免费高清在线观看视频在线观看| 制服人妻中文乱码| 韩国高清视频一区二区三区| 久久亚洲国产成人精品v| 国产极品天堂在线| 狂野欧美激情性xxxx在线观看| 日韩一区二区三区影片| 99re6热这里在线精品视频| 少妇丰满av| 大片免费播放器 马上看| 欧美成人精品欧美一级黄| 免费看光身美女| 婷婷色综合大香蕉| 亚洲精品av麻豆狂野| 丝袜在线中文字幕| 看免费成人av毛片| 五月伊人婷婷丁香| 免费观看性生交大片5| 婷婷色av中文字幕| 亚洲欧美成人精品一区二区| 日产精品乱码卡一卡2卡三| 中文字幕最新亚洲高清| 一级黄片播放器| 欧美老熟妇乱子伦牲交| 亚洲精品久久久久久婷婷小说| 在线观看美女被高潮喷水网站| 国产黄色视频一区二区在线观看| 午夜免费观看性视频| 婷婷色麻豆天堂久久| 久久人人爽av亚洲精品天堂| 飞空精品影院首页| 黄色配什么色好看| 亚洲av成人精品一二三区| 春色校园在线视频观看| videos熟女内射| 欧美三级亚洲精品| 亚洲国产欧美在线一区| 亚洲精品日本国产第一区| 日韩电影二区| 国产精品偷伦视频观看了| 久久影院123| 久久精品久久久久久久性| 成年人免费黄色播放视频| 亚洲天堂av无毛| 国产成人a∨麻豆精品| 国产在线一区二区三区精| 99精国产麻豆久久婷婷| 色婷婷av一区二区三区视频| 久久韩国三级中文字幕| 乱人伦中国视频| 高清午夜精品一区二区三区| 国产成人免费观看mmmm| 亚洲四区av| 国产淫语在线视频| 亚洲激情五月婷婷啪啪| 亚洲国产精品一区三区| 在线亚洲精品国产二区图片欧美 | 日韩成人av中文字幕在线观看| 热re99久久精品国产66热6| 伦理电影大哥的女人| 国产欧美另类精品又又久久亚洲欧美| 女性被躁到高潮视频| 欧美激情 高清一区二区三区| 少妇被粗大猛烈的视频| 热re99久久国产66热| 中文字幕人妻熟人妻熟丝袜美| 色94色欧美一区二区| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区免费毛片| av在线app专区| 制服丝袜香蕉在线| 男女边吃奶边做爰视频| 热re99久久精品国产66热6| 黄色配什么色好看| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久精品久久久| 欧美精品高潮呻吟av久久| 母亲3免费完整高清在线观看 | 香蕉精品网在线| 亚洲精华国产精华液的使用体验| 国产成人精品久久久久久| 国产亚洲最大av| 飞空精品影院首页| 国产有黄有色有爽视频| 国产精品免费大片| 国产日韩欧美视频二区| 免费大片黄手机在线观看| 日日撸夜夜添| 国产老妇伦熟女老妇高清| 少妇被粗大的猛进出69影院 | 99re6热这里在线精品视频| 美女主播在线视频| 王馨瑶露胸无遮挡在线观看| 在线观看一区二区三区激情| 亚洲精品视频女| 一区二区三区四区激情视频| 久久国产精品男人的天堂亚洲 | 久久99精品国语久久久| 精品少妇内射三级| 亚洲欧美一区二区三区国产| 自线自在国产av| 久久热精品热| 9色porny在线观看| 国产亚洲av片在线观看秒播厂| 色5月婷婷丁香| 日韩一本色道免费dvd| 丰满乱子伦码专区| videossex国产| 日韩精品免费视频一区二区三区 | 伦精品一区二区三区| 国产成人免费观看mmmm| 久久久久国产精品人妻一区二区| 精品少妇黑人巨大在线播放| 久久久久久伊人网av| 在线 av 中文字幕| a级片在线免费高清观看视频| 欧美日韩视频精品一区| 黄色配什么色好看| 777米奇影视久久| 国产精品国产三级国产专区5o| 性高湖久久久久久久久免费观看| 777米奇影视久久| 妹子高潮喷水视频| 国产av精品麻豆| 成人影院久久| 中文字幕免费在线视频6| 天天操日日干夜夜撸| 啦啦啦中文免费视频观看日本| 丝袜喷水一区| 91午夜精品亚洲一区二区三区| 在线观看美女被高潮喷水网站| 久久精品国产a三级三级三级| 美女大奶头黄色视频| 午夜免费鲁丝| 又黄又爽又刺激的免费视频.| 22中文网久久字幕| 最近中文字幕高清免费大全6| 免费黄频网站在线观看国产| a级毛片在线看网站| 最近中文字幕高清免费大全6| 亚洲精品一区蜜桃| 成人毛片60女人毛片免费| 精品一区在线观看国产| 91精品国产九色| 插逼视频在线观看| 男女边吃奶边做爰视频| 婷婷色av中文字幕| 大陆偷拍与自拍| 国产国拍精品亚洲av在线观看| 国产视频内射| 亚洲经典国产精华液单| 日本vs欧美在线观看视频| 免费观看的影片在线观看| 成人午夜精彩视频在线观看| 国产黄频视频在线观看| 毛片一级片免费看久久久久| 一级爰片在线观看| 国产黄色免费在线视频| 91久久精品国产一区二区成人| 一级爰片在线观看| freevideosex欧美| 国产av国产精品国产| 免费黄频网站在线观看国产| 一级,二级,三级黄色视频| 免费观看性生交大片5| 新久久久久国产一级毛片| 我的老师免费观看完整版| 精品一区在线观看国产| 熟妇人妻不卡中文字幕| 国产成人精品无人区| 日日爽夜夜爽网站| 精品亚洲成国产av| 天天躁夜夜躁狠狠久久av| 91精品国产国语对白视频| 免费黄频网站在线观看国产| 免费高清在线观看视频在线观看| 男男h啪啪无遮挡| 日韩电影二区| av免费观看日本| 亚洲av免费高清在线观看| 亚洲性久久影院| 99热网站在线观看| 亚洲精品第二区| 欧美日韩综合久久久久久| 久久久精品94久久精品| 狠狠精品人妻久久久久久综合| 天堂中文最新版在线下载| 天美传媒精品一区二区| 国产日韩欧美亚洲二区| 久久久久久久久久人人人人人人| 欧美另类一区| 亚洲精品一区蜜桃| 午夜日本视频在线| 日日摸夜夜添夜夜爱| 美女中出高潮动态图| 亚洲欧美成人综合另类久久久| 99久国产av精品国产电影| 国产探花极品一区二区| 精品久久蜜臀av无| 欧美成人精品欧美一级黄| 亚洲欧美一区二区三区黑人 | 日韩,欧美,国产一区二区三区| 国产国拍精品亚洲av在线观看| 亚洲人成网站在线播| 久久99蜜桃精品久久| 国产精品久久久久成人av| 日韩熟女老妇一区二区性免费视频| 午夜久久久在线观看| 亚洲少妇的诱惑av| 欧美日韩在线观看h| a 毛片基地| 全区人妻精品视频| 国产成人免费无遮挡视频| 亚洲欧美成人精品一区二区| 日韩大片免费观看网站| 国产无遮挡羞羞视频在线观看| 欧美精品一区二区大全| 丝袜美足系列| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久人人爽人人爽人人片va| 特大巨黑吊av在线直播| 男女边摸边吃奶| 99久久综合免费| 久久精品国产亚洲av涩爱| √禁漫天堂资源中文www| 中文字幕人妻熟人妻熟丝袜美| 日韩中字成人| 国产一级毛片在线| 日韩av在线免费看完整版不卡| 亚洲欧美一区二区三区国产| 人妻夜夜爽99麻豆av| 人妻人人澡人人爽人人| 看免费成人av毛片| 18禁在线播放成人免费| 精品视频人人做人人爽| 桃花免费在线播放| 欧美 日韩 精品 国产| 最新的欧美精品一区二区| 99久国产av精品国产电影| 国产精品一区二区三区四区免费观看| 国国产精品蜜臀av免费| 在线观看人妻少妇| 日韩伦理黄色片| 人妻系列 视频| 人成视频在线观看免费观看| 一区二区三区乱码不卡18| 国产精品一区二区三区四区免费观看| 国产视频首页在线观看| 精品一品国产午夜福利视频| 亚洲精品乱久久久久久| 久久久久久久久久久丰满| 精品久久久久久久久亚洲| 99热这里只有是精品在线观看| 午夜福利网站1000一区二区三区| 精品卡一卡二卡四卡免费| 国产一区二区三区综合在线观看 | 久久久久久久大尺度免费视频| 久久99精品国语久久久| 精品一区二区三区视频在线| 亚洲av日韩在线播放| 亚洲av免费高清在线观看| 亚洲第一区二区三区不卡| 卡戴珊不雅视频在线播放| 欧美人与善性xxx| 中文字幕最新亚洲高清| 99久久精品一区二区三区| 一级片'在线观看视频| videos熟女内射| 一区二区av电影网| 中文字幕亚洲精品专区| 成人黄色视频免费在线看| 韩国av在线不卡| 汤姆久久久久久久影院中文字幕| 美女国产高潮福利片在线看| 狠狠婷婷综合久久久久久88av| 狂野欧美白嫩少妇大欣赏| 91成人精品电影| 亚洲av.av天堂| 精品人妻一区二区三区麻豆| 国产精品秋霞免费鲁丝片| 久久久久久人妻| 免费播放大片免费观看视频在线观看| 纯流量卡能插随身wifi吗| 精品久久久久久久久亚洲| 女的被弄到高潮叫床怎么办| 婷婷色综合大香蕉| 国产精品久久久久久久久免| 欧美日韩视频精品一区| 婷婷色av中文字幕| 五月伊人婷婷丁香| 三上悠亚av全集在线观看| 最黄视频免费看| 亚洲美女搞黄在线观看| 亚洲伊人久久精品综合| 中文字幕人妻熟人妻熟丝袜美| 国产成人免费观看mmmm| 丝袜喷水一区| 精品午夜福利在线看| 亚洲第一av免费看| 精品午夜福利在线看| 久久久久久久亚洲中文字幕| 18禁裸乳无遮挡动漫免费视频| 国产日韩欧美亚洲二区| 成人国产av品久久久| 亚洲熟女精品中文字幕| 另类精品久久| av免费观看日本| 飞空精品影院首页| 精品国产乱码久久久久久小说| 如何舔出高潮| 26uuu在线亚洲综合色| 熟妇人妻不卡中文字幕| 精品国产一区二区久久| 97在线人人人人妻| 日韩大片免费观看网站| 亚洲av成人精品一区久久| 高清av免费在线| 两个人免费观看高清视频| 久久影院123| 不卡视频在线观看欧美| 久久国产亚洲av麻豆专区| 国产av精品麻豆| 五月天丁香电影| 亚洲精品视频女| 999精品在线视频| 草草在线视频免费看| 乱码一卡2卡4卡精品| 亚洲av福利一区| 成人午夜精彩视频在线观看| 亚洲av电影在线观看一区二区三区| 亚洲第一av免费看| 男女边吃奶边做爰视频| a级片在线免费高清观看视频| 成人漫画全彩无遮挡| 伦精品一区二区三区| 亚洲欧美日韩卡通动漫| 波野结衣二区三区在线| av不卡在线播放| 亚洲一区二区三区欧美精品| 在线观看www视频免费| 亚洲天堂av无毛| 内地一区二区视频在线|