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

【文章來源】：原子能科學技術. 2020,54(07)北大核心EICSCD

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

【部分圖文】：

AIS方法流程圖

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

DXTRAN球上虛粒子的抽樣方法

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



本文編號：3115421