扩散方程计算反应堆控制棒价值的误差的处理方法研究.pdf

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 扩散方程计算反应堆控制棒价值的 误差的处理方法 高振家李富胡永明 北京清华大学核能与新能源技术研究院， 100084 (+8610nudetector@gmail．com 摘要。本文展示和分析了采用扩散方程相对于输运方程在计算反应堆控制棒价值中产生的误差．同时给出 了控制棒的临界黑度值，当控制棒的黑度值超过临界的时候，由扩散方程计算得到的结果应当予以修正，甚 至只能用输运方程计算．最后，文章给出了处理扩散方程误差的一种方法． 关键词：扩散方程，输运方程，控制棒，黑度，绝对误差，相对误差。 1引言 在实际运行的反应堆中，控制棒是置于反应堆之外，它的位置被水填充，因此，反应堆 cm1，相应的黑 组件热中子的吸收截面非常小，根据下文(1)式的计算方法，大概是0．09 度大概是0．14。即使是停堆，热吸收截面也大概只是0．10cm1，相应的黑度是0．17。因此， 通常情况下扩散计算也能满足精度要求。但是，在一些研究堆中，组件热中子的吸收截面是 非常大的。例如，在含有铪控制棒的反应堆中，组件热中子的吸收截面大概是0．25锄一， 其相应的黑度大于是0．3．在这种情况下，扩散方程计算的结果不能满足精度要求，应该被修 正。 依据简单的堆芯模型，作者分别编写了由扩散方程和输运方程(Sn-16)计算的程序．文 章首先描述了在没有控制棒的反应堆中，扩散方程能够得到和输运方程几乎相同的结果。这 为后面进行控制棒价值的计算提供了基础。然后，在反应堆中插入控制棒，计算其价值，由 扩散方程得到的反应堆控制棒价值的误差将控制棒黑度值的增加而增加。当误差太大以至于 不能被接受时，由扩散方程得到的结果就应当被修正，甚至只能用输运方程计算。因此，本 文计算结果可以为反应堆控制棒计算提供一定参考。 2反应堆模型 2．1反应堆模型的形状 不失一般性，本文采用一维平板模型进行分析，燃料区的宽度是80cm：反射层位于反应 堆的两端，宽度是20cm。根据对称性，选取反应堆的一半进行计算，也就是燃料部分是40cm， 反射层是20cm，如图l所示． ／ 燃料区 。／ 反射层 图一反应堆模型 2．2模型参数 321 第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 IIII■—皇——■——●■■■■■■■■■■曩置_■■■■昌暑昌皇鲁鲁昌蕾皇皇皇薯—暑皇曩■墨皇皇暑詈=詈皇昌昌昌皇=皇詈=詈===昌盲皇篁皇譬曩皇昌暑暑葛篁皇昌昌皇=皇皇暑昌暑詈詈暑暑■—墨覃暑鲁暑暑鲁昌皇暑詈●皇鼻■ 在扩散方程和输运方程计算中的参数都由组件计算程序CASMO产生，见表l。这样就 保证了扩散方程产生的误差主要是由于方法本身，而不是由于不同的截面参数之类。 另外，在输运计算中使用了离散坐标法(S16)。同时，也考虑了能群内部的自散射效应。 在两种方法中，左边界条件是反射边界条件，右边界是空边界条件。 表l反应堆参数 l ∑g_2／(cm’1 I区域 g D／(cm) ∑口／(cm-1)厄，／(cm‘1)Zg·l／(cm‘1) I ) 燃料区三≥2嬲：麓苗％怒嚣3舢．000。0。03％盟淼 l 1．253470．0004991070．0000000．0000000．0555327 反射层 2 0．17756200000000．0000000．000000 g-能群数；D扩散系数；∑口．吸收截面伲，．y·裂