三维离散纵标方法在堆内构件释热率计算中的初步应用.pdfVIP

三维离散纵标方法在堆内构件释热率计算中的初 步应用 1 1 2 1 1 杨寿海 ，陈义学 ，张斌 ，王伟金 ，靳忠敏 (1. 华北电力大学 核科学与工程学院，北京 ，102206) (2. 国家核电技术公司 软件中心，北京 ，100033) 摘要：离散纵标（又称SN ）方法是反应堆堆内构件释热率、RPV 处快中子注量率计算中最常用的方法之一。 计算机的飞速发展和三维离散纵标方法的不断完善，使得三维离散纵标法应用于工程设计成为可能。本文 以美国H.B.Robinson-2 号机组压力容器基准实验模型为研究对象，采用基于三维离散纵标方法程序计算了 压力辐照监督管处中子能谱及典型核素的活度值，计算结果表明：采用三维离散纵标方法得到的辐照监督 管处中子能谱与报告结果吻合较好，典型核素活度的C/M 值(计算值与测量值之比)为1.04±0.04,计算精度满 足工程要求。在此基础上进行释热率分布计算，对三维离散纵标方法应用于堆内释热率计算的工程实践进 行讨论。 关键词：三维离散纵标法；H.B.Robinson-2 基准；堆内构件；释热率 中图分类号：TL315.6 文献标识码：A 1 前言 三维中子与γ射线的空间分布，同时具有进 由于离散纵标方法在解决深穿透问题 行Keff 、固定源、临界浓度及临界尺寸搜索 方面的突出优势，在过去的几十年内它一直 等强大功能,尤其适合于解决“深穿透”辐射 是计算堆内构件释热率及 RPV 快中子注量 屏蔽问题，国际上已经将该程序应用于 率的主流方法。但是受计算机内存容量及计 VENUS-3 的压力容器快中子注量计算的工 算速度的限制，以及三维离散纵标方法本身 程实践中[4] ，同时进行了大量的基准验证[5、 的不完善，前人在计算中大多采用一维、二 6] ，初步证明了该程序应用于核反应堆工程 维离散纵标方法，目前国内屏蔽计算主要采 设计中的可行性与可靠性。 用一维离散纵标法程序ANISN 和二维离散 TORT 程序在我国核工程中的应用分析 纵标法程序 DOT ，并用一维与二维方法合 尚处于起步阶段，本工作针对TORT 程序在 成三维分布，由于在模型建立方面引入了较 大型压水堆中的应用分析，开展探索性尝 多近似，因此合成的三维计算结果的可信度 试。由于三维计算模型输入数据比较复杂， 不高。随着计算机内存容量及计算速度的飞 本工作还使用了三维可视化建模辅助工具， 速发展，同时三维算法及三维辅助工具的开 它不仅能生成 TORT 程序所需要的几何材 发，使三维离散纵标方法应用于反应堆的工 料模型，还能对所构建的几何模型进行可视 程设计计算成为现实的可能[1] 。 化处理，包括源分布的可视化及结果的可视 由于TORT[2]程序在在我国核工程中的 化等。 应用分析尚处于起步阶段，本工作采用三维 本工作堆芯外输运计算采用基于核评 离散纵标程序，针对美国H.B.Robinson-2 号 价数据库ENDF/B-VI[7]开发的MATXS 格式 机组压力容器基准实验[3],得到了辐照监督 的通用多群截面数据库，该数据库是一个 管处中子能谱，并与基准报告的计算结果进 30n×12g、300K、P4 阶展开、包含121 种核 行比较，以验证TORT 程序计算堆内构件释 素的中子、光子及中子光子耦合的多群截面 热率及 RPV 处快中子结果的精度，在此基