Investigations of boron transport in a PWR core with COBAYA3/SUBCHANFLOW inside the NURESIM platform

Abstract Multi-scale, multi-physics problems reveal significant challenges while dealing with coupled neutronic/thermal-hydraulic solutions. Current generations of reactor dynamic codes applied to Light Water Reactors (LWRs), and in the context of this paper, pressurized water reactors (PWRs), are based on 3D neutronic nodal methods coupled with single or two phase flow thermal hydraulic system or sub-channel codes. This paper describes the extension of the coupling scheme between the 3D neutron diffusion codes COBAYA3 and DYN3D, and the sub-channel thermal-hydraulic code SUBCHANFLOW for the simulation of boron dilution transients. This includes the validation of the boron transport model of SUBCHANFLOW. The coupling of COBAYA3 and DYN3D with SUBCHANFLOW is performed inside the NURESIM platform making use of the novel automatic mesh superposition for code coupling. Boron transport models are needed for the simulation of boron dilution transients following a SBLOCA (after loop-seal clearing). In this case, the mixing is a key mechanism determining the positive reactivity insertion in the core. The results obtained with the coupled codes COBAYA3/SUBCHANFLOW, DYN3D/SUBCHANFLOW will be presented and discussed for two transient scenarios; a homogeneous and a heterogeneous boron dilution problem for a mini-core and for a PWR core defined in the NURISP boron dilution Benchmark.

[1]  N. García-Herranz,et al.  The Analytic Coarse-Mesh Finite Difference Method for Multigroup and Multidimensional Diffusion Calculations , 2007 .

[2]  E. Fridman,et al.  Development of the coupled 3D neutron kinetics/thermal-hydraulics code DYN3D-HTR for the simulation of transients in block-type HTGR , 2011 .

[3]  Juan Andrés Lozano Montero.  Desarrollo de un código nodal analítico de difusión neutrónica en multigrupos para reactores nucleares en geometrías 3d rectangulares y hexagonales , 2010 .

[4]  N. García-Herranz,et al.  The analytic nodal diffusion solver ANDES in multigroups for 3D rectangular geometry: Development and performance analysis , 2008 .

[5]  R. Riggs,et al.  Prioritization of generic safety issues , 1983 .

[6]  Soeren Kliem,et al.  Coupling of the neutron-kinetic core model DYN3D with the thermal hydraulic code FLICA-4 within the NURESIM platform , 2010 .

[7]  Soeren Kliem,et al.  Development and verification of the coupled 3D neutron kinetics/thermal-hydraulics code DYN3D-HTR for the simulation of transients in block-type HTGR , 2012 .

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

[9]  Manuel Calleja Reyna.  Improvements in Multi-physics and Multi-scale Methodologies for Safety-related Investigations of Pressurized Water Reactors within the NURESIM Platform , 2013 .

[10]  Sören Kliem.  Ein Modell zur Beschreibung der Kühlmittelvermischung und seine Anwendung auf die Analyse von Borverdünnungstransienten in Druckwasserreaktoren , 2010 .

[11]  Soeren Kliem,et al.  Experiments at the mixing test facility ROCOM for benchmarking of CFD codes , 2008 .

[12]  B. P. Leonard,et al.  The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection , 1991 .

[13]  I. Toumi,et al.  FLICA-4: a three-dimensional two-phase flow computer code with advanced numerical methods for nuclear applications , 2000 .

[14]  Uwe Imke,et al.  Validation of the Subchannel Code SUBCHANFLOW Using the NUPEC PWR Tests (PSBT) , 2012 .

[15]  Robert Stieglitz,et al.  A coupled neutronic/thermal-hydraulic scheme between COBAYA3 and SUBCHANFLOW within the NURESIM simulation platform , 2012 .