A Monte Carlo study on burnup treatment in sodium-cooled reactor with Th fuel

Kerntechnik ◽  
2021 ◽  
Vol 86 (4) ◽  
pp. 302-311
Author(s):  
M. E. Korkmaz ◽  
N. K. Arslan

Abstract Sodium Cooled Reactors is one of the Generation-IV plants selected to manage the long-lived minor actinides and to transmute the long-life radioactive elements. This study presents the comparison between two-designed SFR cores with 600 and 800 MWth total heating power. We have analyzed a conceptual core design and nuclear characteristic of SFR. Monte Carlo depletion calculations have been performed to investigate essential characteristics of the SFR core. The core calculations were performed by using the Serpent Monte Carlo code for determining the burnup behavior of the SFR, the power distribution and the effective multiplication factor. The neutronic and burn-up calculations were done by means of Serpent-2 Code with the ENDF-7 cross-sections library. Sodium Cooled Fast Reactor core was taken as the reference core for Th-232 burnup calculations. The results showed that SFR is an important option to deplete the minor actinides as well as for transmutation from Th-232 to U-233.

2016 ◽  
Vol 6 (2) ◽  
pp. 21-30
Author(s):  
Huu Tiep Nguyen ◽  
Viet Phu Tran ◽  
Tuan Khai Nguyen ◽  
Vinh Thanh Tran ◽  
Minh Tuan Nguyen

This paper presents the results of neutronic calculations using the deterministic and Monte-Carlo methods (the SRAC and MCNP5codes) for the VVER MOX Core Computational Benchmark Specification and the VVER-1000/V392 reactor core. The power distribution and keff value have been calculated for a benchmark problem of VVER core. The results show a good agreement between the SRAC and MCNP5 calculations. Then, neutronic characteristics of VVER-1000/V392 such as power distribution, infinite multiplication factor (k-inf) of the fuel assemblies, effective multiplication factor keff, peaking factor and Doppler coefficient were calculated using the two codes.


2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
Wonkyeong Kim ◽  
Jinsu Park ◽  
Tomasz Kozlowski ◽  
Hyun Chul Lee ◽  
Deokjung Lee

A high-leakage core has been known to be a challenging problem not only for a two-step homogenization approach but also for a direct heterogeneous approach. In this paper the DIMPLE S06 core, which is a small high-leakage core, has been analyzed by a direct heterogeneous modeling approach and by a two-step homogenization modeling approach, using contemporary code systems developed for reactor core analysis. The focus of this work is a comprehensive comparative analysis of the conventional approaches and codes with a small core design, DIMPLE S06 critical experiment. The calculation procedure for the two approaches is explicitly presented in this paper. Comprehensive comparative analysis is performed by neutronics parameters: multiplication factor and assembly power distribution. Comparison of two-group homogenized cross sections from each lattice physics codes shows that the generated transport cross section has significant difference according to the transport approximation to treat anisotropic scattering effect. The necessity of the ADF to correct the discontinuity at the assembly interfaces is clearly presented by the flux distributions and the result of two-step approach. Finally, the two approaches show consistent results for all codes, while the comparison with the reference generated by MCNP shows significant error except for another Monte Carlo code, SERPENT2.


2021 ◽  
Vol 247 ◽  
pp. 04024
Author(s):  
Yurii Bilodid ◽  
Jaakko Leppänen

One of challenges of the Monte Carlo full core simulations is to obtain acceptable statistical variance of local parameters throughout the whole reactor core at a reasonable computation cost. The statistical variance tends to be larger in low-power regions. To tackle this problem, the Uniform-Fission-Site method was implemented in Monte Carlo code MC21 and its effectiveness was demonstrated on NEA Monte Carlo performance benchmark. The very similar method is also implemented in Monte Carlo code Serpent under the name Uniform Fission Source (UFS) method. In this work the effect of UFS method implemented in Serpent is studied on the BEAVRS benchmark which is based on a real PWR core with relatively flat radial power distribution and also on 3x3 PWR mini-core simulated with thermo-hydraulic and thermo-mechanic feedbacks. It is shown that the application of the Uniform Fission Source method has no significant effect on radial power variance but equalizes axial distribution of variance of local power.


Author(s):  
Muhammad Imron ◽  
Donny Hartanto

Abstract This paper presents static and transient solutions for the PWR MOX/UO2 transient benchmark by Serpent 2 Monte Carlo code and open nodal core simulator called ADPRES. The presences of MOX fuels and burn-up variation in the benchmark’s reactor core pose challenges for reactor simulators due to severe flux gradient across fuel assemblies. In this work, the two-step method was used, in which the assembly level two-group constants were generated from single assembly calculations with zero net current boundary conditions using Serpent 2 Monte Carlo code, and later the core calculation was performed using ADPRES open nodal core simulator. Two types of diffusion coefficients were generated: the conventional B1 leakage corrected and Cumulative Migration Method (CMM). Finally, the solutions of Serpent 2/ADPRESS, including multiplication factor, power distribution, control rod worth, and critical boron concentration using both diffusion coefficients were compared against solutions from heterogeneous Serpent 2 calculations where the fuel and cladding are explicitly modeled. The reactor power during transients were also compared qualitatively against other nodal core simulators. The results showed that Serpent 2/ADPRES were able to predict the heterogeneous Monte Carlo solutions very well with reasonable differences. The transient solutions were also quite accurate compared to other nodal core simulators. As for the diffusion coefficients comparison, it was found that the CMM diffusion coefficient provide more accurate solutions for the benchmark compared to the B1 leakage corrected diffusion coefficients.


Author(s):  
Yuxuan Liu ◽  
Ganglin Yu ◽  
Kan Wang

Monte Carlo codes are powerful and accurate tools for reactor core calculation. Most Monte Carlo codes use the point-wise data format, in which the data are given as tables of energy-cross section pairs. When calculating the cross sections at an incident energy value, it should be determined which grid interval the energy falls in. This procedure is repeated so frequently in Monte Carlo codes that its contribution in the overall calculation time can become quite significant. In this paper, the time distribution of Monte Carlo method is analyzed to illustrate the time consuming of cross section calculation. By investigation on searching and calculating cross section data in Monte Carlo code, a new search algorithm called hash table is elaborately designed to substitute the traditional binary search method in locating the energy grid interval. The results indicate that in the criticality calculation, hash table can save 5%∼17% CPU time, depending on the number of nuclides in the material, as well as complexity of geometry for particles tracking.


Author(s):  
Davide Chersola ◽  
Guglielmo Lomonaco ◽  
Guido Mazzini

This paper reports the results of a comparison among JEFF and ENDF/B datasets when used by SERPENT and MONTEBURNS codes on a GFR-like configuration. Particularly, it shows a comparison between the two Monte Carlo based codes, each one adopting three different cross sections dataset, namely JEFF-3.1, JEFF-3.1.2 and ENDF/B-VII.1. Calculations have been carried out on the Allegro reactor, i.e. an experimental GFR-like facility that should be built in EU as GFR demonstrator. Results concern nuclear parameters as effective multiplication factor and fluxes, as well as the atomic densities for some important nuclides versus burnup.


2021 ◽  
Vol 247 ◽  
pp. 06001
Author(s):  
Riku Tuominen ◽  
Ville Valtavirta ◽  
Manuel García ◽  
Diego Ferraro ◽  
Jaakko Leppänen

In coupled calculations with Monte Carlo neutronics and thermal hydraulics the Monte Carlo code is used to produce a power distribution which in practice means tallying the energy deposition. Usually the energy deposition is estimated by making a simple approximation that energy is deposited only in fission reactions. The goal of this work is to study how the accuracy of energy deposition modelling affects the results of steady state coupled calculations. For this task an internal coupling between Monte Carlo transport code Serpent 2 and subchannel code SUBCHANFLOW is used along with a recently implemented energy deposition treatment of Serpent 2. The new treatment offers four energy deposition modes each of which offers a different combination of accuracy and required computational time. As a test case, a 3D PWR fuel assembly is modelled with different energy deposition modes. The resulting effective multiplication factors are within 30 pcm. Differences of up to 100K are observed in the fuel temperatures.


2021 ◽  
Vol 247 ◽  
pp. 06011
Author(s):  
A. Bernal ◽  
M. Pecchia ◽  
D. Rochman ◽  
A. Vasiliev ◽  
H. Ferroukhi

The main goal of this work is to perform pin-by-pin calculations of Swiss LWR fuel assemblies with neutron transport deterministic methods. At Paul Scherrer Institut (PSI), LWR calculations are performed with the core management system CMSYS, which is based on the Studsvik suite of codes. CMSYS includes models for all the Swiss reactors validated against a database of experimental information. Moreover, PSI has improved the pin power calculations by developing models of Swiss fuel assemblies for the Monte Carlo code MCNP, with the isotopic compositions obtained from the In-Core Fuel Management data of the Studsvik suite of codes, by using the SNF code. A step forward is to use a neutron code based on fast deterministic neutron transport methods. The method used in this work is based on a planar Method of Characteristics in which the axial coupling is solved by 1D SP3 method. The neutron code used is nTRACER. Thus, the methodology of this work develops nTRACER models of Swiss PWR fuel assemblies, in which the fuel of each pin and axial level is modelled with the isotopic composition obtained from SNF. This methodology was applied to 2D and 3D calculations of a Swiss PWR fuel assembly. However, this method has two main limitations. First, the cross sections libraries of nTRACER lack some of the isotopes obtained by SNF. Fortunately, this work proves that the missing isotopes do not have a strong effect on keff and the power distribution. Second, the 3D models require high computational memory resources, that is, more than 260 Gb. Thus, the nTRACER code was modified, so now it uses only 8 Gb, without any loss of accuracy. Finally, the keff and power results are compared with Monte Carlo calculations obtained by Serpent.


Sign in / Sign up

Export Citation Format

Share Document