Development and Preliminary Validation of Three-dimensional Subchannel Thermal-hydraulic Analysis Code CorTAF for Numerical Reactor Core
LIU Kai,WANG Mingjun,TIAN Wenxi,QIU Suizheng,SU Guanghui
DOI: https://doi.org/10.7538/yzk.2021.youxian.0860
2022-01-01
Abstract:Reactor core is the pivotal component of nuclear power system, and its integrity is an important prerequisite for safe operation of reactor. Traditional nuclear reactor core thermalhydraulic analysis methods cannot meet the high precision simulation requirements of advanced nuclear power systems in the future. In this paper, based on open source CFD platform OpenFOAM, a coolant flow heat transfer model considering diffusion due to turbulent mixing between adjacent coolant channels and rod bundle structure characteristics of pressurized water reactor (PWR) was established. A fuel rod heat conduction model was built to describe the internal temperature distribution formed by multiple nodes placed along the radial direction of fuel rods, and a coupled heat transfer model were proposed according to the mapping relationship between the external boundary condition of fuel rod heat conduction equation and convective heat transfer energy source term in coolant governing equation. The threedimensional subchannel thermal-hydraulic characteristics analysis code CorTAF for PWR core based on finite volume method was developed. After that, various fuel assembly flow and heat transfer experiments were selected to carry out model validation. For GE3×3 experiment, the calculated results of coolant velocity distribution in each channel of test section agree well with the experimental data, indicating that the CorTAF code can predict the flow characteristics in the fuel assembly effectively against COBRA code with similar computational accuracy. For Weiss experiment, the calculated results of coolant flow rate distribution in each assembly of test section are in good agreement with the experimental data, which illustrates that the CorTAF code can accurately predict the flow characteristics in parallel open assemblies under inlet partially blocked condition. For PNL2×6 experiment, the calculated results of coolant temperature and velocity distribution along center line of experimental section are basically consistent with the experimental data and those calculated by CUPID and MATRA codes, demonstrating that the CorTAF code can obtain the flow and heat transfer characteristics in the rod bundle assembly under condition of transverse nonuniform heating, and the error may be caused by the difference of data acquisition methods. In conclusion, the CorTAF code is proved to be suitable for predicting the flow and heat transfer characteristics in the rod bundle fuel assembly. The work in this paper has reference significance for the development of core safety analysis tools for PWR. Subsequently, further model optimization and thermalhydraulic characteristics analysis of the full PWR core under both steadystate operation and transient accident conditions will be carried out, and research on neutronics and thermalhydraulics coupling will be performed based on OpenFOAM.