自然对流基础上堆顶局部热点分析与设计优化

关键词： 堆顶结构 控制棒驱动机构 流场 温度场 自然对流

在核反应堆中，控制棒驱动机构是反应堆控制和保护系统中重要的伺服机构，控制棒驱动机构是驱动控制棒组件运动、实现准确调节反应堆内反应性变化的重要部件，控制棒驱动机构具有实现反应堆启动、提升功率、保持功率、负荷跟踪、正常停堆和紧急事故停堆等重要的安全功能。

在反应堆的设计分析和设备开发过程中，利用先进的CAE分析技术，既是核电设备国产化的任务和方向，也是技术改进优化，降低成本的重要手段[1]。针对反应堆流场分析研究，常采用计算流体力学（CFD）方法。于浩[2]等通过CFD分析工具对反应堆内的流场分布进行数值模拟，给出反应堆燃料组件入口处流场的速度和流量分配情况。张明[3]等通过三维流场的数值模拟分析，得到了冷却剂在堆内的流动状态。在计算过程中，计算 格的质量决定了计算资源的占用、计算速度以及计算的精度。何培峰[2]等通过HyperMesh软件对堆顶结构进行了高质量六面体 格的划分，并分析了堆顶驱动机构在通风条件下的冷却性能[3]。于浩、余志伟[6][7]等人针对堆顶组件冷却系统进行流场分析，获得了冷却系统的流场分布，为冷却风量的设计计算提供了基础。但堆顶的通风冷却方式，一方面能耗较大，另一方面无法保证绝对的安全可靠，因此基于堆顶结构的自然对流循环冷却方式被提出。邹鹏[8]等通过CFD分析方法，研究了当采用空气自然循环冷却方式时，CRDM群及各线圈的温度分布，探索对CRDM群采用空气自然循环冷却方式的可行性。研究表明，依靠空气的自然对流，可以对CRDM进行有效冷却。

1、模型建立

1.1有限元模型的建立

图1堆顶结构示意图

图2控制棒驱动机构简化模型

利用ANSYSMeshing对堆顶结构进行 格划分，按照计算区域内各类型结构的几何形状特征及其分布来执行相应的块切分操作。如图3、图4所示。其中，对于流体与固体换热区域，流体采取非结构化 格，其余采用结构化 格。流体域 格最大尺寸为6cm，固体域最大 格为2cm，同时对部分典型区域进行 格加密处理。

图3流体域 格划分

图4固体域 格划分

1.2材料物性

(1）结构材料

各种结构材料的物性参数如表1所示。

表1各种结构材料的物性参数

(2）空气物性

如前文所述，在进行空气自然对流流场计算时，空气的密度大小是关键因素。在对流场进行求解时，涉及到稳态和瞬态两种不同计算工况，对流场中流动的空气密度将采用Boussinesq模型进行计算，一方面为了能够使计算更好地收敛，另一方面也确保了计算的准确性。

在动量方程中，为了考虑浮力的影响，则需要添加（ρ-ρ0)g作为浮力项。在Boussinesq模型中可表示为：

式中，β为热膨胀系数。T0为参考温度，ρ0为参考密度。

在计算工况下，空气的比热容为1006.43J/（kg·K），粘度则按照Sutherland公式求解，其表达式为：

式中，μ为粘性，μ0是参考值，单位为kg/（ms),T是流体温度，T0是参考温度，单位为K。S是有效温度，单位是K，也被称为Sutherland常数。

1.3边界条件

(1）控制棒驱动机构（CRDM）承压壳的内部边界条件为第一类边界条件；

(2）流体域与固体域的交界面和不同固体材料之间的交界面均被考虑成可导热的热耦合边界；

(3）对包围CRDM群的空气区域的顶部边界，设为压力边界，其回流温度为环境温度，环境温度为40℃；

(4）考虑围筒内表面与控制棒外表面的热辐射效应，固体内部不考虑热辐射换热的影响；

(5）在运用FLUENT进行模拟分析时，通过对线圈外侧添加虚拟导热厚度的方法，实现了该层壁面边界条件的设置。

2、堆顶局部热点分析及设计优化

2.1堆顶局部热点分析

工况1：全部控制棒驱动机构按照全功率发热进行计算，每台驱动机构按照10kW发热功率计算。

工况2：全部控制棒驱动机构全保持工况进行计算。

图5和图6给出两种稳态工况下，堆顶结构主要构件的温度场分布。两种稳态工况模拟结果对比发现，由于工况1中CRDM功率高于工况2，因此工况1中不论是各线圈的最高温度，还是堆顶各构件的最高温度，包括电缆桥、顶板、电缆托架、抗震板、抗震支撑板、围筒以及电缆通道，均高于工况2。工况1和工况2中围筒温度分布趋势并不相同，由于工况1中线圈功率高于工况2，磁轭围筒的辐射传热量也高于工况2，因此工况1中围筒最高温度分布于围筒进风口周围，工况2中围筒最高温度位于顶部与抗震支撑板接触位置。

图5工况1主要结构温度分布（K)

图6工况2主要结构温度分布（K)

各堆顶结构件局部最高温度见表2所示。

表2部分堆顶构件局部最高热点温度

2.2堆顶局部热点设计优化

图7优化前堆顶结构

图8优化后堆顶结构

为了对比优化方案与原始结构方案的性能，模拟中选取2.1节中工况作为对比工况。

(1）温度场

结构优化后各堆顶结构的温度分布如图9、图10所示。为了进一步分析优化后的方案对整个温度场的影响，分别对不同工况下各堆顶结构的温度进行了分析，如表3所示。对于同一工况下的堆顶结构，由于优化结构后抗震板和抗震支撑板完全封闭，热空气不会如优化结构前直向上流通，因此，有效降低了顶板、电缆桥和电托架的温度，使电缆经过位置温度降低至限值以下。

同时，抗震板和抗震支撑板的封闭，降低了堆顶内的自然循环能力，各线圈温度均有所上升，上升幅度在1-4℃；同时，优化结构后，热空气只能通过出风口流出，围筒上侧空间内热空气发生积聚，使得抗震支撑板、围筒上侧及电缆通道上侧温度升高，但电缆通道温度依然保持在温度限值以下。

图9优化后工况1温度场分布（K)

图10优化后工况2温度场分布（K)

优化改进后各堆顶结构件局部最高温度见表3所示。

表3改进后部分堆顶构件局部最高热点温度

(2）速度场

图11和图12分别展示了不同工况下结构优化前后纵向截面上的速度分布。由图中所示结果可以看出，优化后的速度值要低于相同工况下优化前的值，这也说明优化后的结构降低了堆顶结构局部热点的温度，但也降低了堆顶的自然循环能力，导致堆内其他结构的温度略有升高，但温度的升高值在可接受的范围内。

图11工况1下的纵面速度分布（m/s)

图12工况2下的纵面速度分布（m/s)

图12工况2下的纵面速度分布（m/s)

3、结论

(1）建立了堆顶全尺寸模型，完成了堆顶结构有限元模型的建立于边界条件的设置。

(2）完成了自然对流工况下堆顶局部热点的分析，得出驱动机构在全功率发热时，局部最高热点位于提升线圈处，而全保持工况下，局部最高热点出现在保持线圈处。