蒙特卡罗粒子输运软件JMCT的网格计数功能设计与实现-计算物理.PDFVIP

第３３卷 第５期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３３，Ｎｏ．５ ２０１６年９月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｓｅｐ． ，２０１６ 文章编号：１００１⁃２４６Ｘ（２０１６）０５⁃０５８１⁃０６ 蒙特卡罗粒子输运软件ＪＭＣＴ 的网格计数 功能设计与实现 １，３ ２ ３ １，３ ３ ３ ３ 付元光 ，　 郑建华 ，　 上官丹骅 ，　 李　 瑞 ，　 李　 刚，　 马　 彦 ，　 邓　 力 （１ 中国工程物理研究院高性能数值模拟软件中心，北京　 １０００８８；２ 中国工程物理研究院激光聚变研究中心，四川 绵阳　 ６２１９００； ３ 北京应用物理与计算数学研究所，北京　 １０００９４） 摘　 要：基于蒙特卡罗输运的网格计数能够统计指定区域内密集的分布计数量．为实现ＪＭＣＴ统计局部计数量分 布，设计并实现了ＪＭＣＴ 网格计数功能，该功能支持三种正交几何（ｘｙｚ 直角几何、ｒθｚ 柱几何、ｒθφ球几何）的均匀及 不均匀剖分；介绍ｘｙｚ 直角几何网格计数的算法；基于大亚湾反应堆ｐｉｎ⁃ｂｙ⁃ｐｉｎ模型、Ｖｅｎｕｓ模型以及ＩＴＥＲ次临界 包层模型等初步验证了网格计数功能；应用铀阵列模型，采用若干种网格剖分办法，在单机上对比测试了ｘｙｚ 直角 几何下ＪＭＣＴ和ＭＣＮＰ５ 网格计数功能的串行性能，结果显示ＪＭＣＴ 的计数时间较短，具有较高的效率． 关键词：蒙特卡罗模拟；网格计数；ＪＭＣＴ 中图分类号：ＴＬ３１ 文献标志码：Ａ ０　 引言 基于蒙特卡罗（ＭＣ）方法的粒子输运模拟程序，因具有强大的计数功能，常被用于解决反应堆屏蔽设 计、探测器响应计算等实际问题．一般的ＭＣ程序提供的计数功能主要是基于面几何的流计数、基于几何体 的径迹长度计数和基于点的探测器计数，这些计数的精确度和程序使用者的几何设置有很大关系，得到精确 的局部量计数往往比较难；程序使用者必须通过建立足够密集的几何栅元或探测计数点，才能和确定论程序 达到相同的局部量和分布量计算能力，而密集栅元或大量探测点的建模和储存是比较困难的，以ＭＣＮＰ４Ｃ 为例，有栅元和计数个数的限制，导致其在计算分布量时有所不足． 网格计数（Ｍｅｓｈ Ｔａｌｌｙ）能很好地解决以上问题，其原理为用一套结构网格罩住计算模型中感兴趣的区 域，粒子输运过程仍在原来模型的几何上进行，同时需要统计粒子穿过这套虚拟网格的径迹信息．这种方法 通过虚拟的网格代替实际的栅元或点探测器，在不改变原有几何的情况下，也能得到计数量的分布．同时，在 计数过程中，存储一个虚拟网格只需要存储该网格上径迹长度的累加值，远小于存储一个真实几何栅元所耗 ［１］ 费的内存，因此网格计数的规模可以较大．国际上一些用户群较广的ＭＣ程序（如ＭＣＮＰ５ ）都具备网格计 数功能． ［２］ ＪＭＣＴ是北京应用物理与计算数学研究所（ＩＡＰＣＭ）开发的通用型蒙特卡罗粒子输运模拟程序 ，早期 的ＪＭＣＴ计数功能主要是基于几何栅元的径迹长度计数．使用ＪＭＣＴ计算某些屏蔽问题时，有时不仅仅需要 计算若干局部量，还需要计算全局分布量，而径迹长度计数功能由于几何体个数、建模难度等限制，不能很好 地满足需求，因此开发基于ＪＭＣＴ 的网格计数功能很有必要．本工作设计了ＪＭＣＴ 的网格计数算法，实现了计 数功能，并进行了该功能的数值验证和性能测试． １　 算法 网格计数中的网格一般采用正交坐标系下的结构网格，主要分为笛卡尔坐标系下的ｘｙｚ 型、柱坐标系下 的ｒθｚ 型和球坐标系下的ｒθφ型．在ＪＭＣＴ上实现了三种正交坐标系下的网格计数．下面以ｘｙｚ 型为例介绍计 收稿日期：２０１６－０１－２２；修回日期：２０１６－０３－０２ 基金项目：中国工程物理研究院科学技术发展基金（２０１１Ａ０１