池式铅冷快堆SGTR事故多组分多相流动过程数值模拟研究

Numerical Investigation on Multi-phase Flow of Pool-type Lead-cooled Fast Reactor under SGTR Accident

  • 摘要: 本研究使用欧拉坐标下的多组分多相分析程序ACENA,首先介绍了ACENA程序的基本数学物理模型,然后通过铅铋-氮气两相流动实验HESTIA-2、KYLIN-Ⅱ-S铅铋-水相互作用实验和点堆中子动力学方程解析解,对程序热工水力模块和中子动力学模块进行了验证计算,在此基础上,针对欧洲先进铅冷示范堆ALFRED的设计方案分别开展了热态满功率稳态校核计算和假想无保护蒸汽发生器传热管破裂(SGTR)事故瞬态模拟,重点关注了SGTR事故后铅池内多相流动过程以及包壳最高温度、燃料最高温度、堆芯相对功率以及主容器压力等参数的演变,并分析了断管数量、铅冷却剂循环路径以及所采用的机理模型等影响因素对ACENA程序计算结果的影响。本文研究结果表明,Ishii-Chawla-Suzuki相间曳力系数模型结合Ishii等提出的相间界面面积浓度输运模型能够较好地模拟圆形/环形铅铋流道中上升气泡的扩散迁移特性;通过对KYLIN-Ⅱ-S实验的模拟说明ACENA程序能够较为合理地预测熔融铅基合金-水相互作用过程中,铅池内压力波动和温度瞬变等现象;ACENA程序对ALFRED堆稳态满功率下关键热工参数的计算结果与国际认可的一维系统程序TRACE/FRED的计算结果基本一致,证明了ACENA程序全堆级计算结果的可靠性;对ALFRED堆假想SGTR事故的计算验证了ACENA程序对铅冷快堆SGTR事故下复杂多组分多相流动现象的模拟能立,且计算结果表明合理设计一次侧冷却剂循环路径、尽可能降低管道破损数量均对消减铅冷快堆SGTR事故后果具有重要意义。本工作可为我国池式铅冷快堆SGTR事故安全分析提供技术参考。

     

    Abstract: Steam generator tube rupture (SGTR) accident is one of the accident scenarios that must be considered during the design and safety analysis process of lead-cooled fast reactors (LFR). Focused on this topic, numerous experiments and theoretical research have been conducted in recent decades, which could be generally categorized into two aspects: on the one hand, separation effect experiments and theoretical analysis focusing on the specific phenomena of lead-water interaction process; on the other hand, engineering test facilities and numerical simulations focusing on the overall process of LFR SGTR accident, mainly concentrated on the evolution of macro parameters such as system pressure and lead pool temperature during SGTR accident. In this paper, firstly, the basic mathematical and physical models of ACENA code was introduced; then the code was validated against the LBE-N2 two-phase flow experiment HESTIA-2, the KYILN-Ⅱ-S LBE-H2O interaction experiment and the analytical solution of the point-kinetic neutronic equations; after that the Eulerian multi-phase analysis code ACENA was applied to simulate the postulated unprotected SGTR process of European pool-type lead-cooled fast reactor ALFRED. Attention was paid to the simulated migration process of steam bubbles in the lead pool as well as the evolution of safety parameters such as maximum cladding temperature, fuel centerline temperature, relative nuclear power and reactor vessel pressure. Furthermore, the sensitivity of ruptured tube quantities, lead coolant circulation paths and neutron dynamics parameters on calculation results were analyzed. This paper’s study demonstrates that the Ishii-Chawla-Suzuki interfacial drag coefficient model could accurately predict the diffusion and migration characteristics of rising bubbles in circular/annular LBE flow channels when combined with the interfacial area concentration transport model proposed by Ishii et al; The simulation of KYLIN-Ⅱ-S experiment shows that ACENA code could reasonably reproduce phenomena such as system pressure fluctuations and temperature transients of liquid lead during the interaction process between LBE and water; the calculated results of ALFRED reactor’s key thermal-hydraulic parameters under hot full power condition by ACENA code are generally consistent with results obtained by internationally recognized one-dimensional system code TRACE/FRED, proving the reliability of ACENA code when applied to full reactor system level; The calculation of the hypothetical SGTR accident of the ALFRED reactor verifies the ACENA code’s capability of simulating the complex multi-component multi-phase flow, and the calculation results reveal that designing the primary coolant flow path rationally and reducing the number of ruptured pipes are of great significance to mitigating the potential consequences of LFR SGTR accidents. This work can provide reference for the safety analysis work of China’s pool-type lead-cooled fast reactors.

     

/

返回文章
返回