核反应堆严重事故中燃料与冷却剂相互作用粗混合过程的数值模型验证研究.pdfVIP

核反应堆严重事故中燃料与冷却剂相互作 用粗混合过程的数值模型验证 曹学武1，飞田吉春2 (1上海交通大学核科学与系统工程系，上海市华山路1954号，200030 2日本核燃料循环开发研究所，大洗工学中一tl,，日本) 摘要：熔融燃料与冷却剂相互作用的粗混合过程是核反应堆严重事故机理研究中很重要的物理现象。在粗 混合过程中由于冷却剂的局部蒸发，高温燃料粒子或熔融燃料液滴被一层蒸汽膜包裹，形成了特殊的复合结 构。但是在目前的数值模拟研究中，高温燃料粒子或熔融燃料液滴在冷却剂中的阻力系数大多数是采用Ishii 的模型，并非针对此特殊复合结构开发的。作者基于此特殊复合结构开发出了～个新的阻力系数模型。本文 将通过应用Ishii的模型和新开发的模型分别模拟粗混合实验QUEOS以比较和验证新开发的阻力系数模型。 结果表明新开发的模型较大地改善了粗混合阶段的数值模拟模型中的阻力系数与实验的差异。 关键词：数值模拟；多相流；燃料与冷却剂相互作用；阻力系数 中图分类号：TL364 1 引言 在轻水堆(LWRs)中的堆芯熔化事故CoreMelting 堆(LMFBRs)中的堆芯解体事故CoreDisruptive 象是高温难挥发的熔融燃料与低温易挥发的冷却剂之间的相互作用。从基本过程来看，此相 互作用可分为高温熔融物团块被释放到冷却剂中的破裂过程(Jet Breakup)、与冷却剂的粗混合 过程(Coarse 融物液滴与冷却剂形成一个混合区域，混合区域的状态与范围决定后续爆炸性膨胀的初始条 件与参与可能转换成破坏性机械能的总的热能，近年来各国学者用实验和数值模拟方法对此 过程进行了研究”’。在数值模拟方面，目前国际上正在开发一些基于多相多成分的大型的数 值分析程序”“’51，例如，日本核燃料循环开发研究所开发的SIMMER-111”“1，它是两维、 三速、多相、多成分耦合中子动态模型的核反应堆安全分析计算程序，其中的一个功能就是 模拟燃料与冷却剂相互作用的粗混合过程。该程序模拟粗混合实验”’81结果表明高温金属／氧 化物粒子或熔融物液滴(简称高温粒子)与冷却剂相互作用的混合区域与实验结果有较大的 误差，而且误差随着粒子的温度增加而增大“1。 由于在数值计算中影响混合区域的主要是高温粒子在冷却剂中的阻力，其模型的准确是 基于单个冷粒子在流体中运动开发的…1，大多数程序都是利用此模型。但是在燃料与冷却剂 相互作用时的粗混合阶段，由于高温粒子周围的冷却剂沸腾蒸发，形成一个在冷却剂液体中 被蒸汽膜包裹的特殊复合结构。而且蒸汽不断地从汽液界面上蒸发和在汽膜内流动，从而对 高温粒子产生一个作用力，以至于影响其阻力特性。而在目前开发的数值模型中所用的高温 粒子在冷却剂中的阻力系数模型并没有考虑此特殊复合结构的影响。 本文作者在观察相关实验研究的基础上，开发了一个基于此种特殊复合结构的高温粒子 或熔融物液滴在冷却剂中运动的阻力系数模型…1。本文将通过应用Ishii的阻力系数模型和 323 翌要垄紫学塑系篓模型分别模拟粗混合实验QuEos以比较和验证新开发的阻力系数模型。并 塑誓茎苎塑变塑翌子或熔融物液滴在冷却剂中运动的阻力系数模型嵌入到SIMMER-III程每 中以完善该程序模拟粗混合过程的模型。 2数值模型 燃料与冷却剂相互作用粗混合过程数值模拟的动量方程为m： 一∑k．p吒，．见+日(-r钾．丘_J(1) 口 其中公式右边第四项的动量交换系数为： qq ～·=厶-+曰qq．眈．一咀稍 、q 目}。 (2) y．一矿 对于在冷却剂中运动的高温粒子或熔融物液滴 4 。4是冷却剂与粒子或液滴之间的相对速 度，上述公式中的系数为： 爿*，=三24一丝rd 和 B。．=专咖~CD