Modeling of SFR cores with Serpent–DYN3D codes sequence

Preservation of kinetics parameters generated by Monte Carlo calculations in two-step deterministic calculations

The generation of accurate kinetic parameters such as mean generation time Λ and effective delayed neutron fraction βeff via Monte Carlo codes is established. Employing these in downstream deterministic codes warrants another step to ensure no additional error is introduced by the low-order transport operator when computing forward and adjoint fluxes for bilinear weighting of these parameters. Another complexity stems from applying superhomogenization (SPH) equivalence in non-fundamental mode approximations, where reference and low-order calculations rely on a 3D full core model. In these cases, SPH factors can optionally be computed for only part of the geometry while preserving reaction rates and K-effective, but the impact of such approximations on kinetics parameters has not been thoroughly studied. This paper aims at studying the preservation of bilinearly-weighted quantities in the Serpent–Griffin calculation procedure. Diffusion and transport evaluations of IPEN/MB-01, Godiva, and Flattop were carried out with the Griffin reactor physics code, testing available modeling options using Serpent-generated multigroup cross sections and equivalence data. Verifying Griffin against Serpent indicates sensitivities to multigroup energy grid selection and regional application of SPH equivalence, introducing significant errors; these were demonstrated to be reduced through the use of a transport method together with a finer energy grid.

Superphénix Benchmark Part I: Results of Static Neutronics

In the paper, the specification of a new neutronics benchmark for large sodium cooled fast reactor (SFR) core and results of modeling by different participants are presented. The neutronics benchmark describes the core of the French sodium cooled reactor Superphénix at its startup configuration, which in particular was used for experimental measurement of reactivity characteristics. The benchmark consists of the detailed heterogeneous core specification for neutronic analysis and the results of the reference solution. Different core geometries and thermal conditions from the cold “as fabricated” up to full power were considered. The reference Monte Carlo (MC) solution of serpent 2 includes data on multiplication factor, power distribution, axial and radial reaction rates distribution, reactivity coefficients and safety characteristics, control rods worth, kinetic data. The results of modeling with seven other solutions using deterministic and MC methods are also presented and compared to the reference solution. The comparisons results demonstrate appropriate agreement of evaluated characteristics. The neutronics results will be used in the second phase of the benchmark for the evaluation of transient behavior of the core.

NEUTRONIC MODELLING OF THE FFTF CONTROL ROD WORTH MEASUREMENTS WITH DIFFUSION CODES

This paper presents an assessment of three deterministic core simulators with the focus on the neutronic performance in steady-state calculations of small Sodium cooled Fast Reactor cores. The selected codes are DYN3D, PARCS and the novel multi-physics solver GeN-Foam. By using these codes, the multi-group diffusion solutions are obtained for the selected twenty control rod worth measurements performed during the isothermal physics tests of the Fast Flux Test Facility (FFTF). The identical set of homogenized few-group cross sections applied in the calculations is generated with the Serpent Monte Carlo code. The numerical results are compared with each other as well as with the measured values. The obtained numerical results, such as the multiplication factors and control rod worth values, are in good agreement as compared to the experimental data. Furthermore, a comparison of the radial power distributions is presented between DYN3D, PARCS and GeN-Foam. Ultimately, the power distributions are compared to the full core Serpent solution, demonstrating an adequate performance of the selected deterministic tools. In overall, this study presents a verification and validation of the neutronic solvers applied by DYN3D, PARCS and GeN-Foam to steady-state calculations of SFR cores.

MODELLING ASTRID-LIKE SODIUM-COOLED FAST REACTOR WITH SERPENT-DYN3D CODE SEQUENCE

This study explores the feasibility of applying the Serpent-DYN3D sequence to the analysis of Sodium-cooled Fast Reactors (SFRs) with complex core geometries, such as the ASTRIDlike design. The core is characterised by a highly heterogeneous configuration and was likely to challenge the accuracy of the Serpent-DYN3D sequence. It includes axially heterogeneous fuel assemblies, non-uniform fuel assembly heights and large sodium plena. Consequently, the influence of generation and correction methods of various homogenised, few-group crosssections (XS) on the accuracy of the full-core nodal diffusion DYN3D calculations is presented. An attempt to compare the approximate time effort spent on models preparation against the accuracy of the result is made. Results are compared to reference full-core Serpent MC (Monte Carlo) solutions. Initially, XS data was generated in Serpent using traditional methods (2D single assemblies and 2D super-cells). Full core calculations and MC simulations offered a moderate agreement. Therefore, XS generation with 2D fuel-reflector models and 3D single assembly models was verified. Super-homogenisation (SPH) factors for XS correction were applied. In conclusion, the performed work suggests that Serpent-DYN3D sequence could be used for the analysis of highly heterogeneous SFR designs similar to the studied ASTRID-like, with an only small penalty on the accuracy of the core reactivity and radial power distribution prediction. However, the XS generation route would need to include the correction with SPH factors and generation of XS with various MC models, for different core regions. At a certain point, there are diminishing returns to using more complex XS generation methods, as the accuracy of full-core deterministic calculations improves only slightly, while the time effort required increases significantly.

New version of the reactor dynamics code DYN3D for Sodium cooled Fast Reactor analyses

The reactor dynamics code DYN3D being developed at the Helmholtz-Zentrum Dresden-Rossendorf is currently under extension for Sodium cooled Fast Reactor analyses. This paper provides an overview on the new version of DYN3D to be used for SFR core calculations. The current article shortly describes the newly implemented thermal mechanical models, which can account for thermal expansion effects of the reactor core. Furthermore, the methodology used in Sodium cooled Fast Reactor analyses to generate homogenized few-group cross sections is summarized. The conducted and planned verification and validation studies are briefly presented. Related publications containing more detailed descriptions are outlined for the completeness of this overview.

Assessing reactor physics codes capabilities to simulate fast reactors on the example of the BN-600 Benchmark

This work aims to assess the capabilities of reactor physics codes (initially validated for thermal reactors) to simulate fast sodium cooled reactors. The BFS-62-3A critical experiment from the BN-600 Hybrid Core Benchmark Analyses was chosen for the investigation. Monte-Carlo codes (KENO from SCALE and SERPENT 2.1.23) and the deterministic diffusion code DYN3D-MG are applied to calculate the neutronic parameters. It was found that the multiplication factor and reactivity effects calculated by KENO and SERPENT using the ENDF/B-VII.0 continuous energy library are in a good agreement with each other and with the measured benchmark values. Few-groups macroscopic cross sections, required for DYN3D-MG, were prepared in applying different methods implemented in SCALE and SERPENT. The DYN3D-MG results of a simplified benchmark show reasonable agreement with results from Monte-Carlo calculations and measured values. The former results are used to justify DYN3D-MG implementation for sodium cooled fast reactors coupled deterministic analysis.

The reactor dynamics code DYN3D and its trigonal-geometry nodal diffusion model

The reactor dynamics code DYN3D is a three-dimensional best-estimate tool for simulating steady states and transients of light-water reactors and innovative reactor designs. An overview of the DYN3D features is provided. This paper further focuses on the recently developed trigonal-geometry diffusion model DYN3D-TRIDIF including a description of the underlying nodal approach and the characteristics of trigonal geometries. Via a mesh refinement study by means of a VVER-1000-type core benchmark using a fine-mesh diffusion reference solution, DYN3D-TRIDIF shows spatial convergence. Furthermore, the performance of DYN3D-TRIDIF is verified by means of a single-assembly problem on pin-cell level. Good agreement between DYN3D-TRIDIF and the detailed-geometry transport reference is achieved with an average deviation in power of less than 1%.