溶液堆MUSIG气液两相流动换热数值模拟.pdfVIP

溶液堆MUSIG气液两相流动换热数值模拟 杨立新，巴黎明，李 星 (北京交通大学机电学院，北京100044 摘要：应用欧拉两相流最新CFD模型技术对某溶液堆台架模型堆芯内气液流动传热进行了数值分析。堆芯内离散 相气泡被考虑为不同直径的组分，摒弃了传统均一离散相的假设(如假设所有气泡为同一直径和形状)，把应用群 泡的流动和大小分布。本文将气体离散相采用分为5组不同直径的气泡直径。分别采用LuoandSvendsenc。和P血ce and Blan洲3慎型描述不同直径气泡组之间的破裂和聚合。模拟结果与台架实验结果吻台较好。 关键词：MUSlG；两相流；溶液堆 l引言 医用同位素生产溶液堆是一种利用燃料溶液核反应制造医疗使用的放射性同位素的反应堆。堆芯 内的核反应会产生大量的热量，通过设计台理的盘管式换热器可将热量带走。由于堆芯内的核反应同 时还会产生氢气和氧气，在堆芯内形成气液两项流动传热，使得堆内冷却盘管式换热器的热工水力设 计更加困难。气泡的产生与反应功率分布相关，反应功率与局部温度相关，局部温度又受气液两相流 动影响，构成了流动，传热与核反应复杂的耦台过程。为清楚地了解气泡迁徙引起的流动传热变化， 及对盘管换热器的换热效率的影响，文献[4]对该现象进行了无核反应的热工水力台架模型实验。图l 给出深液堆热工水力模型结构图。该实验模型对整个反应堆模型进行了如下简化： (1)根据溶液堆内功率分布的计算，使用电加热及相应的管路分布模拟功率分布； (2)在水力模型的底部使用气孔板集中进气模拟气泡的产生和流动； (3)使用水代替反应燃料溶液，使用氮气代替气体； (4)使用1：1比例台架模拟溶液堆原型。 奏◆ 侧视图 图1 深液堆热工水力模型结构图 实验结果验证了这种盘管换热器在有气体和无气体通人时的总体换热特性，并测量得到堆芯内局 10l 部温度分布。然而实验结果对进一步设计和改进换热器结构很难提供有效的有理论依据和方法。而且 实际堆运行工况下的流动换热实验很难实现，所以在这种情况下，考虑了实际几何模型、实际流体物 本文应用欧拉两相流最新CFD模型技术对该溶液堆台架模型堆芯内气液流动传热进行了数值分 析。前人在这类气液两项流的数值模拟中，多采取均一离散相的假设，即认为气泡直径和形状不变。 未考虑气泡在迁徙过程中由于聚合和破裂引起的直径和形状变化。而准确模拟气泡的大小分布即意味 着能准确确定两相问传热传质的面积，准确确定气泡在液体中的曳力模型。因此，均一离散相的假设 模拟得到的数值结果往往与实验值误差很大，工程应用的价值低。本文把应用群体平衡理论建立起来 的MUSIG模型考虑到CFD分析中。MUSIG模型用于分析堆芯内气泡的流动和大小分布。分别采用 Luoand Svendsen[21和PrinceandBlanch【31模型描述不同直径气泡组之间的破裂和聚合。模拟结果与台 架实验结果吻台较好。 2数值模型 2．1欧拉多相流模型 本文采用欧拉多相流模型模拟溶液堆堆芯内气液流动传热。欧拉多相流模型中，气液两相均为独 立的连续相共存于堆芯内，流动区域各流体占有份额通过体积分数表示。各相均具有自己的速度场、 温度场和物性参数。通过求解质量，动量和能量守恒方程得到各相的速度温度和相分数分布。本文采 用的欧拉多相流守恒方程如下： 连续方程： a 素(御。)+V(佩Ud)=o (1) 式中：ra代表a相体积分数，m是a密度，％是a相速度，体积分数的总和为l，eP： ∑h=1 (2) 动量方程： 匀。 (3) 素(，巩以)+Vtk(几以。以)】-％Vm+V-【7札(V以+(V乩)7)】+z