Monte Carlo modelling of TRIGA research reactor
Mariam Zoubair
Tarek El Bardouni2010, Radiation Physics and Chemistry
Sign up for access to the world's latest research
checkGet notified about relevant papers
checkSave papers to use in your research
checkJoin the discussion with peers
checkTrack your impact
Abstract
The Moroccan 2 MW TRIGA MARK II research reactor at Centre des Etudes Nuclé aires de la Maâmora (CENM) achieved initial criticality on May 2, 2007. The reactor is designed to effectively implement the various fields of basic nuclear research, manpower training, and production of radioisotopes for their use in agriculture, industry, and medicine. This study deals with the neutronic analysis of the 2-MW TRIGA MARK II research reactor at CENM and validation of the results by comparisons with the experimental, operational, and available final safety analysis report (FSAR) values. The study was prepared in collaboration between the Laboratory of Radiation and Nuclear Systems (ERSN-LMR) from Faculty of Sciences of Tetuan (Morocco) and CENM. The 3-D continuous energy Monte Carlo code MCNP (version 5) was used to develop a versatile and accurate full model of the TRIGA core. The model represents in detailed all components of the core with literally no physical approximation. Continuous energy cross-section data from the more recent nuclear data evaluations (ENDF/B-VI.8, ENDF/B-VII.0, JEFF-3.1, and JENDL-3.3) as well as S(a, b) thermal neutron scattering functions distributed with the MCNP code were used. The cross-section libraries were generated by using the NJOY99 system updated to its more recent patch file ''up259''. The consistency and accuracy of both the Monte Carlo simulation and neutron transport physics were established by benchmarking the TRIGA experiments. Core excess reactivity, total and integral control rods worth as well as power peaking factors were used in the validation process. Results of calculations are analysed and discussed.
Figures (16)
Related papers
Benchmark analysis of the TRIGA MARK II research reactor using Monte Carlo techniques
Annals of Nuclear Energy, 2004
This study deals with the neutronic analysis of the current core configuration of a 3-MW TRIGA MARK II research reactor at Atomic Energy Research Establishment (AERE), Savar, Dhaka, Bangladesh and validation of the results by benchmarking with the experimental, operational and available Final Safety Analysis Report (FSAR) values. The 3-D continuousenergy Monte Carlo code MCNP4C was used to develop a versatile and accurate full-core model of the TRIGA core. The model represents in detail all components of the core with literally no physical approximation. All fresh fuel and control elements as well as the vicinity of the core were precisely described. Continuous energy cross-section data from ENDF/B-VI and ENDF/B-V and S(a; b) scattering functions from the ENDF/B-VI library were used. The consistency and accuracy of both the Monte Carlo simulation and neutron transport physics was established by benchmarking the TRIGA experiments. The effective multiplication factor, power distribution and peaking factors, neutron flux distribution, and reactivity experiments comprising control rod worths, critical rod height, excess reactivity and shutdown margin were used in the validation process. The MCNP predictions and the experimentally determined values are found to be in very good agreement, which indicates that the simulation of TRIGA reactor is treated adequately.
Neutronics analysis of the initial core of the TRIGA Mark II reactor
Nuclear Engineering and Design, 2011
The Atominstitute (ATI) of Vienna University of Technology (VUT) operates a TRIGA Mark II research reactor since March 1962. Its initial criticality was achieved on 7th March 1962 when 57th Fuel Element (FE) was loaded to the core. This paper describes the development of the MCNP model of the TRIGA reactor and its validation through three different experiments i.e. initial criticality, reactivity distribution and a thermal flux mapping experiment in the reactor core. All these experiments were performed on the initial core configuration. The MCNP model includes all necessary core components i.e. FE, Graphite Element GE, neutron Source Element (SE), Central IRradiation channel (CIR) etc. Outside the core, this model simulates the annular grooved graphite reflector, the thermal and thermalizing column, four beam tubes and the reactor water tank up to 100 cm in radial and +60 and −60 cm in axial direction. Each grid position at its exact location is modeled. This model employs the ENDF/B-VI data library except for the Sm-isotopes which are taken from JEFF 3.1 because ENDF/B-VI lacks samarium (Sm) cross sections.
Power and Neutron Flux Calculation for the Puspati Triga Reactor Using MCNP
The Malaysian 1 MW TRIGA MARK II research reactor at Malaysian Nuclear Agency is designed to effectively implement the various fields of basic nuclear research, manpower training, and production of radioisotopes for their use in agriculture, industry, and medicine. This study deals with the calculation of neutron flux and power distribution in PUSPATI TRIGA REACTOR (RTP) 14 th core configuration. The 3-D continuous energy Monte Carlo code MCNP was used to develop a versatile and accurate full model of the TRIGA core and fuels. The model represents in detailed all components of the core with literally no physical approximation. Continuous energy cross-section data from the more recent nuclear data as well as S (α, β) thermal neutron scattering functions distributed with the MCNP code were used. Results of calculations are analyzed and discussed.
A new model with Serpent for the first criticality benchmarks of the TRIGA Mark II reactor
Annals of Nuclear Energy
We present a new model, developed with the Serpent Monte Carlo code, for neutronics simulation of the TRIGA Mark II reactor of Pavia (Italy). The complete 3D geometry of the reactor core is implemented with high accuracy and detail, exploiting all the available information about geometry and materials. The Serpent model of the reactor is validated in the fresh fuel configuration, through a benchmark analysis of the first criticality experiments and control rods calibrations. The accuracy of simulations in reproducing the reactivity difference between the low power (10 W) and full power (250 kW) reactor condition is also tested. Finally, a direct comparison between Serpent and MCNP simulations of the same reactor configurations is presented.
Influence of geometry in TRIGA reactor criticality calculation and reactivity determination using Serpent 2 and MCNPX codes
International Journal of Nuclear Energy Science and Technology, 2017
The IPR-R1 TRIGA Mark I research reactor is located at the Nuclear Technology Development Centre (CDTN), in Belo Horizonte, Brazil. It is operating for more than 50 years and was successfully simulated before. However, new techniques and methods used in nuclear reactors analysis make a further simulation inevitable. In this manuscript, the computational model of an initial core of the IPR-R1 TRIGA reactor was developed employing two different Monte Carlo codes, MCNPX and Serpent 2, to simulate the neutronics behaviour. A new model is suggested, more complete, to improve the simulations results making the model more close the experimental data. This work explores how changes could be inserted in order to make the model closer to reality and if such participation would be noticeable in both codes used. The neutronic parameters obtained from these simulations performed in Serpent 2 are compared to MCNPX simulation results at the same conditions, and the results are compared with previous experimental data.
The RAPID Fission Matrix Approach to Reactor Core Criticality Calculations
Nuclear Science and Engineering, 2018
The Real-time Analysis for Particle transport and In-situ Detection (RAPID) code uses a unique, extremely fast, fission matrix-based methodology to compute the eigenvalue, and three-dimensional, pinwise fission source distribution for reactor, spent fuel pool, and spent fuel cask problems. In this paper, the RAPID fission matrix method is described and analyzed for application to several large pressurized water reactor problems, based on the Organisation for Economic Cooperation and Development/Nuclear Energy Agency Monte Carlo Performance Benchmark problem. In the RAPID methodology, fission matrix coefficients precalculated using the Serpent Monte Carlo code, are then coupled together and solved for different core arrangements. A boundary correction method is used to obtain more accurate fission matrix values near the radial and axial reflectors. Eigenvalues and fission source distributions are compared between RAPID and Serpent reference calculations. In most cases, the eigenvalue differences between methods are less than 10 pcm. For a uniform core model, pinwise fission distributions between the methods differ by a root-mean-square value of , 0:3%, compared to a Serpent uncertainty of , 0:22%. The pinwise, axially dependent (100 axial levels) differences are , 2:2%, compared to a similar Serpent uncertainty of , 2:0%. To achieve these levels of uncertainty, the RAPID calculations are over 2500 times faster than Serpent, not counting the precalculation time.
Neutronic Analysis of Moroccan TRIGA MARK-II Research Reactor using the DRAGON.v5 and TRIVAC.v5 codes
2018
The two-dimensional core neutronic analysis of the Moroccan TRIGA MARK-II research reactor using the DRAGON.v5 code allows us to obtain all databases that will be used in calculations of full TRIGA core in using TRIVAC.v5 code. In this study, the burnup calculations of two-dimensional fuel element and two-dimensional TRIGA core were performed using DRAGON.v5 code, and effective multiplication factor was calculated for full core via TRIVAC.v5 code. The obtained results indicate that the DRAGON.v5 and TRIVAC.v5 codes are capable of neutronic analysis of TRIGA MARK-II research reactor, provided that appropriate treatments are applied in both fuel element and core levels.
Validation of the Monte Carlo Code MVP on the First Criticality of Indonesian Multipurpose Reactor
The validation research works in BATAN are focused using Monte Carlo codes with recent nuclear data on the experimental results. In this paper, the validation results of Monte Carlo code MVP on the first criticality experimental of Indonesia Multipurpose Reactor (RSG GAS reactor) are presented. The MVP code is a continuous energy Monte Carlo code developed by Japan Atomic Energy Research Institute (JAERI). The objective this paper is to show the accuracy of the code using recent nuclear data of JEF-3.0, JENDL-3.3 and ENDF/B-VI.8. The final goal of this research is to use the code as an in-core fuel management code since the code has a module of burn-up calculation (MVP-BURN). The MVP calculations with the three libraries produced k eff values with excellent agreement to experiment data since the maximum differences are less than 0.5%. For the total control rod worth, the maximum difference is 3.6%. Systematically, ENDF/B-VI.8 library gave a maximum difference compared with other libraries. Therefore, the MVP code with recent libraries can be applied for the MTR type reactor with bulky Beryllium reflector.
Measurement and simulation of the neutron flux distribution in the TRIGA Mark II reactor core
Annals of Nuclear Energy, 2015
The neutron flux distribution in the TRIGA Mark II reactor core installed at the Applied Nuclear Energy Laboratory (L.E.N.A.) of the University of Pavia was measured through the neutron activation technique, irradiating Al-Co samples in different positions among the fuel elements and in the Central Thimble. The fast and integral fluxes were simultaneously measured in the same positions analyzing the 27 Al(n,a) 24 Na threshold reaction and the (n,c) activation of 59 Co, respectively. These measurements were then compared with the neutron fluxes evaluated through the MCNP5 reactor model which was developed in the last years, obtaining a good agreement of the results.
Validation of neutronic models and calculation systems by means of experimental results in the IPEN/MB-01 reactor
Annals of Nuclear Energy, 2013
In the year 2008 the National Atomic Energy Commission (CNEA) of Argentina, and the Brazilian Institute of Energetic and Nuclear Research (IPEN-CNEN/SP), under the framework of the Nuclear Energy Argentine-Brazilian Agreement (COBEN), among other projects, included ''Validation and Verification of Calculation Methods used for Research and Experimental Reactors''. At that time, it was established that the validation was to be performed with models implemented in the deterministic codes HUEMUL and PUMA (respectively, cell and reactor codes) developed by CNEA and those implemented in MCNP5. The necessary benchmark data for these validations would correspond to the theoretical-experimental reference cases elaborated in the research reactor IPEN/MB-01 located in the city of São Paulo, Brazil. These benchmarks were previously evaluated and approved for publications in the ICSBEP and IRPhE handbooks. The staff of the Nuclear Design and Analysis Division of the Reactor and Nuclear Power Plant Study Department (ERC) of CNEA, from the argentine side, modeled and performed several calculations with both deterministic (HUEMUL-PUMA) and probabilistic (MCNP5) methods of a great number of physical situations of the reactor, which previously have been studied experimentally and modeled by members of the Nuclear Engineering Center of IPEN, whose results were extensively provided to CNEA. The analyses reveal the great performances of ENDF/B-VII.0 in conjunction with MCNP5 and the HUEMUL and PUMA codes in all benchmark applications.
Related papers
Analysis of neutron flux distribution for the validation of computational methods for the optimization of research reactor utilization
2011
In order to verify and validate the computational methods for neutron flux calculation in TRIGA research reactor calculations, a series of experiments has been performed. The neutron activation method was used to verify the calculated neutron flux distribution in the TRIGA reactor. Aluminium (99.9 wt%)-Gold (0.1 wt%) foils (disks of 5 mm diameter and 0.2 mm thick) were irradiated in 33 locations; 6 in the core and 27 in the carrousel facility in the reflector. The experimental results were compared to the calculations performed with Monte Carlo code MCNP using detailed geometrical model of the reactor. The calculated and experimental normalized reaction rates in the core are in very good agreement for both isotopes indicating that the material and geometrical properties of the reactor core are modelled well. In conclusion one can state that our computational model describes very well the neutron flux and reaction rate distribution in the reactor core. In the reflector however, the accuracy of the epithermal and thermal neutron flux distribution and attenuation is lower, mainly due to lack of information about the material properties of the graphite reflector surrounding the core, but the differences between measurements and calculations are within 10%. Since our computational model properly describes the reactor core it can be used for calculations of reactor core parameters and for optimization of research reactor utilization.
The Method to Predict First Critical Core Loading for Nuclear Reactors
Sakarya University Journal of Science, 2020
Nuclear power plants have an important role in carbon free electricity production in the world. One of the important steps of commissioning a nuclear power plant is the first core loading. This is also called approaching the criticality. Since the number of fuel elements for the criticality is not known, precautions must be taken to prevent safety incidents. Although the procedure is performed on-line such that the neutron counts are measured at each loading of fuel elements to calculate sub-critical multiplication and the number of fuel element to reach criticality were predicted, computer simulations can also be used. In this study, inverse subcritical multiplication method was applied to Istanbul Technical University TRIGA Mark II research reactor first criticality in 1979 by using Monte Carlo simulation code MCNP6.2. Full 3-D model of the reactor was generated for calculations. Both results, experimental and simulation, showed that reactor became critical with 62 fuel elements. The core excess reactivity of 23.1 cents was predicted as 21.7 with the code. The simulation results are in good agreement with experimental results. The methodology and simulations can be used for power reactor analysis as well.
Initial Criticality Analysis of Malaysia’s TRIGA Research Reactor Using TRIGLAV Computer Code
International Journal of Engineering & Technology
Malaysian Nuclear Agency hosts the 1MW Thermal TRIGA MARK II research reactor since 1982. It first initial criticality was achieved on 28th June 1982 with loading of solid fuel elements like Uranium Zirconium Hydride. TRIGA MARK II started its operation of the research reactor on the same year with a 1MW power generation. Training, Research, Isotope Production and General Atomic (TRIGA) is designed to successfully actualize the variety fields of fundamental nuclear research, the manpower training and the production of radioisotopes. This study deals with the initial criticality analysis of the TRIGA research reactor using TRIGLAV reactor physics computer program. For this purpose, a model of its initial core will be developed and simulated using the software and the results will be validated against the experimental result as mentioned in the final safety analysis report (FSAR). The TRIGLAV computer code solves the neutron diffusion equation by using a finite differences method with...
Calculation of LEU Fuel Burn-up and Core Life Time Estimation of BAEC TRIGA Research Reactor Using 2D-TRIGLAV Code
International journal of science & technoledge, 2023
Calculation of LEU Fuel Burn-up and Core Life Time Estimation of BAEC TRIGA Research Reactor Using 2D-TRIGLAV Code 1. Introduction The TRIGA Mark-II research reactor [1] has a maximum thermal neutron flux of 7.46x10 13 n/cm 2 /sec in the middle of the core and is light water-cooled. This graphite-reflected nuclear reactor is intended for continuous operation at a steady state power level of 3 MWt. The TRIGA reactor LEU fuel is made up of burnable poison Erbium, zirconium hydride (primary moderator), and 20 weight percent uranium enriched to 19.7% 235 U. Boron carbide (B4C) serves as the neutron absorber component of the control rods. 100 fuel elements, including five fuelled follower control rods, six control rods, one air follower control rod, 18 graphite dummy elements, one central thimble, one pneumatic transfer system irradiation endpoint, and various light waters make up the BTRR LEU core. All of these components were positioned, supported, and arranged in seven concentric hexagonal rings (A, B, C, D, E, F, and G) of a hexagonal lattice, as shown in figure 1, between the top and bottom grid plates. It was given the go-ahead to carry out a number of nuclear research projects like neutron activation analysis, thermal neutron radiography, and neutron diffraction scattering experiments, as well as to train workers and create radioisotopes for application in medicine, industry, and agriculture. The term 'subcritical' describes a system where the loss of neutrons is greater than the rate of neutron creation [2], and the neutron population gradually declines over time. Criticality is a nuclear term that relates to the balance of neutrons in the core. A system is said to be 'supercritical' when neutron production outpaces neutron loss, increasing the population of neutrons [2, 3]. When the neutron populations are stable, the production and loss of neutrons are perfectly balanced, and the nuclear system is in a critical condition [3]. By comparing the pace at which neutrons are created from fission and other sources to the rate at which they are lost through absorption, scattering, and leakage out of the nuclear reactor core, it is possible to determine the criticality of a system [4]. The neutron diffusion theory code is used to perform this analysis. The configuration of the initial core shape and neutron energy group constants for various homogenized regions of the core, along with the fission spectrum, are required. In this study, neutron group constants were obtained using the well-known 1-D neutron transport code WIMS-D/4 [5] and were utilized for the full core calculations with the TRIGLAV code [6] that is based on a 4-group timeindependent diffusion equation in two-dimensional cylindrical (r,) configuration. The neutron diffusion algorithm is
Criticality qualification of a new Monte Carlo code for reactor core analysis
Annals of Nuclear Energy, 2009
In order to accurately simulate Accelerator Driven Systems (ADS), the utilization of at least two computational tools is necessary (the thermal-hydraulic problem is not considered in the frame of this work), namely: (a) A High Energy Physics (HEP) code system dealing with the ''Accelerator part" of the installation, i.e. the computation of the spectrum, intensity and spatial distribution of the neutrons source created by (p, n) reactions of a proton beam on a target and (b) a neutronics code system, handling the ''Reactor part" of the installation, i.e. criticality calculations, neutron transport, fuel burn-up and fission products evolution. In the present work, a single computational tool, aiming to analyze an ADS in its integrity and also able to perform core analysis for a conventional fission reactor, is proposed. The code is based on the well qualified HEP code GEANT (version 3), transformed to perform criticality calculations. The performance of the code is tested against two qualified neutronics code systems, the diffusion/transport SCALE-CITATION code system and the Monte Carlo TRIPOLI code, in the case of a research reactor core analysis. A satisfactory agreement was exhibited by the three codes.
NIGERIAN RESEARCH REACTOR-1 CORE NEUTRONICS CALCULATIONS - Azande T. S
Feasibility studies for the conversion of the Nigerian Research Reactor-1 (NIRR-1) have been performed using WIMS and CITATION codes at State. In this work, the core neutronics calculation of NIRR-1 concerning mass loading of U-235 in the core, shut down margin (SDM), safety reactivity factor (SRF), control rod worth, and control rod critical depth of insertion were investigated at low enrichment. Two fuel types (UAl 4 and UO 2) were considered and the uranium densities required for the conversion of NIRR-1 core to low enrichment were computed to be 1201g/cc with 20% enrichment, 1144g/cc with 19.75% enrichment, 1274g/cc with 15% enrichment, 1448g/cc with 10% enrichment for UAl 4 fuel type and 1141g/cc with 20% enrichment, 1144g/cc with 19.75% enrichment, 1216g/cc with 15% enrichment, and 1389g/cc with 10% enrichment for UO 2 fuel type. Significantly, higher uranium densities are required to convert NIRR-1 from HEU to LEU – indicating a drastic review of the NIRR-1 core. INTRODUCTION As part of the ongoing global effort to convert HEU reactor cores to Low Enriched Uranium (LEU) cores under the Reduced Enrichment for Research and Test Reactors (RERTR) program, there is need to study the possibilities of converting NIRR-1 core to less than 20% enrichment. Works have been done on Research and Test Reactor core conversion around the globe. In some proposed models, the total number of fuel pins (Khamis and Khattab, 1999) and the core radius/height ratio (Matos and Lell, 2005) has been drastically changed. This brings about noticeable changes in the relative flux values for both inner and outer irradiation sites. In this work, HEU and LEU cores are analyzed using the present UAl 4 fuel and a potential LEU fuel (UO 2 clad in zircalloy) are considered. The existing HEU core was also analyzed to validate the reactor model used. A significant feature of this work is the preservation of the technical and the geometric specification of the reactor so as to maintain the original designed thermal hydraulic of the reactor. A detailed description of the HEU core of NIRR-1 can be found in the Final Safety Analysis Report (FSAR); (Azande, et al, 2010; SAR, 2005).
Neutronics and thermal hydraulic analysis of TRIGA Mark II reactor using MCNPX and COOLOD-N2 computer code
Journal of Physics: Conference Series, 2017
The neutronic analysis of TRIGA Mark II reactor has been performed. A detailed model of the reactor core was conducted including standard fuel elements, fuel follower control rods, and irradiation devices. As the approach to safety nuclear design are based on determining the criticality (keff), reactivity worth, reactivity excess, hot rod power factor and power peaking of the reactor, the MCNPX code had been used to calculate the nuclear parameters for different core configuration designs. The thermal-hydraulic model has been developed using COOLOD-N2 for steady state, using the nuclear parameters and power distribution results from MCNPX calculation. The objective of the thermal-hydraulic model is to determine the thermal safety margin and to ensure that the fuel integrity is maintained during steady state as well as during abnormal condition at full power. The hot channel fuel centerline temperature, fuel surface temperature, cladding surface temperature, the departure from nucleate boiling (DNB) and DNB ratio were determined. The good agreement between experimental data and simulation concerning reactor criticality proves the reliability of the methodology of analysis from neutronic and thermal hydraulic perspective.
Calculation of Control Rod Worth of TRIGA Mark II Reactor Using Evaluated Nuclear Data Library JEFF-3 . 1 . 2
2017
The present work consists of a theoretical study of control rod worth of TRIGA Mark II research reactor using the evaluated nuclear data library JEFF-3.1.2. The data files of the evaluated nuclear data library cannot be used directly as input to neutronic or other applied calculations. These data files are to be converted to preprocessed files and then to post processed into multi-group files, which are then directed into specially formatted working libraries. These libraries are compatible with the neutronic codes. This subject involves computer software using knowledge of ENDF-6 formats, thermal effects, shelf-shielding factors in resolved and unresolved resonance regions, transfer matrices of various Legendre orders etc. The chain codes NJOY99.0, WIMSD-5B and CITATION are used to evaluate the control rod worth of six control rods of TRIGA Mark II research reactor. The obtained results of control rods worth are compared with experimental results of TRIGA at Atomic Energy Research ...
Study of Criticality Safety and Neutronic Performance for a 348-Fuel-Pin Ghana Research Reactor-1 LEU Core Using MCNP Code
World Journal of Nuclear Science and Technology, 2014
The National Nuclear Research Institute of the Ghana Atomic Energy Commission is undertaking steps to convert the Ghana Research Reactor-1 from HEU Core to LEU. The proposed LEU core consists of 12.5% enriched UO 2 fuel elements clad in Zircaloy-4 alloy. This is done in collaboration with Reduced Enrichment for Research and Test Reactor. The versatile MCNP code was used to analyse the neutronics parameters given in the SAR of HEU core, thereby characterizing the core. Subsequently, the LEU core was indentified with necessary changes to the HEU MCNP model. It was ascertained that the reactivity for the LEU core with the same number of fuel pins as the HEU was inadequate, hence the fuel pins were increased from 344 to 348. The neutron flux at the irradiation sites was found to be below the nominal value at full power for the LEU and hence the nominal power was increased to 34 kW for a nominal flux value of 1 × 10 12 n/cm 2 •s. The parameters investigated for the HEU and LEU are shown in this paper.
Neutronic analysis and validation of boiling water reactor core designed by MCNPX code
Annals of Nuclear Energy, 2015
This paper presents a design of boiling water reactor BWR model using MCNPX to develop benchmarks for checking the fuel management computer code packages. MCNPX code based on Monte Carlo method, is used to design a three dimensional model for BWR fuel assembly in typical operating temperature and pressure conditions. A test case was compared with a benchmark problem and good agreement was found. This design is used to study the thermal neutron flux and the pin by pin power distribution through the BWR core assemblies. The fuel used in BWR core is UO 2 with three different types of enrichment (0.711%, 1.76% and 2.19%). This enrichment is distributed in such a way as to flatten the power. The effect of different enrichment values on the radial normalized power distribution is analyzed. The spent fuel in the reactor can be recycled, and plutonium and its isotopes can be extracted.
Loading Preview
Sorry, preview is currently unavailable. You can download the paper by clicking the button above.