蒙特卡罗粒子输运软件JMCT的网格计数功能设计与实现  

Design and Realization of Mesh Tally in General Monte Carlo Particle Transport Code JMCT

作　　者：付元光[1,2] 郑建华[3] 上官丹骅[2] 李瑞[1,2] 李刚[2] 马彦[2] 邓力[2] 

机构地区：[1]中国工程物理研究院高性能数值模拟软件中心,北京100088 [2]北京应用物理与计算数学研究所,北京100094 [3]中国工程物理研究院激光聚变研究中心,四川绵阳621900

出　　处：《计算物理》2016年第5期581-586,共6页Chinese Journal of Computational Physics

基　　金：中国工程物理研究院科学技术发展基金(2011A0103006,2012A0102005);863计划(2012AA01A303)资助项目

摘　　要：基于蒙特卡罗输运的网格计数能够统计指定区域内密集的分布计数量.为实现JMCT统计局部计数量分布,设计并实现了JMCT网格计数功能,该功能支持三种正交几何(xyz直角几何、rθz柱几何、rθφ球几何)的均匀及不均匀剖分;介绍xyz直角几何网格计数的算法;基于大亚湾反应堆pin-by-pin模型、Venus模型以及ITER次临界包层模型等初步验证了网格计数功能;应用铀阵列模型,采用若干种网格剖分办法,在单机上对比测试了xyz直角几何下JMCT和MCNP5网格计数功能的串行性能,结果显示JMCT的计数时间较短,具有较高的效率.Mesh tally function of monte Carlo method can give a detailed and intensive calculation of flux distribution in specific volumes. To realize such function in JMCT mesh tally function are designed and realized. It supports non-uniform mesh in three kinds of orthogonal geometry （ xyz of rectangular coordinates, rθz of cylindrical coordinates, and rθφof spherical coordinates） . Algorithm for rectangular coordinates is discussed. Calculation on DAYAWAN reactor core pin-by-pin model, Venus benchmark model and a 1-D ITER model verifies preliminarily correctness of JMCT mesh tally. Furthermore, U-array benchmark model is used to test serial performance of JMCT mesh tally. Both JMCT and MCNP5 use same xyz mesh grids and run under same condition. It shows that JMCT takes less time consuming and has higher performance dealing with xyz geometry.

关 键 词：蒙特卡罗模拟 网格计数 JMCT 

分 类 号：TL31[核科学技术—核技术及应用]