MCNP程序在多群中子输运计算中应用.pdfVIP

中国核科学技术进展报告（第一卷） 计算物理分卷 Progress Report on China Nuclear Science Technology （Vol.1 ） 2009 年11 月 MCNP 程序在多群中子输运计算中的应用 1 1 2 黄伟兵 ，陈伦寿 ，吴宏春 （1.三门核电有限公司，浙江 台州 317112; 2. 西安交通大学，陕西 西安 710049 ） 摘要：详细介绍了如何利用MCNP 程序进行多群输运计算，并用MCNP 程序计算了多群临界特征值问题和多群固定 源问题的基准算例，计算结果与参考值吻合良好。MCNP 程序在多群输运计算中的成功应用具有重要意义，它不仅 为多群中子输运方程的求解提供了新的解决途径，同时，也为多群输运程序的开发提供了有效的校核工具。 关键词：MCNP ；多群；输运计算 MCNP[1](Monte Carlo N-Particle code) 是由美国Los Alamos 国家实验室X-5 蒙特卡罗小组研制 的用于复杂三维几何粒子输运计算的大型多功能蒙特卡罗程序。它可用于计算中子、光子、电子独 立输运及耦合输运问题，也可用于临界系统本征值问题的计算。在国内，MCNP 程序被广泛应用于 反应堆堆芯物理设计、核临界安全分析、辐射剂量计算、堆物理实验分析、物理保健问题研究及核 辐射损伤研究等各个领域。 MCNP 程序发布时，自身所携带的核数据库为精细的点截面数据库，因此，MCNP 程序的连续 能量输运计算功能深为人们熟知。然而，MCNP 程序强大的多群输运计算能力，在国内却鲜为人知， 至今未见国内有学者发表过关于MCNP 程序在多群计算方面的文章。国内学者基于MCNP 程序自主 开发的蒙特卡罗多群输运计算程序 MCMG ，虽然在多群临界计算上，相比MCNP 程序具有速度上 的优势，但是由于其屏蔽了MCNP 程序的一些功能，MCMG 程序在几何描述及对空间区域进行中 子通量计数上，都存在一些缺陷，并且MCMG 程序不能用于多群固定源问题的计算。 因此，探讨MCNP 程序的多群输运计算功能具有如下重要意义： （1）MCNP 程序是通用的三维粒子输运计算程序，其在几何处理上具有任意性，程序的计算精 度在国际上被广泛认可，因此，MCNP 程序的计算结果可以作为新研制的多群输运程序的参考值， 以便对新研制程序的计算精度进行校核。 （2 ）MCNP 程序的多群输运计算功能，为求解任意几何的多群临界特征值问题及固定源问题提 供了一种新的有效途径。 （3 ）得益于MCNP 程序的广泛应用，其多群输运计算的功能可以在很多领域得以应用。 本文首先介绍了如何利用MCNP 程序进行多群输运计算，然后利用MCNP 程序计算了一些多群 临界问题及固定源问题的基准题，以验证MCNP 程序的多群输运计算能力。 1 多群截面文件格式 在利用MCNP 程序进行多群输运计算之前，首先得提供各材料的多群截面文件。截面文件xsfilen 包含两大部分，分别为信息块和截面数据块。信息块占据截面文件的1 至12 行，主要包含材料的 ZAID 、各反应类型截面在截面数据块中的位置等信息。在此需要强调的是，信息块的填写格式 非常严格，如有些细微的不规范都会导致程序终止运行。如信息块第1 行第1 列必须为空，如果不 作者简介：黄伟兵（1982—），男，浙江省临海市人，助理工程师，硕士，反应堆物理 136 空或者空2 列都会使程序运行出错。另外，信息块中3 至6 行整型和实型数据也须严格遵循规范， 如不规范填写也会导致程序运行中断。截面数据块则给出了材料各能群各反应类型的截面数据。截 面文件的具体格式见文献[1]中的附件F-Ⅲ与F-Ⅷ。编辑截面文件时，可以一种材料对应一个截面文 件，也可以一个截面文件包含全部材料的截面数据。运行程序时，截面文件要求放在与MCNP 程序 的相同目录下。 2 MCNP 程序的多群运行模式 MCNP 程序会根据程序输入卡片的内容自动判断是否进行多群输运计算。多群模式的输入卡片 相对于连续能量模式的输入卡片，只需要增加如下两个额外的卡片： （1）MGOPT 卡 mgopt f ng 其中ng 表示所计算问题