MCNP程序计数分段卡存在的问题及解决办法.pdfVIP

第31卷 第 4期 核 动 力 工 程 V0．1．31．NO．4 2 O l0 年 8月 NuclearPowerEngineering Aug．2 0 10 文章编号：0258-0926(2010)04-0016-05 MCNP程序计数分段卡存在的 问题及解决办法 解家春 ，一，赵守智2，孙 征2，贾宝山 (1．清华大学工程物理系，北京，100084；2．中国原子能科学研究院反应堆工程研究设计所，北京，102314) 摘要：当应用蒙特卡罗程序 (MCNP)的计数分段卡 (Fs卡)配合其他计数卡进行分段计数时，可能给 出错误的计数结果。通过比较 Fs卡对同一几何体的3种模型的分段计数结果发现，填充形式模型的燃料棒 区域的计数结果不正确。其原因在于燃料棒栅元由栅元集合卡 (Universe卡)和填充卡 (Fill卡)填充得到， 而分段计数时MCNP程序认为被填充的栅元仍位于Universe卡定义的几何体处。据此提出了在Universe卡 定义的几何体处对被填充栅元进行分段计数的解决方法。对例题的径向中子注量率分布的计算表明，该方法 能够解决Fs卡对填充形式描述的几何体进行分段计数时出现的错误。 关键词：MCNP程序；FS卡；分段计数；填充；径向中子注量率 中图分类号 ：TL329 文献标识码 ：A 1 引 言 2 FS卡分段计数方法 蒙特卡罗程序 (MCNP)’是国际通用的粒子 本文以Fs卡在中子注量率空间分布计算中 输运计算程序 ，能够计算中子、光子、电子或中 的应用为例，说明Fs卡的分段计数方法。功率密 子．光子．电子耦合输运问题，也可以计算临界系 度分布的计算与之类似。所计算例题为均匀堆， 统的本征值问题t¨。目前，MCNP在反应堆设计 ， 几何结构如图 1所示。堆芯燃料为 UC，日U与 尤其是小型堆的设计中得到了广泛的应用，使用 U的原子数之比为2：8，反射层材料为石墨。所 该程序基本上可以得到所有稳态的物理参数。在 有材料密度见文献 2【】。 这些参数中，与空间位置有关的中子注量率分布 要得到燃料区的径向中子注量率分布(为简单 和功率密度分布等往往需要计数分段卡 (Fs卡) 起见，求全高度的径向分布)有 2种方法：①将 燃料区沿半径方向描述为多个栅元；②燃料区只 配合其他计数卡得到。FS卡是 MCNP程序的一 描述为 1个栅元，使用Fs卡进行径向计数划分。 个重要的辅助计数卡，它可以将栅元或曲面的计 数按照空间位置分段给出，从而避免了仅为计数 目的而将几何体的模型描述为多个栅元或曲面， 减小了编写输入卡的工作量，降低了出错几率。 在实际工作中发现，当应用FS卡进行分段计 三4~50cm 数时可能给出错误的结果，而出现错误的原因与 几何体模型的描述方式密切相关。本文对同一几 图1 均匀堆示意图 Fig．1 ViewofHomogeneousReactor 何体采用不同的描述方式建立了3种计算模型， 并应用Fs卡分别对各模型进行了分段计数。通过 使用 Fs卡计算时的输入文件如下 (只给出 比较计数结果，找出了FS卡在分段计数时出现错 Fs卡相关部分)： 误的原因，并根据这_发现，在不修改源程序的 F4：n1 基础上，提出了在栅元集合卡 (Universe卡)定 FS4—3—4—5—6—7—8—9—10一ll一12一l3 义的几何体处