高等流体力学课程论文气冷堆燃料元件表面流场及温度场模拟

1、气冷堆燃料元件表面流场及温度场模拟1.本文主要研究内容国内外已经有许多学者利用计算流体力学（CFD对模块式高温气冷堆球床内的流场、 温度场做过大量研究，并取得了不 少成果。这方面的研究涉及了燃料元件与壁面或燃料元件之间接触和间隙情况 下的传热和流动的特征，孔隙率对温度分布均匀性的影响等等。本文的主要研 究内容是数值计算获得堆芯元件表面冷却剂氮气的速度及温度分布规律，并研 究不同堆积方式对分布规律的影响。2.球床内流动与换热数学模型实际上，在球床对内，燃料元件的分布是随机不规则分布的， 有些部分可能 是正四棱锥堆积，有些部分是立方体形式堆积；当然，还有可能是其它的堆积方 式。本文主要研究的是正四

2、棱锥堆积于立方体堆积时对温度、速度分布的影响。 因此，建立的物理模型由立方体和正四棱锥两种，采用的计算软件为Fluent。2.1 物理模型在球床对内，球形燃料元件的不同堆积方式对氮气的流动及其温度分布会产 生的影响不同。因此，我们需要建立模型模拟以下情况：（1）球床内温度、压力 和速度分布规律；（2）堆积方式变化时温度分布和流动的改变。基于这两方面， 我们需要建立不同堆积方式的模型。7班4图2.1正四棱锥堆积图2.2和图2.3分别为立方体堆积模型和正四棱锥堆积模型， 其中可以明显 看到，左边为燃料元件的分布，右边为流体区域。在两个模型当中，燃料球的尺 寸均为60mm图2.1中为正四棱锥堆积方式

3、的平面尺寸图。第一层燃料球球心之间的距离为7.002cm,第一层与第二层燃料球球心的距离为3.3912cm。通过对以上两个简单的模型进行数值计算，可以获得球床内冷却剂氮气的速度及温度分 布规律。通过对以上两个模型的计算结果比较和分析不同堆积方式对球床内温度 分布及氮气流动的影响。图2.2立方体堆积模型图2.3正四棱锥堆积模型IJI;图2.4非结构网格2.2 （一模型为了使方程封闭，特别是在高雷诺数下的计算，应采用卜一标准模型，该 模型是在湍动能k方程基础之上引入湍流脉动动能的耗散率 方程形成的，于 1972年提出。其中k为湍动能，为湍流耗散率，分别来自湍流动能方程式和耗 散方程式。它们的计算公

4、式如下：!. 一 二I 、,工(2.4)式中u为流体流速，I为湍流强度。f|gJ=K(2.5)由量纲分析可以得粘性系数 性与湍动能k和耗散率的关系式如下：一 (2.6)其中有：k.q/(2.7)式中q为无量纲常数，通常取0.09, k为湍动能，i为湍流尺度。2.3 物理模型及求解器设置在求解控制方程时采用Fluent软件，选用基于压力的求解器，物理模型采用 非稳态模型，由于本文三维模型网格单元基本为多面体结构，因此选least-quarescell based这种压力梯度来计算控制方程中的导数项。对于近壁面的处理方式采 用标准壁面函数，其中压力-速度解耦一般使用SIMPLE和SIMPLEC,对

5、非定常 流动问题或者包含比平均网倾斜度更高的网格则使用PISO算法，基于本模型网格不规则且是非定常流动模型，因此选用PISO算法，其计算结果会更精确一些； 但是这种算法相比较起SIMPLEC算法，可能收敛较困难，因此需要根据实际需 要调整松弛因子，该因子减小有利于收敛，但是会增加收敛时间，所以需谨慎使 用。由于球床内的燃料元件堆积形成的空隙通道不规则，流体压力变化剧烈， PRESTO!主要用于高旋流，压力急剧变化流。因此压力基分离求解器的压力插 值算法建议使用PRESTO!。这里可以先采用一阶方程进行初始迭代，再进行二 阶迭代计算，若MUSC算法能够满足收敛条件则亦用MUSC算法进行更精确的计

6、2.4 边界条件设定及其实现在本模型当中，由于只取局部进行分析，设定氮气为在高压下的局部模型当 中为不可压缩流体。根据 HTR-10的运行参数知，运行压力为 3MP由于在堆芯 取的局部模型压降相对于运行压力来说比较小，相对压力取0Pa时有利于分析。进出口温度分别为250整个堆芯入口质量流量为4.3kg/s,根据模型入口截面的 大小可以算出正四棱锥模型入口流量为0.03312kg/s,立方体模型的入口流量为0.02454kg/s 。对于进出口的湍动能k和湍流耗散率,可由上述公式计算得出。在流体流动过程当中，球体壁面温度分布不均匀，不能设定球体表面温度，第一 类边界条件不适用，这里采用第二类边界条

7、件。同时选用无滑移光滑表面。从 HTR-10的运行参数可知，其热功率为10MW燃料球总数量为27000个，尺寸为60mm有以下公式可得出:反应堆总热功率燃料球总表确=32747.93W/m:所以，设定壁面的热流量为32747.93W/m23计算结果与分析本章根据第三章所建立的数学模型进行计算， 得到固定床反应堆堆芯内冷却 剂氮气的速度及其温度分布规律，并讨论了燃料球尺寸、孔隙率等对上述规律的 影响。压力分布卡扰动分布“图3.1速度、温度、压力及扰动分布3.1正四棱锥堆积方式从整个模型的分析可以看出，我们已得到温度、速度甚至压力变化的大概规 律，但是要想更好更具体地分析冷却剂氮气的温度、速度以及

8、压力分布，我们还需要从单一燃料元件着手，研究其表面压力以及温度等的分布情况。 在整个球床 堆芯内，除球床周围边界的燃料元件外，几乎所有燃料元件的温度、速度和压力 的分布规律基本相同。因此，取模型当中的单一燃料元件来分析即可， 在此模型 当中，分别取两种情况下其中心的单一燃料元件对比分析。如图3,2所示，在元件间接触点附近出现温度的极大值点。在燃料元件末端 也出现了较高的温度，这是由于在这个地方冷却剂流速很低， 且有涡流形成，使 该区域内冷却剂出现滞留甚至倒流， 导致对流换热很差，从而使元件的局部温度 很高。Portion cnn图3.2元件表面温度分布如图3.3所示，在Z轴-1.6cm和1.6

9、cm的地方都出现了压力的极小值点,这 些地方出现这种现象的原因是由于这些位置间隙通道较小， 且形成类似扩压管结 构，导致压力迅速降低而流速迅速升高。StAtlCPressure1I UhS Kt*：1soot-oi-Z 绐E：1-SOOti即“m (Cim图3.3元件表面压力分布如图3.4所示，该图为元件表面沿冷却剂流动方向的速度分布图，可发现速度最高点与压力的极小值点相对应。 原因如上所述，均是由于元件间狭小缝隙通 道处形成扩压管结构，在此处压力迅速降低，同时速度升高。Vetocty w眄如断SQUt*OCi &QC*R；OEX鹏MW (cm：图3.4元件表面速度分布因此可以得出结论，在燃料

10、元件堆积形式不变，改变尺寸或者增大流速、热 流密度等的情况下，冷却剂在球床内的流动的规律基本一致，分布形式大致相同, 只是在形式或者程度上发生了变化。3 . 2立方体堆积方式当燃料元件堆积形式发生改变时，必然对对球床燃料元件温度、速度和压力 分布和变化的影响。在不改变燃料元件尺寸的前提下，对其堆积形式做一些改变， 分析这一单一变量对上述分布的影响。 把上述正四棱锥堆积方式变为立方体堆积 形式，具分布如下：如图3.5所示，该图是模型中心的单一元件温度分布图。 在坐标-3cm和3cm 处出现了温度的最高值，而在 0点附近也出现了温度分布的较高点。其中坐标 -3cm和3cm处出现温度极值的原因有两个

11、：一是此处位于上下两层燃料元件接 触点；二是此处流速较低，与正四棱锥模型一样，有涡流形成，该区域内冷却剂 出现滞留甚至倒流，对流换热很差。0点附近出现温度较高点的原因单一，就是 燃料元件之间在此处相互接触，接触点附近流体流动较慢，对流换热较差图3.5元件表面温度分布如图3.6所示，沿流动方向，压力在元件表面的变化与正四棱锥堆积形式的 变化有很大区别。在立方体排布形式当中，元件上半部分压降变化很慢，但到下 半部分时，压力会迅速降低，这与上面所述的原理一样，在燃料元件相互形成的 较小空隙处形成了类似扩压管的通道，导致这个位置压力突然降低。两者有一个 共同点就是这个压力迅速降低的位置都是在与周围元件

12、形成的较小空隙处。但是元件堆积方式不同，这个较小间隙位置就处于元件不同位置。图3.6兀件表面压力分布如图3.7所示，速度在0点附近出现了最大值点，在-3cm和3cm处最小， 几乎为零。与图3.13相对应分析可知，在-3cm和3cm处为燃料元件接触点，也 是燃料元件在流动方向的最低和最高位置，在这些位置会出现涡流甚至倒流，速度最低。而在0点附近周围元件形成了类似于扩压管的结构， 压降最快，速度也 增加最快。500*00上523 0Ct-tt1 30-00PCtCKWl 5图3.7元件表面速度分布4 .小结本文不仅计算得到了燃料元件以正四棱锥方式堆积的球床内温度、速度和压力的分布规律，而且得到不同尺寸和不同元件堆积方式模型下的分布规律；讨论了孔隙率对其的影响，得到的主要