基于CFD方法的球床式高温气冷堆稳态热工水力分析.pdfVIP

第36 卷 第12 期 核 技 术 Vol. 36, No.12 2013 年12 月 NUCLEAR TECHNIQUES December 2013 基于CFD 方法的球床式高温气冷堆 稳态热工水力分析 宋士雄1,2 魏 泉1,2 蔡翔舟1,3 郭 威1,3 1 （中国科学院上海应用物理研究所嘉定园区 上海 201800 ） 2 （中国科学院大学 北京 100049 ） 3 （中国科学院核辐射与核能技术重点实验室 上海 201800 ） 摘要 基于CFD 通用计算程序Fluent ，开发了多孔介质流固两相局域非热平衡模型，对PBMR-400 满功率名 义工况下稳态热工水力行为进行了分析。高温气冷堆堆芯的热工水力计算模型通常为带有强内热源宏观的多 孔介质模型，Fluent 程序本身的多孔介质模型为局域热平衡模型，无法恰当描述堆芯热工水力行为。通过开启 Fluent 的自定义标量场(UDS)功能并在动量方程上附加阻力源项，实现了局域非热平衡模型，能够精确地计算 球床堆堆芯热工水力行为。CFD 程序使用了非常精细的网格，可以计算更加精细的流场和温场分布。计算表 明，Fluent 的结果与国际上著名的球床堆热工水力程序THERMIX，TINTE 符合较好，可以作为球床堆热工水 力分析工具。 关键词 CFD (Computational Fluid Dynamics) ，高温气冷堆，多孔介质 中图分类号 TL349 DOI: 10.11889/j.0253-3219.2013.hjs.36.120601 准确可靠的反应堆热工水力分析对于球床式高 界的热耦合，冷却剂边界与反射层边界的热耦合。 温反应堆的设计是非常重要的。目前国际上球床式 PBMR-400 (Pebble Bed Module Reactor)是南非 高温气冷堆热工水力行为分析程序 THERMIX[1−3] 国营电力设计的热功率为 400 MW 、电功率为 和 TINTE[4–5] ，均采用二维轴对称局域非热平衡多 165 MWe 的模块式高温气冷球床式反应堆[14–15] 。反 孔介质模型，流场方程为附加多孔介质阻力源项的 应堆堆芯直接连接氦气轮机，设置了具有闭路直接 准静态方程，冷却剂温场和固相球床的温场分别独 循环特征的能量转换单元，堆芯燃料球采用环形布 立求解，各个模块之间进行迭代并收敛。随着计算 置，具有很强的非能动余热排出能力。反应堆堆芯 机性能的大幅度提高以及算法的改进，在反应堆设 活性区共有大约 452 000 个燃料球，内外半径分别 计及研究中CFD (Computational Fluid Dynamics)方 为1.0 m、1.85 m，活性区高度为11 m。堆芯运行 法开始被采用[6–7]，并且实现了CFD 程序与中子学 压力为9 MPa，入口冷却剂温度约为500 °C，出口 程序、系统热工程序的耦合计算[8−10]。Becker 等[11] 氦气温度约为 900 °C ，堆芯的平均功率密度为 基于 CFD 开发了球床反应堆热工程序，通过与 4.8 MW·m−3 。本文基于Fluent 的自定义标量场及自 SANA 装置的实验结果比较，验证了程序的适用性。 定义函数功能开发了局域非热平衡多孔介质模型， CFD 程序在计算过程中可以使用更加精确的 对PBMR-400 满功率名义工况的热工水力行为进行 流体力学方程与很精细的网格，因此可以对局部流 了分