反应堆堆芯熔化机理研讨.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：056045 反应堆堆芯熔化机理研究 尚智1，陈硕2，赵钧2 1．上海交通大学机械与动力工程学院棱科学与系统工程系，上海200030 2．上海变通大学机械与动力工程学脘航空航天系，上海2删0 精墨：．用分子动力学方法数值模拟了在反应堆发生堆苍熔化严重事故时，熔亿后的熔融金属颗粒从燃 料组件上脱落下，并落到下管板上，随后会继续变形和延展。通过对这个过程的物理机理的计算与描 述，计算出熔融金属颗粒颗粒的形状及其变化过程，从而揭示、研究堆芯融化严重事故下反应堆堆芯 熔化后的力学机理。 关键词：分子动力学，数值模拟，反应堆．严重事故，堆芯熔化 1．引言 无论是何种类型的核反应堆．都要求在安全的环境下运行，并且避免发生任何的运 行事故。为了反应堆的安全运行，世界各嗣的科学家、工程师以及实际运行反应堆的人 员，都进行了大量的研究和预测，保证反应堆的事故发生概率低于10一。由于实际情况 错综复杂，即使在这样低的概率要求下，还是有各种事故相继发生，比如美国三里岛核 泄漏和前苏联切尔诺贝利核电厂事故“1。这些核事故可以归类为核严重事故，针对压 水反应堆有一种堆芯熔化严重事故是人们非常关心的，随着我国核电事业的起步与扩 大，压水堆无疑将成为重点发展堆型，因此对于堆芯熔化严重事故的研究也相应地得到 了重视。 一般研究堆芯熔化严莺事故的学者都将注意力放在发生事故以后的测算和处理上 ”1，或者金属腐蚀“1等事故后情况上。但是很少有人关注在发生堆芯熔化严重事故时， 熔化的金属的运动和变形，以及熔化的金属颗粒落到堆芯下面管板后其运动(流动)是 怎么样的。其中原因除了研究范围超过热工水力分析的范围以外，还有对于这个复杂过 程的研咒手段不完善。实际上，这个过程对于发生堆芯熔化后堆芯内部的瞬态变化以及 事后的处理都是极其重要的。 熔化的金属颗粒落到堆芯下面管板上的流动过程和熔融金属颗粒下落到固体壁面 上的流动过程相类似。实际上，以往的针对熔融金属颗粒下落到固体壁面后流动情形的 研究工作主要集中在实验方面“1，而且还包括了部分浸润和完全浸润等不同的情形。但 是，现有的实验数据与几个互不相同的理论均相符，而且目前还不能够证实哪个理论真 正有效“1。问题的难点在于几乎不可能仅仅改变某一个参数而保持别的参数不变，例如， 347 如果改变了流体种类，那么表面张力以及固液之间的相互作用也因此而改变。因此，十 分有必要发展一种数值模拟技术来详细描述液滴(熔融的金属颗粒)在固体壁面上的延 展，其中，分子动力学不失为一个可行的方法。 本文将用分子动力学(M01ecul”DyrI∞ics)方法来研究堆芯熔化以后，熔化的金属 颗粒下落及其变形过程，对熔化金属颗粒落到下面管板以后的变形和流动过程给出微观 的数值模拟，由此解释熔化金属颗粒的变形及流动的力学机理，并探索为今后的核事故 势能来模拟液滴(熔化的金属颗粒)。以往的研究表明…，影响熔融金属颗粒在固体壁面 延展的因素包括温度、熔融金属颗粒和固体壁面之间的相互作用、固体壁面本身的表面 粗糙情况以及熔融金属颗粒下落的高度和初始速度等等。 2．分子动力学数值方法 在本文中，固体壁面和液滴(熔化的金属颗粒)均由简单的球形分子组成，这些分 子的运动遵循牛顿第二运动定律。一般，分子j的速度i和动量“可表达为 ‘；詈，D，=∑毛+E (1) m 鲥。 其中口为流体分子的质量，如为分子J作用在分子j上的力，L为外场作用力。固体壁 面由两层壁面分子组成(如图1所示)，这些壁面分子处于平面Fc格点上。壁面分子的 质量假定为与流体(熔化的金属颗粒)分子的质量相同。每个壁面分子由非线性弹簧将 其锚在各自的格点上。壁面分子的运动方程与方程(1)相似，只是将外力替换为弹簧力 即可。 图l壁面和熔化的金属颗粒分