基于MC方法地反应堆燃料组件参数非均匀计算研究.pdfVIP

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 基于MC方法的反应堆燃料组件参数 非均匀计算研究水 项骏军周永强吴宏春 (西安交通大学．710049，西安) 摘要：组件均匀化参数计算是反应堆物理设计计算的基础，传统的组件参数计算方法都是经过栅元均匀化 和组件均匀化两步计算。对栅元进行均匀化近似会引入误差．本文研究了基于蒙特卡罗方法的组件参数非均 匀化计算方法，对组件栅元直接进行非均匀计算，从而获得组件参数，提高了组件参数的计算精度． 关键词t蒙特卡罗非均匀组件参数 1引言 精确地计算燃料组件均匀化参数是反应堆物理设计中的一个重要问题，早在70年代到 80年代，国际上研制了许多计算组件均匀化参数的计算软件包，如CASMOllJ， 的组件参数计算程序TPFAPi31。由于当时计算机计算速度和存储容量的限制，进行全组件精 细几何输运计算存在一定的困难，所以这些软件在计算组件均匀化参数时普遍采用了栅元均 匀化方法：首先进行栅元计算，获得各种不同类型栅元的等效均匀化参数；再使用栅元计算 得到的等效均匀化参数进行组件计算。这种均匀化近似会引起栅元内中子通量密度偏离真实 值，从而使计算得到的组件有效增殖系数不够准确141，所以这种“栅元．组件”的两步计算方式 在精度上存在一些不足。 蒙特卡罗方法是一种非确定论方法，与一般确定论方法不同，蒙特卡罗方法具有很强的 几何处理能力，对粒子的跟踪忠实子粒子物理作用过程，计算精度高，缺点是需要跟踪大量 的粒子，计算效率较低。随着当今计算机技术的飞速发展，蒙特卡罗方法在反应堆物理的各 个领域中发挥着越重要的作用。 本文在先进均匀化理论的基础上嘲，结合DRAGON程序和多群蒙特卡罗程序MCMG3B 各自的优点．开发了基于MC方法的输运．燃耗耦合的反应堆燃料组件参数非均匀计算程序 DRAGON．M[C。 新的DRAGON-MC程序同时具备了使用确定论方法和非确定论方法计算的能力；程序结 合了DRAGON和MCMG3B程序各自的优点：(1)DRAGON程序宏观截面的产生考虑了几 (2)MCMG3B能够处理复杂 何、温度及截面的共振自屏效应，这是MCMG3B所不具备的； (3)DRAGON程 的几何问题，同时具有很高的计算精度，这是DRAGON程序所不具备的； (4)DRAGON程序的编辑模块具有很高 序可以进行燃耗计算，这是MCMG3B所不具备的： 的灵活性，能够方便地对计算结果进行并群并区．同时对DRAGON程序和MCMG3B程序做了 —些改进：(1)增加了计算组件功率空间分布的功能；(2)增加了计算组件不连续因子的功能； ’国家自然科学基金资助项目 347 第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 2计算方法简介 序在相同精度下，比MCNP程序计算时间少得多。具备了全组件精细几何非均匀计算的能力。 反应堆的中子输运方程为： 增殖系统的临界问题可归结无外源，定常情况下(1)式的定解问题，也可以归结为下 列方程的本征值问题 ￡◇：二l_朋◇ (2) K啦 !／K够为输运方程的特征值，o为它所对应的特征函数，足∥称为系统的有效增殖系数。 采用源迭代方法，K够的计算可归结为 {‘2：qp箩’(，．!【：j!：2=zX；zg，‘n7，(h-1)(、r、，t’h：。，．．．。：；， I秽’(，)=∑％∑，，窖(，)①≯(，) 当源位置分布未知时，可在裂变区域任意给定一点或多点，初始能群go服从裂变中子能谱分 布以．初始方向Qo服从分布1／4a(即空间各向同性分布)。初始权为