Nuclear data uncertainty propagation to the main physical parameters of NUR research reactor

https://doi.org/10.1016/j.pnucene.2020.103557Get rights and content

Highlights

  • S/U analysis was applied on NUR nuclear research reactor

  • The sensitivity coefficients were calculated and confirmed by the direct perturbation method.

  • The total uncertainty in keff due to cross-section covariance data was obtained.

  • The main nuclide-reaction pairs contributing to the keff uncertainty were identified.

Abstract

An accurate nuclear reactor model with less statistical uncertainty is not enough to have an accurate assessment of the calculated physical parameter value. It is worthwhile to note that another source of uncertainties originating from nuclear data uncertainties must be quantified. To realize such an assessment, sensitivity and uncertainty analysis based on first-order perturbation theory implemented in the SCALE code system was applied to the IV.N (an optimized one in series of configurations) core configuration of the Algerian light water Nuclear Uranium Reactor (NUR). The impact of nuclear data uncertainties on the effective multiplication factor was analyzed. Criticality calculation was performed by using the KENO V.a Monte Carlo code of the SCALE code system and the ENDF/B libraries with 238 and 252 energy group structures. The results were compared to the available measurement and validated. The accuracy of this criticality calculation was judged sufficient to do a sensitivity and uncertainty calculation for the same model.

First, sensitivity and uncertainty calculations were performed using the TSUNAMI3D module of the SCALE code system and ENDF/B-VII.0 library with the covariance data in 44 and 56 energy group structures. Next, the model was updated to the latest library used in SCALE code system version 6.2.3; ENDF/B-VII.1 with 252 energy group and the corresponding 56 group covariance data.

The uncertainty results of NUR effective multiplication factor due to uncertainties in nuclear data were determined. The top contributor is ν (the average number of neutrons per fission) of 235U for all cases. The total uncertainty in keff due to uncertainties in the cross-section data was determined to be around 0.84% using 56 group covariance data. This uncertainty was higher than the uncertainty when using 44 group covariance data.

Introduction

An accurate nuclear reactor model depends on the identification and minimization of the whole errors committed in the modeling approximation, the numerical solution and the input parameters entered (geometries, materials …). Errors in the input parameters, mainly the nuclear data (microscopic cross sections, fission spectra, neutron yield ν , and scattering distributions …) are considered as the most significant source of uncertainties in the reactor physics calculations.

Sensitivity and Uncertainty (S/U) analysis methods were used to quantify the impact of nuclear data uncertainties on effective multiplication factor and reactivity responses. Those methods were first developed for application to fast reactor studies in the 1970s (Weisbin et al., 1976). Several works have revitalized and updated the available S/U computational capabilities for computing sensitivity coefficients such as the adjoint-weighted perturbation capability in MCNP6 (Kiedrowski and Favorite, 2010), a discrete ordinates sensitivity formulation of the first-order perturbation theory in the SUSD3D code (I. Kodeli, 2001) and a fine group calculation from TSUNAMI-3D (Rearden, 2009).

In this work, sensitivity and uncertainty analysis (Broadhead et al., 2004) based on first-order perturbation theory implemented in TSUNAMI-3D of SCALE code system (Oak Ridge National Laboratories, 2018) was applied on IV.N core configuration of NUR research reactor in order to know the nuclear data uncertainty propagation on effective multiplication factor.

First, for the optimized reactor core configuration “IV.N″, the criticality calculation was performed using the KENO V.a Monte Carlo code of the SCALE code system. The main calculated physical parameters out coming from the model were compared to experimental data (Meftah et al., 2006; Madariaga et al., 1989) and show a good agreement (Sellaoui et al., 2014).

Then, sensitivity and uncertainty calculations were performed by TSUNAMI3D using ENDF/B-VII.0 library with the covariance data in 44 and 56 energy groups. Next, the model was updated to the earliest library; ENDF/B-VII.1 with 252 energy group and the corresponding 56 group covariance data (Oak Ridge National Laboratories, 2018).Sensitivity coefficients obtained by sensitivity and uncertainty analysis were confirmed by the direct perturbation sensitivity calculation. Those sensitivities of a computed keff value to cross-section data were coupled with cross-section-covariance data to produce uncertainty in keff.

Section snippets

Code flowchart calculation

In this work, the transport calculation was performed using KENO V.a (the three-dimensional (3D) Monte Carlo criticality transport program) (Busch and Bowman, 2005) of the SCALE 6.2.3 code system. Then sensitivity and uncertainty calculation was performed with TSUNAMI-3D (Tools for Sensitivity and Uncertainty Analysis Methodology Implementation in Three Dimensions) of the same version of the SCALE code system (Oak Ridge National Laboratories, 2018).

TSUNAMI-3D is a control module of the SCALE

NUR reactor model

NUR (Nuclear Uranium Reactor) is a 1 MW nuclear pool-type reactor. Its core is composed of 17 MTR (materials testing reactor) type fuel elements enriched approximately to 19.7%. The reactor uses light water as coolant and moderator, then the graphite blocks as reflector. The reactivity control system of the reactor is made of five Ag–In–Cd absorbing rods: (C1, C2, C3, and C4) and one fine regulating rod (F) (Meftah et al., 2006).

Criticality calculation was performed to the IV-N core

Results and discussions

Criticality calculation followed by sensitivities and uncertainties analysis was performed to the IV.N core configuration model using the number of 10,000 active generations and 30,000 particles per generation for all cases.

Conclusion

Sensitivity and uncertainty analysis of effective multiplication factor for the VI.N configuration with fresh fuel of NUR research reactor was performed using a perturbation method such as implemented in TSUNAMI-3D of SCALE code system.

Results show that keff is most sensitive to the data of the nuclides 1H, 235U, 16O, and C-graphite. These nuclides are the main component in this fresh fuel configuration: Moderator, fissile isotope and reflector. The highest sensitivity of the system is due to (ν

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgements

We thank the scientific and technical staff of the NUR Reactor Division, Commissariat of Atomic Energy and University of Sciences and Technology Houari Boumediene who contributed to this work.

References (14)

  • B. Meftah et al.

    Neutron flux optimization in irradiation channels at NUR research reactor

    Ann. Nucl. Energy

    (2006)
  • B.L. Broadhead

    Uncertainty Analysis Method for S/U Criticality Validation Techniques

    (2000)
  • B.L. Broadhead et al.

    Sensitivity- and uncertainty-based criticality safety validation techniques

    Nucl. Sci. Eng.

    (2004)
  • R.D. Busch et al.

    KENO V.a : Primer: A Primer for Criticality Calculations with SCALE/KENO V.a Using GeeWiz. ORNL/TM-2005/135

    (2005)
  • B.C. Kiedrowski et al.

    Verification of K-Eigenvalue Sensitivity Coefficient Calculations Using Adjoint-Weighted Perturbation Theory in MCNP,” LA-UR-10-04285

    (2010)
  • I. Kodeli

    Multidimensional deterministic nuclear data sensitivity and uncertainty code system, method and application

    Nucl. Sci. Eng.

    (2001)
  • M. Madariaga et al.

    Calculation analysis of the neutronic experimental data coming from the NUR reactor start-up, Berlin, Germany, CPC press

    (1989)
There are more references available in the full text version of this article.

Cited by (0)

View full text