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

研究生课程考核试卷 （适用于课程论文、提交报告） 科 目： 高等流体力学 教 师： 叶建 姓 名： 杨伟，陆祺 ，张鹏青 学 号： 20131002051，20131002021t ，20131002026t 专 业： 动力工程及工程热物理，核科学与技术, 动力工程及工程热物理 类 别： 学术，学术 ,学术 上课时间： 2013 年9月至2013年 11 月 考 生 成 绩： 卷面成绩 平时成绩 课程综合成绩 阅卷评语： 阅卷教师 (签名) 气冷堆燃料元件表面流场及温度场模拟 1.本文主要研究内容 国内外已经有许多学者利用计算流体力学（CFD）对模块式高温气冷堆球床内的流场、温度场做过大量研究，并取得了不少成果。这方面的研究涉及了燃料元件与壁面或燃料元件之间接触和间隙情况下的传热和流动的特征，孔隙率对温度分布均匀性的影响等等。本文的主要研究内容是数值计算获得堆芯元件表面冷却剂氦气的速度及温度分布规律，并研究不同堆积方式对分布规律的影响。 2.球床内流动与换热数学模型 实际上，在球床对内，燃料元件的分布是随机不规则分布的，有些部分可能是正四棱锥堆积，有些部分是立方体形式堆积；当然，还有可能是其它的堆积方式。本文主要研究的是正四棱锥堆积于立方体堆积时对温度、速度分布的影响。因此，建立的物理模型由立方体和正四棱锥两种，采用的计算软件为Fluent。 2.1 物理模型 在球床对内，球形燃料元件的不同堆积方式对氦气的流动及其温度分布会产生的影响不同。因此，我们需要建立模型模拟以下情况：（1）球床内温度、压力和速度分布规律；（2）堆积方式变化时温度分布和流动的改变。基于这两方面，我们需要建立不同堆积方式的模型。 图2.1 正四棱锥堆积 图2.2和图2.3分别为立方体堆积模型和正四棱锥堆积模型，其中可以明显看到，左边为燃料元件的分布，右边为流体区域。在两个模型当中，燃料球的尺寸均为60mm。图2.1中为正四棱锥堆积方式的平面尺寸图。第一层燃料球球心之间的距离为7.002cm,第一层与第二层燃料球球心的距离为3.3912cm。通过对以上两个简单的模型进行数值计算，可以获得球床内冷却剂氦气的速度及温度分布规律。通过对以上两个模型的计算结果比较和分析不同堆积方式对球床内温度分布及氦气流动的影响。 图2.2 立方体堆积模型 图2.3 正四棱锥堆积模型 图2.4非结构网格 2.2 模型 为了使方程封闭，特别是在高雷诺数下的计算，应采用标准模型，该模型是在湍动能k方程基础之上引入湍流脉动动能的耗散率方程形成的，于1972年提出。其中k为湍动能，为湍流耗散率，分别来自湍流动能方程式和耗散方程式。它们的计算公式如下： (2.4) 式中u为流体流速，I为湍流强度。 (2.5) 由量纲分析可以得粘性系数与湍动能k和耗散率的关系式如下： (2.6) 其中有： (2.7) 式中为无量纲常数，通常取0.09，k为湍动能，l为湍流尺度。 2.3物理模型及求解器设置 在求解控制方程时采用Fluent软件,选用基于压力的求解器，物理模型采用非稳态模型，由于本文三维模型网格单元基本为多面体结构，因此选least-quares cell based这种压力梯度来计算控制方程中的导数项。对于近壁面的处理方式采用标准壁面函数，其中压力-速度解耦一般使用SIMPLE和SIMPLEC，对非定常流动问题或者包含比平均网倾斜度更高的网格则使用PISO算法，基于本模型网格不规则且是非定常流动模型，因此选用PISO算法，其计算结果会更精确一些；但是这种算法相比较起SIMPLEC算法，可能收敛较困难，因此需要根据实际需要调整松弛因子，该因子减小有利于收敛，但是会增加收敛时间，所以需谨慎使用。由于球床内的燃料元件堆积形成的空隙通道不规则，流体压力变化剧烈，PRESTO！主要用于高旋流，压力急剧变化流。因此压力基分离求解器的压力插值算法建议使用PRESTO！。这里可以先采用一阶方程进行初始迭代，再进行二阶迭代计算，若MUSCL算法能够满足收敛条件则亦用MUSCL算法进行更精确的计算。 2.4边界条件设定及其实现 在本模型当中，由于只取局部进行分析，设定氦气为在高压下的局部模型当中为不可压缩流体。根据HTR-10的运行参数知，运行压力为3MP，由于在堆芯取的局部模型压降相对于运行压力来说比较小，相对压力