• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.3073-3084
• /
• 2022

The propagation of radiation source uncertainties in spent nuclear fuel (SNF) cask shielding calculations is presented in this paper. The uncertainty propagation employs the depletion and source term outputs of the deterministic code STREAM as input to the transport simulation of the Monte Carlo (MC) codes MCS and MCNP6. The uncertainties of dose rate coming from two sources: nuclear data and modeling parameters, are quantified. The nuclear data uncertainties are obtained from the stochastic sampling of the cross-section covariance and perturbed fission product yields. Uncertainties induced by perturbed modeling parameters consider the design parameters and operating conditions. Uncertainties coming from the two sources result in perturbed depleted nuclide inventories and radiation source terms which are then propagated to the dose rate on the cask surface. The uncertainty analysis results show that the neutron and secondary photon dose have uncertainties which are dominated by the cross section and modeling parameters, while the fission yields have relatively insignificant effect. Besides, the primary photon dose is mostly influenced by the fission yield and modeling parameters, while the cross-section data have a relatively negligible effect. Moreover, the neutron, secondary photon, and primary photon dose can have uncertainties up to about 13%, 14%, and 6%, respectively.

• Nuclear Engineering and Technology
• /
• v.54 no.9
• /
• pp.3260-3271
• /
• 2022

Experience shows that the solution to waste management in any national programme is lengthy and burdened with uncertainties. There are several uncertainties that contribute to the costs associated with spent fuel management. In this work, we have analysed the impact of the current nuclear data on the isotopic composition of the spent fuel and consequently their influence on the main spent fuel observables such as decay heat, activity, neutron multiplication factor, and neutron and photon source terms. Nuclear libraries based on the most general nuclear data ENDF/B-VII.0, ENDF/B-VII.1, ENDF/B-VIII.0 and JEFF-3.3 are considered. A typical NPP Krško fuel assembly is analysed using the Monte Carlo code Serpent 2. The analysis considers burnup of up to 60 GWd/tU and cooling times of up to 100 years. The comparison of results showed significant differences, which should be taken into account when selecting the library and evaluating the uncertainty in determining the characteristics of the spent fuel.

• Nuclear Engineering and Technology
• /
• v.46 no.3
• /
• pp.343-352
• /
• 2014

Uncertainty and sensitivity analyses with respect to nuclear data are performed with depletion calculations for BWR and PWR fuel assemblies specified in the framework of the UAM-LWR Benchmark Phase II. For this, the GRS sampling based tool XSUSA is employed together with the TRITON depletion sequences from the SCALE 6.1 code system. Uncertainties for multiplication factors and nuclide inventories are determined, as well as the main contributors to these result uncertainties by calculating importance indicators. The corresponding neutron transport calculations are performed with the deterministic discrete-ordinates code NEWT. In addition, the Monte Carlo code KENO in multi-group mode is used to demonstrate a method with which the number of neutron histories per calculation run can be substantially reduced as compared to that in a calculation for the nominal case without uncertainties, while uncertainties and sensitivities are obtained with almost the same accuracy.

• Nuclear Engineering and Technology
• /
• v.46 no.3
• /
• pp.281-290
• /
• 2014

In order to properly quantify fission reactor neutronics parameter uncertainties, we have to use covariance data and sensitivity profiles consistently. In the present paper, we establish two consistent methodologies for uncertainty quantification: a self-shielded cross section-based consistent methodology and an infinitely-diluted cross section-based consistent methodology. With these methodologies and the covariance data of uranium-238 nuclear data given in JENDL-3.3, we quantify uncertainties of infinite neutron multiplication factors of light water reactor and fast reactor fuel cells. While an inconsistent methodology gives results which depend on the energy group structure of neutron flux and neutron-nuclide reaction cross section representation, both the consistent methodologies give fair results with no such dependences.

• Nuclear Engineering and Technology
• /
• v.46 no.3
• /
• pp.313-342
• /
• 2014

A Nuclear Energy Agency (NEA), Organization for Economic Co-operation and Development (OECD) benchmark for Uncertainty Analysis in Modeling (UAM) is defined in order to facilitate the development and validation of available uncertainty analysis and sensitivity analysis methods for best-estimate Light water Reactor (LWR) design and safety calculations. The benchmark has been named the OECD/NEA UAM-LWR benchmark, and has been divided into three phases each of which focuses on a different portion of the uncertainty propagation in LWR multi-physics and multi-scale analysis. Several different reactor cases are modeled at various phases of a reactor calculation. This paper discusses Phase I, known as the "Neutronics Phase", which is devoted mostly to the propagation of nuclear data (cross-section) uncertainty throughout steady-state stand-alone neutronics core calculations. Three reactor systems (for which design, operation and measured data are available) are rigorously studied in this benchmark: Peach Bottom Unit 2 BWR, Three Mile Island Unit 1 PWR, and VVER-1000 Kozloduy-6/Kalinin-3. Additional measured data is analyzed such as the KRITZ LEU criticality experiments and the SNEAK-7A and 7B experiments of the Karlsruhe Fast Critical Facility. Analyzed results include the top five neutron-nuclide reactions, which contribute the most to the prediction uncertainty in keff, as well as the uncertainty in key parameters of neutronics analysis such as microscopic and macroscopic cross-sections, six-group decay constants, assembly discontinuity factors, and axial and radial core power distributions. Conclusions are drawn regarding where further studies should be done to reduce uncertainties in key nuclide reaction uncertainties (i.e.: $^{238}U$ radiative capture and inelastic scattering (n, n') as well as the average number of neutrons released per fission event of $^{239}Pu$).

• Nuclear Engineering and Technology
• /
• v.53 no.9
• /
• pp.2830-2838
• /
• 2021

This work deals with the assessment of the burnup capabilities of the Serpent Monte Carlo code to predict spent nuclear fuel (SNF) isotopic concentrations for low-enriched uranium (LEU) fuel at different burnup levels up to 47 MWd/kgU. The irradiation of six UO₂ experimental samples in three different VVER-1000 reactor units has been simulated and the predicted concentrations of actinides up to ²⁴⁴Cm have been compared with the corresponding measured values. The results show a global good agreement between calculated and experimental concentrations, in several cases within the margins of the nuclear data uncertainties and in a few cases even within the reported experimental uncertainties. The differences in the performances of the JEFF3.1.1, ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear data libraries (NDLs) have also been assessed and the use of the newly released ENDF/B-VIII.0 library has shown an increased accuracy in the prediction of the C/E's for some of the actinides considered, particularly for the plutonium isotopes. This work represents a step forward towards the validation of advanced simulation tools against post irradiation experimental data and the obtained results provide an evidence of the capabilities of the Serpent Monte-Carlo code with the associated modern NDLs to accurately compute SNF nuclide inventory concentrations for VVER-1000 type reactors.

• Nuclear Engineering and Technology
• /
• v.52 no.12
• /
• pp.2733-2742
• /
• 2020

This paper presents an assessment of applicability of the multigroup cross sections generated with Monte Carlo tools to the fast reactor analysis based on transport calculations. 33-group cross section sets were generated for simple one- (1-D) and two-dimensional (2-D) sodium-cooled fast reactor problems using the SERPENT code and applied to deterministic steady-state and depletion calculations. Relative to the reference continuous-energy SERPENT results, with the transport corrected P₀ scattering cross section, the k-eff value was overestimated by 506 and 588 pcm for 1-D and 2-D problems, respectively, since anisotropic scattering is important in fast reactors. When the scattering order was increased to P₅, the 1-D and 2-D problem errors were increased to 577 and 643 pcm, respectively. A sensitivity and uncertainty analysis with the PERSENT code indicated that these large k-eff errors cannot be attributed to the statistical uncertainties of cross sections and they are likely due to the approximate anisotropic scattering matrices determined by scalar flux weighting. The anisotropic scattering cross sections were alternatively generated using the MC²-3 code and merged with the SERPENT cross sections. The mixed cross section set consistently reduced the errors in k-eff, assembly powers, and nuclide densities. For example, in the 2-D calculation with P3 scattering order, the k-eff error was reduced from 634 pcm to -223 pcm. The maximum error in assembly power was reduced from 2.8% to 0.8% and the RMS error was reduced from 1.4% to 0.4%. The maximum error in the nuclide densities at the end of 12-month depletion that occurred in ²³⁷Np was reduced from 3.4% to 1.5%. The errors of the other nuclides are also reduced consistently, for example, from 1.1% to 0.1% for ²³⁵U, from 2.2% to 0.7% for ²³⁸Pu, and from 1.6% to 0.2% for ²⁴¹Pu. These results indicate that the scalar flux weighted anisotropic scattering cross sections of SERPENT may not be adequate for application to fast reactors where anisotropic scattering is important.

• Nuclear Engineering and Technology
• /
• v.52 no.9
• /
• pp.1907-1916
• /
• 2020

In this work, the depletion capability implemented in Monte Carlo code MCS is investigated to predict the isotopic compositions of spent nuclear fuel (SNF). By comparison of MCS calculation results to post irradiation examination (PIE) data obtained from one pressurized water reactor (PWR), the validation of this capability is conducted. The depletion analysis is performed with the ENDF/B-VII.1 library and a fuel assembly model. The transmutation equation is solved by the Chebyshev Rational Approximation Method (CRAM) with a depletion chain of 3820 isotopes. 18 actinides and 19 fission products are analyzed in 14 SNF samples. The effect of statistical uncertainties on the calculated number densities is discussed. On average, most of the actinides and fission products analyzed are predicted within ±6% of the experiment. MCS depletion results are also compared to other depletion codes based on publicly reported information in literature. The code-to-code analysis shows comparable accuracy. Overall, it is demonstrated that the depletion capability in MCS can be reliably applied in the prediction of SNF isotopic inventory.

• Journal of Radiation Industry
• /
• v.17 no.4
• /
• pp.471-479
• /
• 2023

In In-situ radioactivity measurement techniques, efficiency calibration models use predefined models to simulate a sample's geometry and radioactivity distribution. However, simplified efficiency calibration models lead to uncertainties in the efficiency curves, which in turn affect the radioactivity concentration results. This study aims to develop an efficiency calibration optimization methodology to improve the accuracy of in-situ gamma radiation measurements for byproducts from industrial facilities. To accomplish the objective, a drive mechanism for rotational measurement of an byproduct simulator and a sample was constructed. Using ISOCS, an efficiency calibration model of the designed object was generated. Then, the sensitivity analysis of the efficiency calibration model was performed, and the efficiency curve of the efficiency calibration model was optimized using the sensitivity analysis results. Finally, the radiation concentration of the simulated subject was estimated, compared, and evaluated with the designed certification value. For the sensitivity assessment of the influencing factors of the efficiency calibration model, the ISOCS Uncertainty Estimator was used for the horizontal and vertical size and density of the measured object. The standard deviation of the measurement efficiency as a function of the longitudinal size and density of the efficiency calibration model decreased with increasing energy region. When using the optimized efficiency calibration model, the measurement efficiency using IUE was improved compared to the measurement efficiency using ISOCS at the energy of ²²⁸Ac (911 keV) for the nuclide under analysis. Using the ISOCS efficiency calibration method, the difference between the measured radiation concentration and the design value for each simulated subject measurement direction was 4.1% (1% to 10%) on average. The difference between the estimated radioactivity concentration and the design value was 3.6% (1~8%) on average when using the ISOCS IUE efficiency calibration method, which was closer to the design value than the efficiency calibration method using ISOCS. In other words, the estimated radioactivity concentration using the optimized efficiency curve was similar to the designed radioactivity concentration. The results of this study can be utilized as the main basis for the development of regulatory technologies for the treatment and disposal of waste generated during the operation, maintenance, and facility replacement of domestic byproduct generation facilities.