scholarly journals Numerical and experimental characterization of the reaction rates in the core of the CNESTEN’s TRIGA Mark II research reactor

2021 ◽  
Vol 253 ◽  
pp. 04006
Author(s):  
H. Ghninou ◽  
A. Gruel ◽  
A. Lyoussi ◽  
C. Reynard-Carette ◽  
C. El Younoussi ◽  
...  

Education, training and isotopes production are the most important uses of the Moroccan 2 MW TRIGA Mark II reactor situated at the National Center for Energy Sciences and Nuclear Techniques (CNESTEN, Morocco). To develop new R&D projects in research reactors, the particular and advanced knowledge of neutron and photon flux distribution, within and around the reactor core, is crucial. In order to precisely preparing the experiments in the CNESTEN’s TRIGA reactor, a detailed model was developed using the 3D continuous energy Monte Carlo code TRIPOLI-4 and the continuous energy cross-section data from the JEFF3.1.1 nuclear data library. This new model was used to carry out preliminary neutron and photon calculations to estimate flux levels in the irradiation channels as well as to calculate kinetic parameters of the reactor, core excess reactivity, integral control rods worth and power peaking factors. As a first step of the validation of the model, the obtained results were compared with the experimental ones available in the Final Safety Analysis Report (FSAR) of the TRIGA reactor. A study is being carried out at the end of which the results will be published as an evaluated benchmark. Furthermore, this work aims at experimentally characterize the reaction rates in various irradiation channels inside and outside the reactor core. The measurements are carried out using the neutron activation technique. To set up the experimental design for the activation experiments a series of preliminary calculations were performed using the TRIPOLI-4 model to calculate the expected gamma flux/intensity levels of various materials after irradiations in different positions in the irradiation facilities. Different activation foils with known characteristics are then irradiated and the activity of several isotopes is measured with the Gamma Spectrometry Method. The measured relative reaction rates are then compared with the calculated ones evaluated through the new TRIPOLI-4 reactor model. Fairly good agreement was found, which indicates that the new computational model is accurate enough to reproduce experiments.

2019 ◽  
Vol 2019 ◽  
pp. 1-6
Author(s):  
Toshio Wakabayashi ◽  
Makoto Takahashi ◽  
Naoyuki Takaki ◽  
Yoshiaki Tachi ◽  
Mari Yano

In a fast reactor, we evaluated a new core concept that prevents severe recriticality after whole-scale molten formation in a severe accident. A core concept in which Duplex pellets including neutron absorber are loaded in the outer core has been proposed. Analysis by the continuous energy model Monte Carlo code MVP using the JENDL-4.0 nuclear data library revealed that this fast reactor core has large negative reactivity due to fuel melting at the time of a severe accident, so that the core prevents recriticality. Regarding the core nuclear and thermal characteristics, the loading of Duplex pellets including neutron absorber in the outer core caused no significant differences from the normal core without Duplex pellets.


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.


Author(s):  
Wankui Yang ◽  
Baoxin Yuan ◽  
Songbao Zhang ◽  
Haibing Guo ◽  
Yaoguang Liu ◽  
...  

Deep penetration problems exist widely in reactor applications, such as SPRR300 (Swimming Pool Research Reactor 300), a light water moderated, enriched uranium fueled research reactor in China. Deterministic transport theory is intrinsically suitable for deep penetration. But there exist some problems when it’s applied in SPRR-300research reactors. First, the reactor core is complicated for geometry description in deterministic theory codes. Monte Carlo method has advantages in complex geometry modeling. And it uses continuous energy cross sections which are independent with specific reactor types and research objections. But usually it’s difficult to converge well enough to deal with deep penetration problems, even though there are a number of variance reduction techniques. Based on the advantages and disadvantages of Monte Carlo and Deterministic method, we proposed a coupled neutron transport calculation method for deep penetration. It combines advantages of these two methods. Firstly, we use Monte Carlo code to finish fine modeling and do the whole reactor core calculation. Domestically developed Monte Carlo code JMCT is used to do the neutron transport calculation. Then homogenized group constants in each mesh are calculated from JMCT output by a self-developed script. Afterwards, we do the whole reactor calculation with deterministic theory code TORT. It directly uses group constants generated by Monte Carlo code. Finally, we can get the deep penetration calculation results from TORT output. Verification is carried out by comparing the group constants of benchmark problem, and by comparing keff calculated by this method with continuous energy Monte Carlo method. Benchmark calculation is conducted with OECD/NEA slab benchmark problem. The comparison shows that group constants generated by this study are in good agreement with results from published references. Then above group constants are applied to 3-dimensional discrete ordinates deterministic theory transport code TORT. But keff calculated by TORT is a little lower than that calculated by Monte Carlo code JMCT. To minimize other influence factors, different Sn/Pn order, and different mesh size in TORT has been tried. Unfortunately the keff difference between these two methods remains. Even though the keff results in this benchmark are less than keff calculated by continuous energy MC method, Benchmark results show that all the group constants generated by this method are in good agreement with existing references. So it can be expected that after further verification and validation, this coupled method can be effectively applied to the deep penetration problem in such kind of research reactors.


2021 ◽  
Vol 247 ◽  
pp. 18001
Author(s):  
Andrej Žohar ◽  
Žiga Štancar ◽  
Paola Batistoni ◽  
Sean Conroy ◽  
Luka Snoj ◽  
...  

Fusion neutronics analysis before and after experiments at JET is traditionally performed using Monte Carlo particle transport code Monte Carlo N-Particle. For redundancy and diversity reasons there is a need of an additional Monte Carlo code, such as Serpent 2, capable of fusion neutronics analysis. In order to validate the Serpent code for fusion applications a detailed model of JET was used. Neutron fluxes and reaction rates were calculated and compared for positions outside the tokamak vacuum vessel, in the vacuum vessel above the plasma and next to a limiter inside the vacuum vessel. For all detector positions with DD and DT neutron sources the difference between neutron fluxes calculated with both Monte Carlo codes were within 2σ statistical uncertainty and for most positions (more than 90 % of all studied positions) even within 1σ uncertainty. Fusion neutronics analysis in the JET tokamak with Serpent took on average 10 % longer but this can be improved by changing the threshold value for determination of the transport method used. With the work presented in this paper the Serpent Monte Carlo code was validated to be a viable alternative to MCNP for fusion neutronics analysis for the JET tokamak.


Author(s):  
Weiyan Yang ◽  
Hongchun Wu ◽  
Liangzhi Cao

The radial power distribution within the fuel rod is important in fuel integrity evaluation. In this paper, the wavelets scaling function expansion method is applied to evaluate the radial power distribution and obtain continuous-energy spectrums in the fuel rod with a temperature distribution. Wavelets scaling function expansion continuous-energy self-shielding method is developed recently. It has been validated and verified by comparison to Monte Carlo calculations. In this method, the continuous-energy cross-sections are processed by NJOY, while the multi-group nuclear data library is jeff31 issued by IAEA. However, for different temperature problems, because of Doppler effect, the continuous-energy data library is different, and should be updated by applying NJOY. The calculation efficiency is a problem needed to be improved. Therefore, in this paper, the interpolation is utilized to obtain the continuous-energy cross-sections for other temperatures between those given temperatures. Also the precision of the temperature interpolation is discussed. Finally, the differences of continuous-energy spectrums, reaction rates and k-inf results are presented and compared with MCNP calculation.


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.


Sign in / Sign up

Export Citation Format

Share Document