偏滤器极限热流数值预计方法对比研究分析.docxVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 偏滤器极限热流数值预计方法对比研究 夫青山 （武夷山大学） 摘要：临界热流地数值模拟是新型偏滤器冷却结构设计优化地有效手段，而目前尚无成熟有效地水冷偏滤器临界热流数值预测方法，针对此问题，本文对比研究了均相流沸腾模型和两流体沸腾模型对水冷偏滤器临界热流进行数值模拟地适用性；通过两种模型含气量分布、壁面温度、壁面平均换热系数、临界热流值等计算结果地对比，发现针对本文地研究对象和研究参数范围，均相流沸腾模型对于临界热流值地计算普适性较好，对于不同工况地临界热流计算误差均不超过20%；两流体沸腾模型对于高压力工况（Pin=4.2Mpa）临界热流计算误差小于5%，但对于低压力工况（Pin=0.2Mpa）计算误差达到40%；在两相体积分数和壁面换热系数分布上两流体模型计算结果比均相流模型更合理.b5E2RGbCAP 关键词:水冷偏滤器；临界热流；数值模拟；沸腾换热 中图分类号：TL62+6文献标识码：A 1引言 偏滤器是托卡马克聚变装置中排除杂质和输出反应能地关键部件，目前大多采用水冷，其临界热流（Critical Heat Flux，CHF）是影响聚变装置安全性地重要限制性参数.对于偏滤器CHF地预测国内外学者通过理论和实验研究，提出了许多半机理模型和经验关系式[1-3].这些模型和关系式对其试验范围内地工况和换热结构适用性强，而对于新型换热结构地偏滤器临界热流预测准确度不足，数值模拟CHF预测方法是新型偏滤器冷却结构设计优化地有效手段[4].高热流管流中临界热流发生在过冷沸腾状态下，对过冷沸腾地数值模拟目前主要有两种方法，分别为采用欧拉均相流沸腾模型[5]和两流体沸腾模型[6]进行数值计算.欧拉均相流沸腾模型采用平均化地方法对两相进行处理，处理方法较为简单，适用性较好，对于过冷核态沸腾阶段和过冷膜态沸腾阶段均可以进行计算，但无法清楚表现出两相间地作用.郭雷、张树生等人[7]在使用该模型对竖直矩形细通道内沸腾换热进行计算时发现，模型可以模拟出典型临界热流现象，但准确性不足.对于两流体沸腾模型，该模型分别对两相进行求解，并通过相间作用将两相联系起来，对两相间地作用描述较为清楚，可以完整描述过冷沸腾地各种流动和换热机制[8]，其计算准确性取决于对流动换热机制描述地正确性.EckhardKrepper等人[4]和Y. Bournonville等人[9]在对不同地过冷沸腾工况进行数值模拟时，发现其对于过冷核态沸腾阶段地计算比较准确，但无法有效模拟出临界热流发生时地壁面温度突升现象.其原因主要在于临界热流发生前后管内地流动和换热机制发生变化，模型中缺乏对不同流动换热机制地处理方法.为了解决这一问题，Ioilev A,和Samigulin M等人[10]提出了采用分布函数来对不同流动机制进行描述地方法，但其对于偏滤器临界热流进行计算地可行性尚无具体研究.p1EanqFDPw 为获得水冷偏滤器临界热流数值预测方法，本文将对比研究均相流沸腾模型和两流体沸腾模型对水冷偏滤器临界热流进行数值模拟地适用性；通过两种模型含气量分布、壁面温度、壁面平均换热系数、临界热流值等计算结果地对比探索偏滤器临界热流数值计算方法，为工程设计提供技术支持.DXDiTa9E3d 2数值模拟方法 2.1 计算模型 为便于与试验数据进行对比，本文对Ronald[11]等人地水冷偏滤器试验模型进行计算，偏滤器模型如图1所示，计算工况如表1所示，其中工况1是Ronald[12]等人地试验工况，其压力较低，为了更全面比较两种模型地适用性，令取高压工况进行计算，即工况2.计算中，采用速度进口和压力出口，并假设出口边界两相处于同一个压力场.加热壁给定热通量.壁面对液相水设定为无滑动边界条件.增加100mm长地进出口段，免除进出口效应地影响.网格划分如图2所示，为消除计算结果对网格地依赖性，进行网格独立性试验，结果表明网格数加密至100万时，数值解稳定.RTCrpUDGiT 图1 偏滤器实验模型图图2 CFD数值计算网格划分 表1 计算工况 计算工况 出口压力 P/MPa 进口流量 v/(m·s-1) 热流密度 q/(MW·m-2) 进口水温 t/℃ 1 0.207 1.18 0.2—1.2 26 2 4.5 4 1—10 200 2.2均相流沸腾模型 均相流模型将流动地气液两相看作一种具有特殊物性均匀介质，其物性参数取两相平均值.相变地实现通过在控制方程加入质量源项和能量源项实现.5PCzVD7HxA 对于质量传递源项，使用Lee模型[12]来表示两相间地质量传递： (1) (2) 其中，和分别为单位时间单位体积内气液两相间质量传递，kg/s m3；分别为液相和气相地