注量展开节块法求解三维六角形几何中子扩散方程 flux expansion nodal method for solving neutron diffusion equations in hexagonal-z geometry.pdfVIP

第 26卷 第 5期 核 动 力 工 程 Vol. 26. No.5 2 0 0 5年 10月 Oct. 2 0 0 5 Nuclear Power Engineering 文章编号 0258-0926(2005)05-0425-06 注量展开节块法求解三维六角形几何 中子扩散方程 夏榜样 谢仲生 西安交通大学核能科学与工程系 710049 摘要 提出了一种在三维六角形几何节块内数值求解中子扩散方程的节块法 该方法把节块内各群中子 注量分布用解析基函数近似展开 为了改善节块耦合关系 提出了一种新的节块边界条件 面平均偏流零次 矩和一次矩同时保持连续 此外 将响应矩阵技术应用于迭代求解过程 使得该方法具有较高的计算效率 基于本文提出的模型 发展了三维六角形组件中子扩散计算程序FEMHEX通过对二维 三维VVER基准问题 校验计算表明 该方法能高效 准确的给出有效增殖系数以及节块功率分布 关键词 六角形几何 中子扩散方程 节块法 解析基函数 中图分类号 TL329.2 文献标识码 A 1 引 言 2 理论模型 在核反应堆堆芯物理设计中 精确地获得反 2.1节块内注量近似展开 应堆内中子注量和功率的空间分布是一个重要的 考虑如图 1所示的反应堆内第n 个均匀化的 课题 在过去20年中 对于正方形燃料组件的反 六角形节块 其三维两群中子扩散方程可写成如 应堆 已经发展了许多完善的先进节块法 而 下的矩阵形式(为了简便 省略去节块n 的编号) 应用于六角形组件反应堆(如VVER型压水堆 快 2 ( ) ( ) F r F1 r 1 中子增殖堆等)的先进节块法相对较少[17] 主要 2 L ( ) 0 (1) F ( ) F r 1 r 2 原因是在正方形几何中行之有效的横向积 上式中所用的符号都是习惯表示方法 其中 分技术应用到六角形几何时出现了奇异项 r=(x,y,z) 矩阵的定义 为 L 问题 从计算结果来看 美国西屋公司的 k 1v v [4] [5] r 1 eff f 1 f 2 ANC-H方法和韩国的 AFEN方法是目前用 D k D L 1 eff 1 (2) 于六角形组件的最为成功的两种方法 二 s r 2 者都有效地克服了横向积分中的奇异项问 D D