Corresponding author: Olga N. Andrianova ( o.n.andrianova@yandex.ru ) Academic editor: Yury Kazansky
© 2021 Olga N. Andrianova, Yury Ye. Golovko, Gleb B. Lomakov, Yevgeniya S. Teplukhina, Gennady M. Zherdev.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Andrianova ON, Golovko YuYe, Lomakov GB, Teplukhina YeS, Zherdev GM (2021) Calculation and experimental analysis of benchmark experiments with a fast neutron spectrum and models of sodium and lead cooled fast reactors using different evaluated nuclear data libraries. Nuclear Energy and Technology 7(2): 103-109. https://doi.org/10.3897/nucet.7.68951
|
The paper presents the results of a comparative analysis of criticality calculations using a Monte-Carlo code with the BNAB-93 and BNAB-RF neutron group constants, as well as with evaluated neutron data files from the Russian ROSFOND evaluated nuclear data library and other evaluated nuclear data libraries (ENDF, JEFF, JENDL) from different years. A set of integral experiments on BFS critical assemblies carried out in different years at the Institute of Physics and Power Engineering (60 different critical configurations) was analyzed. The considered integral experiments are included in the database of evaluated experimental neutronic data used to justify the neutronic performance of sodium and lead cooled fast reactors, to verify codes and nuclear data as well as to estimate uncertainties in neutronic parameters due to the nuclear data uncertainties. It has been shown that the ROSFOND evaluated nuclear data library is a library that minimizes the calculation and experimental discrepancies for the considered set of integral experiments. The paper also presents the results of criticality calculations for models of sodium and lead cooled fast reactors based on different evaluated neutron data libraries and provides estimates for the uncertainty in criticality associated with nuclear data.
Integral experiments, BFS facilities, ENDF, JEFF, JENDL, ROSFOND, BNAB-93, BNAB-RF, sodium cooled fast reactor, lead cooled fast reactor
Refining the neutron data for the isotopes playing a crucial role in achieving the target level of accuracy when predicting analytically the neutronic performance of Generation IV reactors is an issue receiving a great deal of attention. The most prominent efforts were undertaken as part of the international project CIELO (Collaborative International Evaluated Library Organisation) (
ROSFOND, the commonly available official version of the national library of evaluated neutron data, (
One of the factors that hampers the introduction of updated neutron data library versions is the fear of discrepancies to be detected in the reactor performance calculations based on an updated library of neutron data and the predecessor library. At the same time, there is no surprise that such discrepancies can be observed, this being proved by analyzing the world experience in introducing updated versions of evaluated neutron data in the process of which the causes for potential inconsistencies are demonstrated and discussed.
The predictive capability and the efficiency of updated neutron data library versions are compared by matching the results of the performance calculations for critical assemblies with experimental data. In the paper, respective performance indicators were calculated and compared to identify high-priority activities on further enhancing the ROSFOND/BNAB-RF nuclear data libraries. It was used a representative set of experimental data measured on the BFS-1 and BFS-2 critical facilities. Respective calculations were performed using different versions of Russian (ROSFOND, BNAB-RF and BNAB-93) and other evaluated nuclear data libraries (ENDF, JEFF and JENDL) from different years (ENDF/B-V.2 (USA, 1994), ENDF/B-VII.1 (USA, 2011) and ENDF/B-VIII (USA, 2018), JEFF-3.2 (Europe, 2014), JEFF-3.3 (Europe, 2017) and JENDL-4.0u2 (Japan, 2012)) (
Worldwide, the results of integral experiments are used for estimating and improving the accuracy of predictive calculations for reactor facilities and fuel cycle systems under design, planning new experiments, estimating the efficiency of experimental programs, etc. (
Table
Work has been initiated to form the databank for integral experiments based on BFS critical assemblies so that to systematize the accumulated information and enable its use for testing and adjusting evaluated nuclear data files, verifying codes, and estimating the accuracy of predicting the neutronic performance of fast reactors. This databank is designed to supplement the existing information and software tools to support the BFS experimental programs (
The criticality of the BFS critical assemblies shown in Table
Year | Assembly | Critical facility | Materials | Model |
---|---|---|---|---|
Sodium cooled fast reactor | ||||
1990 | 58-1 (1)* | BFS-2 | Pu/UO2/Na/ inert diluter (UO2) | Sodium cooled fast reactors |
1993 | 66-1 (1) | Pu/UO2/Na (UO2/Na 2) | ||
66-B (9) | Pu/U(d)/C/Na/ (UO2/Fe/Cr/Ni) | |||
1996 | 72 (2) | BFS-1 | U(d)/UO2/Pu/Al2O3/Na (UO2) | |
U(90/UO2/ZrH/Na (UO2) | ||||
1997 | 73 (1) | U/Na (UO2/Na) | ||
1998 | 78 (1) | Pu/UO2/Na (UO2/Na 2) | ||
Heavy liquid metal cooled fast reactor | ||||
1991 | 61-1(3) | BFS-1 | Pu/U(d)/C/Pb/Al (Pb/UO2) | Lead cooled fast reactors |
1999 | 77 (2) | Pu/U(d)/UO2/C/Pb (UO2) | ||
2000 | 64-1 (3) | BFS-2 | Pu/U(d)/C/Pb (PbBi) | |
Hydrogen-containing materials (ICSBEP) | ||||
2004–2005 | 97 (4) | BFS-1 | Pu/UO2 (UO2) | Fabrication of MOX fuel |
Pu/UO2/CH2 (UO2) | ||||
2008–2009 | 99 (3) | |||
101 (4) | ||||
Nuclear waste disposal systems (ICSBEP) | ||||
1999 | 79 (5) | BFS-1 | U/Si/CH2 | Waste disposal |
1999 | 81 (6) | Pu/Si/CH2 | ||
Different fuel compositions (ICSBEP) | ||||
1974 | 31(2) | BFS-2 | Pu/UO2 | CNFC |
1976 | 33 (3) | UO2/U(90)O2 | ||
1976 | 35 (3) | U(d)/U(36)/U(90) | ||
1977 | 38 (2) | Pu/U(d) | ||
1980 | 42 (1) | Pu/UO2/CH2 | ||
1985 | 49 (2) | Pu/UO2/CH2 | ||
Thermal reactors (IRPhEP) | ||||
1989 | 57 (1) | BFS-1 | U(36)O2/UO2/CH2/Al (UO2) | Thermal reactor models |
1990 | 59 (1) | Pu/UO2/CH2/Al (UO2) |
Since the configurations of the critical assembly cores have a heterogeneous structure and represent alternating layers of different materials with a thickness of 0.3 to 100 mm, calculations in a group approximation need to take into account the heterogeneous resonant self-shielding of neutron cross-sections. The calculations based on BNAB-93 and BNAB-RF were performed using a subgroup representation of neutron cross-sections in the resolved resonances region for the 238U, 239Pu and Fe isotopes and the neutron cross-section self-shielding factors for the rest of the isotopes and for 238U, 239Pu and Fe in the unresolved resonances region. The corrections for the self-shielding factors were calculated based on the principle of the equivalence of homogeneous and heterogeneous media. The procedures to calculate BFS critical assemblies using Monte Carlo codes in a group approximation are described in (
A statistical analysis was carried out using the calculation and experimental discrepancies for different versions of Russian and other evaluated nuclear data libraries, which makes it possible to conclude on the accuracy of calculations and the efficiency of describing a set of N integral experiments by the given system of constants L. Table
Neutron spectrum | Fast | Intermediate | All | ||||||
---|---|---|---|---|---|---|---|---|---|
Indicator | ? | d | µ | ? | d | µ | ? | d | µ |
ROSFOND | –0.21 | 0.31 | 1.12 | 0.13 | 0.83 | 2.09 | –0.10 | 0.54 | 1.53 |
ENDF/B-VII | –0.27 | 0.53 | 1.51 | 0.23 | 1.00 | 2.51 | –0.09 | 0.72 | 1.92 |
ENDF/B-VIII | –0.22 | 0.40 | 1.31 | 0.37 | 0.88 | 2.70 | –0.02 | 0.60 | 1.91 |
ENDF/B-V | –0.29 | 0.82 | 2.39 | 0.47 | 1.16 | 3.00 | –0.03 | 0.94 | 2.61 |
JEFF 3.2 | –0.18 | 0.41 | 1.47 | 0.30 | 0.85 | 2.46 | –0.01 | 0.59 | 1.87 |
JEFF 3.3 | –0.15 | 0.35 | 1.22 | 0.26 | 0.83 | 2.14 | –0.02 | 0.55 | 1.76 |
JENDL-4.0 | –0.02 | 0.47 | 1.11 | 0.46 | 0.87 | 2.34 | 0.15 | 0.63 | 1.64 |
BNAB-93 | –0.09 | 0.63 | 1.98 | 0.30 | 1.12 | 2.87 | 0.05 | 0.83 | 2.33 |
1. Average value of deviations of estimated multiplication factors from experimental values expressed as a percentage:
,
where ?nL is the calculated value obtained using library L for assembly n assigned to En (experimental value).
2. Root-mean-square deviation of calculation and experimental relations for multiplication factors expressed as a percentage:
,
where ?nL = (CnL/En – 1) is the calculation and experimental deviations in percentage terms obtained using library L for assembly n.
3. Number of standard deviations between the experiment and the criticality calculation, average for the set of N integral experiments, for the system of constants:
,
where de is the relative deviation of the experimental error, %; and dc is the relative value of the calculated (statistical) error, % (
All values of the performance indicators in Table
Average deviations of calculated multiplication factors in transition from ROSFOND to other libraries
Ci/CROSFOND –1, % | ||
---|---|---|
Average | Maximum | |
BNAB-RF | 0.11 | 0.34 |
ENDF/B-VII | 0.23 | 0.86 |
ENDF/B-VIII | 0.33 | 0.66 |
ENDF/B-V | 0.66 | 2.35 |
JEFF 3.2 | 0.26 | 0.70 |
JEFF 3.3 | 0.24 | 0.84 |
JENDL-4.0 | 0.39 | 0.90 |
BNAB-93 | 0.44 | 1.66 |
Table
.
Fig.
Deviations in calculated criticality values caused by the transition from the evaluated nuclear data library’s older version to the new one: 1) magnitude (? BNAB-93/? BNAB-RF –1) in percentage terms; 2) magnitude (? ENDF/B-V/? ENDF/B-VIII –1) in percentage terms. Regions shown by dashed lines: a) sodium cooled fast reactor; b) lead cooled fast reactor; c) fabrication of MOX fuel; d) waste disposal in sand; e) CNFC facilities.
Value of the deviation magnitude in calculations based on different libraries (?i/?j –1) in percentage terms: 1) minimum, 2) maximum. Regions shown by dashed lines: a) sodium cooled fast reactor; b) lead cooled fast reactor; c) fabrication of MOX fuel; d) waste disposal in sand; e) CNFC facilities.
Table
Deviations in keff calculations for fast reactor models involving the use of different evaluated neutron data libraries (C1/C2– 1, %)
Evaluated nuclear data library | Fast reactor 1 (Na+MOX) | Fast reactor 2 (Na+nitride fuel) | Fast reactor 3 (Pb+nitride fuel) | |
---|---|---|---|---|
C1 | C2 | C1/C2– 1, % | ||
BNAB-93 | BNAB-RF | 0.09 | 0.52 | 0.89 |
ENDF/B-V | ENDF/B-VIII | 0.04 | 0.47 | 0.82 |
ENDF/B-VII | 0.18 | 0.21 | -0.05 | |
ENDF/B-VIII | -0.07 | -0.03 | -0.05 | |
JEFF3.3 | 0.95 | 0.99 | 0.78 | |
JEFF3.2 | JEFF3.3 | 0.23 | 0.24 | 0.03 |
Based on covariance matrices and sensitivity coefficients of keff to the variation of the neutron cross-sections, the keff uncertainty were calculated for fast reactor models which have been found to be 1.9% for sodium cooled fast reactors and 2.0% for lead cooled fast reactors.
Based on analyzing the results of the keff calculations for the entire set of the BFS critical assemblies, the following conclusions can be made.
Present-day Russian evaluated nuclear data and constants make it possible to calculate strongly heterogeneous BFS critical assemblies using precise Monte Carlo codes with high accuracy. The transition from ROSFOND to its group version, BNAB-RF, for the considered BFS critical assemblies leads to a deviation of the keff calculations equal to ~ 0.1% (see Table
The transition from BNAB-93 to BNAB-RF leads to an increased accuracy of predicting keff for the BFS critical assemblies both with a fast neutron spectrum and with an intermediate neutron spectrum. On the average, the calculation and experimental discrepancies have decreased by a factor of one and a half to two; the criticality calculation values having changed by not more than 0.44% on the average. The maximum deviation of the keff calculations of up to 2% is observed for assemblies with an intermediate neutron spectrum and plutonium fuel. The maximum deviation value does not exceed the constant uncertainty for systems with an intermediate neutron spectrum, which varies in the limits of 2 to 3% depending on the assembly composition.
The transition from BNAB-93 (1993) to BNAB-RF (2010) leads to deviations of the criticality calculation results for fast reactors with nitride fuel in the limits of 0.5 to 0.9%. The spread in the calculated keff values obtained based on state-of-the-art U.S., European and Japanese evaluated nuclear data libraries from 2017–2018 is not less than ~ 1%.
The transition from JEFF3.3 and ENDF/B-VIII of the 2010 and 2011 versions to the 2017 and 2018 versions did not lead to a more accurate description of the BFS critical assemblies and major deviations in the keff calculations for fast reactor models. However, when comparing the deviations in the BFS critical assembly keff calculations (see Fig.
The observed spread in the calculated keff values for the BFS critical assemblies and fast reactor models as a result of using different evaluated neutron data libraries does not exceed the uncertainty of the calculations caused by the neutron cross-section uncertainties.
The results of the calculation and experimental analysis presented in this paper for a set of experiments on BFS critical assemblies make it possible to conclude the following. The results of the keff calculations based on the ROSFOND evaluated nuclear data library may differ from the results obtained using certain versions of U.S., European and Japanese evaluated nuclear data libraries by about 1%. The differences between the calculated keff values obtained using JEFF 3.3 and ENDF/B-VIII are also at a level of ~ 1%. When comparing different versions of the same library (e.g. ENDF/B-V.2 (1994) and ENDF/B-VIII (2018)) similar deviations are observed in the calculation-experiment results as in the event of the transition from BNAB-93 (1993) to BNAB-RF (2010). It has been shown for the entire set of the considered experiments that the calculation and experimental discrepancies are smaller when updated versions of the evaluated neutron data libraries are used, as compared with the results obtained based on earlier versions of the respective evaluated nuclear data libraries.