Analysis of OECD/NEA medium 1000 MWth sodium-cooled fast reactor using the Monte Carlo Serpent Code and ENDF/B-VIII.0 nuclear data library

Analysis of OECD/NEA medium 1000 MWth sodium-cooled fast reactor using the Monte Carlo Serpent Code and ENDF/B-VIII.0 nuclear data library 1 first-author 1 https://orcid.org/0000-0002-4864-1400 2 3 https://orcid.org/0000-0002-0285-4507 1 corresp Corresponding author, dhartanto@sharjah.ac.ae Corresponding author, dhartanto@sharjah.ac.ae dhartanto@sharjah.ac.ae https://orcid.org/0000-0002-9554-1111 Department of Mechanical and Nuclear Engineering, University of Sharjah, P.O. BOX 27272, Sharjah, United Arab Emirates Department of Mechanical and Nuclear Engineering, University of Sharjah, P.O. BOX 27272, Sharjah, United Arab Emirates Cooperative Major in Nuclear Energy, Graduate School of Engineering, Tokyo City University (TCU), 1-28-1 Tamazutsumi, Setagaya, Tokyo, Japan Cooperative Major in Nuclear Energy, Graduate School of Engineering, Tokyo City University (TCU), 1-28-1 Tamazutsumi, Setagaya, Tokyo, Japan Scientific Computational Division, Nippon Advanced Information Service (NAIS Co., Inc.), 416 Muramatsu, Tokaimura, Ibaraki, Japan Scientific Computational Division, Nippon Advanced Information Service (NAIS Co., Inc.), 416 Muramatsu, Tokaimura, Ibaraki, Japan This study presents a benchmark evaluation of the new ENDF/B-VIII.0 nuclear data library for the Organization for Economic Co-operation and Development/Nuclear Energy Agency Medium 1000 MWth sodium-cooled fast reactor (SFR). The study presented herein covers both SFR core types i.e., metallic fueled (MET-1000) and oxide-fueled (MOX-1000), simulated using the continuous-energy Monte Carlo Serpent2 code. The neutronics performances of the ENDF/B-VIII.0-based simulations were compared mainly to two libraries: ENDF/B-VII.1 and JENDL-4.0. The comparison includes several neutronics parameters evaluated for the beginning and end of the cycle conditions. These parameters include the effective multiplication factor k eff , total effective delayed neutron fraction β eff , sodium void reactivity (∆ρ Na ), Doppler constant (∆ρ Doppler ), and control rod worth (∆ρ CR ). In addition, a sensitivity study was used to reveal the major isotope/reaction pairs contributing to the discrepancy observed in the performance of the three libraries using 33 and 44 energy-group structures. Serpent ENDF/B-VIII.0 Sodium-cooled fast reactor Sensitivity analysis 1 Introduction 1 Under the working party of the Organization for Economic Co-operation and Development/Nuclear Energy Agency (OECD/NEA), a neutronics benchmark study of the Generation-IV sodium-cooled fast reactor (SFR) concepts has been conducted and published using different methods and codes [ 1 1 – 3 3 1–3 ] in conjunction with various nuclear data libraries such as ENDF/B-VII.0 [ 4 4 ], ENDF/B-VII.1 [ 5 5 ], JEFF-3.1 [ 6 6 ], and JENDL-4.0 [ 7 7 ]. The benchmark consists of large 3600 MWth carbide and oxide cores as well as medium 1000 MWth metallic and oxide cores. In this study, the medium SFR cores were calculated using the new version of the nuclear data library, ENDF/B-VIII.0 [ 8 8 ]. The newer library has achieved improvements in many neutron cross-section libraries, including the cross-sections of materials that are commonly used in the SFR, such as U-235 and Pu-239 in the fuel, Fe-56 in the clad, and O-16 in mixed oxide fuel. Because the accuracy of the nuclear data is extreme ly important for the design and safety of a nuclear reactor, the impact of the new ENDF/B-VIII.0 on the neutronics and kinetics parameters of an OECD/NEA Medium 1000 MWth SFR was investigated. Both metallic and oxide cores are considered to distinguish the impact of the neutron spectrum. Several important neutronics parameters are evaluated in the present study, including the effective neutron multiplication factor k , total effective delayed neutron fraction β eff , prompt neutron generation time Λ, Doppler constant K D , coolant void reactivity (CVR), and control rod (CR). Each parameter is calculated at the beginning of cycle (BOC) and end of cycle (EOC) using the continuous-energy Monte Carlo Serpent2 code [ 9 9 ]. Three different modern nuclear data libraries are considered: ENDF/B-VIII.0, ENDF/B-VII.1, and JENDL-4.0. In addition to an inter-library comparison, the calculated results are also compared to the results reported by the working group [ 3 3 ]. It should be noted that the reported result is an average value computed by different calculation methods, including deterministic and stochastic methods, and various nuclear data libraries. 2 SFR: Core Description 1 The layout of the metallic core is illustrated in Fig. 1 Fig. 1 . The core consists of inner and outer fuel regions surrounded by reflectors and shielding. The inner core region has 78 fuel subassemblies, and the outer core region has 102 fuel subassemblies. Meanwhile, there were 114 reflector subassemblies, 66 shielding subassemblies, 15 primary control subassemblies, and 4 secondary control subassemblies. The fuel consists of irradiated U-10Zr metallic fuel, and is contained in HT-9 cladding. The total core height is 480.20 cm, and the active core height is approximately 85.82 cm. In the calculation, the temperature of the fuel was 534 °C, and the temperature of the coolant and structural material was 432.5 °C. In Figure 2 Figure 2 , the radial layout of the oxide core is shown. The core has three fuel regions, which are the inner, middle, and outer regions. The fuel is irradiated UO 2 -(TRU)-O 2 contained in the HT-9 clad. There are 30 inner fuel subassemblies, 90 mid-fuel subassemblies, and 60 outer fuel subassemblies. Meanwhile, the number of non-fuel subassemblies is the same as that in the metallic core. The total core height is also the same as in the metallic core; however, the active core is slightly taller, at approximately 114.94 cm. The fuel temperature is also higher, at approximately 1027 °C, and the other components have the same temperatures as in the metallic core. For detailed information on the benchmark, including the material compositions and assembly dimensions, readers can refer to [ 1 1 ]. 3 ENDF/B-VIII.0 Nuclear Data Library 1 Under the coordination of the Cross-Section Evaluation Working Group (CSEWG), the new ENDF/B-VIII.0 library was released in 2018 [ 8 8 ] and issued in both formats: the traditional ENDF-6 format and the new Generalized Nuclear Database Structure (GNDS) format. This version has represented the most noticeable changes to the ENDF library among all of the previously released versions, which improved the nuclear reaction data library by using the CIELO-project cross-sections, new standards, and thermal scattering data, as a result of full co-operation with the new IAEA standards. The updates in the library involve an updated neutron reaction on the structural materials, minor actinides, dosimetry cross-sections, decay data, fission energy release, light nuclei, charged-particle reactions, and thermal neutron scattering data. The major highlighted and included changes in the ENDF/B-VIII.0 are summarized in Table 1 Table 1 . 1 1 Category 1 1 Changed isotopes 1 1 CIELO evaluation 1 1 1 H, 16 O, 56 Fe, 235 U, 238 U, 239 Pu, including prompt fission neutron spectra (PFNS) and prompt fission gamma spectra (PEGS) 1 1 Light elements 1 1 2 H, 3 He, 6 Li, 9 Be, 10 B, 12,13 C, 35,37 Cl, and 18 O 1 1 Structural materials 1 1 40 Ca, 54,56,57,58 Fe, 58-62,64 Ni, 59 Co, 63,65 Cu, 174-182 Hf, 182-186 W, 105 Rh, and 132 Te 1 1 Nobel gases 1 1 40 Ar, 78 Kr, 124 Xe, and 20−22 Ne 1 1 Minor actinides 1 1 236m Np, 240 Pu, and 241,243 Am 1 1 Misc. materials 1 1 73−75 As, 197 Au, and 190−198 Pt 1 1 Unstable isotopes 1 1 T 1/2 ≥ 1 isotopes 1 1 Primary gammas 1 1 6,7 Li, 11 B, 19 F, 23 Na, 27 Al, 28 Si and 35,37 Cl Moreover, a wider energy range was considered as a standard by adding the integral cross-section of 235 U( n ,f) from 7.8-11 eV, the 30 keV Maxwellian-averaged cross-section of Au( n , γ ), and the higher energy fission reference cross-sections from 200 MeV to 1 GeV for 235 U( n ,f) and 238 U( n ,f), and from 20 MeV to 1 GeV for 209 Bi( n ,f) and nat Pb( n ,f) [ 8 8 ]. The neutron sub-library was expanded to include 557 evaluations, which is 32% of the increment [ 8 8 ]. In return, these changes have a direct impact on the simulations of nuclear criticality. In this study, a comparison of using ENDF/B-VIII.0, ENDF/B-VII.1, and JENDL-4.0 in simulating the nuclear performance of the SFR to see the impact of implementing the new library in contrast with the previous libraries. 4 Methodology 1 4.1 Monte Carlo SERPENT2 Code 2 In this study, a multi-purpose three-dimensional continuous-energy Monte Carlo particle transport code, Serpent2, was used to simulate the neutronics behavior of SFR metallic and oxide cores using the new ENDF/B-VIII.0 library, ENDF/B-VII.1, and JENDL-4.0. The development of Serpent2 was started in 2010 by the VTT Technical Research Centre of Finland [ 9 9 ]. The transport simulation by Serpent2 was conducted based on the k-eigenvalue criticality source method, which limits its applications to a self-sustaining system. The code evaluates the kinetic and delayed neutron parameters automatically, in addition to the entire few group constants required for a coupled nodal diffusion calculation. Based on the user interest, the reaction rate and integral flux detectors for specific regions can be defined and inserted into the code using tallies. Serpent2 is also capable of calculating the sensitivities of various responses to various perturbations using the collision history method and its corresponding uncertainty using the user-provided covariance data. Moreover, the user can construct a 2D or 3D fuel or reactor geometry configuration, where Serpent2 applies the universe based constructive solid geometry model (CSG) [ 10 10 ]. For particle transportation, a collection of the classical surface tracking is used, in addition to the Woodcock delta-tracking method, which is implemented when the given dimensions are smaller than the mean free path of the particles [ 11 11 ]. The interaction physics in Serpent2 depends on the ENDF reaction laws, the conventional collision kinematics, the Doppler-broadening rejection correction (DBRC) method, which treats the free-gas scattering kernel close to resonances [ 12 12 ], and probability table sampling for the unsolved resonance region. In addition, the continuous-energy cross-sections are provided using ACE format libraries. 4.2 Focus of the Current Study 2 In this study, the impact of the changes in the new ENDF/B-VIII.0 on the neutronics and kinetics parameters of an OECD/NEA Medium 1000 MWth SFR as compared to the implementation of older libraries, ENDF/B-VII.1, and JENDL-4.0, was investigated. Several simulations were conducted using the Serpent2 code to evaluate the important neutronics parameters, such as the effective neutron multiplication factor k eff , total effective delayed neutron fraction β eff , prompt neutron generation time Λ, Doppler constant K D , CVR, and CR. A depletion study was conducted, and the BOC and EOC cores were studied. In addition, a sensitivity study was simulated to see the differences between using a 33-energy group and a 44-energy group for each library considered. The 33-energy group is the typical energy group used in the fast reactor analysis [ 13 13 ], and the 44-energy group is the energy group used in the covariance data. 5 Results and Discussion 1 The Monte Carlo Serpent2 code was used to produce the subsequent results. In the calculation, the sub-assembly in the core was explicitly modeled. Each calculation used 100,000 neutron histories and 1,050 total neutron cycles with 50 inactive cycles. Using this calculation condition, the standard deviation of k eff is approximately 7 pcm. 5.1 Multiplication Factor at BOC and EOC 2 The effective neutron multiplication factor k eff values for the metallic and oxide cores at the BOC and EOC are summarized in Tables 2 Tables 2 and 3 3 , respectively. Comparing among the different libraries, JENDL-4.0 provides the highest k eff value of the metallic core at the BOC and EOC than those by ENDF/B-VII.0 and ENDF/B-VII.1. Meanwhile, both JENDL-4.0 and ENDF/B-VII.1 produce a similar k eff , and are higher than those produced by ENDF/VIII.0 for the oxide core. It should be noted that the calculated k eff values are lower than those reported [ 3 3 ]. 1 1 Library 1 1 k eff at BOC 1 1 k eff at EOC 1 1 ENDF/B-VIII.0 1 1 1.02993 ± 0.00007 1 1 1.00580 ± 0.00007 1 1 ENDF/B-VII.1 1 1 1.02808 ± 0.00007 1 1 1.00522 ± 0.00007 1 1 JENDL-4.0 1 1 1.03240 ± 0.00007 1 1 1.00893 ± 0.00007 1 1 Reported result (average) [ 3 3 ] 1 1 1.03578 ± 0.0078 1 1 1.01230 ± 0.0071 1 1 Library 1 1 k eff at BOC 1 1 k eff at EOC 1 1 ENDF/B-VIII.0 1 1 1.02362 ± 0.00006 1 1 1.00567 ± 0.00007 1 1 ENDF/B-VII.1 1 1 1.02610 ± 0.00006 1 1 1.00857 ± 0.00006 1 1 JENDL-4.0 1 1 1.02539 ± 0.00006 1 1 1.00876 ± 0.00007 1 1 Reported result (average) [ 3 3 ] 1 1 1.02860 ± 0.0062 1 1 1.01360 ± 0.0082 The value of k eff at the EOC was taken from the depletion results at 328.5 days. In the burnup calculation, the corresponding decay and neutron fission yield libraries of each nuclear data were used. Figures 3 Figures 3 and 4 4 show k eff as a function of burnup. The change of k eff during burnup is linear for each nuclear data because the cross-sections of the fission products in the fast spectrum are not as large as in the thermal spectrum. The neutron spectra for each library for the metallic and oxide cores are also illustrated in Figs. 5 Figs. 5 and 6 6 . It should be noted that the difference in the neutron spectra can reach approximately 20% within the resonance energy region. 5.2 Sensitivity Coefficient of Multiplication Factor 2 The sensitivity coefficients of k eff for the metallic and oxide cores are shown in Figures 7 Figures 7 and 8 8 , respectively. Serpent2 calculated the sensitivity coefficients based on the collision history-based approach [ 14 14 ]. The sensitivity coefficients are consistent among different nuclear data libraries. It should be noted that Pu-239 ( n,f ), U-238 ( n,f ), Pu-240 ( n,f ), and Pu-241 ( n,f ) give large positive sensitivities for both cores. Meanwhile, U-238 ( n,g ) and ( n,inl ) dominate the negative sensitivities for the metallic core. By contrast, the large negative sensitivities for the oxide core include U-238 (n,g ) and ( n,inl ), Pu-239 ( n,g ), O-16 ( n,el ), Pu-240 ( n,g ), and Fe-56 ( n,inl ). 5.3 k eff Sensitivity 2 An analysis for evaluating the impact of the ENDF/B-VIII.0 library was also conducted by replacing a single isotope from the old library with a new one from the ENDF/B-VIII.0 library at BOC, which is summarized in Tables 4 Tables 4 and 5 5 . It should be noted that there are significant changes in k eff from Fe-56 and U-238 of ENDF/B-VII.1 and Fe-56 and Pu-239 of JENDL-4.0 for both cores. A detailed study was then conducted to evaluate the contribution of each isotope/reaction pair to the change in k eff by comparing the results of the new library ENDF/B-VIII.0 with JENDL-4.0 and ENDF/B-VII.1. If the k eff sensitivity coefficient S σ is defined as follows: 1 1 Nuclear Data 1 1 Isotope Replaced by ENDF-B/VIII.0 1 1 Change (pcm) 1 1 ENDF/B-VII.1 1 1 Na-23 1 1 14 ± 7 1 1 1 1 Fe-56 1 1 -466 ± 7 1 1 1 1 U-235 1 1 6 ± 7 1 1 1 1 U-238 1 1 266 ± 7 1 1 1 1 Pu-239 1 1 86 ± 7 1 1 1 1 Am-241 1 1 8 ± 7 1 1 JENDL-4.0 1 1 Na-23 1 1 68 ± 7 1 1 1 1 Fe-56 1 1 -281 ± 7 1 1 1 1 U-235 1 1 15 ± 7 1 1 1 1 U-238 1 1 32 ± 7 1 1 1 1 Pu-239 1 1 -145 ± 7 1 1 1 1 Am-241 1 1 -94 ± 7 1 1 Nuclear Data 1 1 Isotope Replaced by ENDF-B/VIII.0 1 1 Change (pcm) 1 1 ENDF/B-VII.1 1 1 O-16 1 1 -33 ± 9 1 1 1 1 Na-23 1 1 -8 ± 9 1 1 1 1 Fe-56 1 1 -569 ± 9 1 1 1 1 U-235 1 1 -14 ± 9 1 1 1 1 U-238 1 1 205 ± 9 1 1 1 1 Pu-239 1 1 -42 ± 9 1 1 JENDL-4.0 1 1 Am-241 1 1 -10 ± 9 1 1 1 1 O-16 1 1 16 ± 9 1 1 1 1 Na-23 1 1 96 ± 9 1 1 1 1 Fe-56 1 1 -348 ± 9 1 1 1 1 U-235 1 1 13 ± 9 1 1 1 1 U-238 1 1 78 ± 9 <math display="block" id="M1"><mrow><msub><mi>S</mi><mi>σ</mi></msub><mo>=</mo><mfrac><mi>σ</mi><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></mfrac><mfrac><mrow><mtext>d</mtext><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow><mrow><mtext>d</mtext><mi>σ</mi></mrow></mfrac><mo>,</mo></mrow></math> where <math id="M2"><mi>σ</mi></math> is the multigroup microscopic cross-section obtained from Serpent2. Then, the first-order estimation of the k eff change can be calculated using Eq. (2), where Δ σ is the difference in the cross-section (between different libraries), and <math id="M3"><mrow><msub><mi>S</mi><mi>σ</mi></msub></mrow></math> is the energy-dependent sensitivity coefficient obtained by Serpent2. <math display="block" id="M4"><mrow><mo> </mo><mi>Δ</mi><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow><mi>σ</mi></mfrac><mo>×</mo><msub><mi>S</mi><mi>σ</mi></msub><mo>×</mo><mi>Δ</mi><mi>σ</mi><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo></mrow></math> The results of the major contributors to k eff of Fe-56, U-238, and Pu-239 were obtained for the metallic and oxide SFR core. Moreover, the difference between utilizing a 44-energy group or a 33-energy group of cross-sections and S σ on the k eff was also investigated. Tables 6 Tables 6 through 9 9 summarize the results of each isotope contribution to k eff in total (in the resultant change column). The individual reaction contributions of each isotope are also shown. The resultant change is defined as the summation of the reaction contribution for each isotope to k eff . 1 1 Nuclear Data 1 1 Isotope Replaced by ENDF/B-VIII.0 4 1 Reaction 1 1 Resultant change 1 1 1 1 1 1 ( n,el ) 1 1 ( n,inl ) 1 1 ( n,f ) 1 1 ( n,g ) 1 1 1 1 ENDF/B-VII.1 1 1 Fe-56 1 1 -7.30 1 1 -151.06 1 1 - 1 1 -118.29 1 1 -276.65 1 1 1 1 U-238 1 1 -8.75 1 1 -77.95 1 1 99.06 1 1 123.04 1 1 135.40 1 1 1 1 Pu-239 1 1 0.14 1 1 0.03 1 1 154.15 1 1 -36.08 1 1 118.25 1 1 JENDL-4.0 1 1 Fe-56 1 1 -19.46 1 1 -95.86 1 1 - 1 1 -65.51 1 1 -180.84 1 1 1 1 U-238 1 1 -11.38 1 1 -284.92 1 1 118.37 1 1 185.15 1 1 7.22 1 1 1 1 Pu-239 1 1 -3.01 1 1 -11.92 1 1 26.24 1 1 7.49 1 1 18.80 1 1 Nuclear data 1 1 Isotope by ENDF-B/VIII.0 4 1 Reaction 1 1 Resultant change 1 1 1 1 1 1 ( n , el ) 1 1 ( n , inl ) 1 1 ( n , f ) 1 1 ( n , g ) 1 1 1 1 ENDF/B-VII.1 1 1 Fe-56 1 1 -7.11 1 1 -164.99 1 1 - 1 1 -178.20 1 1 -350.31 1 1 1 1 U-238 1 1 6.64 1 1 -41.05 1 1 82.19 1 1 123.02 1 1 170.81 1 1 1 1 Pu-239 1 1 0.46 1 1 0.04 1 1 53.31 1 1 -44.44 1 1 9.37 1 1 JENDL-4.0 1 1 Fe-56 1 1 -9.95 1 1 -106.99 1 1 - 1 1 -114.10 1 1 -231.03 1 1 1 1 U-238 1 1 -2.41 1 1 -208.27 1 1 92.45 1 1 159.26 1 1 41.02 1 1 1 1 Pu-239 1 1 2.35 1 1 -7.69 1 1 -107.69 1 1 -44.43 1 1 -157.46 The results of the metallic core summarized in Tables 6 Tables 6 and 7 7 illustrate an increase in k eff when using ENDF/B-VII.1 compared to the new library ENDF/B-VIII.0 (for both group structures). The cross-section of U-238 from ENDF/B-VIII.0 increases the k eff by approximately 175 pcm and 135 pcm when analyzed using the 44-energy group and 33-energy group, respectively, and the major reactions are of U-238 ( n,f ) and ( n,g ). Similarly, the cross-section of Pu-239 from ENDF/B-VIII.0 increases the k eff by approximately 121 pcm (44-energy group) and 118 pcm (33-energy-group) and is mostly contributed to by the fission reaction. However, the cross-section of Fe-56 from ENDF/B-VIII.0 decreases the k eff by approximately -278 pcm (44-energy-group) and -276 pcm (33-energy-group). It is also worth mentioning that both ENDF/B-VIII.0 and ENDF/B-VII.1 nuclear data have similar inelastic scattering cross-sections of Pu-239 in this core, which results in a zero discrepancy for the ( n,inl ) reaction. 1 1 Nuclear data 1 1 Isotope replaced by ENDF/B-VIII.0 4 1 Reaction 1 1 Resultant change 1 1 1 1 1 1 ( n,el ) 1 1 ( n,inl ) 1 1 ( n,f ) 1 1 ( n,g ) 1 1 1 1 ENDF/B-VII.1 1 1 Fe-56 1 1 -7.32 1 1 -153.89 1 1 - 1 1 -117.25 1 1 -278.47 1 1 1 1 U-238 1 1 -0.94 1 1 -47.30 1 1 102.18 1 1 121.12 1 1 175.07 1 1 1 1 Pu-239 1 1 -0.23 1 1 0.00 1 1 163.34 1 1 -42.05 1 1 121.05 1 1 JENDL-4.0 1 1 Fe-56 1 1 -19.09 1 1 -108.86 1 1 - 1 1 -59.37 1 1 -187.32 1 1 1 1 U-238 1 1 -7.17 1 1 -277.79 1 1 111.89 1 1 171.87 1 1 -1.21 1 1 1 1 Pu-239 1 1 -5.92 1 1 -14.67 1 1 36.74 1 1 6.39 1 1 22.53 Meanwhile, in the case of the JENDL-4.0 nuclear data library, it is shown that k eff changes by approximately -187 pcm (44-energy group) and -180 pcm (33-energy-group) when Fe-56 from ENDF/B-VIII.0 was used. A minor change of approximately -1.2 pcm (44-energy-group) and 7 pcm (33-energy-group) was noted for U-238 from ENDF/B-VIII.0, which resulted from the balance between the positive contribution of ( n,f ) and ( n,g ) and the negative contribution of elastic and inelastic scattering reactions. Moreover, Pu-239 from ENDF/B-VIII.0 increased k eff by approximately 22 pcm (44-energy group) and 18 pcm (33-energy-group), which was mainly due to the enhancement in the fission reaction. These results for both ENDF/B-VII.1 and JENDL-4.0 are consistent with the trends in Tables 4 Tables 4 and 5 5 . A smaller discrepancy is noted owing to the use of a multi-energy group determine the contribution of each reaction to the k eff . Based on the results for the oxide core in Tables 8 Tables 8 and 9 9 , replacing the Fe-56 cross-section from ENDF/B-VII.1 to the ENDF/B-VIII.0 library has a negative contribution to k eff (-350 pcm for the 44 energy group and -353 pcm for the 33 energy group), whereas the contribution of U-238 was increased (170 pcm for the 44-energy group and 123 pcm for the 33 energy-group). Moreover, the change in keff by Pu-239 was minor (9 pcm for the 44-energy group and -9 pcm for the 33-energy group) owing to the balance between the increment in the fission reaction and the decrement in the ( n , g ) reaction. By contrast, replacing JENDL-4.0 with the new library (ENDF/B-VIII.0) revealed a smaller contribution for U-238 (41 pcm for the 44-energy group and 34 pcm for the 33-energy group). The contribution of Fe-56 from the new nuclear data decreases k eff by approximately -231.03 and -238.61 pcm for the 44-energy-group and 33-energy-group, respectively. Meanwhile, Pu-239 from ENDF/B-VIII.0 showed a decrease in k eff of approximately -157 pcm for the 44-energy group and -169 pcm for the 33-energy group. 1 1 Nuclear data 1 1 Isotope by ENDF-B/VIII.0 4 1 Reaction 1 1 Resultant change 1 1 1 1 1 1 ( n,el ) 1 1 ( n,inl ) 1 1 ( n,f ) 1 1 ( n,g ) 1 1 1 1 ENDF/B-VII.1 1 1 Fe-56 1 1 -8.35 1 1 -160.86 1 1 - 1 1 -183.80 1 1 -353.01 1 1 1 1 U-238 1 1 -6.53 1 1 -64.26 1 1 78.02 1 1 115.99 1 1 123.22 1 1 1 1 Pu-239 1 1 -0.29 1 1 0.06 1 1 35.35 1 1 -44.39 1 1 -9.27 1 1 JENDL-4.0 1 1 Fe-56 1 1 -12.19 1 1 -101.78 1 1 - 1 1 -124.64 1 1 -238.61 1 1 1 1 U-238 1 1 -4.37 1 1 -215.75 1 1 100.91 1 1 154.00 1 1 34.80 1 1 1 1 Pu-239 1 1 0.62 1 1 -7.99 1 1 -118.77 1 1 -43.11 1 1 -169.25 1 1 Tables 10 Tables 10 and 11 11 illustrate the differences between the 44-energy group and 33-energy-group structures. The highest contribution was recorded for U-238 using the ENDF/B-VII.1 library for both the metallic and the oxide cores, which was mainly caused by the inelastic neutron scattering reaction, which is dominant in the fast spectrum region, resulting in more energy groups (44 instead of 33), resulting in a significant difference in this case. However, the difference between using the 44-energy group and the 33-energy group is relatively small and slightly higher for the oxide core than that of the metallic core. 1 1 Nuclear data 1 1 Isotope by ENDF-B/VIII.0 4 1 Reaction 1 1 Resultant change 1 1 1 1 1 1 ( n,el ) 1 1 ( n,inl ) 1 1 ( n,f ) 1 1 ( n , g ) 1 1 1 1 ENDF/B-VII.1 1 1 Fe-56 1 1 -0.02 1 1 -2.83 1 1 - 1 1 1.04 1 1 -1.82 1 1 1 1 U-238 1 1 7.81 1 1 30.65 1 1 3.12 1 1 -1.92 1 1 39.67 1 1 1 1 Pu-239 1 1 -0.37 1 1 -0.03 1 1 9.19 1 1 -5.97 1 1 2.8 1 1 JENDL-4.0 1 1 Fe-56 1 1 0.37 1 1 -13 1 1 - 1 1 6.14 1 1 -6.48 1 1 1 1 U-238 1 1 4.21 1 1 7.13 1 1 -6.48 1 1 -13.28 1 1 -8.43 1 1 1 1 Pu-239 1 1 -2.91 1 1 -2.75 1 1 10.5 1 1 -1.1 1 1 3.73 1 1 Nuclear data 1 1 Isotope by ENDF-B/VIII.0 4 1 Reaction 1 1 Resultant change 1 1 1 1 1 1 ( n , el ) 1 1 ( n , inl ) 1 1 ( n , f ) 1 1 ( n,g ) 1 1 1 1 ENDF/B-VII.1 1 1 Fe-56 1 1 1.24 1 1 -4.13 1 1 - 1 1 5.6 1 1 2.7 1 1 1 1 U-238 1 1 13.17 1 1 23.21 1 1 4.17 1 1 7.03 1 1 47.59 1 1 1 1 Pu-239 1 1 0.75 1 1 -0.02 1 1 17.96 1 1 -0.05 1 1 18.64 1 1 JENDL-4.0 1 1 Fe-56 1 1 2.24 1 1 -5.21 1 1 - 1 1 10.54 1 1 7.58 1 1 1 1 U-238 1 1 1.96 1 1 7.48 1 1 -8.46 1 1 5.26 1 1 6.22 1 1 1 1 Pu-239 1 1 1.73 1 1 0.3 1 1 11.08 1 1 -1.32 1 1 11.79 5.4 Kinetic Parameters 2 The calculated kinetic parameters, i.e., a prompt neutron generation time Λ and the total effective delayed neutron fraction β eff at BOC and EOC, are shown in Figures 9 Figures 9 and 10 10 , respectively. The parameters are adjoint-weighted values calculated using the iterated fission probability method implemented in SERPENT2. The prompt neutron generation time of the metallic core is shorter than that of the oxide core owing to the harder spectrum in the metallic core, whereas the total effective delayed neutron fraction of the metallic core is higher than that of the oxide core because of the lower amount of Pu and the minor actinide in the metallic core. It should be noted that the three different libraries produce similar and consistent kinetic parameters. Meanwhile, the β eff reported [ 3 3 ] is slightly higher than the calculated values. The benchmark results are summarized in Table 12 Table 12 . 1 1 Core 1 1 Condition 1 1 Reported result(average value) [ 3 3 ][pcm] 1 1 ENDF/B-VIII.0[pcm] 1 1 Metallic 1 1 BOC 1 1 345 ± 10 1 1 331 ± 3 1 1 1 1 EOC 1 1 344 ± 12 1 1 329 ± 3 1 1 Oxide 1 1 BOC 1 1 333 ± 15 1 1 320 ± 3 1 1 1 1 EOC 1 1 334 ± 13 1 1 318 ± 3 5.5 Doppler Constant K D 2 Figure 11 Figure 11 shows the Doppler constant K D for both cores at the BOC and EOC. It is defined as the difference in reactivity when the fuel temperature is doubled and at a normal fuel operating temperature. It is clearly shown that the oxide core has a more negative K D than the metallic core owing to the softer spectrum in the oxide core. Moreover, the oxide core contains more U-238 than the metallic core. It is also noted that the three different libraries have similar and consistent K D values. Meanwhile, the reported K D [ 3 3 ] is also slightly higher than the calculated value. Table 13 Table 13 summarizes the benchmark results. 1 1 Core 1 1 Condition 1 1 Reported result(average value) [ 3 3 ](pcm/K) 1 1 ENDF/B-VIII.0(pcm/K) 1 1 Metallic 1 1 BOC 1 1 -346 ± 44 1 1 -337 ± 14 1 1 1 1 EOC 1 1 -348 ± 36 1 1 -340 ± 15 1 1 Oxide 1 1 BOC 1 1 -730 ± 70 1 1 -713 ± 12 1 1 1 1 EOC 1 1 -718 ± 74 1 1 -681 ± 13 5.6 CVR 2 The CVR is summarized in Figure 12 Figure 12 for both cores at the BOC and EOC. It is defined as the difference in reactivity when the coolant is voided and under normal conditions. It is clearly shown that the metallic core is more positive (CVR) than the oxide core because of the harder spectrum shown in the metallic core. It should also be noted that the three different libraries have similar and consistent CVR values. Meanwhile, the reported CVR values [ 3 3 ] are lower than the calculated values except for the oxide core using ENDF/B-VII.1. Table 14 Table 14 summarizes the benchmark results. 1 1 Core 1 1 Condition 1 1 Reported result(average value) [ 3 3 ][pcm] 1 1 ENDF/B-VIII.0[pcm] 1 1 Metallic 1 1 BOC 1 1 2024 ± 407 1 1 2212 ± 10 1 1 1 1 EOC 1 1 2146 ± 435 1 1 2398 ± 10 1 1 Oxide 1 1 BOC 1 1 1831 ± 228 1 1 1897 ± 9 1 1 1 1 EOC 1 1 1922 ± 220 1 1 2039 ± 9 5.7 CR Worth 2 The last result, which is the CR worth, is shown in Fig. 13 Fig. 13 , which is calculated as the reactivity difference when all control subassemblies are withdrawn and inserted. It should also be noted that the three different libraries give similar and consistent CR values. Meanwhile, the reported CR value [ 3 3 ] is higher than the calculated values. Table 15 Table 15 summarizes the benchmark results. 1 1 Core 1 1 Condition 1 1 Reported result(average value) [ 3 3 ][pcm] 1 1 ENDF/B-VIII.0[pcm] 1 1 Metallic 1 1 BOC 1 1 19697 ± 2087 1 1 18547 ± 12 1 1 1 1 EOC 1 1 20497 ± 2228 1 1 19240 ± 14 1 1 Oxide 1 1 BOC 1 1 21605 ± 2021 1 1 20092 ± 12 1 1 1 1 EOC 1 1 22226 ± 2157 1 1 20610 ± 13 6 Summary and Conclusions 1 In this study, the impact of the newly released ENDF/B-VIII.0 nuclear data library on the neutronics and kinetics parameters of two different spectra, namely, metallic and oxide SFR cores, was studied. The newly released Fe-56 cross-section, particularly its ( n,inl ) and ( n,g ) cross-sections, show a significant reduction in k eff compared to those with ENDF/B-VII.1 and JENDL-4.0. In addition, a sensitivity analysis study revealed the major contributors to the resultant change in k eff . Another significant impact includes the U-238 of ENDF/B-VII.1 and Pu-239 of JENDL-4.0. By contrast, the values of the kinetic parameters, Doppler constant, CVR, and CR worth between libraries were found to be consistent. References N.E. Stauff , T.K. Kim , T.A. Taiwo et al ., Benchmark for Neutronic Analysis of Sodium-cooled Fast Reactor Cores with Various Fuel Types and Core Sizes. NEA-NSC-R-2015-9 , Organisation for Economic Co-Operation and Development ( 2016 ). Benchmark for Neutronic Analysis of Sodium-cooled Fast Reactor Cores with Various Fuel Types and Core Sizes. NEA-NSC-R-2015-9 L. Buiron , G. Rimpault , B. Fontaine et al ., Evaluation of large 3600 MWth sodium-cooled fast reactor OECD neutronic benchmarks . In Proc. PHYSOR 2014 ( 2014 ). Evaluation of large 3600 MWth sodium-cooled fast reactor OECD neutronic benchmarks L. Buiron , L., G. Rimpault , B. Fontaine et al ., Evaluation of large 1000 MWth sodium-cooled fast reactor OECD neutronic benchmarks . In Proc. PHYSOR 2014 ( 2014 ) Evaluation of large 1000 MWth sodium-cooled fast reactor OECD neutronic benchmarks M.B. Chadwick , P. Obložinský , M. Herman et al ., ENDF/B-VII. 0: next generation evaluated nuclear data library for nuclear science and technology .  Nuclear Data Sheets ,  107 ( 12 ), 2931 - 3060 ( 2014 ). https://doi.org/10.1016/j.nds.2006.11.001 https://doi.org/10.1016/j.nds.2006.11.001 ENDF/B-VII. 0: next generation evaluated nuclear data library for nuclear science and technology M.B. Chadwick , M. Herman , M., P. Obložinský et al ., ENDF/B-VII. 1 nuclear data for science and technology: cross sections, covariances, fission product yields and decay data .  Nuclear Data Sheets   112 ( 12 ), 2887 - 2996 ( 2011 ). https://doi.org/10.1016/j.nds.2011.11.002 https://doi.org/10.1016/j.nds.2011.11.002 ENDF/B-VII. 1 nuclear data for science and technology: cross sections, covariances, fission product yields and decay data A. Koning , R. Forrest , M. Kellett et al ., The JEFF-3.1 nuclear data library-JEFF report 21. NEA-6190 , Organisation for Economic Co-operation and Development ( 2006 ). The JEFF-3.1 nuclear data library-JEFF report 21. NEA-6190 K. Shibata , O. Iwamoto , T. Nakagawa et al ., JENDL-4.0: a new library for nuclear science and engineering .  Journal of Nuclear Science and Technology   48 ( 1 ), 1 - 30 ( 2011 ). https://doi.org/10.1080/18811248.2011.9711675 https://doi.org/10.1080/18811248.2011.9711675 JENDL-4.0: a new library for nuclear science and engineering D.A. Brown , M.B. Chadwick , R. Capote , et al . ENDF/B-VIII. 0: the 8th major release of the nuclear reaction data library with CIELO-project cross sections, new standards and thermal scattering data .  Nuclear Data Sheets   148 , 1 - 142 ( 2018 ). https://doi.org/10.1016/j.nds.2018.02.001 https://doi.org/10.1016/j.nds.2018.02.001 ENDF/B-VIII. 0: the 8th major release of the nuclear reaction data library with CIELO-project cross sections, new standards and thermal scattering data J. Leppänen , M. Pusa , T. Viitanen , et al ., 2014 . The Serpent Monte Carlo code: Status, development and applications in 2013 . Annals of Nuclear Energy 82 , 142 - 150 ( 2015 ). https://doi.org/10.1016/j.anucene.2014.08.024 https://doi.org/10.1016/j.anucene.2014.08.024 The Serpent Monte Carlo code: Status, development and applications in 2013 E. Woodcock , T. Murphy , P. Hemmings , et al ., Techniques used in the GEM code for Monte Carlo neutronics calculations in reactors and other systems of complex geometry . In  Proc. Conf. Applications of Computing Methods to Reactor Problems, ANL-7050 , 557 , Argonne National Laboratory ( 1965 ) Techniques used in the GEM code for Monte Carlo neutronics calculations in reactors and other systems of complex geometry J. Leppänen , Performance of Woodcock delta-tracking in lattice physics applications using the Serpent Monte Carlo reactor physics burnup calculation code .  Annals of Nuclear Energy 37 ( 5 ), 715 - 722 ( 2010 ). https://doi.org/10.1016/j.anucene.2010.01.011 https://doi.org/10.1016/j.anucene.2010.01.011 Performance of Woodcock delta-tracking in lattice physics applications using the Serpent Monte Carlo reactor physics burnup calculation code B. Becker , R. Dagan , G. Lohnert . Proof and implementation of the stochastic formula for ideal gas, energy dependent scattering kernel .  Annals of Nuclear Energy 36 ( 4 ), 470 - 474 ( 2009 ). https://doi.org/10.1016/j.anucene.2008.12.001 https://doi.org/10.1016/j.anucene.2008.12.001 Proof and implementation of the stochastic formula for ideal gas, energy dependent scattering kernel C.H. Lee , W.S. Yang . MC2-3: multigroup cross section generation code for fast reactor analysis, ANL/NE-11-41 Rev. 2 , Argonne National Laboratory ( 2013 ). MC2-3: multigroup cross section generation code for fast reactor analysis, ANL/NE-11-41 M. Aufiero , A. Bidaud , M. Hursin et al . A collision history-based approach to sensitivity/perturbation calculations in the continuous energy Monte Carlo code Serpent . Annals of Nuclear Energy 85 , 245 - 258 . https://doi.org/10.1016/j.anucene.2015.05.008 https://doi.org/10.1016/j.anucene.2015.05.008 A collision history-based approach to sensitivity/perturbation calculations in the continuous energy Monte Carlo code Serpent 10.1007/s41365-020-00831-5 This work was supported by the Research Institute of Science and Engineering at the University of Sharjah (No. 1802040790-P). 2020-05-19 2020-09-14 2020-10-28 2020-12-04 2020-12-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2020 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-020-00831-5 http://dx.doi.org/10.1007/s41365-020-00831-5 . 2020 Nuclear Science and Techniques 12 31