• Nuclear Engineering and Technology
• /
• v.55 no.3
• /
• pp.1097-1104
• /
• 2023

A preliminary criticality analysis for novel pyrochemical apparatuses for the reprocessing of mixed uranium-plutonium nitride spent nuclear fuel from the BREST-OD-300 reactor was performed. High-temperature processing apparatuses, "metallization" electrolyzer, refinery remelting apparatus, refining electrolyzer, and "soft" chlorination apparatus are considered in this work. Computational models of apparatuses for two neutron radiation transport codes (MCU-FR and MCNP) were developed and calculations for criticality were completed using the Monte Carlo method. The criticality analysis was performed for different loads of fissile material into the apparatuses including overloading conditions. Various emergency situations were considered, in particular, those associated with water ingress into the chamber of the refinery remelting apparatus. It was revealed that for all the considered computational models nuclear safety rules are satisfied.

• Nuclear Engineering and Technology
• /
• v.52 no.11
• /
• pp.2452-2459
• /
• 2020

Neutron spectrum is essential to the safe operation of reactors. Traditional online neutron spectrum measurement methods still have room to improve accuracy for the application cases of wide energy range. From the application of artificial neural network (ANN) algorithm in spectrum unfolding, its accuracy is difficult to be improved for lacking of enough effective training data. In this paper, an adaptive deviation-resistant neutron spectrum unfolding method based on transfer learning was developed. The model of ANN was trained with thousands of neutron spectra generated with Monte Carlo transport calculation to construct a coarse-grained unfolded spectrum. In order to improve the accuracy of the unfolded spectrum, results of the previous ANN model combined with some specific eigenvalues of the current system were put into the dataset for training the deeper ANN model, and fine-grained unfolded spectrum could be achieved through the deeper ANN model. The method could realize accurate spectrum unfolding while maintaining universality, combined with detectors covering wide energy range, it could improve the accuracy of spectrum measurement methods for wide energy range. This method was verified with a fast neutron reactor BN-600. The mean square error (MSE), average relative deviation (ARD) and spectrum quality (Qs) were selected to evaluate the final results and they all demonstrated that the developed method was much more precise than traditional spectrum unfolding methods.

• Journal of Radiation Protection and Research
• /
• v.38 no.1
• /
• pp.37-43
• /
• 2013

In case of neutron dose measurement using TLDs (thermo-luminescence dosimeters), because the neutron energy dependence of the TLD is very high, the calibration of the energy response according to the characteristics of the neutron spectrum of workplace is required. In the present study, the ambient dose equivalent rates inside and around the Long-Counter (neutron detector) with narrow and complex inside in the neutron field of $^{252}Cf$ were evaluated. The calibration factors to account for the neutron energy dependence of TLDs were established for both the bare and $D_2O$ modulated $^{252}Cf$ neutron beams, respectively. The values of the TLD's measurement were compared with the computational results of the MCNPX (Monte Carlo N-Particles transport code). When using the two calibration factors of the TLD than a single calibration factor, the measured and the calculated values at the point of verification outside and inside the Long-Counter were in more good agreement. This results show that TLD should be calibrated in the reference neutron field similar to workplace situation.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3732-3753
• /
• 2023

The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, developed based on the Multi-stage Response-function Transport (MRT) methodology, enables real-time simulation of nuclear systems such as reactor cores, spent nuclear fuel pools and casks, and sub-critical facilities. This paper presents the application of a novel fission matrix-based burnup methodology to the well-characterized JSI TRIGA Mark II research reactor. This methodology allows for calculation of nuclear fuel depletion by combination and interpolation of RAPID's burnup dependent fission matrix (FM) coefficients to take into account core changes due to burnup. The methodology is compared to experimentally validated Serpent-2 Monte Carlo depletion calculations. The results show that the burnup methodology for RAPID (bRAPID) implemented into RAPID is capable of accurately calculating the k_eff burnup changes of the reactor core as the average discrepancies throughout the whole burnup interval are 37 pcm. Furthermore, capability of accurately describing 3D fission source distribution changes with burnup is demonstrated by having less than 1% relative discrepancies compared to Serpent-2. Good agreement is observed for axially and pin-wise dependent fuel burnup and nuclear fuel nuclide composition as a function of burnup. It is demonstrated that bRAPID accurately describes burnup in areas with high gradients of neutron flux (e.g. vicinity of control rods). Observed discrepancies for some isotopes are explained by analyzing the neutron spectrum. This paper presents a powerful depletion calculation tool that is capable of characterization of spent nuclear fuel on the fly while the reactor is in operation.

• Nuclear Engineering and Technology
• /
• v.52 no.10
• /
• pp.2387-2393
• /
• 2020

In this study, a gradient material for shielding neutrons and gamma rays was developed, which consists of epoxy resin, boron carbide (B₄C), lead (Pb) and a little graphene oxide. It aims light weight and compact, which will be applied on the transportable nuclear reactor. The material is made up of sixteen layers, and the thickness and components of each layer were designed by genetic algorithm (GA) combined with Monte Carlo N Particle Transport (MCNP). In the experiment, the viscosities of the epoxy at different temperatures were tested, and the settlement regularity of Pb particles and B₄C particles in the epoxy was simulated by matlab software. The material was manufactured at 25 ℃, the Pb C and O elements of which were also tested, and the result was compared with the outcome of the simulation. Finally, the material's shielding performance was simulated by MCNP and compared with the uniformity material's. The result shows that the shielding performance of gradient material is more effective than that of the uniformity material, and the difference is most noticeable when the materials are 30 cm thick.

• Nuclear Engineering and Technology
• /
• v.53 no.7
• /
• pp.2151-2161
• /
• 2021

Nuclear data (ND) sensitivity and uncertainty (S/U) quantification in shielding applications is performed using deterministic and probabilistic approaches. In this paper the validation of the newly developed deterministic program package ASUSD (ADVANTG + SUSD3D) is presented. ASUSD was developed with the aim of automating the process of ND S/U while retaining the computational efficiency of the deterministic approach to ND S/U analysis. The paper includes a detailed description of each of the programs contained within ASUSD, the computational workflow and validation results. ASUSD was validated on two shielding benchmark experiments from the Shielding Integral Benchmark Archive and Database (SINBAD) - the fission relevant ASPIS Iron 88 experiment and the fusion relevant Frascati Neutron Generator (FNG) Helium Cooled Pebble Bed (HCPB) Test Blanket Module (TBM) mock-up experiment. The validation process was performed in two stages. Firstly, the Denovo discrete ordinates transport solver was validated as a standalone solver. Secondly, the ASUSD program package as a whole was validated as a ND S/U analysis tool. Both stages of the validation process yielded excellent results, with a maximum difference of 17% in final uncertainties due to ND between ASUSD and the stochastic ND S/U approach. Based on these results, ASUSD has proven to be a user friendly and computationally efficient tool for deterministic ND S/U analysis of shielding geometries.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.732-740
• /
• 2021

Scattering resonance of medium mass nuclides leads complex spectrum in the fast reactor, which requires thousands of energy groups in the spectrum calculation. When the broad-group cross sections are collapsed, reaction rate cannot be completely conserved. To eliminate the error from energy collapsing, the Super-homogenization method in energy collapsing (ESPH) was employed in the fast reactor code SARAX. An ESPH factor was derived based on the ESPH-corrected SN transport equation. By applying the factor in problems with reflective boundary condition, both the effective multiplication factor and reaction rate were conserved. The fixed-source iteration was used to ensure the stability of ESPH iteration. However, in the energy collapsing process of SARAX, the vacuum boundary condition was adopted, which was necessary for fast reactors with strong heterogeneity. To further reduce the error caused by leakage, an additional conservation factor was proposed to correct the neutron current in energy collapsing. To evaluate the performance of ESPH with conservation factor, numerical benchmarks of fast reactors were calculated. The results of broad-group calculation agreed well with the direct full-core Monte-Carlo calculation, including the effective multiplication factor, radial power distribution, total control rod worth and sodium void worth.

• Nuclear Engineering and Technology
• /
• v.54 no.12
• /
• pp.4731-4742
• /
• 2022

RBMK fuel assemblies differ from other LWR FA due to a specific arrangement of the fuel rods, the low enrichment, and the used burnable absorber - erbium. Therefore, there is a challenge to adapt modeling tools, developed for other LWR types, to solve RBMK problems. A set of 10 different depletion simulation schemes were tested to estimate the impact on reactivity and spent fuel composition of possible SCALE code options for the neutron transport modelling and the use of different nuclear data libraries. The simulations were performed using cross-section libraries based on both, VII.0 and VII.1, versions of ENDF/B nuclear data, and assuming continuous energy and multigroup simulation modes, standard and user-defined Dancoff factor values, and employing deterministic and Monte Carlo methods. The criticality analysis with burn-up credit was performed for the SFP loaded with RBMK-1500 FA. Spent fuel compositions were taken from each of 10 performed depletion simulations. The criticality of SFP is found to be overestimated by up to 0.08% in simulation cases using user-defined Dancoff factors comparing the results obtained using the continuous energy library (VII.1 version of ENDF/B nuclear data). It was shown that such discrepancy is determined by the higher U-235 and Pu-239 isotopes concentrations calculated.

• Nuclear Engineering and Technology
• /
• v.54 no.2
• /
• pp.523-531
• /
• 2022

In order to study the sensitivity and the uncertainty of the Moroccan research reactor TRIGA Mark II, a model of this reactor has been developed in our ERSN laboratory for use with the N-Particle MCNP Monte Carlo transport codes (version 6). In this article, the sensitivities of the effective multiplication factor of this reactor are evaluated using the ENDF/B-VII.0, ENDF/B-VII.1 and JENDL-4.0 libraries and in 44 energy groups, for the cross sections of the fuel (U-235 and U-238) and the moderator (H-1 and O-16). However, the quantification of the uncertainty of the nuclear data is performed using the nuclear code NJOY99 for the generation and processing of covariance matrices. On the one hand, the highest uncertainty deviations, calculated using the ENDFB-VII.1 and JENDL4.0 evaluations, are 2275, 386 and 330 pcm respectively for the reactions U₂₃₅(n, f), $ U_{235}(n\bar{\nu})$ and H₁(n, γ). On the other hand, these differences are very small for the neutron reactions of O-16 and U-238. Regarding the neutron spectra, in CT-mid plane, they are very close for the three evaluations (ENDF/B-VII.0, ENDF/B-VII.1 and JENDL-4.0). These spectra present two peaks (thermal and fission) around the energies 0.05 eV and 1 MeV.

• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.469-476
• /
• 2020

Research reactors in-core experimental facilities are designed to provide the highest steady state flux for user's irradiation requirements. However, fuel conversion from highly enriched uranium (HEU) to low enriched uranium (LEU) driven by the ongoing effort to diminish proliferation risk, will impact reactor physics parameters. Preserving the reactor capability to produce the needed flux to perform its intended research functions, determines the conversion feasibility. This study investigates the neutron flux in the central experimental facility of two material test reactors (MTR), the IAEA generic10 MW benchmark reactor and the 22 MW s Egyptian Test and Research Reactor (ETRR-2). A 3D full core model with three uranium enrichment of 93%, 45%, and 20% was constructed utilizing the OpenMC particle transport Monte Carlo code. Neutronics calculations were performed for fresh fuel, the beginning of life cycle (BOL) and end of life cycle (EOL) for each of the three enrichments for both the IAEA 10 MW generic reactor and core 1/98 of the ETRR-2 reactor. Criticality calculations of the effective multiplication factor (K_eff) were executed for each of the twelve cases; results show a reasonable agreement with published benchmark values for both reactors. The thermal, epithermal and fast neutron fluxes were tallied across the core, utilizing the mesh tally capability of the code and are presented here. The axial flux in the central experimental facility was tallied at 1 cm intervals, for each of the cases; results for IAEA 10 MW show a maximum reduction of 14.32% in the thermal flux of LEU to that of the HEU, at EOL. The reduction of the thermal flux for fresh fuel was between 5.81% and 9.62%, with an average drop of 8.1%. At the BOL the thermal flux showed a larger reduction range of 6.92%-13.58% with an average drop of 10.73%. Furthermore, the fission reaction rate was calculated, results showed an increase in the peak fission rate of the LEU case compared to the HEU case. Results for the ETRR-2 reactor show an average increase of 62.31% in the thermal flux of LEU to that of the HEU due to the effect of spectrum hardening. The fission rate density increased with enrichment, resulting in 34% maximum increase in the HEU case compared to the LEU case at the assemblies surrounding the flux trap.