两种气液两相流模型的应用和比较研究.pdf

中闺T程热物理学会 传热传质学 学术会汉论义 编蛩：063241 两种气液两相流模型的应用和比较 杨立新 巴黎明 李星 (JL京交通人学机I【l学院北京100044) 联系电话：010 E—mail：Y!苎倒蔓苎x西丛鲤婪：四 摘要：本义心用MUSIG模型和均一直径模型(t!CFD，F命J：对某溶液堆命架模型堆芯内气液iJ}c幼传热进行 丫数值模拟。卉：MuSIG模型中堆芯内离散相气泡被分为为5纽_；川直径的气泡，用于分析堆芯内气泡的 and and 流动和人小分布。采用Luo Svendsen[2]和PrinceBlanch[3]模型描述小M直径气泡组之间 的破裂和聚合。柏：均一直径模型中，堆：卷内的气体被考虑为M一直住的组分，并f1．小考虑』e破裂0聚 合等现象。奉义汁算得到丫两个模型的模拟结果并且对】E进行J，对比圳究。本文得到的计算结果姓，J÷ MUSIG模型对气液两相流的计算：结果’』白架实验结果吻合史好。 关键词：MUSIG均一直径两楣流；!i：f液堆 1．引言 医心同位索生产溶液堆是一种利川燃料溶液核反应制造医疗使用的放射性同位素 的反应堆。7反应堆内会产生人量的热餐和气泡，形成气液两相流动传热，气泡产生速率 和功率分布相互影响，热量由冷却盘管带出堆芯。为消楚地了解气泡迁徙引起的流动传 热变化，及赢管换热器的换热效率，文献【43对该现象进褥了无核反应的热l：水力台架 模型实验。该实验模型对整个反应堆模型进行了如’F简化，如图l所示： ①根据溶液堆内功率分布的计算，使用电加热及相应的管路分布模拟功率分布 ②在水力模型的底部使用气孔板集中进气模拟气泡的产生和流动 ③反应燃料溶液使用水代替，气体使用氮气代替 ④使川1：1比例台架模拟溶液堆原型 实验结果验证了这种盘管换热器在有气体和无气体通入时的总体换热特性，并测量 得剑堆芯内局部温度分布。但是实验结果不能够给出反应堆内气液流动和换热的全部信 息，冈此本文虑刚两种模型分别对该溶液堆台架模型堆芯内气液流动传热进行了数值分 析。本文在数值模拟中，分别采用了均一离散相和把应用群体平衡理论建立起来的MUSIG 模型两种模型，并且对模拟计算结果进行了对比。 2．数值模型 2．1MUSI G模型 Size MUSIG模7鸭全称为MultipleGroup，适合处理离散相具有较人尺度变化的多 项流问题。MuSIG模型将．1卜连续相按尺度分为N纽，并假没所有竹连续相具有一致的速 度，所以简化了动鼙方程数量，可以只求解一个对所有1卜连续相分组都适Ⅲ的动釜方科， 同时保留各分组的近续性方样米描述北连续相的尺度分布。 2．1．1各组连续方程 对1-’1I连续相的第i个5L#t，连续方程为： 昙(剐丹(鹏)一 (1 1162 这里，S；足由Ji破裂、聚合等造成的第i组质苗变化率。显然，所有组之间的质鼙传递 为零，而上L所有组之间的体积分数之和等丁．离散相的体积分数。． 2．1．2群体平衡方程 群体平衡方程考虑了1卜连续相的破裂利聚合效应，是计算+I卜迕续相的尺度分斫i的一个{良 好的方法。群体平衡方程的通j}J表达式为： 昙_+V如吩)=吼一％+最．一球 (2) 这里，_是单伉体积内第i组的数量，吃和反，分别是由丁．破裂和聚合造成的产生率， ∥疗和M’分另0是相应的消火率。 2．1．3破裂模型 在各向同性湍流和概率分布的基础上。该模型假设破裂都是一分为二的， 2．1．4聚合模型 一般假设，两个微粒的聚合通过三个步骤米完成。首先，微粒对其中间的液相进行冲击， 然后液相从该位置流出直到液膜