Static and transient pin-by-pin simulations of a full PWR core with the extended coupled code system DYNSUB

Abstract The evolutionary multi-physics tool developed at the Karlsruhe Institute of Technology is the homogeneous pin-by-pin reactor simulator DYNSUB, an internal coupling of the 3D neutron kinetics code DYN3D developed by Helmholtz Zentrum Dresden Rossendorf and the in-house sub-channel code SUBCHANFLOW. The ultimate goal of the on-going efforts concerning DYNSUB is to provide a cost-effective improved description of light water reactor core behavior with pin-by-pin resolution for both static and transient safety relevant scenarios. A cost-effective computer code is defined to be executable on commodity computing clusters which users/customers commonly have access to. Efforts undertaken to improve DYNSUB’s numerical performance and parallelize the code system are presented in this work. Moreover, the coupled code system has been extended in terms of fuel pin level homogenization corrections and flexible mapping schemes. After optimization and extension DYNSUB is successfully applied to study the OECD/NEA and U.S. NRC PWR MOX/UO 2 core transient benchmark with both fuel assembly/channel and pin level/sub-channel model resolution. Even though further improvements in terms of numerical performance and accuracy of physical models are required, the applicability of DYNSUB pin-by-pin simulations for light water reactor safety analysis is proven in principle in this work.

[1]  Alain Hébert,et al.  Development of a Third-Generation Superhomogénéisation Method for the Homogenization of a Pressurized Water Reactor Assembly , 1993 .

[2]  Tomasz Kozlowski,et al.  Cell Homogenization Method for Pin-by-Pin Neutron Transport Calculations , 2011 .

[3]  K. Smith,et al.  ASSEMBLY HOMOGENIZATION TECHNIQUES FOR LIGHT WATER REACTOR ANALYSIS , 1986 .

[4]  Russian Federation,et al.  VALIDATION OF A PIN-BY-PIN HETEROGENEOUS METHOD AGAINST LWR MOX BENCHMARKS , 2002 .

[5]  Gerhard Wellein,et al.  LIKWID: A Lightweight Performance-Oriented Tool Suite for x86 Multicore Environments , 2010, 2010 39th International Conference on Parallel Processing Workshops.

[6]  Masahiro Tatsumi,et al.  Advanced PWR Core Calculation Based on Multi-group Nodal-transport Method in Three-dimensional Pin-by-Pin Geometry , 2003 .

[7]  Robert Stieglitz,et al.  Investigations of boron transport in a PWR core with COBAYA3/SUBCHANFLOW inside the NURESIM platform , 2014 .

[8]  Kiril Velkov,et al.  Coupled Monte Carlo neutronics and thermal hydraulics for power reactors , 2012 .

[9]  Kord Sterling Smith,et al.  Spatial homogenization methods for light water reactor analysis , 1980 .

[10]  Thomas J. Downar,et al.  High-Fidelity Light Water Reactor Analysis with the Numerical Nuclear Reactor , 2007 .

[11]  R. Stieglitz,et al.  Internal multi-scale multi-physics coupled system for high fidelity simulation of light water reactors , 2014 .

[12]  Benoit Forget,et al.  Benchmark for evaluation and validation of reactor simulations (BEAVRS) , 2013 .

[13]  Thomas J. Downar,et al.  Coupled full core neutron transport/cfd simulations of pressurized water reactors , 2012 .

[14]  Alain Hébert,et al.  A consistent technique for the pin-by-pin homogenization of a pressurized water reactor assembly , 1993 .

[15]  Javier Jiménez Escalante.  Desarrollo e implementación de la descomposición en subdominios del acoplamiento neutrónico-termohidráulico mediante disecciones alternadas con un sistema de cálculo multiescala , 2011 .

[16]  Gerhard Wellein,et al.  Performance Patterns and Hardware Metrics on Modern Multicore Processors: Best Practices for Performance Engineering , 2012, Euro-Par Workshops.

[17]  Tsutomu Ikeno,et al.  Impact of Pin-by-Pin Thermal-Hydraulic Feedback Modeling on Steady-State Core Characteristics , 2005 .

[18]  Tomoyuki Sugiyama,et al.  Clad-to-Coolant Heat Transfer in NSRR Experiments , 2007 .

[19]  Yoshihiro Yamane,et al.  Improvement of the SPH Method for Pin-by-Pin Core Calculations , 2004 .

[20]  D. Lu,et al.  The calculation method for SP3 discontinuity factor and its application , 2014 .

[21]  Kostadin Ivanov,et al.  DYNSUB: A high fidelity coupled code system for the evaluation of local safety parameters – Part II: Comparison of different temporal schemes , 2012 .

[22]  Gerhard Wellein,et al.  LIKWID: Lightweight Performance Tools , 2011, CHPC.

[23]  Edward W. Larsen,et al.  The Simplified P3 Approximation , 2000 .

[24]  N. García-Herranz,et al.  Neighborhood-corrected interface discontinuity factors for multi-group pin-by-pin diffusion calculations for LWR , 2012 .

[25]  M. Avramova,et al.  New developments of the MCNP/CTF/NEM/NJOY code system – Monte Carlo based coupled code for high accuracy modeling , 2013 .

[26]  J. I. Yoon,et al.  Core follow calculation with the nTRACER numerical reactor and verification using power reactor measurement data , 2013 .

[27]  Kan Wang,et al.  The first-principle coupled calculations using TMCC and CFX for the pin-wise simulation of LWR , 2012 .

[28]  Haileyesus Tsige-Tamirat,et al.  Coupled Neutronics Thermal-Hydraulics Analysis using Monte Carlo and Sub-channel Codes , 2012 .

[29]  José Javier Herrero Carrascosa,et al.  Whole Core Pin-by-Pin Coupled Neutronic-Thermal-hydraulic Steady state and Transient Calculations using COBAYA3 code , 2010 .

[30]  Ali R Massih,et al.  Nuclear Fuel Behaviour under Reactivity-initiated Accident (RIA) Condition: State-of-the-art Report , 2010 .

[31]  Richard Sanchez,et al.  Assembly homogenization techniques for core calculations , 2009 .

[32]  Benoit Forget,et al.  Challenges in the development of high-fidelity LWR core neutronics tools , 2013 .

[33]  José Gutiérrez Abascal,et al.  PERFORMANCE OF WHOLE CORE PIN-BY-PIN CALCULATIONS BY DOMAIN DECOMPOSITION THROUGH ALTERNATE DISSECTIONS IN STEADY STATE AND TRANSIENT CALCULATIONS , 2009 .

[34]  Kostadin Ivanov,et al.  DYNSUB: A high fidelity coupled code system for the evaluation of local safety parameters – Part I: Development, implementation and verification , 2012 .

[35]  Kostadin Ivanov,et al.  High fidelity simulation of conventional and innovative LWR with the coupled Monte-Carlo thermal-hydraulic system MCNP-SUBCHANFLOW , 2013 .

[36]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[37]  Vincent Bessiron,et al.  Modelling of Clad-to-Coolant Heat Transfer for RIA Applications , 2007 .

[38]  S. Mittag,et al.  Super-homogenisation factors in pinwise calculations by the reactor dynamics code DYN3D , 2011 .