# Propagation of Cross-Section Uncertainties in Criticality Calculations in the Framework of UAM-Phase I Using MCNPX-2.7e and SCALE-6.1

Authors: C. J. Díez; J. J. Herrero; O. Cabellos; J. S. Martínez
Journal: Science and Technology of Nuclear Installations
 (2013)

---

## Abstract

In the framework of the OECD/NEA project on Benchmark for Uncertainty Analysis in Modeling (UAM) for Design, Operation, and Safety Analysis of LWRs, several approaches and codes are being used to deal with the exercises proposed in Phase I, “Specifications and Support Data for Neutronics Cases.” At UPM, our research group treats these exercises with sensitivity calculations and the “sandwich formula” to propagate cross-section uncertainties. Two different codes are employed to calculate the sensitivity coefficients ofkeff to cross sections in criticality calculations: MCNPX-2.7e and SCALE-6.1. The former uses the Differential Operator Technique and the latter uses the Adjoint-Weighted Technique. In this paper, the main results for exercise I-2 “Lattice Physics” are presented for the criticality calculations of PWR. These criticality calculations are done for a 15×15 TMI fuel assembly at four different states: HZP-Unrodded, HZP-Rodded, HFP-Unrodded, and HFP-Rodded. The results of the two different codes above are presented and compared. The comparison proves a good agreement between SCALE-6.1 and MCNPX-2.7e in keff uncertainty that comes from the sensitivity coefficients calculated by both codes. Differences are found when the sensitivity profiles are analysed, but they do not lead to differences in the keff uncertainty.


---

## Body

## 1. Introduction

As stated in the Introduction of [1], “This benchmark project (…) responds to the needs of estimating confidence bounds for the results from simulations and analysis in real applications.” Phase I deals with the neutronics cases. (i)
Exercise I-1. Cell Physics‌ focused on the derivation of the multigroup microscopic cross section libraries and their uncertainties.(ii)
Exercise I-2. Lattice Physics‌ focused on the derivation of the few-group macroscopic cross section libraries and their uncertainties.(iii)
Exercise I-3. Core Physics‌ focused on the core steady state stand-alone neutronics calculations and their uncertainties. For these calculations, the main source of uncertainty taken into account is the cross section uncertainties which are propagated throughout the different simulation levels.There are mainly two different approaches to propagate uncertainties: The first one is based on a Monte Carlo approach where a large amount of calculations are performed sampling the problem parameters as random variables, and then carrying out a statistical analysis; the second one relies on sensitivity coefficients and the “sandwich formula.” The latter approach is the one employed in this work.The way of obtaining the sensitivity coefficients of the response functions depends on which code is used. In this case, two different techniques are used: the Adjoint-Weighted Technique by SCALE-6.1 [2] and the Differential Operator Technique by MCNPX-2.7e [3]. Previous works [4, 5] have compared both techniques showing their pros and cons.This work is aimed to present how the uncertainty quantification is carried out using the sensitivity approach and how the sensitivity coefficients are calculated with SCALE-6.1 and MCNPX-2.7e. Afterwards, in the framework of Exercise I-2, both codes are used to perform the uncertainty quantification on thekeff of a 15×15 TMI fuel assembly. Finally, their results are compared and conclusions are drawn.


## 2. Uncertainty Quantification Based on Sensitivity Coefficients

The uncertainty quantification based on sensitivity coefficients relies on the “sandwich formula” obtained with the propagation of moments, as presented in [6, Section III.F] and summarized here.BeingR the calculated response function of a system which depends on (α1,…,αn) parameters, R can be approximated by a linear function of (α1,…,αn) using the Taylor series approximation around a nominal value (α10,…,αn0), as presented in (1), where δαi=αi-αi0:
(1)R(α1,…,αn)=R(α10,…,αn0)+∑i=1n‍(∂R∂αi)(α10,…,αn0)δαi.Then, taking the parameters of the system as random variables,R becomes a random variable for which its mean coincides with the response function at the nominal value, (2), and its variance is calculated with (3) known as the “sandwich formula” or the “sandwich rule”:
(2)E(R)=R(α10,…,αn0),(3)var(R)=SVαST,
where (i)
S=(∂R/∂α1,…,∂R/∂αn) is the vector of the sensitivity coefficients.(ii)
Vα is the covariance matrix of the system parameters defined as
(4)(Vα)ij={cov(αi,αj),i≠j,var(αi),i=j. Because usually the sensitivity coefficients are calculated as relative values,
(5)Si′=∂R/R∂αi/αi.
Equation (3) can be rewritten as (6), where the value obtained is the relative standard deviation (rel.std.dev.(R)) and Vα′ is the relative covariance matrix:
(6)(rel.std.dev.(R))2=var(R)R2=∑i=1,j=1n‍(Si′cov(αi,αj)αiαjSj′)=S′Vα′S′T.
Both values, the standard deviation var(R) provided by (3) and the rel.std.dev.(R) calculated in (6), are used as a measure of the uncertainty on the response function. Therefore, the sensitivity coefficients of the response function to the system parameters, Si or Si′, should be calculated and the covariance matrix (Vα) or the relative covariance matrix (Vα′) should be provided.


### 2.1. Calculating the Sensitivity Coefficients

The two codes which perform the criticality calculations in this paper use different methods for calculating the sensitivity coefficients necessary to carry out the uncertainty quantification:(i)
SCALE-6.1 uses the Adjoint-Weighted Technique.(ii)
MCNPX-2.7e uses the Differential Operator Technique.


#### 2.1.1. SCALE-6.1—The Adjoint-Weighted Technique

The Adjoint-Weighted Technique is used in SCALE-6.1 inside the TSUNAMI sequence, and the theory applied is stated in the SAMS module manual [2, Section F22]. Every sensitivity coefficient is calculated as a sum of two terms: explicit and implicit, as presented in
(7)(Sk)complete=(Sk)explicit+(Sk)implicit.
The explicit term, (Sk)explicit, is calculated using the Adjoint-Weighted Technique (called adjoint-based perturbation in the SAMS manual), based on the perturbation δk of the keff given in (8). There, ϕ is the neutron flux of the problem defined by [A-λB]ϕ=0, where A is the operator that represents all of the transport equations except for the fission term, B is the operator that represents the fission term, and λ represents the eigenvalues where the largest one is 1/keff. ϕ† is the adjoint neutron flux of the adjoint problem defined as [A†-λB†]ϕ†=0 where A† and B† are the adjoint operators corresponding to A and B. δA and δB are the perturbed operators due to a perturbation in the cross sections:
(8)δkeffkeff=-〈ϕ†(δA-λδB)ϕ〉〈ϕ†(λδB)ϕ〉.The implicit term,(Sk)implicit, is calculated in order to take into account the effect on keff of perturbing one cross section that affects the resonance-shielded values of other cross sections. This term appears because the transport calculations should be done using groupwise cross sections which require an initial self-shielding adjustment.In any TSUNAMI sequence, the forward and adjoint transport problems are solved in order to calculate the neutron flux and its adjoint using XSDRNPM module for TSUNAMI-1D, NEWT module for TSUNAMI-2D, and KENO-V.a or KENO-VI for TSUNAMI-3D. Then, the SAMS module is applied for calculating the sensitivity coefficients for every energy-group reaction cross section. The sensitivity to the average number of neutrons per fission,ν¯, and to the fission spectrum, χ, is estimated. Also, SAMS module performs the uncertainty quantification on the keff due to the cross section uncertainties.


#### 2.1.2. MCNPX-2.7e—The Differential Operator Technique

The Differential Operator Technique is applied in MCNPX-2.7e to calculate the changeΔkeff due to a perturbation in a cross section, Δσ, by means of the PERT card. It is based on the Taylor series expansion as presented in (9), where σx is the cross section perturbed, σx,0 is its nominal value, and Δσ=σx-σx,0, is the perturbation:
(9)Δkeff=keff(σx)-keff(σx,0)=(∂k∂σ|σx,0Δσ)1st+(12∂2k∂σ2|σx,0(Δσ)2)2nd+⋯=[Δkeff]PERT,1st+[Δkeff]PERT,2nd+⋯.
The first derivative means the derivative of the probability of the random walk occurring. The second one is the derivative of the tally response itself. There is an additional derivative for the changes in the fission source distribution that MCNPX-2.7e does not take into account, but methodologies about how to calculate it have been developed in [4, 7–9] showing also its effect. In [4, 10], it was pointed out that because the scattering affects the fission source spatial distribution more than capture, the sensitivity to scattering cross sections could be wrong. Other limitations are stated in the MCNPX manual [3, Section 5.2.2].For the sensitivity coefficient, only the first derivative is required. Using the first term provided by the PERT card (using the keywordMETHOD=+2), the sensitivity coefficient is calculated with (10), just dividing by keff,0 and px, which is the variation of the cross section from σx=σx,0(1+px). Because [Δkeff]PERT,1st is lineal with respect to the size of the perturbation px, so choosing one value for px is trivial:
(10)Si=1k0px[Δkeff]PERT,1st.
The statistical uncertainty is propagated to the sensitivity coefficient by (11), obtaining the relative standard deviation:
(11)var(Si)Si=var([Δkeff]PERT,1st)[Δkeff]PERT,1st2+var(k0)k02.


## 2.1. Calculating the Sensitivity Coefficients

The two codes which perform the criticality calculations in this paper use different methods for calculating the sensitivity coefficients necessary to carry out the uncertainty quantification:(i)
SCALE-6.1 uses the Adjoint-Weighted Technique.(ii)
MCNPX-2.7e uses the Differential Operator Technique.


### 2.1.1. SCALE-6.1—The Adjoint-Weighted Technique

The Adjoint-Weighted Technique is used in SCALE-6.1 inside the TSUNAMI sequence, and the theory applied is stated in the SAMS module manual [2, Section F22]. Every sensitivity coefficient is calculated as a sum of two terms: explicit and implicit, as presented in
(7)(Sk)complete=(Sk)explicit+(Sk)implicit.
The explicit term, (Sk)explicit, is calculated using the Adjoint-Weighted Technique (called adjoint-based perturbation in the SAMS manual), based on the perturbation δk of the keff given in (8). There, ϕ is the neutron flux of the problem defined by [A-λB]ϕ=0, where A is the operator that represents all of the transport equations except for the fission term, B is the operator that represents the fission term, and λ represents the eigenvalues where the largest one is 1/keff. ϕ† is the adjoint neutron flux of the adjoint problem defined as [A†-λB†]ϕ†=0 where A† and B† are the adjoint operators corresponding to A and B. δA and δB are the perturbed operators due to a perturbation in the cross sections:
(8)δkeffkeff=-〈ϕ†(δA-λδB)ϕ〉〈ϕ†(λδB)ϕ〉.The implicit term,(Sk)implicit, is calculated in order to take into account the effect on keff of perturbing one cross section that affects the resonance-shielded values of other cross sections. This term appears because the transport calculations should be done using groupwise cross sections which require an initial self-shielding adjustment.In any TSUNAMI sequence, the forward and adjoint transport problems are solved in order to calculate the neutron flux and its adjoint using XSDRNPM module for TSUNAMI-1D, NEWT module for TSUNAMI-2D, and KENO-V.a or KENO-VI for TSUNAMI-3D. Then, the SAMS module is applied for calculating the sensitivity coefficients for every energy-group reaction cross section. The sensitivity to the average number of neutrons per fission,ν¯, and to the fission spectrum, χ, is estimated. Also, SAMS module performs the uncertainty quantification on the keff due to the cross section uncertainties.


### 2.1.2. MCNPX-2.7e—The Differential Operator Technique

The Differential Operator Technique is applied in MCNPX-2.7e to calculate the changeΔkeff due to a perturbation in a cross section, Δσ, by means of the PERT card. It is based on the Taylor series expansion as presented in (9), where σx is the cross section perturbed, σx,0 is its nominal value, and Δσ=σx-σx,0, is the perturbation:
(9)Δkeff=keff(σx)-keff(σx,0)=(∂k∂σ|σx,0Δσ)1st+(12∂2k∂σ2|σx,0(Δσ)2)2nd+⋯=[Δkeff]PERT,1st+[Δkeff]PERT,2nd+⋯.
The first derivative means the derivative of the probability of the random walk occurring. The second one is the derivative of the tally response itself. There is an additional derivative for the changes in the fission source distribution that MCNPX-2.7e does not take into account, but methodologies about how to calculate it have been developed in [4, 7–9] showing also its effect. In [4, 10], it was pointed out that because the scattering affects the fission source spatial distribution more than capture, the sensitivity to scattering cross sections could be wrong. Other limitations are stated in the MCNPX manual [3, Section 5.2.2].For the sensitivity coefficient, only the first derivative is required. Using the first term provided by the PERT card (using the keywordMETHOD=+2), the sensitivity coefficient is calculated with (10), just dividing by keff,0 and px, which is the variation of the cross section from σx=σx,0(1+px). Because [Δkeff]PERT,1st is lineal with respect to the size of the perturbation px, so choosing one value for px is trivial:
(10)Si=1k0px[Δkeff]PERT,1st.
The statistical uncertainty is propagated to the sensitivity coefficient by (11), obtaining the relative standard deviation:
(11)var(Si)Si=var([Δkeff]PERT,1st)[Δkeff]PERT,1st2+var(k0)k02.


## 2.1.1. SCALE-6.1—The Adjoint-Weighted Technique

The Adjoint-Weighted Technique is used in SCALE-6.1 inside the TSUNAMI sequence, and the theory applied is stated in the SAMS module manual [2, Section F22]. Every sensitivity coefficient is calculated as a sum of two terms: explicit and implicit, as presented in
(7)(Sk)complete=(Sk)explicit+(Sk)implicit.
The explicit term, (Sk)explicit, is calculated using the Adjoint-Weighted Technique (called adjoint-based perturbation in the SAMS manual), based on the perturbation δk of the keff given in (8). There, ϕ is the neutron flux of the problem defined by [A-λB]ϕ=0, where A is the operator that represents all of the transport equations except for the fission term, B is the operator that represents the fission term, and λ represents the eigenvalues where the largest one is 1/keff. ϕ† is the adjoint neutron flux of the adjoint problem defined as [A†-λB†]ϕ†=0 where A† and B† are the adjoint operators corresponding to A and B. δA and δB are the perturbed operators due to a perturbation in the cross sections:
(8)δkeffkeff=-〈ϕ†(δA-λδB)ϕ〉〈ϕ†(λδB)ϕ〉.The implicit term,(Sk)implicit, is calculated in order to take into account the effect on keff of perturbing one cross section that affects the resonance-shielded values of other cross sections. This term appears because the transport calculations should be done using groupwise cross sections which require an initial self-shielding adjustment.In any TSUNAMI sequence, the forward and adjoint transport problems are solved in order to calculate the neutron flux and its adjoint using XSDRNPM module for TSUNAMI-1D, NEWT module for TSUNAMI-2D, and KENO-V.a or KENO-VI for TSUNAMI-3D. Then, the SAMS module is applied for calculating the sensitivity coefficients for every energy-group reaction cross section. The sensitivity to the average number of neutrons per fission,ν¯, and to the fission spectrum, χ, is estimated. Also, SAMS module performs the uncertainty quantification on the keff due to the cross section uncertainties.


## 2.1.2. MCNPX-2.7e—The Differential Operator Technique

The Differential Operator Technique is applied in MCNPX-2.7e to calculate the changeΔkeff due to a perturbation in a cross section, Δσ, by means of the PERT card. It is based on the Taylor series expansion as presented in (9), where σx is the cross section perturbed, σx,0 is its nominal value, and Δσ=σx-σx,0, is the perturbation:
(9)Δkeff=keff(σx)-keff(σx,0)=(∂k∂σ|σx,0Δσ)1st+(12∂2k∂σ2|σx,0(Δσ)2)2nd+⋯=[Δkeff]PERT,1st+[Δkeff]PERT,2nd+⋯.
The first derivative means the derivative of the probability of the random walk occurring. The second one is the derivative of the tally response itself. There is an additional derivative for the changes in the fission source distribution that MCNPX-2.7e does not take into account, but methodologies about how to calculate it have been developed in [4, 7–9] showing also its effect. In [4, 10], it was pointed out that because the scattering affects the fission source spatial distribution more than capture, the sensitivity to scattering cross sections could be wrong. Other limitations are stated in the MCNPX manual [3, Section 5.2.2].For the sensitivity coefficient, only the first derivative is required. Using the first term provided by the PERT card (using the keywordMETHOD=+2), the sensitivity coefficient is calculated with (10), just dividing by keff,0 and px, which is the variation of the cross section from σx=σx,0(1+px). Because [Δkeff]PERT,1st is lineal with respect to the size of the perturbation px, so choosing one value for px is trivial:
(10)Si=1k0px[Δkeff]PERT,1st.
The statistical uncertainty is propagated to the sensitivity coefficient by (11), obtaining the relative standard deviation:
(11)var(Si)Si=var([Δkeff]PERT,1st)[Δkeff]PERT,1st2+var(k0)k02.


## 3. PWR Calculations on Exercise I-2

### 3.1. Specifications and Modeling Exercise I-2

Exercise I-2 [2, Chapter 3] proposes uncertainty propagation of the input parameters through the lattice physics to output variables for different Light Water Reactors (LWRs). Here, only the first one of three proposed test problems regarding PWR geometry is carried out. It consists in propagating the cross section uncertainties in a 2D fuel assembly model with reflective boundary conditions which will be used later as the standard model for fuel assembly cross-section generation in PWR core analysis.The specifications of the fuel assembly are given in [1, Section 3.2]. It is a 15×15 TMI fuel assembly model which contains gadolinia pins as burnable poison. Also, the specifications of the control rods are found there. The definition of the pin cell is the same as used in Exercise I-1 [1, Section 2.4]. The fuel assembly configuration is presented in Figure 1(a) for the unrodded case and in Figure 1(b) for the rodded case. The fuel pins are in blue, the gadolinia pins in light blue, the control rods in green, cladding in yellow, and water in red.Geometrical model implemented in MCNPX-2.7e for1/8 of the 15×15 TMI fuel assembly.
(a)
Unrodded
(b)
RoddedThe criticality calculations for the15×15 TMI fuel assembly are carried out at Hot Zero Power (HZP) and at Hot Full Power (HFP) conditions for both cases, unrodded and rodded, together with uncertainty quantification on  keff due to the cross section uncertainties.


#### 3.1.1. SCALE-6.1

For the criticality calculations with sensitivity analysis and uncertainty quantification of SCALE-6.1, the TSUNAMI-2D sequence is chosen. There, the NEWT module [2, Section F21] is used for solving the forward and adjoint transport problems. NEWT is an algorithm for two-dimensional analysis in nonorthogonal geometries combining a step characteristics approach plus discrete ordinates. The 238 energy group cross section library of the ENDF/B-VII.0 [2, Section M4] is used with NEWT (using the keyword v7-238). Then, the SAMS module [2, Section F22] is invoked to calculate the sensitivity coefficients in 238 energy groups. Because SCALE-6.1 provides a set of covariance matrices to perform the uncertainty quantification, “SCALE Nuclear Data Covariance Library” [2, Section M19] that is given in 44 energy group [2, Table M4.2.1], the sensitivity profiles are collapsed from 238 to 44 energy group structure inside SAMS, but the latter profile is not supplied to the user. But collapsing to 44 energy groups can be done with the VIBE tool [11].


#### 3.1.2. MCNPX-2.7e

For the criticality calculations of MCNPX-2.7e, the KCODE card [3, Section 5.2.2] is applied and activated only for neutrons. The ENDF/B-VII.0 cross section library is used, taking the processed cross section files for MCNPX from [12]. Because the cross section library used for 107Ag in [12] comes from JEFF-3.1.1, the 107Ag ENDF/B-VII.0 cross section file is processed and included in the set of cross section libraries used by MCNPX-2.7e to substitute the previous one.After preparing the inputs for the criticality calculations, PERT cards are added in order to calculate the sensitivity coefficients ofkeff to the reaction cross sections given in Table 1.Table 1
Reaction cross sections perturbed using PERT cards in the MCNPX-2.7e calculations.
Reaction
Isotopes
σ
(
n
,
n
)
152,154,155,156,157,158,160
Gd, 238U, 107,109Ag, 113Cd, 115In
σ
(
n
,
n
′
)
238
U
σ
(
n
,
f
)
235,238
U
σ
(
n
,
γ
)
152,154,155,156,157,158,160
Gd, 238U, 107,109Ag, 113Cd, 115InTo calculate the same sensitivity coefficients as SCALE-6.1, the keyword that sets to which reaction cross section is assigned the perturbation on thekeff, RXN, is specified taking into account the equivalences between RXN and the identifiers of the SCALE reaction cross sections presented in [4, Table V]. In order to compare the sensitivity profiles provided by SCALE-6.1 and MCNPX-2.7e, the PERT cards should use the same 44 energy group structure used by SCALE-6.1. px is set to 100 to calculate the sensitivity coefficients by (10).Once the sensitivity coefficients are calculated, the “SCALE Nuclear Data Covariance Library” is processed by the VIEWCVX code (provided as a module of ERRORJ code [13]) to obtain the relative covariance matrix of the reactions given in Table 1. Then, the “sandwich formula” given in (6) is calculated, obtaining the uncertainty on keff.


### 3.2. Comparison ofkeff Values

Thekeff values are presented in Table 2. MCNPX-2.7e provides larger values than SCALE-6.1 for all cases, but the difference between codes is reduced when going from unrodded to rodded case. These differences appear due to the use of multigroup cross-sections or continuous energy cross-sections, the different methodologies used in each code to solve the neutron transport equations, and the different processing codes used for preparing the cross section data files: NJOY v99.364 for MCNPX-2.7e and AMPX for SCALE-6.1.Table 2
keff  values for the different states of the 15×15 TMI fuel assembly obtained by MCNPX-2.7e and SCALE-6.1, taking the SCALE results as the reference for the difference in pcm.
Unrodded
Rodded
HZP
HFP
HZP
HFP
MCNPX-2.7.e
1.41768
±
16 pcm
1.40441
±
16 pcm
1.07386
±
22 pcm
1.06146
±
25 pcm
SCALE-6.1
1.41227
1.39802
1.07160
1.05834
Difference
383 pcm
457 pcm
211 pcm
295 pcm


### 3.3.   Comparison of thekeff Uncertainty

The uncertainty results obtained by SCALE-6.1 and MCNPX-2.7e are presented in the tables: for the unrodded fuel assembly at HZP (Table3(a)), and at HFP (Table 3(b)); for the rodded case at HZP (Table 4(a)), and at HFP (Table 4(b)). These tables present in the first two columns the reaction pair of the variance-covariance matrix that contributes to the keff uncertainty. The contribution is given in rel.std.dev.(%) in the third and fourth columns. Only the contributions with ||rel.std.dev.(%)||>0.001 are presented, sorted in descending order. The square root of the sum of the square rel.std.dev. provides the total uncertainty in keff as rel.std.dev.(%). The ratio between MCNPX-2.7e and SCALE-6.1 values is given for each contribution and for the total values. Two more rows are added to show the importance of the 235U  ν¯ and 235U  χ uncertainty contributions to keff, calculated by SCALE-6.1. The last row shows the total uncertainty given by SCALE-6.1, taking into account all the possible reactions for which there is uncertainty information, and the ratio to the total uncertainty provided by MCNPX-2.7e.Table 3
Comparison of  the uncertainty contribution tokeff of each variance-covariance matrix of the reaction pair presented in the first two columns given in rel.std.dev.(%) calculated by SCALE-6.1 and MCNPX-2.7e for unrodded case.
(a)
HZP
Reaction
Reaction
MCNP-2.7e
SCALE-6.1
Ratio
238U
(
n
,
γ
)
238U
(
n
,
γ
)
0.23145
0.24582
0.942
235U
(
n
,
γ
)
235U
(
n
,
γ
)
0.19674
0.19970
0.985
235U
(
n
,
f
)
235U
(
n
,
γ
)
0.10733
0.10895
0.985
238U
(
n
,
n
′
)
238U
(
n
,
n
′
)
0.09078
0.09455
0.960
235U
(
n
,
f
)
235U
(
n
,
f
)
0.08368
0.08500
0.985
238U
(
n
,
f
)
238U
(
n
,
f
)
0.01363
0.01366
0.998
238
U
(n, n)
238
U
(n,
γ
)
0.01215
−0.00474
−2.566
238
U
(n, n)
238
U
(n, n)
0.01185
0.02265
0.523
157Gd
(
n
,
γ
)
157Gd
(
n
,
γ
)
0.00645
0.00683
0.945
155Gd
(
n
,
γ
)
155Gd
(
n
,
γ
)
0.00476
0.00488
0.977
238U
(
n
,
f
)
238U
(
n
,
γ
)
0.00246
0.00257
0.957
156Gd
(
n
,
γ
)
156Gd
(
n
,
γ
)
0.00129
0.00128
1.006
238
U
(n, n)
238
U
(n, f)
−0.00116
0.00065
−1.787
238U
(
n
,
n
)
238U
(
n
,
n
′
)
−0.00476
−0.00504
0.944
Total
0.34577
0.35940
0.962
235U
ν
¯
235U
ν
¯
—
0.26937
235U
χ
235U
χ
—
0.08531
Total SCALE
0.46441
0.774
(b)
HFP
Reaction
Reaction
MCNP-2.7e
SCALE-6.1
Ratio
238U
(
n
,
γ
)
238U
(
n
,
γ
)
0.23793
0.25293
0.941
235U
(
n
,
γ
)
235U
(
n
,
γ
)
0.19738
0.20025
0.986
235U
(
n
,
f
)
235U
(
n
,
γ
)
0.10697
0.10860
0.985
238U
(
n
,
n
′
)
238U
(
n
,
n
′
)
0.09808
0.09962
0.985
235U
(
n
,
f
)
235U
(
n
,
f
)
0.08333
0.08470
0.984
238
U
(n, n)
238
U
(n,
n

′
)
0.02302
−0.00562
−4.093
238U
(
n
,
n
)
238U
(
n
,
γ
)
0.01634
0.00943
1.732
238U
(
n
,
f
)
238U
(
n
,
f
)
0.01407
0.01412
0.997
238
U
(n, n)
238
U
(n, n)
0.01358
0.02520
0.539
157Gd
(
n
,
γ
)
157Gd
(
n
,
γ
)
0.00653
0.00693
0.943
155Gd
(
n
,
γ
)
155Gd
(
n
,
γ
)
0.00480
0.00496
0.967
238U
(
n
,
f
)
238U
(
n
,
γ
)
0.00259
0.00272
0.953
156Gd
(
n
,
γ
)
156Gd
(
n
,
γ
)
0.00135
0.00132
1.023
238
U
(n, n)
238
U
(n, f)
−0.00137
0.00070
−1.959
Total
0.35331
0.36607
0.965
235U
ν
¯
235U
ν
¯
—
0.26834
235U
χ
235U
χ
—
0.08823
Total SCALE
0.46984
0.779Table 4
Comparison of  the uncertainty contribution tokeff of each variance-covariance matrix of the reaction pair presented in the first two columns given in rel.std.dev.(%) calculated by SCALE-6.1 and MCNPX-2.7e for rodded case.
(a)
HZP
Reaction
Reaction
MCNP-2.7e
SCALE-6.1
Ratio
238U
(
n
,
γ
)
238U
(
n
,
γ
)
0.21670
0.23086
0.939
235U
(
n
,
γ
)
235U
(
n
,
γ
)
0.17403
0.17602
0.989
238U
(
n
,
n
′
)
238U
(
n
,
n
′
)
0.16286
0.16089
1.012
235U
(
n
,
f
)
235U
(
n
,
γ
)
0.11027
0.11189
0.985
235U
(
n
,
f
)
235U
(
n
,
f
)
0.10708
0.10936
0.979
109Ag
(
n
,
γ
)
109Ag
(
n
,
γ
)
0.07569
0.07831
0.967
107Ag
(
n
,
γ
)
107Ag
(
n
,
γ
)
0.06786
0.06965
0.974
115In
(
n
,
γ
)
115In
(
n
,
γ
)
0.03640
0.03627
1.003
238U
(
n
,
f
)
238U
(
n
,
f
)
0.02242
0.02237
1.002
238
U
(n, n)
238
U
(n, γ)
0.01812
0.01185
1.529
238
U
(n,n)
238
U
(n,n)
0.01357
0.02340
0.580
157Gd
(
n
,
γ
)
157Gd
(
n
,
γ
)
0.00747
0.00817
0.915
113Cd
(
n
,
γ
)
113Cd
(
n
,
γ
)
0.00744
0.00738
1.008
155Gd
(
n
,
γ
)
155Gd
(
n
,
γ
)
0.00556
0.00584
0.952
107
A
g
(n, n)
107
A
g
(n, n)
0.00347
0.00042
8.193
238U
(
n
,
f
)
238U
(
n
,
γ
)
0.00322
0.00337
0.955
115
I
n
(n, n)
115
I
n
(n, n)
0.00310
0.00164
1.888
156Gd
(
n
,
γ
)
156Gd
(
n
,
γ
)
0.00162
0.00154
1.050
109
A
g
(n, n)
109
A
g
(n, n)
0.00145
0.00031
4.695
238
U
(n, n)
238
U
(n,  f)
−0.00139
0.00088
−1.579
238
U
(n, n)
238
U
(n, n

′
)
−0.03164
−0.00757
4.178
Total
0.37315
0.38518
0.969
235U
ν
¯
235U
ν
¯
—
0.25594
235U
χ
235U
χ
—
0.13360
Total SCALE
0.49199
0.783
(b)
HFP
Reaction
Reaction
MCNP-2.7e
SCALE-6.1
Ratio
238U
(
n
,
γ
)
238U
(
n
,
γ
)
0.22281
0.23783
0.937
235U
(
n
,
γ
)
235U
(
n
,
γ
)
0.17426
0.17639
0.988
238U
(
n
,
n
′
)
238U
(
n
,
n
′
)
0.16419
0.16909
0.971
235U
(
n
,
f
)
235U
(
n
,
γ
)
0.10985
0.11143
0.986
235U
(
n
,
f
)
235U
(
n
,
f
)
0.10703
0.10925
0.980
109Ag
(
n
,
γ
)
109Ag
(
n
,
γ
)
0.07689
0.07969
0.965
107Ag
(
n
,
γ
)
107Ag
(
n
,
γ
)
0.06925
0.07100
0.975
115In
(
n
,
γ
)
115In
(
n
,
γ
)
0.03707
0.03690
1.005
238U
(
n
,
f
)
238U
(
n
,
f
)
0.02313
0.02310
1.001
238
U
(n, n)
238
U
(n, γ)
0.02141
0.01592
1.345
238
U
(n, n)
238
U
(n, n

′
)
0.01567
−0.00854
−1.835
238
U
(n, n)
238
U
(n, n)
0.01276
0.02589
0.493
157Gd
(
n
,
γ
)
157Gd
(
n
,
γ
)
0.00770
0.00825
0.934
113Cd
(
n
,
γ
)
113Cd
(
n
,
γ
)
0.00762
0.00748
1.019
107
A
g
(n, n)
107
A
g
(n, n)
0.00587
0.00043
13.758
155Gd
(
n
,
γ
)
155Gd
(
n
,
γ
)
0.00569
0.00592
0.961
109
A
g
(n, n)
109
A
g
(n, n)
0.00450
0.00031
14.391
238U
(
n
,
f
)
238U
(
n
,
γ
)
0.00341
0.00356
0.958
115
I
n
(n, n)
115
I
n
(n, n)
0.00262
0.00167
1.569
156Gd
(
n
,
γ
)
156Gd
(
n
,
γ
)
0.00166
0.00159
1.048
238
U
(n, n)
238
U
(n, f)
−0.00151
0.00094
−1.599
Total
0.37973
0.39380
0.964
235U
ν
¯
235U
ν
¯
—
0.25451
235U
χ
235U
χ
—
0.13788
Total SCALE
0.49966
0.788There is a good agreement between MCNPX-2.7e and SCALE-6.1 results except when the(n,n) reaction cross section is involved. There is a contribution in the HZP unrodded case in which the (n,n) reaction is involved and both codes remain in good agreement: the 238U (n,n)-238U (n,n′) covariance matrix contribution. But this result is obtained by chance, because as showed later, there are large differences in the 238U (n,n) sensitivity profile between both codes, while for the 238U (n,n′) sensitivity profile there is almost no difference. Therefore, after applying the “sandwich formula,” the contribution given by each code is almost the same due to the combination of the differences in the 238U (n,n) sensitivity profile.However, the(n,n) reactions are included in the total values given by MCNPX; they have no impact because their contribution to the total is negligible.It is necessary to implement more PERT cards in MCNPX-2.7e, because only with the ones calculated, the total uncertainty in thekeff is being underestimated around a 23%. This 23% comes mainly because MCNPX-2.7e is not able to calculate the sensitivities to 235U  ν¯ and 235U  χ yet. But this capability could be implemented as presented in [14].The differences between HZP and HFP are explained later through the comparison of the sensitivity profiles, because the variance-covariance matrices used in the “sandwich formula” do not change between cases.


### 3.4. Comparison of Sensitivity Profiles

The sensitivity profiles of238U reaction cross sections calculated by MCNPX-2.7e and SCALE-6.1 are presented in Figures 2 and 3, and for 235U reaction cross sections in Figure 4, both for unrodded and rodded cases at HZP. MCNPX values are given in red, with dashed lines for their statistical uncertainty because of the Monte Carlo approach of solving the transport problem, and SCALE ones are in blue.Figure 2
Sensitivity profiles of238U reaction cross sections calculated by MCNPX-2.7e and SCALE-6.1 for unrodded case at HZP.Figure 3
Sensitivity profiles of238U reaction cross sections calculated by MCNPX-2.7e and SCALE-6.1 for rodded case at HZP.Sensitivity profiles of235U reaction cross sections calculated by MCNPX-2.7e and SCALE-6.1 for unrodded and rodded cases at HZP.
(a)
Unrodded
(b)
RoddedThere is good agreement between all reaction cross sections except for(n,n). Special remark should be done for the large statistical uncertainty on the (n,n′) reaction while its mean value remains in good agreement with the SCALE-6.1 values. For the other isotope-reaction cross sections, the same trend is observed: good agreement for all reaction cross section except for the (n,n) reaction. These kinds of differences in (n,n) reaction were found before in [4], where the scattering reactions did not agree well with the SCALE results because of the deficiencies of the Differences Operator Technique used in the PERT card.The temperature effect on the sensitivity profiles is analysed. Only for the unrodded case there is a noticeable change in the238U (n,γ) sensitivity profile; the sensitivity coefficients of the low energy groups are increased when going from HZP to HFP due to the Doppler broadening.The effect of the control rods can be observed comparing Figures2 (unrodded) and 3 (rodded). The 238U (n,f) and (n,n′) reactions are affected: the sensitivity coefficients of the low energy groups are increased when going from unrodded to rodded case. The 157Gd (n,γ) sensitivity coefficients follow the same trend: they are increased at low energies. The effect on 235U reactions can be seen comparing Figures 4(a) and 4(b). The 235U (n,f) sensitivity coefficients are increased in the whole range while the 235U (n,γ) ones are decreased only at low energies, when changing from unrodded to rodded configuration. These changes provoke the variations observed in the contributions of these reactions to the keff uncertainty.Another source of difference, apart from the methodology used by each code to calculate the sensitivity coefficients, is that SCALE-6.1 provides the sensitivity profiles in 238 energy groups, while for MCNPX-2.7e the sensitivity profiles are calculated in 44 energy groups. Therefore, the energy group structure could lead to differences in between both codes.


### 3.5. Ranking the Importance of the Reactions

SCALE-6.1 provides a useful result: the integrated sensitivity coefficients of every reaction involved in the criticality calculations. It is the sum of the sensitivity coefficients over all energy groups of the same reaction. This value reflects how sensitive is thekeff to this reaction. Then, the reactions can be ranked as a function of their relevance to keff with their integrated sensitivity coefficient.With the sensitivity profiles provided by MCNPX-2.7e, the integrated sensitivity values can be calculated. Because comparing the integrated sensitivity coefficients is equivalent to compare sensitivity profiles, only the SCALE-6.1 results are shown in this section.Table5 presents the integrated sensitivity values for the unrodded case at HZP and HFP, while Table 6 shows the results for the rodded case at HZP and HFP. The most relevant reactions are the ones for 235U, 238U, and 1H. In the rodded case, the (n,γ) reactions of 107Ag, 109Ag, and 115In arise as important.Table 5
Largest integrated sensitivity coefficients of the reaction cross sections for the unrodded case at HZP and HFP sorted in descending order provided by SCALE-6.1.
(a)
HZP
Reaction
Int.sen.coef.
235U
ν
¯
9.4533
E
-
01
235U
(
n
,
f
)
2.7342
E
-
01
238U
(
n
,
γ
)
-
1.9175
E
-
01
1H
(
n
,
n
)
1.5701
E
-
01
235U
(
n
,
γ
)
-
1.4756
E
-
01
238U
ν
¯
5.4659
E
-
02
1H
(
n
,
γ
)
-
4.5765
E
-
02
238U
(
n
,
f
)
2.6307
E
-
02
238U
(
n
,
n
)
1.4825
E
-
02
(b)
HFP
Reaction
Int.sen.coef.
235U
ν
¯
9.4404
E
-
01
235U
(
n
,
f
)
2.7389
E
-
01
238U
(
n
,
γ
)
-
1.9659
E
-
01
1H
(
n
,
n
)
1.6365
E
-
01
235U
(
n
,
γ
)
-
1.4806
E
-
01
238U
ν
¯
5.5948
E
-
02
1H
(
n
,
γ
)
-
4.4564
E
-
02
238U
(
n
,
f
)
2.7192
E
-
02
238U
(
n
,
n
)
1.5870
E
-
02Table 6
Largest integrated sensitivity coefficients of the reaction cross sections for the rodded case at HZP and HFP sorted in descending order provided by SCALE-6.1.
(a)
Rodded at HZP
Reaction
Int.sen.coef.
235U
ν
¯
9.2720
E
-
01
235U
(
n
,
f
)
3.6753
E
-
01
1H
(
n
,
n
)
2.3743
E
-
01
238U
(
n
,
γ
)
-
1.8097
E
-
01
235U
(
n
,
γ
)
-
1.2798
E
-
01
238U
ν
¯
7.1728
E
-
02
238U
(
n
,
f
)
4.3054
E
-
02
1H
(
n
,
γ
)
-
3.3060
E
-
02
109Ag
(
n
,
γ
)
-
2.1940
E
-
02
238U
(
n
,
n
)
1.6606
E
-
02
107Ag
(
n
,
γ
)
-
1.3914
E
-
02
115In
(
n
,
γ
)
-
1.2778
E
-
02
(b)
Rodded at HFP
Reaction
Int.sen.coef.
235U
ν
¯
9.2540
E
-
01
235U
(
n
,
f
)
3.6931
E
-
01
1H
(
n
,
n
)
2.4283
E
-
01
238U
(
n
,
γ
)
-
1.8549
E
-
01
235U
(
n
,
γ
)
-
1.2804
E
-
01
238U
ν
¯
7.3537
E
-
02
238U
(
n
,
f
)
4.4462
E
-
02
1H
(
n
,
γ
)
-
3.2066
E
-
02
109Ag
(
n
,
γ
)
-
2.2270
E
-
02
238U
(
n
,
n
)
1.7615
E
-
02
107Ag
(
n
,
γ
)
-
1.4158
E
-
02
115In
(
n
,
γ
)
-
1.2976
E
-
02


## 3.1. Specifications and Modeling Exercise I-2

Exercise I-2 [2, Chapter 3] proposes uncertainty propagation of the input parameters through the lattice physics to output variables for different Light Water Reactors (LWRs). Here, only the first one of three proposed test problems regarding PWR geometry is carried out. It consists in propagating the cross section uncertainties in a 2D fuel assembly model with reflective boundary conditions which will be used later as the standard model for fuel assembly cross-section generation in PWR core analysis.The specifications of the fuel assembly are given in [1, Section 3.2]. It is a 15×15 TMI fuel assembly model which contains gadolinia pins as burnable poison. Also, the specifications of the control rods are found there. The definition of the pin cell is the same as used in Exercise I-1 [1, Section 2.4]. The fuel assembly configuration is presented in Figure 1(a) for the unrodded case and in Figure 1(b) for the rodded case. The fuel pins are in blue, the gadolinia pins in light blue, the control rods in green, cladding in yellow, and water in red.Geometrical model implemented in MCNPX-2.7e for1/8 of the 15×15 TMI fuel assembly.
(a)
Unrodded
(b)
RoddedThe criticality calculations for the15×15 TMI fuel assembly are carried out at Hot Zero Power (HZP) and at Hot Full Power (HFP) conditions for both cases, unrodded and rodded, together with uncertainty quantification on  keff due to the cross section uncertainties.


### 3.1.1. SCALE-6.1

For the criticality calculations with sensitivity analysis and uncertainty quantification of SCALE-6.1, the TSUNAMI-2D sequence is chosen. There, the NEWT module [2, Section F21] is used for solving the forward and adjoint transport problems. NEWT is an algorithm for two-dimensional analysis in nonorthogonal geometries combining a step characteristics approach plus discrete ordinates. The 238 energy group cross section library of the ENDF/B-VII.0 [2, Section M4] is used with NEWT (using the keyword v7-238). Then, the SAMS module [2, Section F22] is invoked to calculate the sensitivity coefficients in 238 energy groups. Because SCALE-6.1 provides a set of covariance matrices to perform the uncertainty quantification, “SCALE Nuclear Data Covariance Library” [2, Section M19] that is given in 44 energy group [2, Table M4.2.1], the sensitivity profiles are collapsed from 238 to 44 energy group structure inside SAMS, but the latter profile is not supplied to the user. But collapsing to 44 energy groups can be done with the VIBE tool [11].


### 3.1.2. MCNPX-2.7e

For the criticality calculations of MCNPX-2.7e, the KCODE card [3, Section 5.2.2] is applied and activated only for neutrons. The ENDF/B-VII.0 cross section library is used, taking the processed cross section files for MCNPX from [12]. Because the cross section library used for 107Ag in [12] comes from JEFF-3.1.1, the 107Ag ENDF/B-VII.0 cross section file is processed and included in the set of cross section libraries used by MCNPX-2.7e to substitute the previous one.After preparing the inputs for the criticality calculations, PERT cards are added in order to calculate the sensitivity coefficients ofkeff to the reaction cross sections given in Table 1.Table 1
Reaction cross sections perturbed using PERT cards in the MCNPX-2.7e calculations.
Reaction
Isotopes
σ
(
n
,
n
)
152,154,155,156,157,158,160
Gd, 238U, 107,109Ag, 113Cd, 115In
σ
(
n
,
n
′
)
238
U
σ
(
n
,
f
)
235,238
U
σ
(
n
,
γ
)
152,154,155,156,157,158,160
Gd, 238U, 107,109Ag, 113Cd, 115InTo calculate the same sensitivity coefficients as SCALE-6.1, the keyword that sets to which reaction cross section is assigned the perturbation on thekeff, RXN, is specified taking into account the equivalences between RXN and the identifiers of the SCALE reaction cross sections presented in [4, Table V]. In order to compare the sensitivity profiles provided by SCALE-6.1 and MCNPX-2.7e, the PERT cards should use the same 44 energy group structure used by SCALE-6.1. px is set to 100 to calculate the sensitivity coefficients by (10).Once the sensitivity coefficients are calculated, the “SCALE Nuclear Data Covariance Library” is processed by the VIEWCVX code (provided as a module of ERRORJ code [13]) to obtain the relative covariance matrix of the reactions given in Table 1. Then, the “sandwich formula” given in (6) is calculated, obtaining the uncertainty on keff.


## 3.1.1. SCALE-6.1

For the criticality calculations with sensitivity analysis and uncertainty quantification of SCALE-6.1, the TSUNAMI-2D sequence is chosen. There, the NEWT module [2, Section F21] is used for solving the forward and adjoint transport problems. NEWT is an algorithm for two-dimensional analysis in nonorthogonal geometries combining a step characteristics approach plus discrete ordinates. The 238 energy group cross section library of the ENDF/B-VII.0 [2, Section M4] is used with NEWT (using the keyword v7-238). Then, the SAMS module [2, Section F22] is invoked to calculate the sensitivity coefficients in 238 energy groups. Because SCALE-6.1 provides a set of covariance matrices to perform the uncertainty quantification, “SCALE Nuclear Data Covariance Library” [2, Section M19] that is given in 44 energy group [2, Table M4.2.1], the sensitivity profiles are collapsed from 238 to 44 energy group structure inside SAMS, but the latter profile is not supplied to the user. But collapsing to 44 energy groups can be done with the VIBE tool [11].


## 3.1.2. MCNPX-2.7e

For the criticality calculations of MCNPX-2.7e, the KCODE card [3, Section 5.2.2] is applied and activated only for neutrons. The ENDF/B-VII.0 cross section library is used, taking the processed cross section files for MCNPX from [12]. Because the cross section library used for 107Ag in [12] comes from JEFF-3.1.1, the 107Ag ENDF/B-VII.0 cross section file is processed and included in the set of cross section libraries used by MCNPX-2.7e to substitute the previous one.After preparing the inputs for the criticality calculations, PERT cards are added in order to calculate the sensitivity coefficients ofkeff to the reaction cross sections given in Table 1.Table 1
Reaction cross sections perturbed using PERT cards in the MCNPX-2.7e calculations.
Reaction
Isotopes
σ
(
n
,
n
)
152,154,155,156,157,158,160
Gd, 238U, 107,109Ag, 113Cd, 115In
σ
(
n
,
n
′
)
238
U
σ
(
n
,
f
)
235,238
U
σ
(
n
,
γ
)
152,154,155,156,157,158,160
Gd, 238U, 107,109Ag, 113Cd, 115InTo calculate the same sensitivity coefficients as SCALE-6.1, the keyword that sets to which reaction cross section is assigned the perturbation on thekeff, RXN, is specified taking into account the equivalences between RXN and the identifiers of the SCALE reaction cross sections presented in [4, Table V]. In order to compare the sensitivity profiles provided by SCALE-6.1 and MCNPX-2.7e, the PERT cards should use the same 44 energy group structure used by SCALE-6.1. px is set to 100 to calculate the sensitivity coefficients by (10).Once the sensitivity coefficients are calculated, the “SCALE Nuclear Data Covariance Library” is processed by the VIEWCVX code (provided as a module of ERRORJ code [13]) to obtain the relative covariance matrix of the reactions given in Table 1. Then, the “sandwich formula” given in (6) is calculated, obtaining the uncertainty on keff.


## 3.2. Comparison ofkeff Values

Thekeff values are presented in Table 2. MCNPX-2.7e provides larger values than SCALE-6.1 for all cases, but the difference between codes is reduced when going from unrodded to rodded case. These differences appear due to the use of multigroup cross-sections or continuous energy cross-sections, the different methodologies used in each code to solve the neutron transport equations, and the different processing codes used for preparing the cross section data files: NJOY v99.364 for MCNPX-2.7e and AMPX for SCALE-6.1.Table 2
keff  values for the different states of the 15×15 TMI fuel assembly obtained by MCNPX-2.7e and SCALE-6.1, taking the SCALE results as the reference for the difference in pcm.
Unrodded
Rodded
HZP
HFP
HZP
HFP
MCNPX-2.7.e
1.41768
±
16 pcm
1.40441
±
16 pcm
1.07386
±
22 pcm
1.06146
±
25 pcm
SCALE-6.1
1.41227
1.39802
1.07160
1.05834
Difference
383 pcm
457 pcm
211 pcm
295 pcm


## 3.3.   Comparison of thekeff Uncertainty

The uncertainty results obtained by SCALE-6.1 and MCNPX-2.7e are presented in the tables: for the unrodded fuel assembly at HZP (Table3(a)), and at HFP (Table 3(b)); for the rodded case at HZP (Table 4(a)), and at HFP (Table 4(b)). These tables present in the first two columns the reaction pair of the variance-covariance matrix that contributes to the keff uncertainty. The contribution is given in rel.std.dev.(%) in the third and fourth columns. Only the contributions with ||rel.std.dev.(%)||>0.001 are presented, sorted in descending order. The square root of the sum of the square rel.std.dev. provides the total uncertainty in keff as rel.std.dev.(%). The ratio between MCNPX-2.7e and SCALE-6.1 values is given for each contribution and for the total values. Two more rows are added to show the importance of the 235U  ν¯ and 235U  χ uncertainty contributions to keff, calculated by SCALE-6.1. The last row shows the total uncertainty given by SCALE-6.1, taking into account all the possible reactions for which there is uncertainty information, and the ratio to the total uncertainty provided by MCNPX-2.7e.Table 3
Comparison of  the uncertainty contribution tokeff of each variance-covariance matrix of the reaction pair presented in the first two columns given in rel.std.dev.(%) calculated by SCALE-6.1 and MCNPX-2.7e for unrodded case.
(a)
HZP
Reaction
Reaction
MCNP-2.7e
SCALE-6.1
Ratio
238U
(
n
,
γ
)
238U
(
n
,
γ
)
0.23145
0.24582
0.942
235U
(
n
,
γ
)
235U
(
n
,
γ
)
0.19674
0.19970
0.985
235U
(
n
,
f
)
235U
(
n
,
γ
)
0.10733
0.10895
0.985
238U
(
n
,
n
′
)
238U
(
n
,
n
′
)
0.09078
0.09455
0.960
235U
(
n
,
f
)
235U
(
n
,
f
)
0.08368
0.08500
0.985
238U
(
n
,
f
)
238U
(
n
,
f
)
0.01363
0.01366
0.998
238
U
(n, n)
238
U
(n,
γ
)
0.01215
−0.00474
−2.566
238
U
(n, n)
238
U
(n, n)
0.01185
0.02265
0.523
157Gd
(
n
,
γ
)
157Gd
(
n
,
γ
)
0.00645
0.00683
0.945
155Gd
(
n
,
γ
)
155Gd
(
n
,
γ
)
0.00476
0.00488
0.977
238U
(
n
,
f
)
238U
(
n
,
γ
)
0.00246
0.00257
0.957
156Gd
(
n
,
γ
)
156Gd
(
n
,
γ
)
0.00129
0.00128
1.006
238
U
(n, n)
238
U
(n, f)
−0.00116
0.00065
−1.787
238U
(
n
,
n
)
238U
(
n
,
n
′
)
−0.00476
−0.00504
0.944
Total
0.34577
0.35940
0.962
235U
ν
¯
235U
ν
¯
—
0.26937
235U
χ
235U
χ
—
0.08531
Total SCALE
0.46441
0.774
(b)
HFP
Reaction
Reaction
MCNP-2.7e
SCALE-6.1
Ratio
238U
(
n
,
γ
)
238U
(
n
,
γ
)
0.23793
0.25293
0.941
235U
(
n
,
γ
)
235U
(
n
,
γ
)
0.19738
0.20025
0.986
235U
(
n
,
f
)
235U
(
n
,
γ
)
0.10697
0.10860
0.985
238U
(
n
,
n
′
)
238U
(
n
,
n
′
)
0.09808
0.09962
0.985
235U
(
n
,
f
)
235U
(
n
,
f
)
0.08333
0.08470
0.984
238
U
(n, n)
238
U
(n,
n

′
)
0.02302
−0.00562
−4.093
238U
(
n
,
n
)
238U
(
n
,
γ
)
0.01634
0.00943
1.732
238U
(
n
,
f
)
238U
(
n
,
f
)
0.01407
0.01412
0.997
238
U
(n, n)
238
U
(n, n)
0.01358
0.02520
0.539
157Gd
(
n
,
γ
)
157Gd
(
n
,
γ
)
0.00653
0.00693
0.943
155Gd
(
n
,
γ
)
155Gd
(
n
,
γ
)
0.00480
0.00496
0.967
238U
(
n
,
f
)
238U
(
n
,
γ
)
0.00259
0.00272
0.953
156Gd
(
n
,
γ
)
156Gd
(
n
,
γ
)
0.00135
0.00132
1.023
238
U
(n, n)
238
U
(n, f)
−0.00137
0.00070
−1.959
Total
0.35331
0.36607
0.965
235U
ν
¯
235U
ν
¯
—
0.26834
235U
χ
235U
χ
—
0.08823
Total SCALE
0.46984
0.779Table 4
Comparison of  the uncertainty contribution tokeff of each variance-covariance matrix of the reaction pair presented in the first two columns given in rel.std.dev.(%) calculated by SCALE-6.1 and MCNPX-2.7e for rodded case.
(a)
HZP
Reaction
Reaction
MCNP-2.7e
SCALE-6.1
Ratio
238U
(
n
,
γ
)
238U
(
n
,
γ
)
0.21670
0.23086
0.939
235U
(
n
,
γ
)
235U
(
n
,
γ
)
0.17403
0.17602
0.989
238U
(
n
,
n
′
)
238U
(
n
,
n
′
)
0.16286
0.16089
1.012
235U
(
n
,
f
)
235U
(
n
,
γ
)
0.11027
0.11189
0.985
235U
(
n
,
f
)
235U
(
n
,
f
)
0.10708
0.10936
0.979
109Ag
(
n
,
γ
)
109Ag
(
n
,
γ
)
0.07569
0.07831
0.967
107Ag
(
n
,
γ
)
107Ag
(
n
,
γ
)
0.06786
0.06965
0.974
115In
(
n
,
γ
)
115In
(
n
,
γ
)
0.03640
0.03627
1.003
238U
(
n
,
f
)
238U
(
n
,
f
)
0.02242
0.02237
1.002
238
U
(n, n)
238
U
(n, γ)
0.01812
0.01185
1.529
238
U
(n,n)
238
U
(n,n)
0.01357
0.02340
0.580
157Gd
(
n
,
γ
)
157Gd
(
n
,
γ
)
0.00747
0.00817
0.915
113Cd
(
n
,
γ
)
113Cd
(
n
,
γ
)
0.00744
0.00738
1.008
155Gd
(
n
,
γ
)
155Gd
(
n
,
γ
)
0.00556
0.00584
0.952
107
A
g
(n, n)
107
A
g
(n, n)
0.00347
0.00042
8.193
238U
(
n
,
f
)
238U
(
n
,
γ
)
0.00322
0.00337
0.955
115
I
n
(n, n)
115
I
n
(n, n)
0.00310
0.00164
1.888
156Gd
(
n
,
γ
)
156Gd
(
n
,
γ
)
0.00162
0.00154
1.050
109
A
g
(n, n)
109
A
g
(n, n)
0.00145
0.00031
4.695
238
U
(n, n)
238
U
(n,  f)
−0.00139
0.00088
−1.579
238
U
(n, n)
238
U
(n, n

′
)
−0.03164
−0.00757
4.178
Total
0.37315
0.38518
0.969
235U
ν
¯
235U
ν
¯
—
0.25594
235U
χ
235U
χ
—
0.13360
Total SCALE
0.49199
0.783
(b)
HFP
Reaction
Reaction
MCNP-2.7e
SCALE-6.1
Ratio
238U
(
n
,
γ
)
238U
(
n
,
γ
)
0.22281
0.23783
0.937
235U
(
n
,
γ
)
235U
(
n
,
γ
)
0.17426
0.17639
0.988
238U
(
n
,
n
′
)
238U
(
n
,
n
′
)
0.16419
0.16909
0.971
235U
(
n
,
f
)
235U
(
n
,
γ
)
0.10985
0.11143
0.986
235U
(
n
,
f
)
235U
(
n
,
f
)
0.10703
0.10925
0.980
109Ag
(
n
,
γ
)
109Ag
(
n
,
γ
)
0.07689
0.07969
0.965
107Ag
(
n
,
γ
)
107Ag
(
n
,
γ
)
0.06925
0.07100
0.975
115In
(
n
,
γ
)
115In
(
n
,
γ
)
0.03707
0.03690
1.005
238U
(
n
,
f
)
238U
(
n
,
f
)
0.02313
0.02310
1.001
238
U
(n, n)
238
U
(n, γ)
0.02141
0.01592
1.345
238
U
(n, n)
238
U
(n, n

′
)
0.01567
−0.00854
−1.835
238
U
(n, n)
238
U
(n, n)
0.01276
0.02589
0.493
157Gd
(
n
,
γ
)
157Gd
(
n
,
γ
)
0.00770
0.00825
0.934
113Cd
(
n
,
γ
)
113Cd
(
n
,
γ
)
0.00762
0.00748
1.019
107
A
g
(n, n)
107
A
g
(n, n)
0.00587
0.00043
13.758
155Gd
(
n
,
γ
)
155Gd
(
n
,
γ
)
0.00569
0.00592
0.961
109
A
g
(n, n)
109
A
g
(n, n)
0.00450
0.00031
14.391
238U
(
n
,
f
)
238U
(
n
,
γ
)
0.00341
0.00356
0.958
115
I
n
(n, n)
115
I
n
(n, n)
0.00262
0.00167
1.569
156Gd
(
n
,
γ
)
156Gd
(
n
,
γ
)
0.00166
0.00159
1.048
238
U
(n, n)
238
U
(n, f)
−0.00151
0.00094
−1.599
Total
0.37973
0.39380
0.964
235U
ν
¯
235U
ν
¯
—
0.25451
235U
χ
235U
χ
—
0.13788
Total SCALE
0.49966
0.788There is a good agreement between MCNPX-2.7e and SCALE-6.1 results except when the(n,n) reaction cross section is involved. There is a contribution in the HZP unrodded case in which the (n,n) reaction is involved and both codes remain in good agreement: the 238U (n,n)-238U (n,n′) covariance matrix contribution. But this result is obtained by chance, because as showed later, there are large differences in the 238U (n,n) sensitivity profile between both codes, while for the 238U (n,n′) sensitivity profile there is almost no difference. Therefore, after applying the “sandwich formula,” the contribution given by each code is almost the same due to the combination of the differences in the 238U (n,n) sensitivity profile.However, the(n,n) reactions are included in the total values given by MCNPX; they have no impact because their contribution to the total is negligible.It is necessary to implement more PERT cards in MCNPX-2.7e, because only with the ones calculated, the total uncertainty in thekeff is being underestimated around a 23%. This 23% comes mainly because MCNPX-2.7e is not able to calculate the sensitivities to 235U  ν¯ and 235U  χ yet. But this capability could be implemented as presented in [14].The differences between HZP and HFP are explained later through the comparison of the sensitivity profiles, because the variance-covariance matrices used in the “sandwich formula” do not change between cases.


## 3.4. Comparison of Sensitivity Profiles

The sensitivity profiles of238U reaction cross sections calculated by MCNPX-2.7e and SCALE-6.1 are presented in Figures 2 and 3, and for 235U reaction cross sections in Figure 4, both for unrodded and rodded cases at HZP. MCNPX values are given in red, with dashed lines for their statistical uncertainty because of the Monte Carlo approach of solving the transport problem, and SCALE ones are in blue.Figure 2
Sensitivity profiles of238U reaction cross sections calculated by MCNPX-2.7e and SCALE-6.1 for unrodded case at HZP.Figure 3
Sensitivity profiles of238U reaction cross sections calculated by MCNPX-2.7e and SCALE-6.1 for rodded case at HZP.Sensitivity profiles of235U reaction cross sections calculated by MCNPX-2.7e and SCALE-6.1 for unrodded and rodded cases at HZP.
(a)
Unrodded
(b)
RoddedThere is good agreement between all reaction cross sections except for(n,n). Special remark should be done for the large statistical uncertainty on the (n,n′) reaction while its mean value remains in good agreement with the SCALE-6.1 values. For the other isotope-reaction cross sections, the same trend is observed: good agreement for all reaction cross section except for the (n,n) reaction. These kinds of differences in (n,n) reaction were found before in [4], where the scattering reactions did not agree well with the SCALE results because of the deficiencies of the Differences Operator Technique used in the PERT card.The temperature effect on the sensitivity profiles is analysed. Only for the unrodded case there is a noticeable change in the238U (n,γ) sensitivity profile; the sensitivity coefficients of the low energy groups are increased when going from HZP to HFP due to the Doppler broadening.The effect of the control rods can be observed comparing Figures2 (unrodded) and 3 (rodded). The 238U (n,f) and (n,n′) reactions are affected: the sensitivity coefficients of the low energy groups are increased when going from unrodded to rodded case. The 157Gd (n,γ) sensitivity coefficients follow the same trend: they are increased at low energies. The effect on 235U reactions can be seen comparing Figures 4(a) and 4(b). The 235U (n,f) sensitivity coefficients are increased in the whole range while the 235U (n,γ) ones are decreased only at low energies, when changing from unrodded to rodded configuration. These changes provoke the variations observed in the contributions of these reactions to the keff uncertainty.Another source of difference, apart from the methodology used by each code to calculate the sensitivity coefficients, is that SCALE-6.1 provides the sensitivity profiles in 238 energy groups, while for MCNPX-2.7e the sensitivity profiles are calculated in 44 energy groups. Therefore, the energy group structure could lead to differences in between both codes.


## 3.5. Ranking the Importance of the Reactions

SCALE-6.1 provides a useful result: the integrated sensitivity coefficients of every reaction involved in the criticality calculations. It is the sum of the sensitivity coefficients over all energy groups of the same reaction. This value reflects how sensitive is thekeff to this reaction. Then, the reactions can be ranked as a function of their relevance to keff with their integrated sensitivity coefficient.With the sensitivity profiles provided by MCNPX-2.7e, the integrated sensitivity values can be calculated. Because comparing the integrated sensitivity coefficients is equivalent to compare sensitivity profiles, only the SCALE-6.1 results are shown in this section.Table5 presents the integrated sensitivity values for the unrodded case at HZP and HFP, while Table 6 shows the results for the rodded case at HZP and HFP. The most relevant reactions are the ones for 235U, 238U, and 1H. In the rodded case, the (n,γ) reactions of 107Ag, 109Ag, and 115In arise as important.Table 5
Largest integrated sensitivity coefficients of the reaction cross sections for the unrodded case at HZP and HFP sorted in descending order provided by SCALE-6.1.
(a)
HZP
Reaction
Int.sen.coef.
235U
ν
¯
9.4533
E
-
01
235U
(
n
,
f
)
2.7342
E
-
01
238U
(
n
,
γ
)
-
1.9175
E
-
01
1H
(
n
,
n
)
1.5701
E
-
01
235U
(
n
,
γ
)
-
1.4756
E
-
01
238U
ν
¯
5.4659
E
-
02
1H
(
n
,
γ
)
-
4.5765
E
-
02
238U
(
n
,
f
)
2.6307
E
-
02
238U
(
n
,
n
)
1.4825
E
-
02
(b)
HFP
Reaction
Int.sen.coef.
235U
ν
¯
9.4404
E
-
01
235U
(
n
,
f
)
2.7389
E
-
01
238U
(
n
,
γ
)
-
1.9659
E
-
01
1H
(
n
,
n
)
1.6365
E
-
01
235U
(
n
,
γ
)
-
1.4806
E
-
01
238U
ν
¯
5.5948
E
-
02
1H
(
n
,
γ
)
-
4.4564
E
-
02
238U
(
n
,
f
)
2.7192
E
-
02
238U
(
n
,
n
)
1.5870
E
-
02Table 6
Largest integrated sensitivity coefficients of the reaction cross sections for the rodded case at HZP and HFP sorted in descending order provided by SCALE-6.1.
(a)
Rodded at HZP
Reaction
Int.sen.coef.
235U
ν
¯
9.2720
E
-
01
235U
(
n
,
f
)
3.6753
E
-
01
1H
(
n
,
n
)
2.3743
E
-
01
238U
(
n
,
γ
)
-
1.8097
E
-
01
235U
(
n
,
γ
)
-
1.2798
E
-
01
238U
ν
¯
7.1728
E
-
02
238U
(
n
,
f
)
4.3054
E
-
02
1H
(
n
,
γ
)
-
3.3060
E
-
02
109Ag
(
n
,
γ
)
-
2.1940
E
-
02
238U
(
n
,
n
)
1.6606
E
-
02
107Ag
(
n
,
γ
)
-
1.3914
E
-
02
115In
(
n
,
γ
)
-
1.2778
E
-
02
(b)
Rodded at HFP
Reaction
Int.sen.coef.
235U
ν
¯
9.2540
E
-
01
235U
(
n
,
f
)
3.6931
E
-
01
1H
(
n
,
n
)
2.4283
E
-
01
238U
(
n
,
γ
)
-
1.8549
E
-
01
235U
(
n
,
γ
)
-
1.2804
E
-
01
238U
ν
¯
7.3537
E
-
02
238U
(
n
,
f
)
4.4462
E
-
02
1H
(
n
,
γ
)
-
3.2066
E
-
02
109Ag
(
n
,
γ
)
-
2.2270
E
-
02
238U
(
n
,
n
)
1.7615
E
-
02
107Ag
(
n
,
γ
)
-
1.4158
E
-
02
115In
(
n
,
γ
)
-
1.2976
E
-
02


## 4.  Conclusions

The propagation of cross section uncertainties in criticality calculations for a15×15 TMI fuel assembly in two different configurations, unrodded and rodded, at two different states, HZP and HFP, has been performed using two different approaches carried out by each of the codes used: (i) MCNPX-2.7e—Differential Operator Technique and (ii) SCALE-6.1—Adjoint-Weighted Technique.The two approaches have been presented and compared through this exercise. Thekeff and its uncertainty due to cross-section uncertainties have been calculated, distinguishing between its different contributors. The most important ones for this exercise are 235U-ν¯, χ, (n,f),(n,γ),238U-(n,n′),(n,f),(n,γ), and the covariance 235U (n,f)-235U (n,γ). Additionally, the sensitivity coefficients and the integrated sensitivity coefficients have been calculated and compared.In general, both approaches/codes are in good agreement, in spite of the differences in thekeff values. There are only two exceptions observed: one is for the contribution of the (n,n) reaction cross sections which are not properly calculated by MCNPX-2.7e because the deficiencies of the PERT card to calculate the (n,n) sensitivity coefficients, and the another is for the contribution of 238U(n,n)-(n,n′) covariance whose agreement between codes comes by chance because of the combination of the differences in the (n,n) sensitivity profile. The good agreements between codes are corroborated by the comparison of the sensitivity profiles given by each code.The inability of calculating the contribution due toν¯ and χ values makes MCNPX-2.7e to underestimate the uncertainty in keff at least in a 23% compared with SCALE-6.1. It is important to remark that the most relevant reactions are not always the main source of the uncertainty, such as the case of 1H which  is in the top list of reactions for which keff is sensitive to, but it is not in the list of the most important contributions to the keff uncertainty because of the low uncertainty of its reactions.Thus, MCNPX-2.7e can deal with the uncertainty quantification problem as SCALE-6.1 does, but improvements should be done in the PERT card capabilities such as the proper calculation of the sensitivity coefficients of(n,n) reaction cross sections, the inclusion of the perturbation due to the change in the fission source distribution, and the implementation of the perturbation of the ν¯ and χ values.



---

*Source: 380284-2013-03-27.xml*