發(fā)布時(shí)間：2021-04-02 14:18

　　反應(yīng)堆屏蔽計(jì)算中經(jīng)常出現(xiàn)厚屏蔽、小探測(cè)器問(wèn)題,常規(guī)蒙特卡羅方法難以有效解決。基于自動(dòng)重要抽樣（AIS）方法,本文提出了小探測(cè)器自動(dòng)重要抽樣（SDAIS）方法,并針對(duì)小探測(cè)器問(wèn)題,優(yōu)化了AIS方法的虛粒子賭分裂算法。該方法在自主開(kāi)發(fā)的蒙特卡羅屏蔽程序MCShield上進(jìn)行了實(shí)現(xiàn)。使用NUREG/CR-6115 PWR基準(zhǔn)題驗(yàn)證該方法的正確性和計(jì)算效率。結(jié)果表明,SDAIS方法可有效地解決厚屏蔽小探測(cè)器問(wèn)題,相比AIS方法及傳統(tǒng)的重要性方法,計(jì)算效率提升1～2個(gè)量級(jí)。 

【文章來(lái)源】：原子能科學(xué)技術(shù). 2020,54(07)北大核心EICSCD

【文章頁(yè)數(shù)】：7 頁(yè)

【部分圖文】：

AIS方法流程圖

DXTRAN球上的虛粒子通過(guò)抽樣產(chǎn)生,并被賦予相應(yīng)的權(quán)重,以保證無(wú)偏性。圖3為DXTRAN球上的虛粒子抽樣方法。設(shè)粒子在P1發(fā)生碰撞,DXTRAN球心為P0,DXTRAN外球半徑為RO,內(nèi)球半徑為RI,P1P0長(zhǎng)度為L(zhǎng)。DXTRAN外球、內(nèi)球在P1P0方向上極角最大值分別為θO、θI,其余弦值分別為ηO、ηI。虛粒子位置PS在DXTRAN外球上,具體坐標(biāo)通過(guò)抽樣以P1P0方向?yàn)榛鶞?zhǔn)的極角和方位角確定。極角η余弦值按式(2)的概率密度函數(shù)p(η)抽樣。圖3 DXTRAN球上虛粒子的抽樣方法

DXTRAN球上虛粒子的抽樣方法

【參考文獻(xiàn)】：
期刊論文
[1]Improved algorithms and coupled neutron-photon transport for auto-importance sampling method[J]. 王鑫,李君利,武禎,邱睿,李春艷,梁漫春,張輝,剛直,徐紅.  Chinese Physics C. 2017(01)
[2]粒子輸運(yùn)蒙特卡羅模擬現(xiàn)狀概述[J]. 鄧力,李剛.  計(jì)算物理. 2010(06)
[3]內(nèi)照射小器官劑量計(jì)算中減方差技巧的比較和應(yīng)用[J]. 張嶄,李君利,武禎,范佳錦,曾志.  清華大學(xué)學(xué)報(bào)(自然科學(xué)版). 2007(S1)
[4]一種改進(jìn)的計(jì)算探測(cè)器校正因子的相關(guān)抽樣方法[J]. 武禎,李君利,程建平.  高能物理與核物理. 2006(08)



本文編號(hào)：3115421