Validating the Serpent-Ants Calculation Chain Using BEAVRS Fresh Core HZP Data

Author(s):  
Ville Valtavirta ◽  
Antti Rintala ◽  
Unna Lauranto

Abstract The Serpent Monte Carlo code and the Serpent-Ants two step calculation chain are used to model the hot zero power physics tests described in the BEAVRS benchmark. The predicted critical boron concentrations, control rod group worths and isothermal temperature coefficients are compared between Serpent and Serpent-Ants as well as against the experimental measurements. Furthermore, radial power distributions in the unrodded and rodded core configurations are compared between Serpent and Serpent-Ants. In addition to providing results using a best practices calculation chain, the effects of several simplifications or omissions in the group constant generation process on the results are estimated. Both the direct and two-step neutronics solutions provide results close to the measured values. Comparison between the measured data and the direct Serpent Monte Carlo solution yields RMS differences of 12.1 mg/kg, 25.1 × 10-5 and 0.67 × 10-5 K-1 for boron, control rod worths and temperature coefficients respectively. The two-step Serpent-Ants solution reaches a similar level of accuracy with RMS differences of 17.4 mg/kg, 23.6 × 10-5 and 0.29 × 10-5 K-1. The match in the radial power distribution between Serpent and Serpent-Ants was very good with the RMS and maximum for pin power errors being 1.31 % and 4.99 % respectively in the unrodded core and 1.67 %(RMS) and 8.39 % (MAX) in the rodded core.

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.


2020 ◽  
Vol 225 ◽  
pp. 03007
Author(s):  
Tanja Goričanec ◽  
Domen Kotnik ◽  
Žiga Štancar ◽  
Luka Snoj ◽  
Marjan Kromar

An approach for calculating ex-core detector response using Monte Carlo code MCNP was developed. As a first step towards ex-core detector response prediction a detailed MCNP model of the reactor core was made. A script called McCord was developed as a link between deterministic program package CORD-2 and Monte Carlo code MCNP. It automatically generates an MCNP input from the CORD-2 data. A detailed MCNP core model was used to calculate 3D power distributions inside the core. Calculated power distributions were verified by comparison to the CORD-2 calculations, which is currently used for core design calculation verification of the Krško nuclea power plant. For the hot zero power configuration, the deviations are within 3 % for majority of fuel assemblies and slightly higher for fuel assemblies located at the core periphery. The computational model was further verified by comparing the calculated control rod worth to the CORD-2 results. The deviations were within 50 pcm and considered acceptable. The research will in future be supplemented with the in-core and ex-core detector signal calculations and neutron transport outside the reactor core.


2016 ◽  
Vol 96 ◽  
pp. 332-343 ◽  
Author(s):  
M.V. Shchurovskaya ◽  
V.P. Alferov ◽  
N.I. Geraskin ◽  
A.I. Radaev ◽  
A.G. Naymushin ◽  
...  

2016 ◽  
Vol 2 (2) ◽  
Author(s):  
Haykel Raouafi ◽  
Guy Marleau

The Canadian-SCWR is a heavy-water moderated supercritical light-water-cooled pressure tube reactor. It is fueled with CANada deuterium uranium (CANDU)-type bundles (62 elements) containing a mixture of thorium and plutonium oxides. Because the pressure tubes are vertical, the upper region of the core is occupied by the inlet and outlet headers render it nearly impossible to insert vertical control rods in the core from the top. Insertion of solid control devices from the bottom of the core is possible, but this option was initially rejected because it was judged impractical. The option that is proposed here is to use inclined control rods that are inserted from the side of the reactor and benefit from the gravitational pull exerted on them. The objective of this paper is to evaluate the neutronic performance of the proposed inclined control rods. To achieve this goal, we first develop a three-dimensional (3D) supercell model to simulate an inclined rod located between four vertical fuel cells. Simulations are performed with the SERPENT Monte Carlo code at five axial positions in the reactor to evaluate the effect of coolant temperature and density, which varies substantially with core height, on the reactivity worth of the control rods. The effect of modifying the inclination and spatial position of the control rod inside the supercell is then analyzed. Finally, we evaluate how boron poisoning of the moderator affects their effectiveness.


2021 ◽  
Vol 247 ◽  
pp. 04021
Author(s):  
Marton Szogradi

In order to meet modern industrial and scientific demands the Kraken multi-physics platform’s development was recently launched at VTT Technical Research Centre of Finland. The neutronic solver of the framework consists of two calculation chains, providing full core solutions by the Serpent high fidelity code (1) and the AFEN/FENM-based reduced-order diffusion solver called Ants (2) capable of handling square and hexagonal geometries in steady-state. Present work introduces the simulation of a large 3600 MWth Sodium-cooled Fast Reactor (SFR) described within the activities of the Working Party on Scientific Issues of Reactor Systems (WPRS) of OECD. Full-core 3D results were obtained by Serpent for carbide- and oxide-fuel cores, moreover group constants were generated for Ants utilizing 2D super-cell and single assembly infinite lattice models of Serpent. The continuous-energy Monte Carlo method provided the reference results for the verification of the reduced-order method. Implementing the spatially homogenized properties, 3D solutions were obtained by Ants as well for both core configurations. Comparison was made between the various core designs and codes based on reactivity feedbacks (Doppler constant, sodium voiding, control rod worth) considering power distributions. Regarding reactivity sensitivity on geometry, axial fuel- and radial core expansion coefficients were evaluated as well.


Author(s):  
Hao Luo ◽  
Mancang Li ◽  
Shanfang Huang ◽  
Minyun Liu ◽  
Kan Wang

2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
A. Rais ◽  
D. Siefman ◽  
G. Girardin ◽  
M. Hursin ◽  
A. Pautz

In order to analyze the steady state and transient behavior of the CROCUS reactor, several methods and models need to be developed in the areas of reactor physics, thermal-hydraulics, and multiphysics coupling. The long-term objectives of this project are to work towards the development of a modern method for the safety analysis of research reactors and to update the Final Safety Analysis Report of the CROCUS reactor. A first part of the paper deals with generation of a core simulator nuclear data library for the CROCUS reactor using the Serpent 2 Monte Carlo code and also with reactor core modeling using the PARCS code. PARCS eigenvalue, radial power distribution, and control rod reactivity worth results were benchmarked against Serpent 2 full-core model results. Using the Serpent 2 model as reference, PARCS eigenvalue predictions were within 240 pcm, radial power was within 3% in the central region of the core, and control rod reactivity worth was within 2%. A second part reviews the current methodology used for the safety analysis of the CROCUS reactor and presents the envisioned approach for the multiphysics modeling of the reactor.


2021 ◽  
Vol 247 ◽  
pp. 07019 ◽  
Author(s):  
Margaux Faucher ◽  
Davide Mancusi ◽  
Andrea Zoia

In this work, we present the first dynamic calculations performed with the Monte Carlo neutron transport code TRIPOLI-4R with thermal-hydraulics feedback. For this purpose, the Monte Carlo code was extended for multi-physics capabilities and coupled to the thermal-hydraulics subchannel code SUBCHANFLOW. As a test case for the verification of transient simulation capabilities, a 3x3-assembly mini-core benchmark based on the TMI-1 reactor is considered with a pin-by-pin description. Two reactivity excursion scenarios initiated by control-rod movement are simulated starting from a critical state and compared to analogous simulations performed using the Serpent 2 Monte-Carlo code. The time evolution of the neutron power, fuel temperature, coolant temperature and coolant density are analysed to assess the multi-physics capabilities of TRIPOLI-4. The stabilizing e_ects of thermal-hydraulics on the neutron power appear to be well taken into account. The computational requirements for massively parallel calculations are also discussed.


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.


Sign in / Sign up

Export Citation Format

Share Document