1. 引言
钠冷快堆在完全丧失电源的情况下,紧急停堆后要依靠一回路和事故余热排放系统进行冷却。其中经过堆芯建立的一回路自然循环是堆芯燃料冷却的唯一方式。由于钠冷快堆组件采用绕丝缠绕棒束的固定方式,使得冷却剂通道变得狭窄紧凑。这种布置方式加强了冷却剂的横向交混,增强了燃料棒与冷却剂之间的换热,但是这种布置方式增加了流道的阻力。燃料组件棒束的阻力占了整个循环回路的90%以上,所以对燃料棒束阻力的研究对于自然循环过程中换热的关键。很多的学者进行了一系列的试验研究,提出了带绕丝燃料组件的摩擦阻力系数关联式;如Cheng与Todreas、Rehme、Engel根据实验得出计算带绕丝棒束摩擦阻力系数的经验关联式 [1] [2] [3]。早期由于计算机技术不发达,对快堆的研究只能用实验来完成,这使得完成一项工作需要耗费很大的成本。随着计算机技术的发展,可以使用计算流体力学工具来完成对快堆的研究。如Kurt完整地演示了使用CFD软件模拟快堆燃料棒束的全过程 [4];Govindha等利用CFD软件对快堆带绕丝燃料组件棒束进行了模拟研究 [5]。但是大多的CFD模拟研究都是在高雷诺数条件下进行的,很少有学者在低雷诺数条件下进行研究 [6] [7] [8]。近些年来我国火力发电总比重持续下降,非化石能源在电力生产能源中的占比在持续提升,能源结构不断得到优化。但是风力发电、太阳能和水力发电对自然资源都有着较高的要求,导致其电力生产环节稳定性较差,而核电因为其用核资源作为燃料供给,运行功率和电力供给水平稳定具有较强优势 [9] [10] [11]。低雷诺数的流动工况对于快堆运行安全也是很重要的,层流和过渡流在快堆停堆后会出现。本文利用CFD软件STAR-CCM +对快堆燃料棒束进行低雷诺数下的水力特性分析。
2. 几何建模及网格划分
2.1. 几何模型
在钠冷快堆燃料组件中棒束是按照三角形格栅的形式排列,利用金属绕丝固定,组装在一个六角形套管中。复杂的绕丝结构给建模带来很大的困难,主要体现在绕丝与其接触的棒束是线接触,且与棒束相切,这就很难做出适合这种结构的网格。通常情况下都要对绕丝结构作简化,使线接触变成面接触,简化前后的结构如图1(a)所示。本文选取37棒燃料组件作为研究对象,单棒总长度为2700 mm,燃料棒直径7 mm、节距8.6 mm;绕丝直径1.5 mm、绕丝螺距H为150 mm;由于计算资源的限制,本次CFD模拟选择37棒的加热段作为研究对象,长度为1000 mm,图1(b)~(d)为37棒带绕丝燃料组件,图1(d)为计算流体体积模型。
2.2. 网格划分
STAR-CCM+生成的网格结构被归类为自动非结构化网格,生成的网格如图2(a)所示,从图可以看到在绕丝靠近组件壁处的网格密集;图2(b)显示的是绕丝与棒束相交之处的网格,获得这样高质量的网格是很有挑战的。以进出口压降作为选择网格的衡量指标,选择网格数在3407万与13,110万之间作一个网格无关性研究,从图3可以发现在网格数为9636万之后的压降变化很小,考虑到现有的有限计算资源,所以选择网格数为9636万的燃料组件作为研究对象来进行计算研究。
(a)
(b)
(c)
(d)
Figure 1. Schematic diagram of the calculation model
图1. 计算模型的示意图
(a)
(b)
Figure 2. Schematic diagram of the grid
图2. 网格的示意图
2.3. 边界条件和计算模型设置
本次模拟采用液态钠为冷却剂,在计算中使用表1中所示的钠物性参数。
![](Images/Table_Tmp.jpg)
Table 1. Physical properties of liquid sodium (T in K)
表1. 液态钠物性参数(T为开尔文温度)
在计算等温流动时,冷却剂为液态钠,入口温度为280℃,不同雷诺数对应不同的流量;出口压力设置为0 Pa。计算时采用层流和湍流模型分别计算,计算结果作比较。自然循环情况下,组件内的冷却剂流动为层流,但是在绕丝存在的情况下,会形成局部的湍流或者过渡流。所以选择层流和湍流模型分别进行计算。
3. 计算结果及分析
计算结果与以公开的Cheng经验公式、Rehme经验公式和Engel经验公式作对比,从而得到与实际情况比较接近的模型。
简化的Cheng与Todreas经验关系式如下:
层流:
(1)
湍流:
(2)
过渡流动:
(3)
其中:
式中:
分别为层流和湍流关联式的系数;
分别为层流临界雷诺数和湍流临界雷诺数;
分别为棒束直径和绕丝直径,单位:m;
为棒-束栅格距。
Rehme经验关系式如下:
(4)
其中:
式中:
分别为棒束直径和绕丝直径,单位:m;
为棒束栅格距,单位:m;
为燃料棒数;
为湿润周长,单位:m;H为螺距,单位:m。
Engel经验关联式如下:
层流:
(5)
湍流:
(6)
过渡流:
(7)
其中:
3.1. 等温流动计算结果分析
37棒带绕丝棒束组件摩擦阻力系数计算结果示于图4。从图4中可以看出,利用层流模型计算的结果和利用k-ω湍流模型计算的结果非常的接近,两者的差距都小于0.05%以内。在Re小于600时,模拟的结果与Cheng经验关系式的结果比较接近,偏差小于13.21%;当Re在600~3000时,模拟结果与Rehme关系式的计算结果较接近,偏差小于13%;而Engel关系式计算结果整体与模拟的结果偏差较大。所以在计算类似CFR600燃料组件的带绕丝燃料棒束组件时,使用Rehme经验关系式来进行低Re条件下的阻力系数计算准确度较高。图5显示了轴向摩擦阻力系数的分布情况。从图5可以看出,摩擦阻力系数在入口区域的值要远远大于发展段的值,并且在距入口150 mm处时的摩擦阻力系数已经与充分发展段的值相差无几,说明在一个螺距之后摩擦阻力系数回复正常值,入口段效应消除;在出口处,摩擦阻力系数有明显的下降,但是变化的数值不是很大,是微小的变化,所以出口段对于摩擦阻力系数的影响是小的。同时也可以看出,在Re增大的过程中,入口段摩擦阻力系数的差距也越来越小。
![](//html.hanspub.org/file/3-1270608x46_hanspub.png)
Figure 4. Comparison of the CFD calculation results of the friction resistance coefficient of the 37-bar-tape-wire-wound assembly with various relational expressions
图4. 37棒带绕丝组件摩擦阻力系数CFD计算结果与各关系式的比较
![](//html.hanspub.org/file/3-1270608x47_hanspub.png)
Figure 5. The distribution of the axial friction coefficient of the 37-rod fuel assembly with wound wires
图5. 带绕丝37棒燃料组件轴向摩擦系数分布
3.2. 棒束内部流场分析
以Re = 435为例进行内部流场分析。图6显示了在加热条件下绕丝棒束在前两个螺距范围内,每隔1/6螺距的距离取得截面的速度分布云图。观察图中速度的分布情况,发现在第一段绕丝距离内,速度的峰值在0.0793~0.0866 m/s之间;而在第二段绕丝距离内,峰值速度在0.0796~0.0818 m/s之间,速度的波动变化较小。在75 mm之后的速度分布图,整体的分布都可以认为是前一个速度分布图顺时针旋转30˚得到的,可以认为入口距离是要小于一个螺距的距离,即从150 mm处开始冷却剂进入充分发展段。从图中可以很明显的看出,出现最大速度的位置与绕丝和六角形壁面相对位置有关,绕丝与壁面的距离较近,则速度较大,相反速度就较小。边子通道的流速要远高于中间子通道的速度值,这主要由于边子通道比中间子通道的流通面积较大,流动阻力较小,流体流入组件后因个子通道的流动阻力不同产生再流
量分配,进而导致边子通道的流体流速高于中间子通道。此外,在绕丝的横向扰流的作用下,组件内的流体流动将会随着绕丝的缠绕方向呈周期性螺旋上升,靠近绕丝的子通道流动会被增强,因此六个边子通道中离螺旋绕丝最近的子通道流速最大。从图中可以看出每一截面的整体速度分布都是由低到高逐渐递增的,方向是跟随绕丝缠绕方向。峰值速度出现在外围的7根棒周围,而相对侧的棒则出现低流速。
4. 结论
本研究主要进行了钠冷快堆单盒组件的数值模拟,对钠冷快堆在自然循环条件下的流动特性有了初步的认识,得出如下结果:
1) 通过不同雷诺数计算得出的摩擦阻力系数与现存的经验关系式进行对比发现,在Re小于600时,模拟的结果与Cheng经验关系式的结果比较接近,偏差小于13.21%;当Re在600~3000时,模拟结果与Rehme关系式的计算结果较接近,偏差小于13%;而Engel关系式计算结果整体与模拟的结果偏差较大。所以在计算类似CFR600燃料组件的带绕丝燃料棒束组件时,使用Rehme经验关系式来进行低Re条件下的阻力系数计算准确度较高。
2) 通过分析燃料组件轴向摩擦阻力系数分布情况发现,摩擦阻力系数在入口区域的值要远远大于发展段,在150 mm之后也就是第一个螺距之后摩擦阻力系数恢复了正常值;出口段对于摩擦阻力系数的影响较小。
3) 通过分析不同位置的速度分布云图发现,入口距离是要小于一个螺距的距离;出现最大速度的位置与绕丝和六角形壁面相对位置有关,绕丝与壁面的距离较近,则速度较大,相反速度就较小。