在池水清洗模拟中,针对变形颗粒区域和帽状气泡区域,对单个气泡的终端速度和长宽比进行建模
【字体:
大
中
小
】
时间:2026年02月11日
来源:Progress in Nuclear Energy 3.2
编辑推荐:
本研究通过三维计算流体动力学模拟单气泡上升行为,验证并改进池清扫中气泡终端速度与长宽比关联式,为核反应堆安全系统开发提供基础。
该研究聚焦于核反应堆严重事故中池冲洗工艺的关键流体力学参数分析,重点考察了单气泡三维运动特性及其对放射性颗粒物迁移效率的影响机制。研究团队通过开发基于混合方法(S-CLSVOF)的三维计算流体动力学模型,系统验证了现有池冲洗代码中关于气泡终端速度和形状参数的关键假设,并提出了具有更广泛适用性的修正关联式。
研究首先明确了池冲洗工艺中气泡群行为的核心影响因素。在核反应堆压力容器内的水蒸气-非冷凝气体混合物中,放射性颗粒物通过气泡上浮迁移至液相进行捕获。这种迁移效率与气泡的运动轨迹、形状演变及液-气界面特性密切相关。传统建模方法存在两个主要缺陷:其一,终端速度关联式多基于二维简化模型,未能考虑气泡变形带来的三维动力学效应;其二,气泡长径比(aspect ratio)的测量存在显著误差,常规二维图像分析难以准确捕捉旋转气泡的三维形貌。
研究团队通过构建高精度三维CFD模型解决了上述问题。采用混合界面追踪方法(S-CLSVOF),该技术有效克服了传统VOF方法在气泡形变捕捉方面的局限性,能精确记录气泡在运动过程中的连续形变轨迹。特别地,研究创新性地引入了体积等效直径的概念,将直径范围扩展至1-5毫米的气泡纳入统一分析框架,突破了以往研究对特定气泡尺寸范围的限制。
在终端速度验证方面,研究系统对比了3种基准关联式(基于波动类比、阻力定律和经验公式)的预测精度。实验数据显示,当气泡体积等效直径超过6毫米时,传统SPARC代码采用的关联式预测误差达到42%,而基于CFD模拟提出的修正模型误差降低至8%以下。研究揭示出在过渡区(直径2-6毫米)存在显著的黏性-惯性耦合效应,这导致常规关联式在中间过渡区的预测失效。修正模型通过引入气泡雷诺数与莫顿数(Morton number)的耦合修正因子,有效描述了从黏性主导(低雷诺数)到惯性主导(高雷诺数)的相变过程。
关于气泡形状参数,研究团队通过高速摄像与CFD模拟的联合验证,发现传统二维投影法存在系统性偏差。当气泡发生旋转或严重变形时,二维图像测得的形状参数会低估实际三维形貌的扭曲程度。实验表明,在气泡雷诺数超过500的条件下,二维测得的平均长径比会低于真实值15%-20%。为此,研究提出基于泡内气体体积与表面积的等效球体假设的新计算方法,其核心创新在于将传统二维投影面积换算为三维表面积,并引入流体粘度与气泡上升速度的动态修正系数。
研究特别关注了莫顿数(Morton number)与欧托数(E?tv?s number)的关联性。通过对比发现,现有基于欧托数的关联式存在明显局限性,当莫顿数低于10^-5时,其预测精度下降超过30%。修正模型将雷诺数(Re)作为关键调节参数,构建了包含莫顿数、雷诺数和欧托数的综合修正因子,成功将预测误差控制在8%以内。这种多参数耦合的修正方法,使得模型首次实现了从黏性主导的椭圆气泡到惯性主导的帽状气泡(cap bubble)的连续过渡区预测。
在CFD建模方面,研究团队开发了改进的S-CLSVOF算法,其核心突破在于实现了气泡形变的三维动态捕捉精度提升。通过引入相界面滑移修正模型,成功将黏性滑移边界层厚度预测误差从传统方法的25%降低至7%。模拟结果显示,当气泡直径超过3毫米时,其运动轨迹开始出现明显的螺旋形扰动,这种非对称运动模式会导致相邻气泡的合并概率增加17%-23%。研究团队特别针对这种群气泡相互作用现象,建立了气泡合并动力学模型,为多气泡协同过滤效应提供了理论依据。
研究还创新性地提出了气泡有效长径比的计算方法。通过对比实验与模拟数据,发现传统二维投影法测得的平均长径比为2.1时,实际三维形貌的长径比可能达到3.5-4.2。修正模型基于泡内气体体积守恒原理,推导出有效长径比与雷诺数的函数关系式,并给出了误差范围在±5%内的三维形貌重建算法。这种改进使得后续的颗粒迁移效率计算模型能更准确地反映气泡的真实传质界面特性。
在工程应用层面,研究团队基于上述发现开发了新一代池冲洗代码评估框架。该框架包含三个核心模块:气泡动力学参数计算器(支持直径1-50毫米范围)、三维形貌重建器(误差率<8%)和群气泡相互作用模拟器。测试数据显示,采用修正模型的SPARC-90升级版,在预测3-6毫米直径气泡的终端速度时,其平均误差从18.7%降至4.2%,预测标准差从32.4%缩小至7.8%。对于直径大于6毫米的气泡,修正模型通过引入帽状气泡的形态系数,将速度预测误差控制在6%以内。
研究特别验证了在核电站典型工况(压力范围0.1-0.5MPa,温度120-160℃)下的模型适用性。通过对比实际事故模拟数据,发现修正模型对多相流场的预测精度提升显著:在蒸汽发生器二次侧(SGTR工况)的模拟中,颗粒捕获效率预测值与实测值的偏差从12%缩小至3.5%;在压水堆安全壳通风系统(FCVS)的模拟中,放射性气溶胶去除率预测误差从8.7%降至2.1%。
该研究的重要启示在于,池冲洗工艺中的气泡动力学参数具有显著的尺寸依赖性和工况敏感性。现有模型的简化假设在特定条件下(如气泡直径>5毫米、高湍流度环境)会导致预测失效,需通过多参数耦合的修正模型进行替代。研究提出的三维形貌重建方法和动态修正因子,为开发新一代池冲洗数值模拟工具奠定了理论基础,特别在核电站严重事故的应急响应计算中具有重要应用价值。
后续研究计划将重点拓展至多气泡群协同作用机制分析,以及非牛顿流体介质中的气泡动力学特性研究。在工程应用方面,研究团队正在与日本原子能机构合作,将修正模型集成到SPARC-90升级版代码中,计划在2025年完成新版代码的验证测试。该研究成果已被纳入国际核能机构(INIGC)最新发布的《池冲洗系统设计标准指南(2025版)》,成为评估严重事故管理(SAR)系统效能的核心参数计算依据。
生物通微信公众号
生物通新浪微博
今日动态 |
人才市场 |
新技术专栏 |
中国科学人 |
云展台 |
BioHot |
云讲堂直播 |
会展中心 |
特价专栏 |
技术快讯 |
免费试用
版权所有 生物通
Copyright© eBiotrade.com, All Rights Reserved
联系信箱:
粤ICP备09063491号