Search

Recent developments of the reactor dynamics code DYN3D have introduced the micro-depletion model which allows for explicit calculation of radioactive decay heat. Such a unique combination of nodal diffusion, thermal hydraulic (T/H) and depletion solvers allows DYN3D to perform fuel cycle depletion and obtain detailed core isotopic concentration and decay heat distributions. More importantly, this capability allows for the generation of time-dependent spatial decay heat distribution, which is essential to predict the behavior of the reactor under off-nominal scenarios. The new sequence utilizes considerably less computational resources than coupled Monte Carlo-T/H-depletion systems, but with comparable accuracy. This capability was recently tested on a limited number of simple unit cell models. The main objective of this work is to further verify the decay heat calculation capabilities of DYN3D by applying it to a considerably more realistic and detailed full core model. For the purpose of the current analysis a 3D full core model of Advanced Burner Reactor (ABR) was adopted from the OECD/NEA Benchmark for Neutronic Analysis of Sodium-Cooled Fast Reactor Cores with Various Fuel Types and Core Sizes. In this work, the Monte Carlo code Serpent was used to generate macro- and microscopic parameters, and the neutron diffusion code, DYN3D, was used to perform neutronic and depletion analyses. Detailed spatial isotopic and decay heat distributions obtained with DYN3D were verified against the equivalent Serpent reference 3D full core solution. Results indicate very good agreement between the Serpent-DYN3D code sequence and the reference Serpent solutions, with a discrepancy in total decay heat on the order of 0.5%.

In order to compensate the excess reactivity over the fuel cycle in pressurized water reactors, boric acid is added to the reactor coolant. During certain accident conditions, the formation and the subsequent transportation towards the reactor core of a slug of coolant with reduced boron concentration can lead to a positive reactivity insertion. In this paper, a new approach dealing with such heterogeneous boron dilution scenarios is presented. This approach is based on the use of computational fluid dynamics methods for the whole reactor pressure vessel with direct coupling to a 3D neutron kinetic model of the reactor core. The application of this approach is demonstrated on a main coolant pump start scenario for hot-shutdown conditions using the coupled ANSYS CFX/DYN3D code. The calculation results show the enormous influence of the coolant mixing on the reactivity behavior of the reactor core.

This paper describes the derivation of the nodal flux expansion method HEXNEM3, its implementation into the nodal diffusion code DYN3D and the corresponding testing versus benchmarks. As in the earlier versions of expansion method HEXNEM1 and HEXNEM2, the neutron flux in a hexagonal node is expanded into superposition of orthogonal polynomials and exponential functions. The main difference of the HEXNEM3 method is the additional use of tangentially weighted exponential functions and the coupling of neighboring nodes by tangentially weighted fluxes and currents on node surfaces. The HEXNEM3 method was tested in several benchmark problems, including numerical benchmarks with given cross sections set and reference solutions by fine-mesh finite difference diffusion and a real plant benchmark with Monte Carlo reference solution. The test results demonstrate good agreement with reference solutions and improvement of method accuracy in comparison with HEXNEM1 and HEXNEM2.

The residual radioactive decay heat plays an important role in some accident scenarios and, therefore, needs to be accurately calculated when performing accident analyses. The current reactor simulation codes used for accident analysis account for the residual decay heat by means of simplified models. Typically, these models rely on semi-empirical correlations which are defined over a limited range of burnup and fuel types. Therefore, the applicability of such correlations is limited and any deviation from the definition range may lead to high uncertainties, which is detrimental in evaluating safety margins. Reactor dynamic code DYN3D was originally developed for transient and accident analysis of LWRs. In DYN3D, the residual radioactive decay heat calculation is based on the German national standard DIN Norm 25463 model. The applicability of this model is limited to a low enriched uranium dioxide fuel for light water reactors. This paper describes a new general decay heat calculation model implemented in DYN3D. The radioactive decay rate of each nuclide in each spatial node is calculated by recently implemented depletion module and the cumulative released heat is used to obtain the spatial distribution of the decay power for every time step. Such explicit approach is based on first principles and is free from approximations and, thus, can be applied to any reactor system (e.g. thermal and fast) and fuel type. The proposed method is verified through code-to-code comparison with the Serpent 2 Monte Carlo code results.

High conversion LWRs concepts typically rely on a heterogeneous core configuration, where fissile zones are interspersed with fertile blanket zones in order to achieve a high conversion ratio. Modeling such a heterogeneous structure of these cores represents a significant challenge to the conventional reactor analysis methods. It was recently suggested to overcome such difficulties, in particular, for the case of axially heterogeneous reduced moderation BWRs, by introducing an additional set of discontinuity factors in axial direction at the interfaces between fissile and fertile fuel assembly zones. However, none of the existing nodal diffusion core simulators have the capability of accounting for discontinuity of homogeneous nodal fluxes in axial direction since the fuel composition of conventional LWRs is much more axially uniform. In this work, we modified the nodal diffusion code DYN3D by introducing such a capability. The new version of the code was tested on a series of reduced moderation BWR cases with Th–U233 and U–Pu–MA fuel. The library of few-group homogenized cross sections and the data required for the calculation of discontinuity factors were generated using the Monte Carlo transport code Serpent. The results obtained with the modified version of DYN3D were compared with the reference Monte Carlo solutions and were found to be in good agreement. The current analysis demonstrates that high conversion LWRs can in principle be modeled using existing nodal diffusion core simulators.

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%.

The few-group cross section libraries, used by reactor dynamics codes, are affected by the spectral history effect – a dependence of fuel cross sections not only on burnup, but also on local spectral conditions during burnup. Neglecting this effect leads to an additional component of error in neutron-physical characteristics. Two solution approaches to this problem implemented in the reactor dynamic code DYN3D are described and compared in this paper: a cross section correction method based on239Pu concentration and separate cross sections treatment for each axial layer of reactor core. Steady-state and burnup characteristics of a PWR and a WWER-1000 cores, calculated by DYN3D with and without cross section corrections, are compared. An impact of the correction on transient calculations is studied for a control rod ejection example. Studies have shown a significant influence of spectral history on axial power and burnup distributions as well as on calculated cycle length. Two different correction methods have shown similar major effects.