Burn up calculations for the Iranian miniature reactor: A reliable and safe research reactor

Presenting neutronic calculations pertaining to the Iranian miniature research reactor is the main goal of this article. This is a key to maintaining safe and reliable core operation. The following reactor core neutronic parameters were calculated: clean cold core excess reactivity (ρ_ex), control rod and shim worth, shut down margin (SDM), neutron flux distribution of the reactor core components, and reactivity feedback coefficients. Calculations for the fuel burnup and radionuclide inventory of the Iranian miniature neutron source reactor (MNSR), after 13 years of operational time, are carried out. Moreover, the amount of uranium burnup and produced plutonium, the concentrations and activities of the most important fission products, the actinide radionuclides accumulated, and the total radioactivity of the core are estimated. Flux distribution for both water and fuel temperature increases are calculated and changes of the central control rod position are investigated as well. Standard neutronic simulation codes WIMS-D4 and CITATION are employed for these studies. The input model was validated by the experimental data according to the final safety analysis report (FSAR) of the reactor. The total activity of the MNSR core is calculated including all radionuclides at the end of the core life and it is found to be equal to 1.3 × 10³Ci. Our investigation shows that the reactor is operating under safe and reliable conditions.     

Investigation of BeO as a reflector for the low power research reactor

Calculations of the fuel burnup, core excess reactivity, and the reactivity worths of the top beryllium shim plates for two reflector types (beryllium and beryllium oxide (BeO)) in the Miniature Neutron Source Reactor (MNSR) have been presented in this paper using the GETERA and MCNP4C codes. The results showed that the reactor infinity multiplication factors were 1.7030 and 1.6824, the core unadjusted excess reactivities were 31.9 and 5.0 mk, and the reactivity worths of the top beryllium shim plates were 22 and 19 mk for the BeO and Be reflectors respectively. Finally, using the beryllium oxide instead of the existing Be reflector in the MNSR reactor increased the core excess reactivity and reactor operation time.     

Time dependent burn-up and fission products inventory calculations in the discharged fuel of the Syrian MNSR

The Syrian Miniature Neutron Source Reactor (MNSR), a 30 kW, 89.8% HEU fueled (U-Al), went critical in March, 1996. By operating the reactor at nominal power for 2.5 h/day, the estimated core life is 10 years. This paper presents the results of fuel burn-up and depletion analysis of the MNSR fuel lattice using the ORIGEN 2 code. A one-group cross-section data base for the ORIGEN 2 computer code was developed for the Syrian MNSR research reactor. The ORIGEN 2 predicted burn-up dependent actinide compositions of MNSR spent fuel using the newly developed data base show a good agreement with the published results in the literature. In addition, the burn-up characteristics of MNSR spent fuel was analyzed with the new data base. Finally, to study the effect of burn-up on the reactivity, the microscopic cross-sections of the fission products calculated by the WlMS code (using the number densities of fission products generated by the ORIGEN 2 code as a function of burn-up time), were used as an input for the CITATION code calculations. The results contained in this paper could be used in performing criticality safety analysis and shielding calculations for the design of a spent fuel storage cask for the MNSR core.     

Calculations of fuel burn-up and radionuclide inventory in the syrian miniature neutron source reactor using the WIMSD4 code

Calculations of the fuel burn up and radionuclide inventory in the Miniature Neutron Source Reactor after 10 years (the reactor core expected life) of the reactor operating time are presented in this paper. The WIMSD4 code is used to generate the fuel group constants and the infinite multiplication factor versus the reactor operating time for 10, 20, and 30 kW operating power levels. The amounts of uranium burnt up and plutonium produced in the reactor core, the concentrations and radioactivities of the most important fission product and actinide radionuclides accumulated in the reactor core, and the total radioactivity of the reactor core are calculated using the WIMSD4 code as well.     

Burnup analysis and in-core fuel management study of the 3 MW TRIGA MARK II research reactor

The principal objective of this study is to formulate an effective optimal fuel management strategy for the TRIGA MARK II research reactor at AERE, Savar. The core management study has been performed by utilizing four basic types of information calculated for the reactor: criticality, power peaking, neutron flux and burnup calculation. This paper presents the results of the burnup calculations for TRIGA LEU fuel elements. The fuel element burnup for approximately 20 years of operation was calculated using the TRIGAP compute code. The calculation is performed in one-dimensional radial geometry in TRIGAP. Inter-comparison of TRIGAP results with other two calculations performed by MVP-BURN and MCNP4C-ORIGEN2.1 show very good agreement. Reshuffling at 20,000 MWh step provides the highest core lifetime of the reactor, which is 64,500 MWh. Besides, the study gives valuable insight into the behaviour of the reactor and will ensure better utilization and operation of the reactor in future.     

Assessment of fuel conversion from HEU to LEU in the Syrian MNSR reactor using the MCNP code

Assessment of fuel conversion from high enriched uranium (HEU) to low enriched uranium (LEU) fuel in the Syrian MNSR reactor was conducted in this paper. Three 3-D neutronic models for the Syrian MNSR reactor using the MCNP-4C code were developed to assess the possibility of fuel conversion from 89.87% HEU fuel (UAl₄–Al) to 19.75% LEU fuel (UO₂). The first model showed that 347 fuel rods with HEU fuel were required to obtain a reactor core with 5.17 mk unadjusted excess reactivity. The second model showed that only 200 LEU fuel rods distributed in the reactor core like the David star figure were required to obtain a reactor core with 4.85 mk unadjusted excess reactivity. The control rod worth using the LEU fuel was enhanced. Finally, the third model showed that distribution of 200 LEU fuel rods isotropically in the 10 circles of the reactor core failed to convert the fuel since the calculated core unadjusted excess reactivity for this model was 10.45 mk. This value was far beyond the reactor operation limits and highly exceeded the current MNSR core unadjusted excess reactivity (5.17 mk).     

Calculation of the power distribution in the fuel rods of the low power research reactor using the MCNP4C code

The Monte Carlo method, using the MCNP4C code, was used in this paper to calculate the power distribution in 3-D geometry in the fuel rods of the Syrian Miniature Neutron Source Reactor (MNSR). To normalize the MCNP4C result to the steady state nominal thermal power, the appropriate scaling factor was defined to calculate the power distribution precisely. The maximum power of the individual rod was found in the fuel ring number 2 and was found to be 105 W. The minimum power was found in the fuel ring number 9 and was 79.9 W. The total power in the total fuel rods was 30.9 kW. This result agrees very well with nominal power reported in the reactor safety analysis report which equals 30 kW. Finally, the peak power factors, which are defined as the ratios between the maximum to the average and the maximum to the minimum powers were calculated to be 1.18 and 1.31 respectively.     

Analysis of relative release rates of radionuclides from the RBMK-1500 reactor fuel elements

Specific activities (concentrations) of fission products (FP) and activation products in spent fuel elements of the RBMK-1500 reactor were calculated using SCALE 5 computer code. Different burnup (5.1–21.0 MWd/kg) fuel assemblies were experimentally investigated. Activities of radionuclides present in the coolant water of storage cases of defective fuel elements were experimentally measured and analyzed. Experimental results provide a basis for a quantitative analysis of radionuclide release from spent fuel of the RBMK-1500 reactor. Relative release rates of radionuclides from the fuel matrix were assessed based on a comparison of experimental results with theoretical calculations. On the basis of analysis results released fission and activation products can be divided into several groups according to their release rates from fuel; this can be generalized for radionuclides with similar chemical properties.     

Minimization of uranium load in Miniature Neutron Source Reactors

The uranium load in the Syrian MNSR is minimized by reducing the clad thickness to a standard value of 0.38 mm instead of 0.60 mm based on a 3-D model of the reactor that included all reactor components. More than 31 fuel rods are saved. The effects of the reduction of the fuel load in the core on both the reactor safety and performance, in relation to its use as a tool for Neutron Activation Analyses, are analyzed.     

An investigation on recycling the recovered uranium from electro-refining process in a CANDU reactor

Feasibility studies for recycling the recovered uranium from electro-refining process of pyroprocessing into a Canada Deuterium Uranium (CANDU) reactor have been carried out with a source term analysis code ORIGEN-S, a reactor lattice analysis code WIMS-AECL, and a Monte Carlo analysis code MCNPX. The uranium metal can be recovered in a solid cathode during an electro-refining process and has a form of a dendrite phase with about 99.99% expecting recovery purity. Considering some impurities of transuranic (TRU) elements and fission products in the recovered uranium, sensitivity calculations were also performed for the compositions of impurities. For a typical spent PWR fuel of 3.0 wt.% of uranium enrichment, 30 GWD/tU burnup and 10 years cooling, the recovered uranium exhibited an extended burnup up to 14 GWD/tU. And among the several safety parameters, the void reactivity at the equilibrium state was estimated 15 mk. Additionally, a simple sphere model was constructed to analyze surface dose rates with the Monte Carlo calculations. It was found that the recovered uranium from the spent PWR fuel by electro-refining process has a significant radioactivity depending on the impurities such as fission products.     

Analysis of fuel rod behavior under normal operating conditions in Super Fast Reactor

A Super Fast Reactor is a pressure-vessel type, fast spectrum supercritical water-cooled reactor (SCWR) that is presently researched in a Japanese project. A preliminary core has been designed with 1.59E+06 W/m³ of power density [1]. In order to ensure the fuel rod integrity, the fuel rod behaviors under the normal operating conditions are analyzed using FEMAXI-6 code. Three types of the limiting fuel rods, with the maximum cladding surface temperature (MCST), maximum power peak (MPP) and maximum discharge burnup (MDB), are chosen to cover all the fuel rods in the core. The power histories of these fuel rods are taken from the neutronics calculation results in the core design. The available design range of the fuel rod design parameters, such as the initial gas plenum pressure, gas plenum length, grain size and pellet-cladding gap size, are found out in order to satisfy the following design criteria: (1) Maximum fuel centerline temperature should be less than 1900 °C. (2) Maximum cladding stress in circumstance direction should be less than 100 MPa. (3) Pressure difference on the cladding should be less than 1/3 of buckling collapse pressure. (4) Compressive stress to yield strength ratio should be less than 0.2. (5) Cumulative damage fraction (CDF) on the cladding should be less than 1.0. Finally the improved fuel rod design is proposed.     

Effects of core excess reactivity and coolant average temperature on maximum operable time of NIRR-1 miniature neutron source reactor

We appraised in this study the effects of core excess reactivity and average coolant temperature on the operable time of the Nigeria Research Reactor-1 (NIRR-1), which is a miniature neutron source reactor (MNSR). The duration of the reactor operating time and fluence depletion under different operation mode as well as change in core excess reactivity with temperature coefficient was investigated over a period of five years. Our result shows that there is a strong dependence of reactor operating time on core excess reactivity and temperature coefficient. It was observed in 2004 that with a cold core excess reactivity of 3.77 mk, at full-power flux of 1.0 × 10¹² n cm⁻² s⁻¹ the reactor operated for 5 continues hours. At half-power flux of 0.5 × 10¹² n cm⁻² s⁻¹ and under the same excess reactivity condition, the reactor reaches 8 h of operation. However, re-measurements done in 2009 shows that excess reactivity of the reactor has reduced to 2.80 mk, the operable time at full flux dropped to 3.5 h while that of half-power became 7 h. We also investigated the reactor's energy consumption within the period under study and found to be much more in 2008 compared to the previous years. We infer that the amount of fluence consumed and the excessive reactor usage in 2008 has contributed significantly to the reduction of the reactor's excess reactivity in that year. The results obtained here revels that for an MNSR with a clean core excess reactivity between 3.5 mk and 4.0 mk, 5 and 8 h are the maximum operable times under full and half-power flux conditions, respectively. Negative deviation from these optimum times is therefore an indication of a drop in excess reactivity and the need for beryllium shims addition.     

Effect of burnup on the neutronic parameters of ITU TRIGA Mark II research reactor

This study aims to estimate burnup of the fuel elements for the Istanbul Technical University TRIGA Mark II Research and Training Reactor using a Monte Carlo-based burnup-depletion code. Effect of burnup on the core neutronic parameters, effective core multiplication factor, fast/epithermal/thermal neutron fluxes, and core-average neutron spectrum, and incoming neutron spectrum of the piercing beam port (PBP), is investigated at the Beginning of Life (BOL) and End of Life (EOL). Operational data peculiar to a selected operation sequence, which contains positions of CRs, power level of the reactor, material temperatures and latest core map, are used to determine the current fuel burnup of fuel elements at the time under consideration. A specific operation sequence is selected for the analysis. Furthermore, all control rods are considered fully withdrawn to assess the excess reactivity. Results are obtained using MONTEBURNS2 with ENDFB/V-II.1 neutron/photon library for a full power of 250 kW. Neutron cross-section libraries at the full-power operating temperatures are generated using NJOY. From the results, the calculated burnup values of the core at the sequence considered and EOL are found to be 420 MWh and 560 MWh, respectively. Remaining excess reactivity is calculated to be less than 0.3 $. It is observed that core average thermal neutron flux reduces by 1 % while the fast and epithermal neutron fluxes remain almost unchanged.     

Temperature and burnup reactivities and operational lifetime for the submersion-subcritical, safe space (S4) reactor

The He–Xe gas-cooled, S4 reactor has a sectored, Mo–14%Re solid core for avoidance of single point failures in reactor cooling and Closed Brayton Cycle (CBC) energy conversion. The reactor core is loaded with UN fuel and each of its three sectors is thermal-hydraulically coupled to a separate CBC loop and radiator panels. The solid core minimizes voids, and the BeO reflectors are designed to easily disassemble upon impact, ensuring that the bare S4 reactor is sufficiently subcriticial when submerged in wet sand or seawater and flooded with seawater, following a launch abort accident. Spectral shift absorber (SSA) additives in the core and thin SSA coatings on the outer surface of the core can also be used to ensure subcriticality in such an accident. This paper investigates the effects of various SSAs (Re, Ir, Eu-151, B-10 and Gd-155) on the temperature and burnup reactivity coefficients and the operating lifetime of the S4 reactor at a steady thermal power of 550 kW. The calculations of the burnup, reactivity feedback coefficient used a mixture of the top 10 light and top 10 heavy fission products plus Sm-149 and are performed for isothermal reactor core and reflector temperatures of 1200 and 900 K. In this fast spectrum space reactor, SSAs markedly increase fuel enrichment and decrease the burnup reactivity coefficient, but only slightly decrease the temperature, reactivity feedback coefficient. With no SSAs, the UN fuel enrichment is lowest (58.5 wt.%), the temperature and burnup reactivity coefficients are the highest (−0.2709 ¢/K and −1.3470 $/at.%), and the estimated operating lifetime is the shortest (7.6 years). The temperature and burnup reactivity coefficients decrease to −0.2649 ¢/K and −1.0230 $/at.%, and the operating lifetime increases to 8.3 years when rhenium additives are used. With europium-151 and gadolinium-155 additions, fuel enrichment (91.5 and 94 wt.%) and operating lifetime (9.9 and 9.8 years) are the highest and both the temperature reactivity feedback coefficient (−0.2382 and −0.2447 ¢/K) and the burnup reactivity coefficient (−0.9073 and −0.8502 $/at.%) are the lowest.     

Fuel management of the HTR-10 including the equilibrium state and the running-in phase

The mode of fuel management of the HTR-10 was studied, including the simulation of the fuel shuffling process and the measurement of the burnup of a fuel element. The prior consideration was the design of the equilibrium state. Based on this the fuel loading of the initial core and the fuel shuffling mode from the initial core through the running-in phase into the equilibrium state were studied. The code system VSOP was used for the physical layout of the HTR-10 at the equilibrium state and in the running-in phase. For the equilibrium state, in order to lessen the difference between the peak and the average burnup, 5-fuel-passage-through-the-core was chosen for the fuel management. The average burnup of the spent fuel for the equilibrium core is 80 000 MWd t⁻¹, and the peak value of it is less than 100 000 MWd t⁻¹ when the burnup of the recycled fuel element is under 72 000 MWd t⁻¹. The mixture of fuel element and graphite element was used for the initial core loading, the volume fractions of the fuel and the graphite elements were 0.57 and 0.43, respectively. During the running-in phase, the volume fraction of graphite will decrease with the fresh fuel elements being loaded from the top of the core and the graphite elements discharged from the bottom of the core. The fuel shuffling mode is similar to that of the equilibrium state. The burnup limit of recycled fuel element is also 72 000 MWd t⁻¹ and the peak burnup is less than 100 000 MWd t⁻¹. Finally the core will be full of fuel elements with a certain profile of burnup and reaches the equilibrium state. According to the characteristics of the pebble-bed high temperature gas-cooled reactor, a calibrating method of concentration of ¹³⁷Cs was proposed for the measurement of fuel burnup.     

Development and validation of a fast reactor core burnup code – FARCOB

A large fast breeder reactor requires the accurate estimation of power produced in different parts of the reactor core and blanket during any operating condition for a safe and economic operation through out reactor life time. A fast reactor core simulation code FARCOB based on multigroup diffusion theory has been developed in IGCAR for core simulation of PFBR reactor under construction. FARCOB uses centre mesh differencing scheme with triangular meshes in the X–Y plane. Steady state solution results match exactly with those of other reputed codes DIF3D and VENTURE for SNR-300 benchmarks. For burnup simulation, core is divided into radial and axial burnup zones and burnup equations are solved at constant power. Burnable fuel and blanket number densities are found and stored for each mesh, so that the user can alter burnup zones and core geometry after a burnup step. For validation, results of FARCOB has been compared with results of other institutes in two burnup benchmarks (ANL 1000 MWe benchmark and BN-600 hybrid core benchmark). It is found that FARCOB results match well with those of the other institutes.     

Design of a fuel element for a lead-cooled fast reactor

The options of a lead-cooled fast reactor (LFR) of the fourth generation (GEN-IV) reactor with the electric power of 600 MW are investigated in the ELSY Project. The fuel selection, design and optimization are important steps of the project. Three types of fuel are considered as candidates: highly enriched Pu-U mixed oxide (MOX) fuel for the first core, the MOX containing between 2.5% and 5.0% of the minor actinides (MA) for next core and Pu-U-MA nitride fuel as an advanced option. Reference fuel rods with claddings made of T91 ferrite-martensitic steel and two alternative fuel assembly designs (one uses a closed hexagonal wrapper and the other is an open square variant without wrapper) have been assessed. This study focuses on the core variant with the closed hexagonal fuel assemblies. Based on the neutronic parameters provided by Monte-Carlo modeling with MCNP5 and ALEPH codes, simulations have been carried out to assess the long-term thermal-mechanical behaviour of the hottest fuel rods. A modified version of the fuel performance code FEMAXI-SCK-1, adapted for fast neutron spectrum, new fuels, cladding materials and coolant, was utilized for these calculations. The obtained results show that the fuel rods can withstand more than four effective full power years under the normal operation conditions without pellet-cladding mechanical interaction (PCMI). In a variant with solid fuel pellets, a mild PCMI can appear during the fifth year, however, it remains at an acceptable level up to the end of operation when the peak fuel pellet burnup ∼80 MW d kg⁻¹ of heavy metal (HM) and the maximum clad damage of about 82 displacements per atom (dpa) are reached. Annular pellets permit to delay PCMI for about 1 year. Based on the results of this simulation, further steps are envisioned for the optimization of the fuel rod design, aiming at achieving the fuel burnup of 100 MW d kg⁻¹ of HM.     

Flow of kinetic parameters in a typical swimming pool type research reactor

Calculations were performed to estimate the variation in kinetic parameters (delayed neutron fraction and prompt neutron generation time) in different core configurations of a typical swimming pool type research reactor. Pakistan research Reactor-1 (PARR-1) was employed for this study. The effect due to burnup of the core was also studied. Calculations were performed with the help of computer codes WIMSD/4 and CITATION. Precursors yield was modified according to the neutron flux averaging only. This is the simple way to calculate the precursor yield for a particular core. The kinetic parameters are different for different core configurations. The β_eff decreases with 1.33 × 10⁻⁶/% burnup whereas prompt neutron generation time increases with 6.42 × 10⁻⁸ s/% burnup. The results were compared with safety analysis report and with published values and were found in good agreement. This study provides the confidence to understand the change in the kinetic parameters of research reactors with core change and also with burnup of the core.     

Monte Carlo simulation of a conceptual thermal column in the Syrian miniature neutron source reactor using MCNP-4C

Based on probabilistic approach, the MCNP-4C code has been used effectively to simulate the Syrian MNSR reactor core and all its surrounding components in three dimensions, including a preliminary conceptual design of a thermal column to be installed later. For verification and validation purposes, reactor calculations include: criticality and control rod worth. Values of these parameters are 1.00517 and 6.54 mk, respectively. The thermal column is to be installed in the water of the reactor pool. Optimal conditions for this thermal column were tested using the already developed model. Optimization focused on the most suitable position for placement of the column in the water pool, dimensions, and material. The aim was to have a thermal neutron flux of 1 × 10⁹ n cm⁻² s⁻¹ in the center of thermal column, and resonant and fast neutron fluxes to be as low as possible as well.     

Determination of major kinetic parameters of the Syrian MNSR for different fuel loading using Monte Carlo technique

A comprehensive 3-D model of the Syrian MNSR reactor has been developed using the MCNP-4C code aiming at accurate predicting of key core physics parameters. For the currently utilized HEU fuel (89.87% UAl4-Al) and two possible alternative LEU fuels (UO2 12%, and UO2 20%) the main core kinetics parameters like prompt neutron generation time, effective delayed neutron fraction, clean cold core excess reactivity and reactivity feedback coefficients of moderator temperature have been calculated. In this regard the role of particle weight loss on capture, fission and escape in determining the temperature effect of reactivity has been evaluated. The calculated results for the HEU fuel agree well with experimental values. The evaluated kinetics parameters are being used in accomplishing necessarily safety analyses related to the conversion of MNSR reactor to low enriched uranium.