孙贤祖 刘兴民 刘林顶

(中国原子能科学研究院 北京 102413)

堆外核测量系统，又被称为堆外中子注量率监测系统，其功能是为了准确地探测核功率和功率变化，是确保核安全的重要系统。堆外核测量系统运行时有通风系统对其进行冷却，以保证内部设备处于正常工作温度。

国内外学者对反应堆堆坑内部的数值模拟研究已有很多。何培峰等人[1]对华龙一号反应堆的一体化堆顶通风结构建立三维模型，研究了在不同通风量下CRDM线圈组件之间的气流速度和压降分布情况。张志强等人[2]采用理论分析和数值模拟的方法，对堆顶流场和局部热点进行了分析，并对堆顶结构进行设计优化。陈宇清等人[3]通过实验与数值模拟的方法，研究了环形夹层入口水温、环形夹层外部风速等因素对环形夹层内自然对流传热的影响。习蒙蒙等人[4]利用三维数值模拟方法，研究了某三代压水堆堆坑不同进风方向和不同进风参数条件下，反应堆堆坑通风系统内的流场与温度场分布情况。马崇扬等人[5]计算了丧失外电源事故下中国实验快堆堆坑内部的温度分布。王予烨等人[6]利用CFD软件对正常运行工况下的示范快堆堆坑空气流域进行三维数值模拟。HUNG T C 等人[7]以SAFR为设计参考，对反应堆非能动空冷系统进行了三维数值模拟。CRAH A 等人[8]使用CFD 方法对欧洲钠冷快堆反应堆冷却系统在堆坑内的余热排出能力进行了分析。SONG M S等人[9]利用CFD方法计算了冷却剂在压力容器间隙中的流速和压力分布，并研究了热工水力方面的瞬态变化。目前以堆外核测量系统为研究对象的研究尚有不足。

为验证示范快堆堆外核测量系统已有的通风系统，为一体化对系统设计提供参考，文章通过建立真实堆外核测量系统计算模型，计算和分析了不同进风口位置、进风量、进风温度对系统内部温度场与流场的影响。

1 堆外核测量系统结构与通风系统

示范快堆堆外核测量系统位于堆容器外侧堆坑混凝土环墙中的封闭方腔中，方腔前部有石墨-铅屏蔽装置，系统本体由多个石墨块体组成，外部覆有钢覆面，底部由多个混凝土柱支撑。方腔上方有14 个堆外电离室通道倾斜插入系统，通道中用于放置探测器。堆外核测量系统结构如图1所示。

图1 堆外核测量系统结构

堆外核测量系统冷却空气由方腔下前处4个送风口送入，送风口半径为125 mm。冷却空气一部分由堆外电离室通道流出，对系统本体及电离室通道内探测器进行冷却，另一部分则在系统和方腔的间隙中进行冷却，并由方腔后上方的5个排风口流出，排风口半径为55 mm。为确保电离室通道内探测器正常运行，须保证测量精度，堆外核测量通风系统设计要求则要保证系统探测器正常工作环境即电离室通道内温度不高于45 ℃。

2 数值计算模型

2.1 系统模型

三维几何模型采用SpaceClaim 软件进行建模，模型以石墨-铅屏蔽板为前边界，其余面以混凝土表面上的钢覆面为边界。模型中固体域包括混凝土支柱、系统内部组件、电离室通道等，环腔和电离室内部区域为流体域。建模时对模型中不影响计算精度的结构进行了适当的简化处理。计算模型如图2所示。

图2 计算模型

2.2 网格模型、边界条件及计算方法

本文采用ANSYS mesh模块进行网格划分，对形状规则的几何模型，如系统内部的石墨块体，采取六面体结构化网格进行划分；对流体域和固体域交界的复杂区域采用四面体非结构化网格，并对进出分口等区域进行网格加密。

为进行网格无关性验证，分别建立5 组不同网格进行计算，计算结果具体见图3。由图3 可以看出，网格数量从233万增长到348万时，电离室通道和后出口的平均温度基本没有变化。基于上述网格无关性验证，最终网格总体数量选定为292万。网格模型如图4所示。

图3 网格无关性验证

图4 网格模型

边界条件的设置参照已有的相关模拟研究[4,6-7]。其中进口边界条件中，空气设置为理想气体，进风温度参照通风手册设定为25 ℃，进风速度根据通风手册规定进风量7 000 m3/h 与进风口面积计算为9.9 m/s。壁面温度的设置上，远离堆坑方向的后壁面温度选定为40 ℃，其余壁面的温度选混凝土壁面最高温度50 ℃，以保证计算的保守性，其中石墨-铅屏蔽壁面额外添加来自堆容器保温层的80 ℃辐射热源。石墨块发热量为1×10-4W，可忽略不计。

空气进口为速度入口，出口为压力出口，指定为静压，湍流模型采用理想k-epsilon模型，近壁面处理采用enhanced wall treatment 壁面函数，压力速度耦合选用Coupled方法。并在竖直方向加入重力作用。

3 数值计算与结果分析

3.1 标准工况下流场与温度场

以进风温度25 ℃、进风速度9.9 m/s为基准参数进行模拟计算，以探究堆外核测量系统内部流畅及温度场特征。图5为石墨-铅屏蔽壁面温度云图，屏蔽板底部因靠近送风口，与冷却空气接触最为直接，呈现出相对低温的区域。冷却空气进而流入上方屏蔽板与箱体间隙，对屏蔽板进行冷却，但随着冷却空气吸热升温，冷却效果逐渐降低，屏蔽板中间区域的温度从下至上呈现升高的趋势。屏蔽板两侧区域因与较大空间的两侧间隙相连，冷却空气不断流出，故散热效果最差，温度最高。

图5 石墨-铅屏蔽壁面温度云图

图6 为系统流线图和后出口处速度矢量图，由图可以看出，冷却空气自下方进风口流入后，一部分从三排共14根堆外电离室通道中流出，其余部分则流经各个面间隙，最终汇流并从后出风口流出。利用Fluent计算各出风口质量流量，算得堆外电离室通道出口质量流量总计1.52 kg/s，后出风口质量流量总计0.74 kg/s。堆外电离室通道通风量约为后出风口通风量的一倍，原因是电离室通道与后出风口相比距离送风口更为接近，且电离室通道入口均匀分布在下间隙流道中，导致通道分流的流量更大。

图6 流线图和后出口处速度矢量图

表1为堆外核测量系统各部位温度值，由表1数据可以看出，在标准工况下，电离室通道出风口、系统后出风口和系统箱体平均温度均小于设计要求规定的45 ℃。电离室出风口平均温度低于后出风口平均温度，原因是电离室通道中空气流量要高于外部流量，且电离室通道流道简单，流程短；箱体外部流道复杂，流程长，冷却空气需流经各个方向间隙，并与高温混凝土墙体直接接触进行换热，故后出风口温度要更高。

表1 堆外核测量系统各部位温度 （单位：℃）

3.2 不同进风口位置计算结果

保持进风温度，进风速度不变，通过改变进风口位置进行数值模拟，以探究进风口位置对通风冷却的影响。进风口位置调整如图7所示。其中，位置1为现有方案中下送风口位置，对比方案将送风口向远离堆坑方向布置，即位置2、位置3处，

图7 底视图及进风口位置调整示意图

对比结果如表2所示，其计算结果表明，随着送风口位置不断远离堆容器方向，堆外电离室通道上出口平均温度呈上升趋势，后出口平均温度呈下降趋势。比较各组数据侧切面速度矢量图可知，流程变长，流道变复杂是电离室通道出口温度上升的原因。系统箱体温度及各出口质量流量并无明显变化。

表2 各组数据结果

3.3 不同进风参数计算结果

保持送风位置不变，通过改变送风温度、送风量以探究通风系统参数变化对堆外和测量系统工作环境的影响。

第一组数据保持进风速度为9.9 m/s，依次将进风温度设置为20 ℃、30 ℃、35 ℃、40 ℃，并记录各组数据的电离室出口平均温度、后出口平均温度及箱体平均温度，计算结果如图8 所示。由图8 可以看出，各项平均温度基本随进风温度线性增长。当送风温度设定为45 ℃时，电离室通道出口平均温度已上升至46 ℃，箱体平均温度上升至47.3 ℃，超过设计要求规定的45 ℃温度限值，因此，堆外核测量系统通风系统在正常工作时应确保送风温度不高于45 ℃。

第二组数据保持进风温度为25 ℃，分别将进风速度设置为6 m/s、8 m/s、12 m/s，并记录各组数据的电离室出口平均温度、后出口平均温度及箱体平均温度，计算结果如图9 所示。从图中可以看出，随着进风速度的增加，各项平均温度均逐步降低，但降低的幅度有限，送风量的影响不及进风温度的影响深。在对堆外核测量系统通风系统进行调节时，应更注重对送风温度的调节。

4 结论

本文利用CFD 软件对不同送风参数与进风位置的情况下对示范快堆堆外核测量系统进行了三维数值模拟，得出以下结论。

（1）计算了标准工况下，即进风温度25 ℃、进风量7 000 m3/h，堆外核测量系统内部各部位温度值，结果表明：标准工况下堆外核测量系统能够在设计要求的温度限制内工作。

（2）进风口位置远离堆容器方向会使得堆外电离室通道出口平均温度上升，方腔后出口平均温度下降，但对核测量系统整体温度变化影响不大，后续反应堆可根据设备设计要求对进风口位置进行相应调整。

（3）与进风量相比，进风温度对核测量系统的温度变化影响更为显著，在设计及运行过程中，应更多地对进风温度做出调整。

本文计算所得的数据，流场与温度场分布特征可为后续一体化快堆相应系统建设提供参考。