Explicit decay heat calculation in the nodal diffusion code DYN3D

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

Hybrid microscopic depletion model in nodal code DYN3D

The paper presents a general hybrid method that combines the micro-depletion technique with correction of micro- and macro- diffusion parameters to account for the spectral history effects. The fuel in a core is subjected to time- and space-dependent operational conditions (e.g. coolant density), which cannot be predicted in advance. However, lattice codes assume some average conditions to generate cross sections (XS) for nodal diffusion codes such as DYN3D. Deviation of local operational history from average conditions leads to accumulation of errors in XS, which is referred as spectral history effects. Various methods to account for the spectral history effects, such as spectral index, burnup-averaged operational parameters and micro-depletion, were implemented in some nodal codes. Recently, an alternative method, which characterizes fuel depletion state by burnup and ²³⁹Pu concentration (denoted as Pu-correction) was proposed, implemented in nodal code DYN3D and verified for a wide range of history effects. The method is computationally efficient, however, it has applicability limitations. The current study seeks to improve the accuracy and applicability range of Pu-correction method. The proposed hybrid method combines the micro-depletion method with a XS characterization technique similar to the Pu-correction method. The method was implemented in DYN3D and verified on multiple test cases. The results obtained with DYN3D were compared to those obtained with Monte Carlo code Serpent, which was also used to generate the XS. The observed differences are within the statistical uncertainties.This is the author accepted manuscript. The final version is available from Elsevier via http://dx.doi.org/10.1016/j.anucene.2016.02.01

Experience and Perspective of Best-Estimate Approach Application for RIA Analysis

The use of best-estimate approach for WWER safety analysis in RIA is considered. The relevance of this problem is concerned with small margin to acceptance criteria under the conservative approach and becomes stronger under power uprate of nuclear power plants. Previous experience in this area for WWER-1000 reactor types is overviewed. The necessity to extend these activities for successful implementation of the best-estimate approach is noticed and areas of further work are discussed.Рассмотрены вопросы использования и реализации подходов наилучшей оценки для анализа безопасности ВВЭР в реактивностных авариях. Актуальность проблемы связана с малыми запасами до критериев приемлемости при реализации консервативного подхода, особенно усиливается в условиях повышения номинального уровня мощности реакторной установки. Представлен краткий обзор предыдущего опыта в этой области для реакторов ВВЭР-1000. Указана необходимость расширения работ для успешной реализации подходов наилучшей оценки, обсуждены направления дальнейшей деятельности.Розглянуто питання використання та реалізації підходів найкращої оцінки для аналізу безпеки ВВЕР у реактивнісних аваріях. Актуальність проблеми пов’язана з малими запасами до критеріїв прийнятності при реалізації консервативного підходу, що особливо посилюється в умовах підвищення номінального рівня потужності реакторної установки. Представлено короткий огляд попереднього досвіду в цій галузі для реакторів ВВЕР-1000. Зазначено необхідність розширення робіт для успішної реалізації підходів найкращої оцінки, обговорено напрями подальшої діяльності

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

NEUTRON NOISE PATTERNS FROM COUPLED FUEL-ASSEMBLY VIBRATIONS

The neutron flux fluctuation magnitude of KWU-built PWRs shows a hitherto unexplained correlation with the types of loaded fuel assemblies. Also, certain measured long-range neutron flux fluctuation patterns in neighboring core quadrants still lack a closed understanding of their origin. The explanation of these phenomena has recently revived a new interest in neutron noise research. The contribution at hand investigates the idea that a synchronized coolant-driven vibration of major parts of the fuel-assembly ensemble leads to these phenomena. Starting with an assumed mode of such collective vibration, the resulting effects on the time-dependent neutron-flux distribution are analyzed via a DYN3D simulation. A three-dimensional representation of the time-dependent bow of all fuel assemblies is taken into account as a nodal DYN3D feedback parameter by time-dependent variations of the fuel-assembly pitch. The impact of its variation on the cross sections is quantified using a cross-section library that is generated from the output of corresponding CASMO5 calculations. The DYN3D simulation qualitatively reproduces the measured neutron-flux fluctuation patterns. The magnitude of the fluctuations and its radial dependence are comparable to the measured details. The results imply that collective fuel-assembly vibrations are a promising candidate for being the key to understand long-known fluctuation patterns in KWU built PWRs. Further research should elaborate on possible excitation mechanisms of the assumed vibration modes

Axial discontinuity factors for the nodal diffusion analysis of high conversion BWR cores

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. © 2013 Elsevier Ltd. All rights reserved

NEUTRON NOISE PATTERNS FROM COUPLED FUEL-ASSEMBLY VIBRATIONS

The neutron flux fluctuation magnitude of KWU-built PWRs shows a hitherto unexplained correlation with the types of loaded fuel assemblies. Also, certain measured long-range neutron flux fluctuation patterns in neighboring core quadrants still lack a closed understanding of their origin. The explanation of these phenomena has recently revived a new interest in neutron noise research. The contribution at hand investigates the idea that a synchronized coolant-driven vibration of major parts of the fuel-assembly ensemble leads to these phenomena. Starting with an assumed mode of such collective vibration, the resulting effects on the time-dependent neutron-flux distribution are analyzed via a DYN3D simulation. A three-dimensional representation of the time-dependent bow of all fuel assemblies is taken into account as a nodal DYN3D feedback parameter by time-dependent variations of the fuel-assembly pitch. The impact of its variation on the cross sections is quantified using a cross-section library that is generated from the output of corresponding CASMO5 calculations. The DYN3D simulation qualitatively reproduces the measured neutron-flux fluctuation patterns. The magnitude of the fluctuations and its radial dependence are comparable to the measured details. The results imply that collective fuel-assembly vibrations are a promising candidate for being the key to understand long-known fluctuation patterns in KWU built PWRs. Further research should elaborate on possible excitation mechanisms of the assumed vibration modes

Verification of coupled 3D fuel cycle analysis with stable Monte Carlo based code, BGCore, against the nodal diffusion DYN3D code

Previous studies suggested that different schemes for coupling Monte Carlo (MC) neutron transport with burnup and thermal hydraulic feedbacks may potentially be numerically unstable. This issue can be resolved by application of new Stochastic Implicit Mid-Point (SIMP) methods. In order to assure numerical stability, the new methods do require additional computational effort. The instability issue however, is problem-dependent and does not necessarily occur in all cases. Therefore, blind application of the unconditionally stable coupling schemes, and thus incurring extra computational costs, may not always be necessary. In this paper, we attempt to develop an intelligent diagnostic mechanism, which will monitor numerical stability of the calculations and, if necessary, switch from simple and fast explicit coupling scheme to more computationally expensive but unconditionally stable one. To illustrate this diagnostic mechanism, we performed a coupled burnup and TH analysis of a single BWR fuel assembly. The reference solution was obtained by the state of the art nodal diffusion code - DYN3D. Very good agreement was observed in all neutronic and TH parameters. The results indicate that the developed algorithm can be easily implemented in any MC based code for monitoring of numerical instabilities. The proposed monitoring method has negligible impact on the calculation time even for realistic 3D multi-region full core calculations

Spectral history model in DYN3D: Verification against coupled Monte-Carlo thermal-hydraulic code BGCore

This research focuses on the verification of a recently developed methodology accounting for spectral history effects in 3D full core nodal simulations. The traditional deterministic core simulation procedure includes two stages: (1) generation of homogenized macroscopic cross section sets and (2) application of these sets to obtain a full 3D core solution with nodal codes. The standard approach adopts the branch methodology in which the branches represent all expected combinations of operational conditions as a function of burnup (main branch). The main branch is produced for constant, usually averaged, operating conditions (e.g. coolant density). As a result, the spectral history effects that associated with coolant density variation are not taken into account properly. Number of methods to solve this problem (such as micro-depletion and spectral indexes) were developed and implemented in modern nodal codes. Recently, we proposed a new and robust method to account for history effects. The methodology was implemented in DYN3D and involves modification of the few-group cross section sets. The method utilizes the local Pu-239 concentration as an indicator of spectral history. The method was verified for PWR and VVER applications. However, the spectrum variation in BWR core is more pronounced due to the stronger coolant density change. The purpose of the current work is investigating the applicability of the method to BWR analysis. The proposed methodology was verified against recently developed BGCore system, which couples Monte Carlo neutron transport with depletion and thermal-hydraulic solvers and thus capable of providing a reference solution for 3D simulations. The results clearly show that neglecting the spectral history effects leads to a very large deviation (e.g. 1700 pcm in multiplication factor) from the reference solution. Application of the Pu-correction method results in a very good agreement between DYN3D and BGCore on the order of 200 pcm in kinf