scholarly journals Methods and Models for the Coupled Neutronics and Thermal-Hydraulics Analysis of the CROCUS Reactor at EFPL

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.


2020 ◽  
Vol 239 ◽  
pp. 22011 ◽  
Author(s):  
Peng Hong Liem ◽  
Zuhair ◽  
Donny Hartanto

The results of criticality, sensitivity and uncertainty (S\U) analyses on the first core criticality of the Indonesian 30 MWth Multipurpose Reactor RSG GAS (MPR-30) using the recent nuclear data libraries (ENDF/B-VII.1 and JENDL-4.0) and analytical tools available at present (WHISPER-1.1) are presented. Two groups of criticality benchmark cases were carefully selected from the experiments conducted during the first criticality approach and control rod calibrations. The C/E values of effective neutron multiplication factor (k) for the worst case was found around 1.005. Large negative sensitivities were found in (n,e-mail:γ) reaction of H-1, U-235, Al-27, U-238 and Be-9 while large positive sensitivities were found in U-235 (total nu and fission), H-1 (elastic), Be-9 (free gas, elastic) and H-1 S(α,β) (lwtr.20t, inelastic). The S\U analysis results concluded that the uncertainties of k originated from the nuclear data were found around 0.6% which covered well the [C/E-1] values. Differences in the sensitivities amongst the two nuclear data libraries were also identified, and recommendation for improving the nuclear data library was given.



Author(s):  
Yuta Maruyama ◽  
Satoshi Imura ◽  
Junto Ogawa ◽  
Shuhei Miyake

Mitsubishi Heavy Industries (MHI) has developed the SPARKLE code, which is a PWR plant system transient analysis code that includes a three-dimensional (3D) neutronics module coupled with a thermal-hydraulics module. MHI has performed a study of the applicability of the SPARKLE code to the events which are associated with dynamic changes in power distribution, such as the rod ejection event or the steam line break event. In this paper, MHI has applied the SPARKLE code to the control rod drop event (drop of multiple rods), which features such a power distribution change. In addition, the neutron flux detection is dependent on the location of the dropped rods in this event, which can be dynamically calculated in the SPARKLE code. By applying the SPARKLE code to the control rod drop event, it was confirmed that the safety margin for this event is sufficiently larger than the margin calculated using the current safety analysis method, even if the appropriate conservative assumptions are made.



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):  
Masao Yamanaka

AbstractExcess reactivity and control rod worth are generally considered important reactor physics parameters for experimentally examining the neutron characteristics of criticality in a core, and for maintaining safe operation of the reactor core in terms of neutron multiplication in the core. For excess reactivity and control rod worth at KUCA, as well as at the Fast Critical Assembly in the Japan Atomic Energy Agency, special attention is given to analyzing the uncertainty induced by nuclear data libraries based on experimental data of criticality in representative cores (EE1 and E3 cores). Also, the effect of decreasing uncertainty on the accuracy of criticality is discussed in this study. At KUCA, experimental results are accumulated by measurements of excess reactivity and control rod worth. To evaluate the accuracy of experiments for benchmarks, the uncertainty originated from modeling of the core configuration should be discussed in addition to uncertainty induced by nuclear data, since the uncertainty from modeling has a potential to cover the eigenvalue bias more than uncertainty by nuclear data. Here, to investigate the uncertainty of criticality depending on the neutron spectrum of cores, it is very useful to analyze the reactivity of a large number of measurements in typical hard (EE1) and soft (E3) spectrum cores at KUCA.



2020 ◽  
Vol 239 ◽  
pp. 22006
Author(s):  
Donny Hartanto ◽  
Bassam Khuwaileh ◽  
Peng Hong Liem

This paper presents the benchmark evaluation of the new ENDF/B-VIII.0 nuclear library for the OECD/NEA Medium 1000 MWth Sodium-cooled Fast Reactor (SFR). There are 2 SFR cores: metallic fueled (MET-1000) and oxide fueled (MOX-1000). The continuous-energy Monte Carlo Serpent2 code was used as the calculation tool. Various nuclear libraries such as ENDF/B-VII.1 and JENDL-4.0 were included to be compared with the newest ENDF/B-VIII.0. The evaluated parameters are k,βeff, sodium void reactivity (∆ρNa), Doppler constant (∆ρDoppler), and control rod worth (∆ρCR).



2021 ◽  
Vol 2048 (1) ◽  
pp. 012029
Author(s):  
Suwoto ◽  
H Adrial ◽  
T Setiadipura ◽  
Zuhair ◽  
S Bakhri

Abstract One of the main critical issues on a nuclear reactor is safety and control system. The control rod worth plays an important role in the safety and control of nuclear reactors. The control rods worth calculation is used to specify the safety margin of the reactor. The main objective of this work is to investigate impact of different nuclear data libraries on calculating the control rod reactivity worth on small pebble bed reactor. Calculation of the control rod reactivity worth in small high temperature gas cooled reactor has been conducted using the Monte Carlo N-Particle 6 (MCNP6) code coupled with a different nuclear data library. Famous evaluated nuclear data libraries such as JENDL-40u, ENDF/B-VII.1 and JEFF-3.2 continuous cross section-energy data libraries were used. The overall calculation results of integral control rod worth show that the ENDF/B-VII.1, JENDL-40u and JEFF-3.2 files give values of - 17.814%☐k/k, -18.0204 %☐k/k and -18.0267%☐k/k, respectively. Calculations using ENDF/B-VII.1 give a slightly lower value than the others, while the JENDL-4.0u file gives results that are close to JEFF-3.2 file. The different nuclear data libraries have a relatively small impact on the control rod worth of small pebble bed reactor. Accurate prediction by simulation of control rod worth is very important for the safety operation of all reactor types, especially for new reactor designs.



Author(s):  
Tomáš Czakoj ◽  
Evžen Losa

Three-dimensional Monte Carlo code KENO-VI of SCALE-6.2.2 code system was applied for criticality calculation of the LR-0 reactor core. A central module placed in the center of the core was filled by graphite, lithium fluoride-beryllium fluoride (FLIBE), and lithium fluoride-sodium fluoride (FLINA) compounds. The multiplication factor was obtained for all cases using both ENDF/B-VII.0 and ENDF/B-VII.1 nuclear data libraries. Obtained results were compared with benchmark calculations in the MCNP6 using ENDF/B-VII.0 library. The results of KENO-VI calculations are found to be in good agreement with results obtained by the MCNP6. The discrepancies are typically within tens of pcm excluding the case with the FLINA filling. Sensitivities and uncertainties of the reference case with no filling were determined by a continuos-energy version of the TSUNAMI sequence of SCALE-6.2.2. The obtained uncertainty in multiplication factor due to the uncertainties in nuclear data is about 650 pcm with ENDF/B-VII.1.



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.



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.



Sign in / Sign up

Export Citation Format

Share Document