棒束间复杂流动传热模拟计算.pdfVIP

　第31卷第2期 原 子 能 科 学 技 术 . 31, . 2 V o l N o 　1997年3月 A tom ic Energy Science and T echno logy M ar. 1997 棒束间复杂流动传热模拟计算 杨志民 ( 中国原子能科学研究院反应堆工程研究设计所, 北京, 102413) 采用多孔体模型对流体在棒束间的复杂流动传热进行了模拟计算, 并将模拟计算结果与实验 数据进行了比较。对该问题的研究结果可应用于快堆安全分析以及工业换热设备内流场和温场的 模拟计算。 关键词　多孔体　S IM PL E 算法　棒束　流动传热 快堆燃料组件内有 61 根燃料棒, 燃料棒直径为 6 mm , 棒间距为 7 mm 。为了给安全分析 提供依据, 对如此稠密的棒束局部堵流时形成的速度和温度分布状况必须作出预测。在快堆中 间热交换器、蒸汽发生器、冷凝器及冷阱等换热设备中, 为了提高换热性能, 管束中设置的挡板 导致出现流体在复杂管束中流动换热问题, 也需预测其内部的速度和温度分布。流体沿棒束流 动过程中遇到挡板的情况是上述情况中的 1 个基本问题, 对该问题的研究可以推广应用到其 它复杂棒束流动问题中去。 采用多孔体模型能较好的解决上述复杂的流动传热计算问题[ 1- 3 ]。多孔体模型把流体流 动空间的障碍物引入被计算单元中, 用多孔度、穿透率及分布阻力和分布热源等参数来考虑障 碍物对流体在该微元中流动的影响, 分别以质量守恒、动量守恒和能量守恒的形式给出。在流 体区域中, 固体的存在一是减小了流动面积, 从而影响到流速及其相关量; 二是改变了能量和 动量传递。前者可以通过引入体积多孔度和表面穿透率来修正, 后者可以通过在能量方程中引 入分布热源和在动量方程中引入分布阻力来处理。准连续区域的计算, 实际上就是用多孔体取 代实际区域中含有的固体进行计算的。计算模型采用均一化方法, 在计算模型中把实际区域中 各处的固体和流体, 按相同的体积多孔度, 相同的表面穿透率, 和在界面上具有相同的能量和 动量传递来处理。 本文应用 S IM PL E 方法对多孔体模型编制程序, 对在棒束中存在挡板的流体流动进行模 拟计算。 1　模拟对象 ( ) 模拟对象为文献[4 ] 中的实验装置 图1 , 实验工质为钠, 实验段由实验段壳体、内六边形 收稿日期:　收到修改稿日期: © 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved. 156 原子能科学技术　　第31卷 管、冷冻密封头和 19 根带绕丝定位的电加热 模拟棒组成。棒直径 6 mm , 棒中心距 7 mm , 绕丝直径 1 mm , 绕丝螺距 144 mm 。中心 24 个子通道堵塞, 挡板厚度 8 mm 。每根元件额 定发热功率 113 kW , 发热段长度 400 mm 。 阻塞板位于模拟元件发热段起始位置处, 距 流道进口 185 , 钠额定流量为 84 3, mm m h 入口温度分别为 250 ℃和 400 ℃。 2　计算模型[2, 3 ] 多孔体模型在质量、动量、能量守恒方程 中考虑了障碍物的影响, 引进了多孔度、穿透 率、分布阻力及分布热源。 连续性方程为 + d iv ( U ) = 0 t 　　动量方程为 ( U ) + d iv ( U U ) =