基于CFD的棱柱型模块式高温气冷堆典型事故研究

CFD-based Investigation on Typical Accident for Prismatic-type Modular High Temperature Gas-cooled Reactor

  • 摘要: 为验证和评估棱柱型模块式高温气冷堆设计的固有安全性,需针对代表性事故工况开展计算分析。目前针对棱柱型堆芯的模块式高温气冷堆尚缺少专用的事故分析程序。本研究基于通用CFD程序COMSOL针对堆芯活性区域和压力容器建立三维模型,包括燃料和冷却剂通道、石墨慢化剂、侧反射层以及压力容器;非能动余热排出系统采用对流边界条件简化模拟。采用C++编写点堆模块求解中子动力学,并通过动态链接库(DLL)与COMSOL实现耦合。首先计算了正常运行工况下的稳定状态;然后以该结果作为初始条件,选取3个典型事故瞬态工况开展了数值模拟,包括未失压丧失强迫流动冷却(PLOFC)事故、未失压丧失强迫流动冷却且未能停堆(PLOFC+ATWS)事故以及反应性引入且未能停堆(RIA+ATWS)事故;最后针对压力容器壁与非能动余热排出系统的辐射发射率开展了敏感性分析。计算结果表明:在本文分析的事故条件下,燃料最高温度均低于安全限值(1 620 ℃)且具有较大的裕量,因此均能保证堆芯燃料结构的完整性。对于PLOFC事故,提高非能动余热排出系统的换热能力能显著缓解事故后果,但对于ATWS类事故影响趋势则正好相反,需进一步开展综合分析和模型验证。

     

    Abstract: Analyses of representative accident scenarios are imperative to assess and evaluate the inherent safety features of a prismatic modular high temperature gascooled reactor (HTGR). The accident progression of a HTGR is typically characterized by the coupling of multiphysics phenomena. However, computational codes available for the HTGR accident analyses are still relatively rare at present, especially those dedicated to the prismatictype core design. In addition, many of these codes are either lumpedparameter that cannot give the detailed threedimension (3D) distributions of key variables (e.g. fuel temperature), or are incapable of capturing the coupling of major phenomena. As a consequence, a generic CFD code, i.e. COMSOL multiphysics, was employed in this work to create a detailed 3D geometry representing the reactor, including fuel assemblies containing fuel and coolant channels, graphite moderator blocks, side reflector as well as the reactor pressure vessel. Specifically, the passive residual heat removal system (RCCS) was simply treated as a convective boundary condition due to the lack of detailed design scheme right now. A point kinetics module was prepared in this work using C++ language to account for the transient neutron kinetics. The Gear algorithm was applied to numerically solve the system of equations of the point kinetics, given that the stiffness is relatively large. The C++ was then compiled into dynamic link library (DLL), which was accessed by the COMSOL simulation in the form of an external function, with the fuel and moderator temperatures evaluated by COMSOL as input arguments. The steady state under normal operation condition was first calculated, which served as the initial condition for the subsequent transient analyses for three selected typical accidents, including pressurized loss of forced cooling (PLOFC), PLOFC without scram (PLOFC+ATWS), and reactivity insertion without scram (RIA+ATWS) accidents. Finally, sensitivity analysis was carried out on the effect of radiative emissivity between reactor pressure wall and RCCS. The simulation results show that the peak fuel temperature throughout the whole accident progression is remained far below its upper limit value (1 620 ℃) with a considerably large margin, implying that the integrity of fuel assemblies can be successfully guaranteed under the accident scenarios concerned in the present work. Moreover, enhancement of the RCCS heat removal capacity can significantly relief the accident consequence for PLOFC accident but lead to an inverse trend for ATWS accident, which underlines the importance of further comprehensive analysis and modeling validation. To conclude, the computational model developed in this work along with the derived results can serve as the starting point for future development of the accident analysis tool, which is applicable for the incoming indepth safety analysis and assessment works.

     

/

返回文章
返回