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.     

Investigation on spent fuel characteristics of reduced-moderation water reactor (RMWR)

The spent fuel characteristics of the reduced-moderation water reactor (RMWR) have been investigated using the SWAT and ORIGEN codes. RMWR is an advanced LWR concept for plutonium recycling by using the MOX fuel. In the code calculation, the ORIGEN libraries such as one-group cross-section data prepared for RMWR were necessary. Since there were no open libraries for RMWR, they were produced in this study by using the SWAT code. New libraries based on the heterogeneous core modeling in the axial direction and with the variable actinide cross-section (VXSEC) option were produced and selected as the representative ORIGEN libraries for RMWR. In order to investigate the characteristics of the RMWR spent fuel, the decay heat, the radioactivity and the content of each nuclide were evaluated with ORIGEN using these libraries. In this study, the spent fuel characteristics of other types of reactors, such as PWR, BWR, high burn-up PWR, full-MOX-PWR, full-MOX-BWR and FBR, were also evaluated with ORIGEN.

It has been found that about a half of the decay heat of the RMWR spent fuel comes from the actinides nuclides. It is the same with the radioactivity. The decay heat and the radioactivity of the RMWR spent fuel are lower than those of full-MOX-LWRs and FBR, and are the same level as those of the high burn-up PWR. The decay heat and the radioactivity from the fission products (FPs) in the spent fuel mainly depend on the burn-up and the burn-up time rather than the reactor type. Therefore, the decay heat and the radioactivity from FPs in the RMWR spent fuel are smaller, reflecting its relatively long burn-up time resulted from its core characteristics with the high conversion ratio. The radioactivity from the actinides in the spent fuel mainly depends on the ²⁴¹Pu content in the initial fuel, and the decay heat mainly depends on ²³⁸Pu and ²⁴⁴Cm. The contribution of ²⁴⁴Cm is much smaller in RMWR than in MOX-LWRs because of the difference in the spectrum. In addition, from the waste disposal point of view, the characteristics of the heat generation FP elements, the platinum group metals, Mo and the long-lived FPs (LLFPs) were also investigated.     

TRIGA fuel burn-up calculations and its confirmation

The Cesium (Cs-137) isotopic concentration due to irradiation of TRIGA Fuel Elements FE(s) is calculated and measured at the Atominstitute (ATI) of Vienna University of Technology (VUT). The Cs-137 isotope, as proved burn-up indicator, was applied to determine the burn-up of the TRIGA Mark II research reactor FE. This article presents the calculations and measurements of the Cs-137 isotope and its relevant burn-up of six selected Spent Fuel Elements SPE(s). High-resolution gamma-ray spectroscopy based non-destructive method is employed to measure spent fuel parameters. By the employment of this method, the axial distribution of Cesium-137 for six SPE(s) is measured, resulting in the axial burn-up profiles. Knowing the exact irradiation history and material isotopic inventory of an irradiated FE, six SPE(s) are selected for on-site gamma scanning using a special shielded scanning device developed at the ATI. This unique fuel inspection unit allows to scan each millimeter of the FE. For this purpose, each selected FE was transferred to the fuel inspection unit using the standard fuel transfer cask. Each FE was scanned at a scale of 1 cm of its active length and the Cs-137 activity was determined as proved burn-up indicator. The measuring system consists of a high-purity germanium detector (HPGe) together with suitable fast electronics and on-line PC data acquisition module. The absolute activity of each centimeter of the FE was measured and compared with reactor physics calculations. The ORIGEN2, a one-group depletion and radioactive decay computer code, was applied to calculate the activity of the Cs-137 and the burn-up of selected SPE. The deviation between calculations and measurements was in range from 0.82% to 12.64%.     

A model for fission gas release and gaseous swelling of the uranium dioxide fuel coupled with the FALCON code

A comprehensive model GRSW-A was developed to analyse the processes of fission gas release, gaseous swelling and microstructural evolutions in the uranium dioxide fuel during base irradiation and under transient conditions. The GRSW-A analysis incorporates a number of models published in open literature, as well as some original models that were already published by the authors elsewhere. Consequently, only the most prominent aspects of GRSW-A and its coupling with the FALCON fuel behaviour analysis and licensing code are described in this paper. The analysis of fuel behaviour in the REGATE experiment is presented, which includes the base irradiation of the fuel segment in a PWR to a burn-up of about 50 MWd/kgU, which was followed by a power ramp in the SILOE research reactor. Besides, the generalized data on fission gas release (FGR) in PWR fuel during the base irradiation up to a burn-up of about 70 MWd/kgU is interpreted using coupled FALCON and GRSW-A. Moreover, a mechanistic interpretation of the published data for pellet swelling during the base irradiation up to a burn-up of 100 MWd/kgU is put forward. In all the cases, the coupled FALCON/GRSW-A analysis has shown the improved prediction capability compared to the original FALCON MOD01, which is achieved due to the account for the mutual effect of thermal and, in particular, high-burn-up-assisted mechanisms of fission gas release and swelling under steady-state and transient conditions.     

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.     

Post-shutdown decay power and radionuclide inventories in the discharged fuels of HEU and potential LEU miniature neutron source reactors

Assessment of fission product and actinide content along with the time variation of decay power of discharged fuels of both HEU and LEU cores of MNSRs have been carried out for once-through cycle using the ORIGEN2 computer code. The results for the LEU core have been compared with the corresponding values for the current HEU core of MNSRs. For the HEU and the potential LEU UO₂, U-9Mo discharged fuels, the ORIGEN2 computed isotopic and total activity values have been found in good agreement with the corresponding results obtained by using the WIMSD4 code. All three MNSR fuels show fission product dominated activity behavior for post-shutdown periods up to about 10³ years during which, the total activity decreases by as much as 10⁶ times. The residual actinide activity shows smaller variations as the three discharged fuels decay thru 10⁶ years. The time variation of the decay power follows the same behavior as the corresponding total activity values during the fission product dominated period. A decrease from initial values of 154.76, 162.6,160.39 W to the final values 9.35 × 10⁻⁵, 2.1 × 10⁻³, 1.7 × 10⁻³ W has been found for the standard HEU, and potential UO₂, U-9Mo LEU fuels correspondingly during this time. The standard HEU fuel shows smallest decay power values while the UO₂ and U-9Mo LEU fuels have comparable values for time spans from 10³ to about 10⁶ years.     

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

Thermal hydraulic and safety analysis for core conversion (HEU–LEU) of Syrian Miniature Neutron Source Reactor

The paper presents the behavior and properties analysis of the low enriched uranium fuel compared with the original high enriched uranium fuel. The MNSR reactor core was modeled with both fuel materials and the reactor behavior was studied during the steady state and abnormal conditions. The MERSAT code was used in the analysis. The steady state thermal hydraulic analysis results were compared with that obtained from the experimental results hold during commissioning the Syrian MNSR. Comparison with experimental data shows that the steady-state behavior of the HEU core was accurately predicted by the MERSAT code calculations. The validated model was then used to analyze LEU cores with two proposed UO₂ fuel pin designs. With each LEU core, the steady state and 3.77 mk rod withdrawal transient were run and the results were compared with the available published data in the literatures for the low enriched uranium fuel core. The results reveal that the low enriched uranium fuel showed a good behavior and the peak clad temperatures remain well below the clad melting temperature during reactivity insertion accident.     

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.     

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.     

Effect of fission products yield library on the calculation of the nuclide composition of the fuel in sodium-cooled fast reactors

The calculation of the composition of irradiated fuel for different degrees of burnup is a basic problem in the analysis of nuclear-radiological safety of objects holding spent fuel assemblies. The yield of fission products is one of the important initial indicators in burnup calculations. Methods for compiling libraries of fission products yield on the basis of the ENDF/B up-to-date evaluated nuclear data files are described. The nuclide composition of uranium oxide and uranium-plutonium-zirconium metal fuel in sodium-cooled fast reactors is analyzed by means of high-precision calculations performed with different fission product yields libraries using different computer codes MONTEBURNS–MCNP5–ORIGEN2 and the results are presented.     

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.     

基于蒙特卡罗方法的三维燃耗计算研究

采用通过编写连接MCNP程序和ORIGEN2程序的接口处理程序的方法进行快中子系统的燃耗计算。由MCNP、ORIGEN2、接口处理程序和截面文件组成的软件系统可用于燃料或堆芯非均匀布置快中子系统的燃料同位素成分和燃耗反应性损失计算,在燃耗反应性损失计算中采用了伪裂变产物的方法。介绍程序系统的研制情况,并给出用该软件系统计算中国实验快堆首炉堆芯和OECD/NEAMOX燃料快堆基准题的燃耗计算结果。     

Calculation of the activity of fission products in the primary coolant of the eastern-type pressurized water reactor (VVER1000-V446) of the Bushehr Nuclear Power Plant at normal full power operational condition

Calculation of the primary circuit's coolant activation due to fission products (FPs) has been investigated for the eastern-type pressurized water reactor (VVER1000-V446). The reactor has been considered under normal full power operational condition for the first fuel cycle. Determination of the reactor coolant activity is based on time-dependent fission product core inventories. ORIGEN2.1 code has been used to determine the time-dependent fission product core inventories. The fission products activity in the primary coolant is calculated using a set of ordinary differential equations (ODEs) which governs the FPs concentration in the primary coolant. Results for 87 FPs have been calculated. The results of these calculations have been found to agree well with the corresponding available values found in the Final Safety Analysis Report (FSAR) of the Bushehr Nuclear Power Plant (BNPP).     

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.     

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.     

Determination of 79Se and 135Cs in Spent Nuclear Fuel for Inventory Estimation of High-Level Radioactive Wastes

⁷⁹Se and ¹³⁵Cs are long-lived fission products and are found in high-level radioactive waste (HLW). The estimation of their inventories in HLW is essential for the safety assessment of geological disposal, owing to their mobility in the strata. In this study, the amounts of ⁷⁹Se and ¹³⁵Cs in a spent nuclear fuel solution were measured. About 5 g of irradiated UO₂ fuel discharged from a commercial Japanese pressurized water reactor (PWR) with a burn-up of 44.9 GWd/t was sampled and dissolved with 50mL of 4M nitric acid in a hot cell for 2 h. After Se and Cs were chemically separated, the amounts of ⁷⁹Se and ¹³⁵Cs in the spent nuclear fuel solution were measured by inductively coupled plasma quadrupole mass spectrometry (ICP-QMS). The amounts of ⁷⁹Se and 135Cs were 5:2 ± 1:5 and 447 ± 40 g/MTU, respectively. The results presented in this study, which are the first postirradiation experimental data in Japan, showed good agreement with those obtained by the ORIGEN2 code using the data library of JENDL-3.3.     

Modeling of WWER-440 fuel pin behavior at extended burn-up

Currently, there is an ongoing effort to increase fuel discharge burn-up of all LWRs fuel including WWERs as much as possible in order to decrease power production cost. Therefore, burn-up is expected to be increased from 60 to 70 MWd/kg U. The change in the fuel radial power distribution as a function of fuel burn-up can affect the radial fuel temperature distribution as well as the fuel microstructure in the fuel pellet rim. Both of these features, commonly termed the “rim effect.” High burn-up phenomena in WWER-440 UO₂ fuel pin, which are important for fission gas release (FGR) were modeled. The radial burn-up as a function of the pellet radius and enrichment has to be known to determine the local thermal conductivity.In this paper, the radial burn-up and fissile products distributions of WWER-440 UO₂ fuel pin were evaluated using MCNP4B and ORIGEN2 codes. The impact of the thermal conductivity on predicted FGR calculations is needed. For the analysis, a typical WWER-440 fuel pin and surrounding water moderator are considered in a hexagonal pin well. The thermal release and the athermal release from the pellet rim were modeled separately. The fraction of the rim structure and the excessive porosity in the rim structure in isothermal irradiation as a function of the fuel burn-up was predicted. A computer program; RIMSC-01, is developed to perform the required FGR calculations. Finally, the relevant phenomena and the corresponding models together with their validation are presented.     

Extension of the TRANSURANUS burn-up model

The validation range of the model in the TRANSURANUS fuel performance code for calculating the radial power density and burn-up in UO₂ fuel has been extended from 64 MWd/kgHM up to 102 MWd/kgHM, thereby improving also its precision. In addition, the first verification of calculations with post-irradiation examination data is reported for LWR-MOX fuel with a rod average burn-up up to 45 MWd/kgHM. The extension covers the inclusion of new isotopes in order to account for the production of ²³⁸Pu. The corresponding one-group cross-sections used in the equations rely on results obtained with ALEPH, a new Monte Carlo burn-up code. The experimental verification is based on electron probe microanalysis (EPMA) and on secondary ion mass spectrometry (SIMS) as well as radiochemical data of fuel irradiated in commercial power plants. The deviations are quantified in terms of frequency distributions of the relative errors. The relative errors on the burn-up distributions in both fuel types remain below 12%, corresponding to the experimental scatter.