蒸汽发生器二回路过冷沸腾的数值模拟.pdfVIP

华 南 理 工 大 学 学 报 （自然科 学版） 第４１卷 第５期 ＪｏｕｒｎａｌｏｆＳｏｕｔｈＣｈｉｎａＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ Ｖｏｌ．４１　Ｎｏ．５ ２０１３年５月 （ＮａｔｕｒａｌＳｃｉｅｎｃｅＥｄｉｔｉｏｎ） Ｍａｙ　２０１３ 文章编号：１０００５６５Ｘ（２０１３）０５００８７０６ 　 蒸汽发生器二回路过冷沸腾的数值模拟 张小英　丁斐　陈佳跃 （华南理工大学 电力学院，广东 广州５１０６４０） 摘　要：为深入研究核电蒸汽发生器二回路侧汽液两相的沸腾传热和流动特性，采用 ＲＰＩ模型对过冷沸腾区域壁面的热流分配进行划分，以此修正 ＣＦＤ程序中的两流体模 型，并利用文献中的实验结果验证了修正后模型的适定性．最后以大亚湾压水堆核电站为 例，采用该模型对蒸汽发生器内二回路预热段单元通道内的过冷沸腾进行计算，获得了通 道内流体空泡份额、速度、温度、热流量分配等的分布情况． 关键词：蒸汽发生器；二回路；过冷沸腾；两流体模型；数值模拟 ＋ 中图分类号：ＴＬ３５３．１３　　 ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１０００５６５Ｘ．２０１３．０５．０１４ 　　蒸汽发生器是压水堆核电站中连接一、二回路 于过冷沸腾区内对壁面热量的划分［１０１２］． 的关键设备，也是核电厂工作运行中故障发生最多 目前ＣＦＤ程序已经被广泛地用于两相领域的研 的设备之一．在蒸汽发生器二回路的预热段易出现 究［１２１４］．为了充分认识蒸汽发生器二次侧预热段低流 过冷沸腾．由于过冷沸腾区域内不稳定的两相流动 速情况下过冷沸腾现象和机理，文中基于ＣＦＤ程序 径向存在较大的温差，导致两相之间出现热力学与 现有的两流体模型，嵌入ＲＰＩ过冷沸腾模型对其修正 动力学不平衡，有可能引发自然循环系统不稳定性． ［１５］ 和简化，并将该模型的模拟结果与Ｂａｒｔｏｌｏｍｅｊ等 的 文献［１］通过实验研究了自然循环系统中过冷沸腾 实验数据对比，验证了ＲＰＩ模型用于低流速过冷沸 流动的不稳定性，文献［２］基于 ＲＥＬＡＰ５程序分析 腾模拟问题的适应性；采用修正后的两流体模型对 了单通道内自然循环流动的不稳定性．然而，当前已 大亚湾９００ＭＷ压水堆核电站蒸汽发生器二回路的 公布的关于蒸汽发生器二次侧流体的热工水力分析 过冷沸腾进行预测，得出蒸汽发生器二回路管束通 ［３６］ 一般不考虑预热段的过冷沸腾 ，或简化为单相 道内过冷沸腾的流动参数特征． ［７］ 换热处理 ．为优化蒸汽发生器运行稳定性，对蒸 汽发生器二回路的过冷沸腾的研究是很有必要的． １　理论模型 Ｓａｈａ模型和Ｌｅｖｙ模型是目前国际上应用最为 ［８９］ １．１　相间热量与质量的交换 广泛的过冷沸腾模型．但孙奇等 通过实验证明， 在低流速过冷沸腾下，Ｓａｈａ模型和Ｌｅｖｙ模型预测结 在过冷沸腾起始点，主流温度仍低于流体的饱 果均明显高于实验值，因此两种模型不能准确预测 和温度，流道壁面开始出现气泡时，液相向气相转 低