A comparative neutronic study of the standard HEU core and various potential LEU alternatives for a typical MNSR system

A comparative study has been performed for neutronic analysis of highly enriched in uranium (HEU) and potential low enriched in uranium (LEU) cores for the Pakistan Research Reactor-2 (PARR-2) taken as a typical miniature neutron source reactor (MNSR) system. The group constant generation has been carried out using transport theory code WIMS-D4 and a detailed five-group RZ-model has been used in the CITATION code for multigroup diffusion theory analysis. The neutronic analysis of the 90% HEU reference and potential LEU alternative: UO₂, U₃Si₂ and U9Mo, cores has been carried out yielding 11%, 20.7% and 14.25% enrichments with corresponding values of excess reactivity: 4.33, 4.30 and 4.07 mk. These results have been found in good agreement with recently reported Monte Carlo-based transport theory calculations. The diffusion theory-based calculated values of thermal flux profiles for axial as well as for radial directions have been found to agree well with the corresponding experimental measurements. The UO₂-based LEU core has been found having flux spectrum closest to the reference core while U9Mo core has significantly harder flux spectrum at irradiation site.     

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.     

Monte Carlo simulation of core physics parameters of the Syrian MNSR reactor

A 3-D neutronic model for the Syrian Miniature Neutron Source Reactor (MNSR) was developed earlier to conduct the reactor neutronic analysis using the MCNP-4C code. The continuous energy neutron cross sections were evaluated from the ENDF/B-VI library. This model is used in this paper to calculate the following reactor core physics parameters: the clean cold core excess reactivity, calibration of the control rod and calculation its shut down margin, calibration of the top beryllium shim plate reflector, the axial neutron flux distributions in the inner and outer irradiation positions and calculations of the prompt neutron life time (_lp$l_p$) and the effective delayed neutron fraction (_βeff${\beta }_{\mathit{eff}}$). Good agreements are noticed between the calculated and the measured results. These agreements indicate that the established model is an accurate representation of Syrian MNSR core and will be used for other calculations in the future.     

原型微堆低浓化初步研究

利用蒙特卡罗计算程序,对高浓铀为燃料的原型微堆的有效增殖因数、控制棒价值、上铍反射层价值以及辐照座内的中子注量率等参数进行了计算。将计算值与实验结果进行了比较,两者基本相符。在原型微堆堆芯尺寸保持不变的情况下,将堆芯燃料元件芯体用富集度为12.5%UO₂替换UAl和用锆包壳替换铝包壳,对堆芯燃料低浓化方案进行了计算,给出了方案的计算结果。并利用RELAP5程序计算了原型微堆低浓铀堆芯阶跃引入4.0 mk反应性情况下反应堆的相关参数。     

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).     

The effect of temperature and control rod position on the spatial neutron flux distribution in the Syrian miniature neutron source reactor

The effect on the spatial neutron flux distribution for both of water and fuel temperature increase as well as the change in the control rod position are presented in the Syrian miniature neutron source reactor (MNSR). The cross-sections of all the reactor components at different temperatures are generated using the WIMSD4 code. These group constants are used then in the CITATION code to calculate the spatial neutron flux distribution at different water and fuel temperatures and different control rod positions using four energy groups. This work shows that the increase in water and fuel temperatures during the reactor daily operating time does not affect the spatial neutron flux distribution in the reactor. The change in the control rod position does not affect as well the spatial neutron flux distribution in the reactor except in the region around the control rod position.     

Calculation of fuel burnup and radionuclide inventory in the Syrian miniature neutron source reactor using the GETERA code

Calculations of the fuel burnup and radionuclide inventory in the Syrian miniature neutron source reactor (MNSR) after 10 years (the reactor core expected life) of the reactor operation time are presented in this paper using the GETERA code. The code is used to calculate 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 burntup and plutonium produced in the reactor core, the concentrations and radionuclides of the most important fission products and actinide radionuclides accumulated in the reactor core, and the total radioactivity of the reactor core were calculated using the GETERA code as well. It is found that the GETERA code is better than the WIMSD4 code for the fuel burnup calculation in the MNSR reactor since it is newer, has a bigger library of isotopes, and is more accurate.     

Computational analysis of neutronic parameters for TRIGA Mark-II research reactor using evaluated nuclear data libraries

The aim of this study is to analyze the neutronic parameters of TRIGA Mark-II research reactor using the chain of NJOY-WIMS-CITATION computer codes based on evaluated nuclear data libraries CENDL-2.2 and JEFF-3.1.1. The nuclear data processing code NJOY99.0 has been employed to generate the 69 group WIMS library for the isotopes of TRIGA core. The cell code WIMSD-5B was used to generate the cross sections in CITATION format and then 3-dimensional diffusion code CITTATION was used to calculate the neutronic parameters of the TRIGA Mark-II research reactor. All the analyses were performed using the 7-group macroscopic cross section library. The CITATION test-runs using different cross section sets based on different models applied in WIMS calculations have shown a strong influence of those models on the final integral parameters. Some of the cells were specially treated with PRIZE options available in WIMSD-5B to take into account the fine structure of the flux gradient in the fuel-reflector interface region. It was observed that two basic parameters, the effective multiplication factor, k_eff and the thermal neutron flux, were in good agreement among the calculated results with each other as well as the measured values. The maximum power densities at the hot spot were 1.0446E02 W/cc and 1.0426E02 W/cc for the libraries CENDL-2.2 and JEFF-3.1.1 respectively. The calculated total peaking factors 5.793 and 5.745 were compared to the original SAR value of 5.6325 as well as MCNP result. Consequently, this analysis will be helpful to enhance the neutronic calculations and also be used for the further thermal–hydraulics study of the TRIGA core.     

The behavior of reactor power and flux resulting from changes in core-coolant temperature for a miniature neutron source reactor

In this work, measurements were performed to verify the theoretical predictions of reactor power and flux parameters that result from changes in core inlet temperature (T_in) and the temperature difference between the coolant inlet and outlet (ΔT) in the Nigeria Research Reactor-1 (NIRR-1), which is a Miniature Neutron Source Reactor (MNSR). The measured data shows that there is a strong dependence of the reactor power on coolant temperature in agreement with the design of MNSR. The experimental parameters were found to be in good agreement with data obtained using a semi-empirical relationship between the reactor power, flux parameters, core inlet temperature, and the coolant temperature rise. The relationship was therefore used to predict the power level of NIRR-1 from its neutron flux parameters to which it has been found to be proportional. The variation of T_in and ΔT with the reactor power and flux was also investigated and the results obtained are hereby discussed.     

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.     

Design of a permanent epithermal neutron irradiation site in the Syrian MNSR

A permanent epithermal neutron irradiation site was designed in the Syrian Miniature Neutron Source Reactor (MNSR) by using cadmium as a thermal neutron shielding material. This site was designed by Cd-shielding the internal surface of the outer aluminum tube of the FOIS (First Outer Irradiation Site) in the MNSR. The MCNP-4C calculations showed that, to have a permanent epithermal neutron irradiation site for the ENAA using the cadmium, it is necessary to add the top beryllium shims of the reactor to compensate for the reactivity losses due to the neutrons absorption in the cylindrical cadmium shell. The activation detectors were used to measure the thermal and epithermal neutron fluxes in the FOIS. Distribution of the thermal neutron flux along the vertical direction of the outer irradiation capsule used in the FOIS has been found using MCNP-4C code, and experimentally by irradiating five copper wires. Good agreements were obtained between the calculated and the measured results.     

Neutronic analysis for core conversion (HEU–LEU) of Pakistan research reactor-2 (PARR-2)

Neutronic analyses for the core conversion of Pakistan research reactor-2 (PARR-2) from high enriched uranium (HEU) fuel to low enriched uranium (LEU) fuel has been performed. Neutronic model has been verified for 90.2% enriched HEU fuel (UAl₄–Al). For core conversion, UO₂ fuel was chosen as an appropriate fuel option because of higher uranium density. Clad has been changed from aluminum to zircalloy-4. Uranium enrichment of 12.6% has been optimized based on the design basis criterion of excess reactivity 4 mk in miniature neutron source reactor (MNSR). Lattice calculations for cross-section generation have been performed utilizing WIMS while core modeling was carried out employing three dimensions option of CITATION. Calculated neutronic parameters were compared for HEU and LEU fuels. Comparison shows that to get same thermal neutron flux at inner irradiation sites, reactor power has to be increased from 30 to 33 kW for LEU fuel. Reactivity coefficients calculations show that doppler and void coefficient values of LEU fuel are higher while moderator coefficient of HEU fuel is higher. It is concluded that from neutronic point of view LEU fuel UO₂ of 12.6% enrichment with zircalloy-4 clad is suitable to replace the existing HEU fuel provided that dimensions of fuel pin and total number of fuel pins are kept same as for HEU fuel.     

Design of the thermal neutron beam for neutron radiography at the Syrian MNSR

Neutron beam design was studied at the Syrian reactor (MNSR, 30 kW) with a view to generating thermal neutron beam in the vertical irradiation sites for neutron radiography. The design of the neutron collimator was performed using MCNP4C and the ENDF/B-V cross-section library. Thermal, epithermal and fast neutron energy ranges were selected as <0.4 eV, 0.4 eV–10 keV, >10 keV, respectively. To produce a good neutron beam quality, bismuth was used as photon filter. In this design, the L/D ratio of this facility had the value of 125. The thermal neutron flux at the beam exit was about 2.548 × 10⁵ n/cm² s. If such neutron beam were built into the Syrian MNSR many scientific applications would be available using the neutron radiography.     

Design calculation of an epithermal neutronic beam for BNCT at the Syrian MNSR using the MCNP4C code

This article describes the design calculation of an epithermal neutronic beam for the boron neutron capture therapy at the Syrian MNSR by using the MCNP4C code and ENDF/B-V cross-section library. To produce a high flux of epithermal neutrons at the beam exit, the moderator/filter from Al, Cd, Fluental and Bi was used with Pb as reflector for neutrons along the beam. In addition, the Bi lined collimator with Li₂CO₃-PE and Pb at the end. The calculated beam parameters under 30.0 kW of reactor power at the beam exit are Ф_epi = 2.83 × 10⁸ n/cm² s, D_f/Ф_epi = 7.98 × 10⁻¹¹ cGy cm²/n, D_γ/Ф_epi = 1.70 × 10⁻¹¹ cGy cm²/n, Φ_epi/Φ_the = 0.05 and J_n⁺/Ф_n = 0.77. As well as, the calculated values of the advantage depth and advantage ratio are 7.51 cm and 3.49, respectively. If such beam was built into the Syrian MNSR the scientific applications of the reactor would increase.     

3D neutron diffusion computational code based on GFEM with unstructured tetrahedron elements: A comparative study for linear and quadratic approximations

In the present study, the comparison between the results obtained from the linear and quadratic approximations of the Galerkin Finite Element Method (GFEM) for neutronic reactor core calculation was reported. The sensitivity analysis of the calculated neutron multiplication factor, neutron flux and power distributions in the reactor core vs. the number of the unstructured tetrahedron elements and order of the considered shape function was performed. The cost of the performed calculation using linear and quadratic approximation was compared through the calculation of the FOM. The neutronic core calculation was performed for both rectangular and hexagonal geometries. Both the criticality and fixed source calculations were done using the developed GFEM-3D computational code. An acceptable accuracy with low computational cost is the main advantage of applying the unstructured tetrahedron elements. The generated unstructured tetrahedron elements with Gambit software were used in the GFEM-3D computational code via a developed interface. The criticality calculation was benchmarked against the valid data for IAEA-3D and VVER-1000 benchmark problems. Also, the neutron fixed source calculation was validated through the comparison with the similar computational code. The results show that the accuracy of the calculation for the both linear and quadratic approximations improves vs. the number of elements. Quadratic approximation gives acceptable results for almost all considered number of the elements, while the results obtained from the linear approximation have good accuracy for only high number of the elements.     

Rejuvenation of spent fuel in the force-free helical reactor (FFHR)

This study considers spent fuel rejuvenation potential of the force-free helical reactor (FFHR), which is a relevant heliotron-type D-T fusion reactor. For this purpose, three different spent fuels were selected, i.e.: (1) CANDU, (2) PWR-UO₂, and (3) PWR-MOX spent fuels. The spent fuel (volume fraction of 60%), spherically prepared and cladded with SiC (volume fraction of 10%), was located in the fuel zone (FZ) in the blanket of the FFHR. The FZ was cooled with high-pressured helium gas (volume fraction of 30%) for the nuclear heat transfer. The neutronic calculations were performed by solving the Boltzmann transport equation with the help of the neutron transport code XSDRNPM-S/SCALE4.3. The calculations of the time dependent atomic densities of the isotopes were performed for an operation period (OP) of up to 4 years by 75% plant factor (η) under a first-wall neutron load (P) of 1.5 MW/m². The temporal variations of the atomic densities of the isotopes in the spent fuel composition and other physical parameters were calculated for a discrete time interval (Δt) of 1/12 year (one month), by using the interface program (code). In all investigated spent fuel cases, the tritium self-sufficiency is maintained for DT fusion driver along the OP. The CANDU spent fuel becomes usable in a conventional CANDU reactor after a regeneration time of ∼6 months. The CFFE value approaches 3.5% in the blanket fuelled with the PWR-UO₂ and -MOX spent fuels after 40 and 34 months, respectively. The plutonium component can never reach a nuclear weapon grade quality during spent fuel rejuvenation. Consequently, the blanket has high neutronic performance for the rejuvenation of the spent fuels.     

Tokamak D-T neutron source models for different plasma physics confinement modes

Neutronic studies of European demonstration fusion power plant (DEMO) have been so far based on plasma physics low confinement mode (L-mode). Future tokamaks, nevertheless, may likely use alternative confinement modes such as high or advanced confinement modes (H&A-mode). Based on analytical formulae used in plasma physics, H&A-modes D-T neutron sources formulae are proposed in this paper. For that purpose, a tokamak random neutron source generator, TRANSGEN, has been built generating bidimensional (radial and poloidal) neutron source maps to be used as input for neutronics Monte-Carlo codes (TRIPOLI-4 and MCNP5). The impact of such a source on the neutronic behavior of the European DEMO-2007 Helium-cooled lithium–lead reactor concept has been assessed and compared with previous results obtained using a L-mode neutron source. An A-mode neutron source map from TRANSGEN has been used with the code TRIPOLI-4. Assuming the same fusion power, results show that main reactor global neutronic parameters, e.g. tritium breeding ratio and neutron multiplication factor, evolved slightly when compared to present uncertainties margin. However, local parameters, such as the neutron wall loading (NWL), change significantly compared to L-mode shape: from ?22% to +11% for NWL.     

Calculation of the effective delayed neutron fraction and the prompt neutron lifetime with the indigenously modified CEBIS code: MCEBIS

The CEBIS code has been modified to enable the calcination of both the effective delayed neutron fraction and prompt neutron generation lifetime in any nuclear thermal reactor, especially reactor types such as TRIGA, SLOWPOKE, and MNSR. The new version, called MCEBIS, includes sonie special subroutines which will be called up as part of the input to calculate the above two dynamic parameters. In addition, some control flags have been added to recognize any important reactor components such as beryllium as a reflector or heavy water as moderator and pence calculate their photo-neutron fractions.

The MCEBIS code has been tested using two reactor models: TRIGA and MNSR. These models were developed mainly to verify the modified code. Each model represents a 1-D neutronics model of the reactor. Calculated results for the effective delayed neutron fraction and prompt neutron generation lifetime in both reactors have been compared with published data. Good agreement with published results has been established.     

Monte Carlo calculation of the core of the heavy water zero power reactor (HWZPR) using MCNP4c

The MCNP4c code, based on the probabilistic approach, was used to simulate 3D configuration of the core of the heavy water zero power reactor (HWZPR). In present work, first, all of the constituents of the core such as fuel pellets, fuel element, moderator (D₂O) and annular graphite reflector were modeled using MCNP4c code. Then calculations of axial and radial neutron fluxes were performed in three energy groups such as thermal (0-0.625 eV), epithermal (0.625-550 eV), and fast (0.550-20 MeV). The cadmium ratio was calculated as well and the neutron flux parameters such as extrapolated height (H_e), extrapolated radius (R_e) and physical center of the core (z₀) were computed using cadmium ratio. Comparison of the neutron flux parameters with the experimental data showed that the MCNP4c model of the HWZPR was validated.