HFETR与MCNP程序的接口程序的开发研究.pdf

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 邓才玉邱立青傅蓉 (中国核动力研究设计院，成都610005) 的堆芯描述．该程序根据HFETR堆芯装载布置、燃料元件功率分布、燃耗深度及控制棒高度等参数，自动生 ，成跹NP程序的输入文件．结果表明：该接口程序对HFETR堆芯的材料及结构描述是准确的，我们所开发的接 口程序是有效的，从而提高了HFETR蒙特卡罗计算效率及准确性． MCNPIhP 关键词：tlFLrrR蒙特卡罗方法HCNP接口程序 1前言 利用MCNP程序对物理问题进行计算，首先必须对问题进行详细准确的描述。这关系到问 题计算的可靠性和准确性。对核反应堆而言，就必须对反应堆堆芯进行准确的描述，包括堆 芯内几何结构、各种不同的材料以及源的描述Ⅲ。 I-IFETR堆芯有三百多个六角形栅元埋1，各栅元位置布置的组件不同，其描述各不相同，如 燃料元件占据一个栅元，要描述这个燃料元件栅元的几何结构则必须描述34个曲面和29个 几何体及其材料(如图1)。要描述整个堆芯结构，一般需要描述1万多个曲面和3千多个几何 体及其材料；并且堆芯装载稍有变动，则相应的许多曲面和几何体都要重新描述。HFETR堆芯 内现在一般布置有70--80盒燃料元件，每盒燃料元件的燃耗和功率都不一样，则须描述420 480个栅元的材料和源分布，且材料的具体组成也不宜确定。 因此HFETR堆芯的MCNP程序描述非常困难、复杂，靠人工对HFETR堆芯进行精细描述 所花的时间很长而且容易出错。 、八U汐 粕、，．7，ll。I氟19’器，Tt、翟舞表I 曩喀．L霉．∞。堰．“、lS．1嚣，∞、嚣．拍．蕊。驾为菇 曩‘“碍、lT、Zl、薹身慧#元舞■萏—k 。～，， 图1燃科元件29层几何体及材料示意图 647 第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 2MCNPl NP程序 考虑到IIF-ETR堆芯的跹NP程序描述的难度和复杂性，采用自行研制开发的前处理程序 置、燃料元件功率分布和燃耗、控制棒高度等参数，来生成跹NP程序输入文件所要求的堆芯 描述，包括曲面描述、栅元几何描述、栅元材料描述和源的描述，从而自动完成需由用户填 写的托NP程序输入文件，对堆芯进行蒙特卡罗计算。 2．1程序编写 由于HFETR堆芯描述非常复杂，如果采用面向过程的计算机程序语言来描述将十分烦琐 和困难¨1的，因此采用面向对象的计算机程序语言c++来描述。源程序包括MCNPINP．h和 Visual MCNPINP．cpp两个文件，约2200行。采用Microsoftc++6．0H1编译通过。 硬件要求：intelpentitun200及其兼容或以上，32M以上内存； 9x、Windows XF。 操作系统：Windows 2000或Windows 采用C++语言创建一个DXSY类，该类自己处理几何描述和材料描述，堆芯内每一个实际 栅元就是DXSY类的一个具体实例：对象。 DxSY类如下： classDXSY { private： int type∥栅元类型 float shape√／栅元形状 float dense；／／栅元材料密度 double x：／／栅元中心位置x坐标 double