Research Article  Open Access
Propagation of CrossSection Uncertainties in Criticality Calculations in the Framework of UAMPhase I Using MCNPX2.7e and SCALE6.1
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 crosssection uncertainties. Two different codes are employed to calculate the sensitivity coefficients of to cross sections in criticality calculations: MCNPX2.7e and SCALE6.1. The former uses the Differential Operator Technique and the latter uses the AdjointWeighted Technique. In this paper, the main results for exercise I2 “Lattice Physics” are presented for the criticality calculations of PWR. These criticality calculations are done for a TMI fuel assembly at four different states: HZPUnrodded, HZPRodded, HFPUnrodded, and HFPRodded. The results of the two different codes above are presented and compared. The comparison proves a good agreement between SCALE6.1 and MCNPX2.7e in 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 uncertainty.
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 I1. Cell Physics focused on the derivation of the multigroup microscopic cross section libraries and their uncertainties. (ii)Exercise I2. Lattice Physics focused on the derivation of the fewgroup macroscopic cross section libraries and their uncertainties. (iii)Exercise I3. Core Physics focused on the core steady state standalone 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 AdjointWeighted Technique by SCALE6.1 [2] and the Differential Operator Technique by MCNPX2.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 SCALE6.1 and MCNPX2.7e. Afterwards, in the framework of Exercise I2, both codes are used to perform the uncertainty quantification on the of a 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.
Being the calculated response function of a system which depends on parameters, can be approximated by a linear function of using the Taylor series approximation around a nominal value , as presented in (1), where :
Then, taking the parameters of the system as random variables, 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”: where (i) is the vector of the sensitivity coefficients. (ii) is the covariance matrix of the system parameters defined as Because usually the sensitivity coefficients are calculated as relative values, Equation (3) can be rewritten as (6), where the value obtained is the relative standard deviation (rel.std.dev.()) and is the relative covariance matrix: Both values, the standard deviation provided by (3) and the rel.std.dev.() 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, or , should be calculated and the covariance matrix () or the relative covariance matrix () 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)SCALE6.1 uses the AdjointWeighted Technique. (ii)MCNPX2.7e uses the Differential Operator Technique.
2.1.1. SCALE6.1—The AdjointWeighted Technique
The AdjointWeighted Technique is used in SCALE6.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 The explicit term, , is calculated using the AdjointWeighted Technique (called adjointbased perturbation in the SAMS manual), based on the perturbation of the given in (8). There, is the neutron flux of the problem defined by , where is the operator that represents all of the transport equations except for the fission term, is the operator that represents the fission term, and represents the eigenvalues where the largest one is . is the adjoint neutron flux of the adjoint problem defined as where and are the adjoint operators corresponding to and . and are the perturbed operators due to a perturbation in the cross sections:
The implicit term, , is calculated in order to take into account the effect on of perturbing one cross section that affects the resonanceshielded values of other cross sections. This term appears because the transport calculations should be done using groupwise cross sections which require an initial selfshielding 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 TSUNAMI1D, NEWT module for TSUNAMI2D, and KENOV.a or KENOVI for TSUNAMI3D. Then, the SAMS module is applied for calculating the sensitivity coefficients for every energygroup 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 due to the cross section uncertainties.
2.1.2. MCNPX2.7e—The Differential Operator Technique
The Differential Operator Technique is applied in MCNPX2.7e to calculate the change 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 is the cross section perturbed, is its nominal value, and , is the perturbation: 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 MCNPX2.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 ].
For the sensitivity coefficient, only the first derivative is required. Using the first term provided by the PERT card (using the keyword METHOD=+2), the sensitivity coefficient is calculated with (10), just dividing by and , which is the variation of the cross section from . Because is lineal with respect to the size of the perturbation , so choosing one value for is trivial: The statistical uncertainty is propagated to the sensitivity coefficient by (11), obtaining the relative standard deviation:
3. PWR Calculations on Exercise I2
3.1. Specifications and Modeling Exercise I2
Exercise I2 [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 crosssection generation in PWR core analysis.
The specifications of the fuel assembly are given in [1, Section 3.2]. It is a 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 I1 [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.
(a) Unrodded
(b) Rodded
The criticality calculations for the 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 due to the cross section uncertainties.
3.1.1. SCALE6.1
For the criticality calculations with sensitivity analysis and uncertainty quantification of SCALE6.1, the TSUNAMI2D 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 twodimensional analysis in nonorthogonal geometries combining a step characteristics approach plus discrete ordinates. The 238 energy group cross section library of the ENDF/BVII.0 [2, Section M4] is used with NEWT (using the keyword v7238). Then, the SAMS module [2, Section F22] is invoked to calculate the sensitivity coefficients in 238 energy groups. Because SCALE6.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. MCNPX2.7e
For the criticality calculations of MCNPX2.7e, the KCODE card [3, Section ] is applied and activated only for neutrons. The ENDF/BVII.0 cross section library is used, taking the processed cross section files for MCNPX from [12]. Because the cross section library used for ^{107}Ag in [12] comes from JEFF3.1.1, the ^{107}Ag ENDF/BVII.0 cross section file is processed and included in the set of cross section libraries used by MCNPX2.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 of to the reaction cross sections given in Table 1.

To calculate the same sensitivity coefficients as SCALE6.1, the keyword that sets to which reaction cross section is assigned the perturbation on the , 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 SCALE6.1 and MCNPX2.7e, the PERT cards should use the same 44 energy group structure used by SCALE6.1. 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 .
3.2. Comparison of Values
The values are presented in Table 2. MCNPX2.7e provides larger values than SCALE6.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 crosssections or continuous energy crosssections, 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 MCNPX2.7e and AMPX for SCALE6.1.

3.3. Comparison of the Uncertainty
The uncertainty results obtained by SCALE6.1 and MCNPX2.7e are presented in the tables: for the unrodded fuel assembly at HZP (Table 3(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 variancecovariance matrix that contributes to the uncertainty. The contribution is given in rel.std.dev.(%) in the third and fourth columns. Only the contributions with rel.std.dev.(%) are presented, sorted in descending order. The square root of the sum of the square rel.std.dev. provides the total uncertainty in as rel.std.dev.(%). The ratio between MCNPX2.7e and SCALE6.1 values is given for each contribution and for the total values. Two more rows are added to show the importance of the ^{235}U and ^{235}U uncertainty contributions to , calculated by SCALE6.1. The last row shows the total uncertainty given by SCALE6.1, taking into account all the possible reactions for which there is uncertainty information, and the ratio to the total uncertainty provided by MCNPX2.7e.
(a) HZP  
 
(b) HFP  

(a) HZP  
 
(b) HFP  

There is a good agreement between MCNPX2.7e and SCALE6.1 results except when the reaction cross section is involved. There is a contribution in the HZP unrodded case in which the reaction is involved and both codes remain in good agreement: the ^{238}U ^{238}U covariance matrix contribution. But this result is obtained by chance, because as showed later, there are large differences in the ^{238}U sensitivity profile between both codes, while for the ^{238}U 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 ^{238}U sensitivity profile.
However, the 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 MCNPX2.7e, because only with the ones calculated, the total uncertainty in the is being underestimated around a 23%. This 23% comes mainly because MCNPX2.7e is not able to calculate the sensitivities to ^{235}U and ^{235}U 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 variancecovariance matrices used in the “sandwich formula” do not change between cases.
3.4. Comparison of Sensitivity Profiles
The sensitivity profiles of ^{238}U reaction cross sections calculated by MCNPX2.7e and SCALE6.1 are presented in Figures 2 and 3, and for ^{235}U 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.
(a) Unrodded
(b) Rodded
There is good agreement between all reaction cross sections except for . Special remark should be done for the large statistical uncertainty on the reaction while its mean value remains in good agreement with the SCALE6.1 values. For the other isotopereaction cross sections, the same trend is observed: good agreement for all reaction cross section except for the reaction. These kinds of differences in 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 the ^{238}U 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 Figures 2 (unrodded) and 3 (rodded). The ^{238}U and reactions are affected: the sensitivity coefficients of the low energy groups are increased when going from unrodded to rodded case. The Gd sensitivity coefficients follow the same trend: they are increased at low energies. The effect on ^{235}U reactions can be seen comparing Figures 4(a) and 4(b). The ^{235}U sensitivity coefficients are increased in the whole range while the ^{235}U 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 uncertainty.
Another source of difference, apart from the methodology used by each code to calculate the sensitivity coefficients, is that SCALE6.1 provides the sensitivity profiles in 238 energy groups, while for MCNPX2.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
SCALE6.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 the to this reaction. Then, the reactions can be ranked as a function of their relevance to with their integrated sensitivity coefficient.
With the sensitivity profiles provided by MCNPX2.7e, the integrated sensitivity values can be calculated. Because comparing the integrated sensitivity coefficients is equivalent to compare sensitivity profiles, only the SCALE6.1 results are shown in this section.
Table 5 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 ^{235}U, ^{238}U, and ^{1}H. In the rodded case, the reactions of Ag, ^{109}Ag, and ^{115}In arise as important.
(a) HZP  
 
(b) HFP  

(a) Rodded at HZP  
 
(b) Rodded at HFP  
