高温气冷堆控制棒区不连续因子的计算与应用 calculation and application of discontinuity factors in control rod region for high temperature gas-cooled reactor.pdfVIP

第29卷第6期 核动力工程 、，01．29．No．6 0 2OO8 年12月 NuclearPower Dec．2O8 Engineering 文章编号：0258·0926(2008)06-Oool-06 高温气冷堆控制棒区不连续因子的计算与应用 周旭华，李富，王登营，颜见秋，吕炜枫，韩仁余 (清华大学核能与新能源技术研究设计院．北京．Io0084) 摘要：应用不连续因子理论修正的扩散方程，对高温气冷堆侧反射层中控制棒区域的强吸收体和空洞区 的局部均匀化进行处理。解决了均匀化区域不存在裂变中子源、边界存在强中子流所带来的的难题，并以输 运计算的解为基准进行验算。结果表明：对控制棒区域可进行局部均匀化以及采用简化的边界条件计算不连 续因子；用不连续因子修正的扩散方程可准确地处理控制棒。采用不连续因子修正的扩散方程计算高温气冷 堆控制棒在计算精度、计算时间上均有优势。 关键词：高温气冷堆；控制棒；均匀化；不连续因子；扩散方程 中图分类号：TL329文献标识码：A 1 引 言 究控制棒价值带来误差。为此，有必要对每根控 球床式高温气冷堆的控制棒布置在堆芯外的 制棒插入形成的强吸收体区和控制棒拔出形成的 侧反射层中。由于反射层是中子吸收截面小的石 空洞区进行局部均匀化，在全堆模型中模拟实际 墨，而控制棒区是中子吸收截面很大的碳化硼， 的控制棒位置，且为了保证计算精度，需采用不 同时，控制棒的抽出在这一区域又形成空洞，使 连续因子修正的扩散方法进行计算。 这一区域物质的非均匀性非常强烈；可以预见， 侧反射层中的中子注量率将急剧地变化；而且在 2不连续因子的需求 运行过程中还会出现部分控制棒插入、部分控制 如图1所示的反应堆模型中，侧反射层均匀 棒拔出的复杂情况。因此，直接使用中子扩散方 布置有24根控制棒，控制棒的插入可以形成强吸 法进行物理分析计算已经不合适。另一方面，进 收体区域，控制棒的拔出可以形成空洞区域。各 行全堆细网输运计算虽是较准确的方法，但是计 种控制棒组合状态都可以由强吸收体区域和空洞 算量过大(尤其是针对三维模型)，且收敛速度很 区域排列组合而成。分别选取包括强吸收体区域 慢，并不实用。对控制棒区域进行均匀化，然后 或空洞区域的l／24堆芯扇区作为均匀化模型对其 利用扩散方法进行分析计算更为可行，现已在工 进行均匀化(图2)，并假设模型上下边界处为反射 程设计上广泛采用。 边界条件，求出强吸收体区域和空洞区域的均匀 通常的均匀化策略是将控制棒所在的整个圆 化截面以及边界处的不连续因子，然后将均匀化 环区域进行整体均匀化【l‘2J，得到的均匀化参数 参数和不连续因子应用于全堆扩散计算。 在整个控制棒圆环区域是相同的。然而，采用该 均匀化参数计算得到控制棒附近(包括堆芯)的中 3不连续因子的原理 子注量率分布值与真实值存在较大的误差；特别 对某个非均匀化区域进行均匀化处理的目标 对部分控制棒插入、部分控制棒拔出的组合情况 是：保证均匀化前后区域的总反应率守恒和边界 更难处理。这对于研究堆芯功率分布、燃耗变化、 处的中子流守恒‘31。均匀化后的群参数计算如式 探测器读数等情况是不可接受的。同时，中子注 (1)、式(2)。 量率分布误差还会