Issue |
EPJ Nuclear Sci. Technol.
Volume 8, 2022
|
|
---|---|---|
Article Number | 3 | |
Number of page(s) | 7 | |
DOI | https://doi.org/10.1051/epjn/2021028 | |
Published online | 17 January 2022 |
https://doi.org/10.1051/epjn/2021028
Regular Article
Impact of H in H2O thermal scattering data on criticality calculation: uncertainty and adjustment
1
Reactor Physics and Thermal hydraulic Laboratory, Paul Scherrer Institut,
Villigen, Switzerland
2
International Atomic Energy Agency, Vienna International Centre,
PO Box 100,
1400
Vienna, Austria
* e-mail: dimitri-alexandre.rochman@psi.ch
Received:
20
September
2021
Received in final form:
29
November
2021
Accepted:
1
December
2021
Published online: 17 January 2022
In this paper, the impact of the thermal scattering data for H in H20 is estimated on criticality benchmarks, based on the variations of the CAB model parameters. The Total Monte Carlo method for uncertainty propagation is applied for 63 keff criticality cases, sensitive to H in H20. It is found that their impact is of a few tenth of pcm, up to 300 pcm maximum, and showing highly non-linear distributions. In a second step, an adjustment is proposed for these thermal scattering data, leading to a better agreement between calculated and experimental keff values, following an increase of scattering contribution. This work falls into the global approach of combining advanced theoretical modelling of nuclear data, followed by possible adjustment in order to improve the performances of a nuclear data library.
© D. Rochman et al., Published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The evaluation of neutron-induced nuclear data, e.g. cross sections, fission yields, spectra, is a continuous process, mainly because the nuclear reaction theory is still improving, new measurements are frequently available, and because there is a need from the user community. Examples of recently evaluated nuclear data are the Nuclear Energy Agency JEFF-3.3 library [1], the US ENDF/B-VIII.0 library [2], the Japanese JENDL-4.0 library [3] and finally the TENDL libraries [4]. The need from users might concern better estimations of current reactor quantities (boron letdown curves, power maps) [5–7], spent fuel quantities (decay heat, source terms) [8–12], or quantities for advanced systems [13–15]. A commonly-used neutronics parameter is the neutron multiplication factor, or keff, which plays for instance a crucial role in spent fuel management and in burnup credit [16–18]. The validation of existing simulation codes and the assessment of the uncertainties from various sources is a necessity, often imposed by national and international regulations (see Refs. [19–21] and references within). In this context, the impact of the nuclear data uncertainties needs to be evaluated in a regular time manner, with respect to the release of new nuclear data libraries and their associated covariance information.
This is precisely the first goal of this paper: providing uncertainties for a number of criticality-safety benchmarks [22], due to the H in H2O thermal scattering data (also called Thermal Scattering Laws, or TSL). These TSL are from the latest evaluations based on the CAB model [23], and included in the ENDF/B-VIII.0 and TENDL-2021 libraries. Regarding their uncertainties, reference [24] is used as a reference with a number of modifications (see Sect. 2). The uncertainty propagation method is the Total Monte Carlo method [25], leading to keff uncertainties (later called Δkeff), corresponding to one standard deviation of the observed probability density functions (pdf), as well as higher moments. Similar works on uncertainty propagation and covariance generation for TSL can be found in references [26,27]. In the first reference, the TMC method (and fast TMC) was applied, whereas the second reference defines TSL covariance matrices, with an application to UO2 and MOX fuels.
The second goal of this paper is to propose adjusted TSL, or at least to demonstrate the feasibility of the adjustment, based on the result of the TMC method. Such possibility was already explored almost ten years ago and presented in reference [28] for the TSL based on a different evaluation model. The method, called the Petten method, simply proposes the best performing sampled (random) TSL, based on the uncertainty propagation for criticality benchmarks, as the new recommended evaluation. Such method was also applied to Cu and 239Pu [29,30]. In this approach, the best performing sampled TSL is the one leading to the smallest χ2. Other methods of adjustments have been proposed since, apart from the traditional GLLS method (Generalized Linear Least Squares), based on non-linear Monte Carlo parameter search, see for instance references [31–36], the Unified Monte Carlo method [37] and references [38,39] for a comparison and review between GLLS and the Monte Carlo-based methods. Naturally, based on the first part of this study (uncertainty propagation), we are calculating χ2 values for each sampled TSL, and then proposing the best performing one as an example of adjustment.
Details of the study are provided in the next sections: first a recall on the sampling method of the TSL, followed by uncertainty values for a number of criticality benchmarks. Additionally, other moments are obtained, such as the skewness and kurtosis, indicating a non-linear propagation of uncertainties from the TSL to keff. Finally, the last section will analyze the χ2 distribution and draw conclusions and a potential adjustment of the TSL.
2 Variation of thermal scattering data
The same approach and sampled TSL files were presented and used in reference [40]. It is basically based on the sampling of the parameters of the CAB model [23]. The nominal (not perturbed) parameters of this model can be found in reference [41] and as mentioned, such TSL were included in some of the modern nuclear data libraries. The variations of the model parameters are derived from reference [42], with a number of modification, following the recommendations of the authors. Such variations are as follows:
σs (elastic cross section): ± 0.2%
Δ (scaling factor): ± 10%
C (diffusion coefficient): ± 0.5%
ωt (translational weight): ± 15%
ωc (continuous spectrum weight), ω1 and ω2 (1st and 2nd oscillator weights: correlated with ωt based on equation (1)
E1 (first oscillator energy): ± 5%, and
E2 (second oscillator energy): ± 30%.
For the oscillator weights ωc, ω1 and ω2, a full correlation is considered with the following equation: (1)
Such approach follows the recommendation of reference [24]. All variations are a Normal distributions. Uncertainties correspond to one standard deviation (1σ). As in the case of reference [40], 3000 sampled TSL are produced, and processed with NJOY into ACE files [43].
3 Results
The impact of the H in H2O thermal scattering data are here assessed on criticality-safety benchmarks. Such benchmarks are obtained from the ICSBEP handbook [22]. Such handbook contains a large number of benchmarks (simulation input files and experimental quantities such as keff), but only some of them are sensitive to TSL. Those benchmarks are mostly thermal ones (the neutron fluence is mainly in the thermal energy region), but not only. In the following, a total of 63 benchmarks are considered, among which only one is considered as fast (called pmf11-1) and one intermediate (usi1-1); the other ones are thermal (61 benchmarks). The nomenclature used here is similar to the handbook: all benchmark names consist of three letters: the first one denotes the type of fuel (“l” for low-enriched 235U, “h” for high-enriched 235U, “u” for 233U fuel, “m” for a mixture of 235U and 239Pu, and “p” for 239Pu fuel), the second one defines the physical form (“s” for solution and “c” for compound), and finally the last letter defines the neutron spectrum (“t” for thermal, “i” for intermediate and “f” for fast). All benchmark names are continued with an additional information (after a dash), indicating its subnumber, referring to the ICSBEP cases.
The benchmark calculations are performed with MCNP6; apart from the sampled TSL ACE files, the other nuclear data files correspond to the beta version of the JEFF-4.0 library. The method of uncertainty propagation is the Total Monte Carlo (or TMC) [25]; it was already presented in many publications and it is summarized here as follows: it consists in repeating i times the same calculation, each time with a different sampled quantity of interest as input. In other words, each random file differs slightly from the previous one and averaging over all of them would return a credible standard deviation and covariance matrix. The averaging however takes place after the applied calculations. In our case, the quantity of interest corresponds to the TSL data, as explained in Section 2. In this work, i is equal to 3000 and it allows to achieve an acceptable convergence for various moments of the observed keff probability density functions. In all uncertainty values, the statistical uncertainties from the MCNP6 is already removed, assuming that the observed variance is the sum of the effect of the TSL and statistics. Finally, together with the uncertainty Δ keff (or one standard deviation), the quantity MAD (for Median Absolute Deviation) is also provided. This is another measure of the width of a distribution, less affected by extreme values (both the standard deviation and MAD values will significantly differ for skewed distributions; for a Normal distribution, the MAD is equal to 0.67 of the standard deviation). The definition of the MAD is Median(|xi-), with xi the observed realization for the sample i, and the average over all i. Results for the average, uncertainty, skewness and kurtosis for each considered benchmark are presented in Table 1.
Uncertainties (standard deviation and MAD), skewness and kurtosis on criticality benchmarks due to thermal scattering data of H in H20 obtained from MCNP6 calculations. The base library (except the thermal scattering data) is JEFF-4.0β. The statistical uncertainties are already removed.
3.1 Uncertainties and other moments
Compared to uncertainties calculated in reference [28], the present values are in general smaller. This is not surprising as the model parameters in reference [28] were modified in an ad-hoc way, being one of the first trials of the TMC method on TSL. As observed, the maximum uncertainties are less than 300 pcm (and in general less than 100 pcm). It is worth noticing that for skewed distributions with a strong kurtosis, as the majority of cases presented here, one standard deviation does not refer to 68% of the population: for lct8.5, 68.3% of the cases are within the mean ± 1σ, but for lst2.1, 73% are within the mean ± 1σ. One can also compare the present results for the lct and pst benchmarks with the values presented in reference [27]. In this reference, uncertainties due to TSL for the MISTRAL-1 experiment, similar to the lct benchmarks, range from 70 to 150 pcm. Such uncertainties were obtained by a direct perturbation approach, based on different model uncertainties, and are of the same order of magnitude compared to the present values. In the case of the MOX cases, reference [27] reports uncertainties between 110 and 200 pcm for the MISTRAL-2 system, which can be compared to the present values for the pst benchmarks. In this case too, uncertainties are relatively similar.
Generally, the uncertainties (as well as skewness and kurtosis) are similar within types of benchmarks: for instance all lct benchmarks have negative skewnesses, all lct7 cases present uncertainties close to the same values. Examples of some pdf are presented in Figure 1. These distributions are sensibly different, with strong different skewnesses (for the pst and uct cases), zero skewness (the lct case), or zero kurtosis (corresponding to a flat distribution, for mct2). Given that the model parameters for the TSL production were sampled from Normal distributions, such a variety of keff pdf indicates non linear cases of uncertainty propagation.
Similarities in pdf within families of benchmarks are certainly related to the fact that these benchmarks share similar geometry, enrichment, and possibly reflectors. This can be observed for the majority of cases, but with a number of exceptions, as for the uct1 benchmarks and the cases 3 and 7: uncertainties strongly differ, and skewnesses are far from zero with opposite signs. Both benchmarks are also strongly anti-correlated, with a Pearson correlation factor of −0.8. These two cases nevertheless present very different physical systems (e.g. rectangular and hexagonal cores with different reflectors).
The values from the lct benchmarks (being clusters of UO2 rods in water) can be compared to the k∞ uncertainties from reference [40] (being for a UO2 pincell, 3.5% enriched, burned up to 100 MWd/kgU; same sampled TSL files were used), where the maximum uncertainty was about 30 pcm. In the case of the lct benchmarks, the maximum is about 80 pcm. These two values differ, but are also relatively close, given that a direct comparison is not possible (different keff values, different enrichments, and different geometry).
Based on these random keff values (3000for each benchmark), the adjustment step can be performed, as presented in the next section.
Fig. 1 Probability density function for a number of criticality benchmarks, obtained by randomly changing the H in H2O thermal scattering data. |
3.2 Adjustment
The previous results were obtained based on 3000 sampled TSL; they all differ and lead to different benchmark keff values. Let’s i denotes a specific sampled TSL, and n the benchmark case (n goes from 1 to N = 63). For each i, a χ2 value can be calculated as follows: (2)
ki,n,calc and kn,exp correspond to the calculated and experimental keff values for the benchmark n, respectively. For instance, for i = 0, which corresponds to the beta version of the future JEFF-4.0 library, called JEFF-4T0, . By considering all i sampled cases, the obtained χ2 distribution can be calculated and is presented in Figure 2. The shape of this distribution is typical of χ2 values, as seen in references [28–30]. As it can be easily understood, some sampled TSL do not perform well (high χ2), whereas some of them are performing better than the original evaluated TSL. Based on such distribution, the Petten method is excessively simple: select the sampled TSL with the lowest χ2. Naturally, as the TSL have a limited effect on the benchmark calculated values, the low tail of the χ2 distribution does not reach χ2 = 0, but as observed, noticeable improvements in terms of C/E are possible. Out of the 3000 sampled TSL, the one with the smallest χ2 is number 1871. The calculated keff from this specific TSL are presented in Figure 3, compared with the nominal (unperturbed) TSL evaluation. The changes on keff values from the nominal calculations and from the TSL file number 1871 are not that striking, but the χ2 values were definitively reduced. One has to keep in mind that the TSL induces low uncertainties on keff, meaning that their impact is nevertheless limited. One can also repeat this work and double the uncertainties on the CAB parameters; this way, it is likely to find an even better solution. Naturally, this would imply that the TSL are “responsible” for more than they should be, given that other nuclear data can also be adjusted. Users who are not able or prepared to do the full uncertainty quantification with all 3000 random files, can then at least use random file number 1871 to get the best result.
This specific TSL run number 1871 was obtained with different CAB model parameters, see Table 2 (parameters ωc, ω1 and ω2 are fully correlated with ωt using Eq. (1)). The parameters for the file number 1871 are one example of adjustments with integral benchmarks. Table 2 is also showing parameters for random files performing almost as good as number 1871. As observed, their parameters are sensibly different, and it is difficult to extract a pattern. The E2 for run 1871 might not be suitable, as such value is not valid for the second oscillator, which was closer to 400 meV. To avoid such result, a rejection condition could be used, based on water physics. This indicates that the adjustment solely based on this number of integral benchmarks might not be enough to obtain a unique solution, either because there is not a unique solution, or because the number of benchmarks is not enough. Additionally, the χ2 definition can take into account differential data (measured cross sections), thus leading to a compromise between integral and differential agreements. More advanced method can also be applied, based on Bayesian Monte Carlo, or methods borrowed from Machine Learning, e.g. separating integral and differential data between a training set and a validation set. Such more advanced methods can be applied in future work.
Fig. 2 Distribution of the χ2 per degree offreedom considering sampled (random) H in H2O thermal scattering data. The vertical black line indicates the position of χ2∕N = 1. |
Fig. 3 Calculated over Experimental values for the considered criticality benchmarks. The case “Best” corresponds to a selected H in H2O thermal scattering data case. |
Nominal (unperturbed) CAB parameters (file number 0) and sampled parameters leading to the minimum χ2 (file number 1871).
4 Conclusion
This work presents an estimation of the impact of the thermal scattering data for H in H2O (or TSL) on keff of criticality benchmarks. The method used is the sampling of the CAB model parameters, leading to sampled thermal scattering data in the form of ACER files, and finally to sampled keff values. It was found that the impact of such TSL is modest compared to other nuclear data (less than 300 pcm). It was also found that the effect of TSL is highly non linear. Based on these results, a selection of “best performing” sampled TSL is proposed, using a simple χ2 definition. In conclusion, this paper is an application of the TMC and Petten methods towards nuclear data adjustment, more precisely thermal scattering data adjustment. If used in combination of traditional nuclear data evaluation methods, it can help to reach a better agreement between simulated and experimental integral values.
Author contribution statement
The authors have equally contributed to the various steps of technical work and to the preparation of the paper.
Acknowledgements
We would like to thank J.I. Marquez Damian from the European Spallation Source and G. Noguere from the CEA Cadarache for their support and efforts in helping for the production of random TSL data. This work would not have been possible without their help.
It was also partly funded by the European Union’s Horizon 2020 Research and Innovation Programme with the SANDA project, under grant agreement No 847552.
References
- A.J.M. Plompen, et al., The joint evaluated fission and fusion nuclear data library, JEFF-3.3, Eur. Phys. J. A 56, 181 (2020) [Google Scholar]
- D.A. Brown 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, Nucl. Data Sheets 148, 1 (2018) [CrossRef] [Google Scholar]
- K. Shibata et al., JENDL-4.0: a new library for nuclear science and engineering, J. Nucl. Sci. Technol. 48, 1 (2011) [Google Scholar]
- A. Koning, D. Rochman, J. Sublet, N. Dzysiuk, M. Fleming, S. van der Marck, TENDL: complete nuclear data library for innovative nuclear science and technology, Nucl. Data Sheets 155, 1 (2019) [Google Scholar]
- A. Aures, W. Bernnat, F. Bostelmann, J. Bousquet, B. Krzykacz-Hausmann, A. Pautz, K. Velkov, W. Zwermann, Reactor simulations with nuclear data uncertainties, proceedings of the ANS Best Estimate Plus Uncertainty International Conference (BEPU 2018), BEPU2018-KN-E7, Real Collegio, Lucca, Italy, May 13-19, 2018 [Google Scholar]
- D. Rochman, A. Dokhane, A. Vasiliev, H. Ferroukhi, M. Hursin, Nuclear data uncertainties for core parameters based on Swiss BWR operated cycles, Ann. Nucl. Energy 148, 107727 (2020) [CrossRef] [Google Scholar]
- O. Cabellos, E. Castro, C. Ahnet, C. Holgado, Propagation of nuclear data uncertainties for PWR core analysis, Nucl. Eng. Technol. 46, 299 (2014) [CrossRef] [Google Scholar]
- I.C. Gauld, M.L. Williams, F. Michel-Sendis, J.S. Martinez, Integral nuclear data validation using experimental spent nuclear fuel compositions, Nucl. Eng. Technol. 49, 1226 (2017) [CrossRef] [Google Scholar]
- M. Frankl, M. Hursin, D. Rochman, A. Vasiliev, H. Ferroukhi, Nuclear data uncertainty quantification in criticality safety evaluations for spent nuclear fuel geological disposal, Appl. Sci. 11, 6499 (2021) [CrossRef] [Google Scholar]
- D. Rochman, A. Vasiliev, A. Dokhane, H. Ferroukhi, Uncertainties for Swiss LWR spent nuclear fuels due to nuclear data, Eur. Phys. J. Nuclear Sciences & Technologies 4, 6 (2018) [CrossRef] [EDP Sciences] [Google Scholar]
- G. Ilas, H. Liljenfeldt, Decay heat uncertainty for BWR used fuel due to modeling and nuclear data uncertainties, Nucl. Eng. Des. 319, 176 (2017) [CrossRef] [Google Scholar]
- A. Rintala, Evaluating the Effect of Decay and Fission Yield Data Uncertainty on BWR Spent Nuclear Fuel Source Term, in Proceedings of the 29th international conference on nuclear energy for the new Europe, NENE 2020, Sept. 7–10, Portoroz, Slovenia (2020) [Google Scholar]
- I. Trivedi, J. Hou, G. Grasso, K. Ivanov, F. Franceschini, Nuclear data uncertainty quantification and propagation for safety analysis of lead-cooled fast reactors, Sci. Technol. Nucl. Instal. 3961095 (2020) [Google Scholar]
- L.A. Bernstein, D.A. Brown, A.J. Koning, B.T. Rearden, C.E. Romano, A.A. Sonzogni, A.S. Voyles, W. Younes, Our future nuclear data needs, Annu. Rev. Nucl. Part. Sci. 69, 109 (2019) [CrossRef] [Google Scholar]
- F. Bostelmann, G. Ilas, W.A. Wieselquist, Key Nuclear Data Impacting Reactivity in Advanced Reactors, Oak Ridge National Laboratory, Report ORNL/TM-2020/1557 (2020) [CrossRef] [Google Scholar]
- W.A. Wieselquist, K.S. Kim, G. Ilas, I.C. Gault, Comparison of burnup credit uncertainty quantification methods, in Proceedings of the ANS NCSD 2013 - Criticality Safety in the Modern Era: Raising the Bar, Wilmington, NC, September 29-October 3, 2013, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2013) [Google Scholar]
- A. Vasiliev, J. Herrero, M. Pecchia, D. Rochman, H. Ferroukhi, S. Caruso, Preliminary assessment of criticality safety constraints for Swiss spent nuclear fuel loading in disposal canisters, Materials 12, 494 (2019) [CrossRef] [Google Scholar]
- I.C. Gauld, U. Mertyurek, Validation of BWR spent nuclear fuel isotopic predictions with applications to burnup credit, Nucl. Eng. Des. 345, 110 (2019) [CrossRef] [Google Scholar]
- Safety assessment and verification for nuclear power plants, IAEA Safety Standards Series NS-G-1.2 (International Atomic Energy Agency, Vienna, Austria, 2009) [Google Scholar]
- Best estimate safety analysis for nuclear power plants: Uncertainty evaluation, IAEA Safety Report Series No. 52 (International Atomic Energy Agency, Vienna, Austria, 2008) [Google Scholar]
- E.J. Bonano, J.E. Meacham, G.J. Appel, Radioactive waste management: it’s not all science and engineering, in Proceedings of the conference on Management of Spent Fuel from Nuclear Power Reactors: Learning from the Past, Enabling the Future, International Atomic Energy Agency, Vienna, Austria, 2020, on CD-ROM, IAEA-CN-272/65, ID65 (2020) [Google Scholar]
- J.B. Briggs Ed., International Handbook of evaluated Criticality Safety, Benchmark Experiments, NEA/NSC/DOC(95)03/I (Organization for Economic Co-operation and Development, Nuclear Energy Agency, 2004) [Google Scholar]
- J.I. Marquez Damian, J.R. Granada, D.C. Malaspina, CAB models for water: a new evaluation of the thermal neutron scattering laws for light and heavy water in ENDF-6 format, Ann. Nucl. Energy 65, 280 (2014) [CrossRef] [Google Scholar]
- G. Noguere, J.P. Scotta, C. De Saint Jean, P. Archier, Covariance matrices of the hydrogen neutron cross sections bound in light water for the JEFF-3.1.1 neutron library, Ann. Nucl. Energy 104, 132 (2017) [CrossRef] [Google Scholar]
- A.J. Koning, D. Rochman, Towards sustainable nuclear energy: putting nuclear physics to work, Ann. Nucl. Energy 35, 2024 (2008) [CrossRef] [Google Scholar]
- L. Maul, Thermal Scattering Law Uncertainties and Propagation into Small Thermal Fission Reactors, PhD Thesis, University of New South Wales, Sydney, Australia (2018) [Google Scholar]
- J.P. Scotta, G. Noguère, J.I. Marquez Damian, Generation of the 1H in H2O neutron thermal scattering law covariance matrix of the CAB model, Eur. Phys. J. Nuclear Sciences & Technologies 4, 32 (2018) [CrossRef] [EDP Sciences] [Google Scholar]
- D. Rochman, A. Koning, Improving the H in H2O thermal scattering data using the Petten method, Nucl. Sci. Eng. 172, 287 (2012) [CrossRef] [Google Scholar]
- D. Rochman, A. Koning, Evaluation and adjustement of the neutron-induced reactions of 63,65Cu, Nucl. Sci. Eng. 170, 265 (2012) [CrossRef] [Google Scholar]
- D. Rochman, A. Koning, How to randomly evaluate nuclear data: a new method applied to 239Pu, Nucl. Sci. Eng. 169, 68 (2011) [CrossRef] [Google Scholar]
- E. Bauge, M. Dupuis, S. Hilaire, S. Peru, A.J. Koning, D. Rochman, S. Goriely, Connecting the dots, or nuclear data in the age of supercomputing, Nucl. Data Sheets 118, 32 (2014) [CrossRef] [Google Scholar]
- D. Rochman, E. Bauge, A. Vasiliev, H. Ferroukhi, Correlation nu-sigma-chi in the fast neutron range via integral information, Eur. J. Phys. Nuclear Sciences & Technologies 3, 14 (2017) [CrossRef] [EDP Sciences] [Google Scholar]
- D. Rochman, E. Bauge, A. Vasiliev, H. Ferroukhi, S. Pelloni, J.Ch. Koning, Sublet, Monte Carlo nuclear data adjustment via integral information, Eur. Phys. J. Plus 133, 537 (2018) [CrossRef] [Google Scholar]
- D. Rochman, A. Vasiliev, H. Ferroukhi, S. Pelloni, E. Bauge, A.J. Koning, Correlation nu-sigma for U-Pu in the thermal and resonance neutron range via integral information, Eur. Phys. J. Plus 134, 453 (2019) [Google Scholar]
- D. Kumar, S.B. Alam, H. Sjöstrand, J.M. Palau, C. De Saint Jean, Nuclear data adjustment using Bayesian inference, diagnostics for model fit and influence of model parameters, Eur. Phys. J. Web Conf. 239, 13003 (2020) [CrossRef] [EDP Sciences] [Google Scholar]
- C. De Saint Jean, P. Archier, E. Privas, G. Noguere, B. Habert, P. Tamagno, Evaluation of neutron-induced cross sections and their related covariances with physical constraints, Nucl. Data Sheets 148, 383 (2018) [CrossRef] [Google Scholar]
- C.W. Chapman, Thermal neutron scattering evaluation framework, PhD Thesis, Georgia Institute of Technology, USA, August 2017 [Google Scholar]
- D. Siefman, M. Hursin, D. Rochman, S. Pelloni, A. Pautz, Stochastic vs. sensitivity-based integral parameter and nuclear data adjustments, Eur. Phys. J. Plus 133, 429 (2018) [Google Scholar]
- E. Ivanov, C. de Saint Jean, V. Sobes, Nuclear data assimilation, scientific basis and current status, Eur. J. Phys. NuclearSciences & Technologies 7, 9 (2021) [Google Scholar]
- D. Rochman, M. Hursin, A. Vasiliev, H. Ferroukhi, Impact of H in H2O thermalscattering data on depletion calculation: keff, nuclide inventory and decay heat, submitted to Eur. Phys. J. Nuclear Sciences & Technologies (2021) [Google Scholar]
- tsl-HinH2O, J.I. Marquez Damian, J.R. Granada, D. Roubtsov, https://github.com/marquezj/tsl-HinH2O, September 2020 [Google Scholar]
- J.P. Scotta, G. Noguere, D. Bernard, J.I. Marquez Damian, A. Santamarina, Impact of the thermal scattering law of H in H2O on the isothermal temperature reactivity coefficients for UOX and MOX fuel lattices in cold operating conditions, Eur. Phys. J. Nuclear Sciences & Technologies 2, 28 (2016) [CrossRef] [EDP Sciences] [Google Scholar]
- R.E. MacFarlane, A.C. Kahler, Methods for processing ENDF/B-VII with NJOY, Nucl. Data Sheets 111, 2739 (2010) [CrossRef] [Google Scholar]
Cite this article as: D. Rochman, A. Vasiliev, H. Ferroukhi, A. Koning, J.Ch. Sublet, Impact of H in H2O thermal scattering data on criticality calculation: uncertainty and adjustment, EPJ Nuclear Sci. Technol. 8, 3 (2022)
All Tables
Uncertainties (standard deviation and MAD), skewness and kurtosis on criticality benchmarks due to thermal scattering data of H in H20 obtained from MCNP6 calculations. The base library (except the thermal scattering data) is JEFF-4.0β. The statistical uncertainties are already removed.
Nominal (unperturbed) CAB parameters (file number 0) and sampled parameters leading to the minimum χ2 (file number 1871).
All Figures
Fig. 1 Probability density function for a number of criticality benchmarks, obtained by randomly changing the H in H2O thermal scattering data. |
|
In the text |
Fig. 2 Distribution of the χ2 per degree offreedom considering sampled (random) H in H2O thermal scattering data. The vertical black line indicates the position of χ2∕N = 1. |
|
In the text |
Fig. 3 Calculated over Experimental values for the considered criticality benchmarks. The case “Best” corresponds to a selected H in H2O thermal scattering data case. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.