Coupling a CFD code with neutron kinetics and pin thermal models for nuclear reactor safety analyses

Abstract Most system codes are based on the one-dimensional lumped-parameter method, which is unsuitable to simulate multi-dimensional thermal-hydraulics problems. CFD method is a good tool to simulate multi-dimensional thermal-hydraulics phenomena in the nuclear reactor, which can increase the accuracy of analysis results. However, since there is no neutron kinetics model and pin thermal model in current CFD codes, the application of the CFD method in the area of nuclear reactor safety analyses is still limited. Coupling a CFD code with the neutron kinetics model (PKM) and the pin thermal model (PTM) is a good way to use CFD code to simulate multi-dimensional thermal-hydraulics problems of nuclear reactors. The motivation for this work is to develop a CFD/neutron kinetics coupled code named FLUENT/PK for nuclear reactor safety analyses by coupling the commercial CFD code named FLUENT with the point kinetics model (PKM) and the pin thermal model (PTM). The mathematical models and the coupling method are described and the unprotected transient overpower (UTOP) accident of a liquid metal cooled fast reactor (LMFR) is chosen as an application case. As a general validation, the calculated results are used to compare with that of another multi-physics coupled code named SIMMER-III and good agreements are achieved for various characteristic parameters.

[1]  Sergio Chiva,et al.  A New Coupled CFD/Neutron Kinetics System for High Fidelity Simulations of LWR Core Phenomena: Proof of Concept , 2014 .

[2]  V. I. Subbotin,et al.  Investigation of heat transfer in the turbulent flow of liquid metals in tubes , 1961 .

[3]  Nam-il Tak,et al.  Investigation on turbulent heat transfer to lead–bismuth eutectic flows in circular tubes for nuclear applications , 2006 .

[4]  Pavel V. Tsvetkov,et al.  Fast spectrum reactors , 2012 .

[5]  M. Aufiero DEVELOPMENT OF ADVANCED SIMULATION TOOLS FOR CIRCULATING-FUEL NUCLEAR REACTORS , 2014 .

[6]  P. Savva,et al.  Multi-physics and multi-scale methods used in nuclear reactor analysis , 2014 .

[7]  Francesco Saverio D'Auria,et al.  Use of Computational fluid dynamics codes for safety analysis of nuclear reactor systems , 2003 .

[8]  Werner Maschek,et al.  Transient analyses for accelerator driven system PDS-XADS using the extended SIMMER-III code , 2005 .

[9]  K. Nishihara,et al.  Transient analyses for lead–bismuth cooled accelerator-driven system , 2013 .

[10]  Hongli Chen,et al.  Development and application of a thermal-hydraulics analysis code for small modular natural circulation lead or lead-alloy cooled fast reactors , 2014 .

[11]  Jia Bao-shan Preliminary Research on Nuclear Reactor Thermal-Hydraulic Multi-scale Coupled Simulation , 2010 .

[12]  Klas Jareteg,et al.  Development of an integrated deterministic neutronic/thermal-hydraulic model using a CFD solver , 2012 .

[13]  Zhao Chen,et al.  CFD Analysis of Thermal Stratification of China Lead Alloy Cooled Research Reactor (CLEAR-I) , 2013 .

[14]  Barry D Ganapol,et al.  A highly accurate algorithm for the solution of the point kinetics equations , 2013 .

[15]  J. Vierendeels,et al.  CFD analysis of the MYRRHA primary cooling system , 2011 .

[16]  Hongli Chen,et al.  Study of core flow distribution for small modular natural circulation lead or lead-alloy cooled fast reactors , 2014 .

[17]  Werner Maschek,et al.  Application of safety analysis code SIMMER-IV to blockage accidents in FASTEF subcritical core , 2014 .

[18]  Weimin Han,et al.  Numerical Solution of Ordinary Differential Equations: Atkinson/Numerical , 2009 .

[19]  A. Nahla,et al.  An efficient technique for the point reactor kinetics equations with Newtonian temperature feedback effects , 2011 .