# SCALE Modeling of Selected Neutronics Test Problems within the OECD UAM LWR’s Benchmark

Authors: Luigi Mercatali; Kostadin Ivanov; Victor Hugo Sanchez
Journal: Science and Technology of Nuclear Installations
 (2013)

---

## Abstract

The OECD UAM Benchmark was launched in 2005 with the objective of determining the uncertainty in the simulation of Light Water Reactors (LWRs) system calculations at all the stages of the coupled reactor physics—thermal hydraulics modeling. Within the framework of the “Neutronics Phase” of the Benchmark the solutions of some selected test cases at the cell physics and lattice physics levels are presented. The SCALE 6.1 code package has been used for the neutronics modeling of the selected exercises. Sensitivity and Uncertainty analysis (S/U) based on the generalized perturbation theory has been performed in order to assess the uncertainty of the computation of some selected reactor integral parameters due to the uncertainty in the basic nuclear data. As a general trend, it has been found that the main sources of uncertainty are the238U (n,γ) and the 239Pu nubar for the UOX- and the MOX-fuelled test cases, respectively. Moreover, the reference solutions for the test cases obtained using Monte Carlo methodologies together with a comparison between deterministic and stochastic solutions are presented.


---

## Body

## 1. Introduction

In recent years there has been an increasing demand from nuclear research, industry, safety, and regulation bodies for best estimate predictions of Light Water Reactors (LWRs) performances to be provided with their confidence bounds. In addition to the establishment of LWRs best-estimate calculations for design and safety analysis, understanding uncertainties of evaluated reactor parameters is important for introducing appropriate design margins and deciding where additional efforts should be undertaken to reduce those uncertainties. In order to address those issues, an in-depth discussion on “Uncertainty Analysis in Modeling” started to take place in 2005 within the OECD/NEA Nuclear Science Committee, which led to the creation of a dedicated Expert Group and to the launching of a Benchmark exercise, the OECD UAM (Uncertainty Analysis in Modeling) LWR Benchmark [1]. The proposed technical approach is to establish a benchmark for uncertainty analysis in best-estimate modeling and coupled multiphysics and multiscale LWR analysis, using as bases a series of well-defined problems with complete sets of input specifications and reference experimental data. The objective is to determine the uncertainty in LWR system calculations at all stages of coupled reactor physics/thermal hydraulics calculation. The UAM benchmark has been conceived to be structured in three different phases, being Phase I the “Neutronics Phase,” Phase II the “Core Phase,” and Phase III the “System Phase.” Additionally, each benchmark phase is subdivided in a number of different Exercises in order to propagate the full chain of uncertainty in the modeling across different scales (multi-scale) and physics phenomena (multi-physics). The present paper is devoted to the solutions of some selected test problems within the Exercises I-1 and I-2 of Phase I. The Exercise I-1 is entitled “Cell Physics” and is focused on derivation of the multigroup microscopic cross-section libraries. Its objective is to address the uncertainties due to the basic nuclear data as well as the impact of processing the nuclear and covariance data, selection of multi-group structure, and self-shielding treatment. Within Exercise I-1 the uncertainties in the evaluated Nuclear Data Libraries (NDLs) are propagated into multigroup microscopic cross-sections. In Exercise I-2 multi-group cross-section uncertainties are input uncertainties which are then propagated through the lattice physics calculations to few-group cross-section uncertainties.


## 2. Description of the Test Cases

Within the framework of Exercise I-1 different fuel pin-cell test problems have been defined representing both square and triangular pitches. The two types of basic geometries for the unit cells are schematized in Figure1.Figure 1
Types of geometries for the fuel pin-cell test cases within Exercise I-1.In this paper the following test cases have been considered.(a)
Two-dimensional fuel pin-cell problems representative of Boiling Water Reactor (BWR) Peach Bottom 2 (PB-2) [2], Pressurized Water Reactor (PWR) Three Mile Island 1 (TMI-1) [3], and Koyloduz-6 VVER-1000 [4]. Each pin-cell model has to be analyzed at Hot Full Power conditions (HFP) as well as at Hot Zero Power conditions (HZP). To enhance the differences between the three cases (PWR, BWR, and VVER) the HFP case of the BWR is defined to be calculated at 40% void fraction (with a corresponding moderator density (ρ) of 460.72 kg/m3) instead of 0%. Hence the PWR and BWR cases are for square pitch but with different spectra, while the VVER case is for triangular pitch.(b)
Fuel pin-cell test problems from the KRITZ-2 LEU critical experiments [5].(c)
PWR MOX (MOX 9.8% Pu) pin-cell case representative of Generation 3 PWR designs, which in the following text of this paper will be referred to as GEN-III [6].Within Exercise I-2, different stand-alone neutronics single Fuel Assembly (FA) and minicore test problems have been proposed. In this paper we will present the solutions for the following test cases.(a)
BWR PB-2 assembly model [2]: the assembly “type 2” of the initial loading of the Peach Bottom 2 nuclear power plant is chosen for this exercise.(b)
PWR TMI-1 assembly model [3].(c)
GEN-III assembly models [6]: one MOX and three UOX FAs types with different 235U enrichment and Gd content are available for this exercise.The parameter specifications as well as the operating conditions for all the test cases analyzed in the present paper are summarized in Table1. The six types of FAs considered in our analysis are shown in Figure 2.Table 1
Pin-cell data for the test cases of Exercises I-1 and Exercise I-2.
Parameter
BWR
PWR
VVER
KRITZ-2 : 1
GEN-III
FA geometry
7
×
7
15
×
15
—
—
17
×
17
FA pitch (mm)
152.4
218.11
—
—
216.1
Fuel rods per assembly
49
208
—
—
265
Number of guide tubes per FA
—
16
—
—
24
Number of instrumentation tubes per FA
—
1
—
—
—
Number of Gd pins per FA
—
4
—
—
—
Guide tube outside diameter (mm)
—
13.462
—
—
12.07
Guide tube inside diameter (mm)
—
12.649
—
—
11.27
Instrumentation tube outside diameter (mm)
—
12.522
—
—
—
Instrumentation tube inside diameter (mm)
—
11.201
—
—
—
Unit cell (mm)
18.75
14.427
12.75
14.85
12.62
Fuel pellet diameter (mm)
12.1158
9.391
7.56
10.58
8.253
Fuel pellet material
UO2
UO2
UO2
UO2
MOX
Fuel density (g/cm3)
10.42
10.283
10.4
—
—
Fuel enrichment (w/o)
2.93
4.85
3.3
1.86
9.8 (Pu)
Central void diameter (mm)
—
—
1.4
—
—
Central void material
—
—
Dry air
—
—
Cladding outside diameter (mm)
14.3002
10.928
9.1
12.25
9.487
Cladding thickness (mm)
0.9398
0.673
0.69
0.74
0.578
Cladding material
Zircaloy-2
Zircaloy-4
Zr + 1% Nb
Zircaloy-2
Zircaloy-4
Cladding density (g/cm3)
6.55
6.55
—
—
—
Gap material
He
He
He
He
He
Moderator material
H2O
H2O
H2O
H2O
H2O
Fuel temperature (K)
HZP
552.833
551
552.15
292.7
—
HFP
900
900
900
521.5
900
Cladding temperature (K)
HZP
552.833
551
552.15
—
—
HFP
600
600
600
—
610
Moderator temperature (K)
HZP
552.833
551
552.15
—
—
HFP
557
562
560
—
584
Moderator density (kg/m3)
HZP
753.978
766
767
—
—
HFP
460.72
748.4
752.5
—
—FAs for the test cases of Exercise I-2: BWR (a), PWR (b), GEN-III Type 1 (c), GEN-III Type 2 (d), GEN-III Type 3 (e), and GEN-III Type 4 (f).
(a)
(b)
(c)
(d)
(e)
(f)


## 3. Theoretical Approach and Computational Method

The basic problem of the neutronics is the solution of the integral-differential Boltzmann equation for the neutron transport, which is a linear equation requiring the treatment of seven independent variables: three in space, two in angle, one in energy of the incident neutrons, and time. As a consequence of such a complexity, one has to keep in mind that even if the accuracy in the predictions of the modern transport codes (both Monte Carlo and deterministic) is continuously improving, there will be always approximations introduced in the calculational procedure. Examples of uncertainties are the ones originated from the basic nuclear reaction data, from the geometrical description of the problem, and from the material compositions. The knowledge of the approximations used in the analysis and of the overall calculational uncertainties is therefore essential to gain confidence in the results obtained, and sensitivity analysis and uncertainty evaluation (S/U) are the main instruments for dealing with the sometimes scarce knowledge of the input parameters used in the simulation tools [8]. For sensitivity analysis, sensitivity coefficients are the key quantities that have to be evaluated. They are determined and assembled, using different methodologies, in a way that when multiplied by the variation of the corresponding input parameter, they will quantify the impact on the targeted quantities whose sensitivity is referred to. There are two main methodologies developed for sensitivity and uncertainty analysis. One is the forward (direct) calculation method based either on the numerical differentiation or on a stochastic method, and the other is the adjoint method based on the perturbation theory [9]. In general, the forward approach is preferable when there are few input parameters that can vary and many output parameters of interest. The contrary is true for the adjoint methodology, which is the one mainly adopted in rector physics, as the source of uncertainty is mainly related to the neutron cross-sections that can represent a very notable number of variables (up to several hundred thousand). Moreover, the linear property of the Boltzmann equation makes the adjoint approach even more attractive. Since all the analysis for the benchmark cases presented in this paper has been carried out using perturbation methodologies, let us briefly recall the theoretical background of these techniques.From a general point of view one can represent a generic integral reactor parameterQ (i.e., the keff, a reactivity coefficient, a reaction rate, etc.) as a function of cross-sections:
(1)Q=f(σ1,σ2,…,σJ),
where σ1,σ2,…,σJ represent cross sections by isotope, type of reaction, and energy range (or energy group in a multi-group representation).The variation ofQ due to variations of cross-sections σ can be expressed using perturbation theories to evaluate the sensitivity coefficients S as follows [8–12]:
(2)δQQ=∑jSjδσjσj,
where the sensitivity coefficients are formally given by
(3)Sj=∂Q∂σj·σjQ.For practical purposes, one can consider the sensitivity coefficient as divided into two components as follows:(4)δQQ=∑jSjδσjσj+(∂Q∂σe·σeQ)·δσeσe=I+D,
where the terms I and D are generally referred to as “indirect” and “direct” effect, respectively. The D term in (4) reflects the hypothesis of a direct dependence of the parameter Q on only the energy dependent detector cross-section σe. The I term in (4) is the response perturbation due to flux perturbations. The indirect term of (4) consists also of two components, namely, the explicit and implicit ones [13]. The explicit component of the indirect effect comes from the flux perturbation caused by perturbing any multi-group cross-section appearing explicitly in the transport equation. The implicit component of the indirect effect is associated with self-shielding perturbations; in other words perturbing the cross section of one nuclide may change the self-shielded cross section of another nuclide, which causes additional flux perturbations. As an example if one considers the hydrogen, perturbing the H elastic value has an explicit effect because the flux is perturbed due to changes in H moderation. However there is also an implicit effect because changing the H data perturbs the self-shielded 238U absorption cross section, which causes another flux perturbation.Let us now consider a ratio responseR characterized by the macroscopic cross-sections Σ1 and Σ2 as follows:
(5)R=〈Σ1Φ〉〈Σ2Φ〉,
where in (5) the brackets 〈,〉 indicate the integration over the phase space and Φ is the homogeneous flux. In this case the sensitivity coefficients are given by
(6)SjR=〈Ψ_,σjΦ_〉,
where Ψ_R* is the solution of
(7)M*Ψ_R*=dRdΦ,
where M* is the adjoint Boltzmann operator. The uncertainties associated to the cross-section can be represented in the form of a variance-covariance matrix:
(8)Dσ=[d11⋯d1J⋮⋱⋮d1J⋯dJJ],
where the elements dij represent the variances and covariances of the nuclear data. Once the sensitivity coefficients and the Dσ matrix are available, the variance (i.e., the uncertainty) of the generic integral parameter Q can be expressed as
(9)var(Q)=∑j,iJSjSidij.All the calculations presented in this paper have been performed by means of the SCALE 6.1 code system [14] using ENDF/B-VII.0 nuclear data [15]. SCALE (Standardized Computer Analysis for Licensing Evaluations) is a modular code system developed at Oak Ridge National Laboratory to perform analysis for criticality safety, reactor physics and radiation shielding applications. SCALE calculations typically use sequences that execute a predefined series of executable modules to compute particle fluxes and responses (multiplication factor, reaction rates, etc.). SCALE also includes modules for sensitivity and uncertainty analysis (S/U) of calculated responses. The S/U codes in SCALE are collectively referred to as TSUNAMI (Tools for Sensitivity and Uncertainty Analysis Methodology Implementation) [7, 16]. The techniques used in TSUNAMI to generate sensitivity information are based on the widely used adjoint-based perturbation theory approach described above. The flow diagram of the TSUNAMI calculations is shown in Figure 3.Figure 3
General flow diagram of the TSUNAMI code [7].The calculation procedure for the multi-group cross-section processing is based on a rigorous mechanism using the continuous energy solvers BONAMI and CENTRM [17] for self-shielding in the unresolved and resolved resonance regions, respectively, for appropriately weighting multi-group cross-sections using a continuous energy spectrum. The CENTRM module performs transport calculation using ENDF-based point data on an ultrafine energy grid (typically 30.000–70.000 energy points) to generate effectively continuous energy neutron flux solutions in the resonance and thermal ranges. This is used to weight the multi-group cross-sections to be utilized in the subsequent transport calculations. After the cross-sections are processed, the TSUNAMI-1D sequence performs two criticality calculations, solving the forward and adjoint forms of the Boltzmann equation, respectively, using the XSDRNPM discrete ordinate code [18]. In this step an energy discretization based on a 238-group structure is adopted. The sequence then calls the SAMS module, specifically SAMS5 [19], in order to compute the sensitivity coefficients. Once sensitivities are available, the uncertainty on the integral parameters of interest due to the uncertainty in the basic nuclear data are evaluated according to (9) using the so-called 44GROUPCOV covariance matrix [20]. The 44GROUPCOV matrix comprehends a total of 401 materials in a 44-group energy structure. The library includes evaluated covariances obtained from ENDF/B-VII, ENDF/B-VI, and JENDL3.3 for more than 50 materials. It is assumed [20] that covariances taken from one data evaluation, such as ENDF/B-VI or JENDL-3.3, can also be applied to other evaluations of the same data, such as ENDF/B-VII. All other nuclear data uncertainties have been estimated from approximations in which the uncertainty assessment is decoupled from the original evaluation procedure.


## 4. Results

Results for thekeff values and associated uncertainties related to the benchmark test cases of Exercises I-1 and I-2 are summarized in Table 2. As expected, because of the negative fuel Doppler coefficient, the reactivities computed for all the test cases at HFP conditions are consistently lower than those at HZP conditions. The total uncertainties of the keff have been evaluated to be ~0.5%–0.6% for all the test cases with the exception of the GEN-III Type 4 case within Exercise I-2 where the computed uncertainty is higher by around a factor of two because of the presence of the plutonium isotopes in the fuel (MOX fuel). In Figure 4 the five reaction cross-sections which contribute the most to the keff uncertainty for the test cases of Exercise I-1 are shown. While for the UOX-fuelled test cases the main contribution to the total uncertainty is due to the 238U (n,γ) followed by the 235U nubar (average number of neutrons per fission reaction—ν) and 235U (n,γ); for the MOX-fuelled test case considered (GEN-III) the predominant component to the total uncertainty comes from the 239Pu nubar followed by the 238U (n,n′) and 239Pu (n,fission). By definition the reason for these main contributions to the uncertainty can be due to the highest sensitivities associated to such reactions, to the highest value of the associated covariances, or to a combination of both. As an example, in the case of the 238U (n,γ), on one hand, the keff is quite sensitive to its value (especially in the unresolved resonance regions), and on the other hand its evaluation is still quite “uncertain”, and evaluated cross-sections from various sources differ by more than their assigned uncertainties [21]. In Table 3 the explicit and implicit contributions to the total sensitivity coefficient of the 238U absorption cross-sections are given. One can observe that for some oxygen and uranium isotopes cross-sections the implicit part computed by the TSUNAMI code is not negligible.Table 2
Exercises I-1 and I-2 results:keff.
Test cases
k
eff
Uncertainty (%ΔR/R)
Exercise I-1
BWR
HZP
1.34050
5.23
E
-
01
HFP
1.22270
6.16
E
-
01
PWR
HZP
1.42290
4.82
E
-
01
HFP
1.40424
4.89
E
-
01
VVER
HZP
1.34498
5.13
E
-
01
HFP
1.32725
5.20
E
-
01
KRITZ-2 : 1
Cold
1.23394
5.87
E
-
01
Hot
1.18584
6.31
E
-
01
GEN-III
HFP
1.09591
5.20
E
-
01
Exercise I-2
PWR
HZP
1.11029
5.00
E
-
01
HFP
1.07736
5.56
E
-
01
BWR
HZP
1.41009
4.64
E
-
01
HFP
1.39351
4.71
E
-
01
GEN-III Type 1
1.25325
4.87
E
-
01
GEN-III Type 2
1.12304
4.94
E
-
01
GEN-III Type 3
1.04501
5.03
E
-
01
GEN-III Type 4
1.07008
9.68
E
-
01Table 3
Exercise I-1: PWR test case at HFP: energy integrated sensitivity coefficients of the238U absorption cross-section with respect to different cross-section isotopes.
Nuclide
Reaction
Explicit
Implicit
Sensitivity
16O
Total
3.00
E
-
02
1.19
E
-
02
4.19
E
-
02
16O
Scattering
2.99
E
-
02
1.19
E
-
02
4.18
E
-
02
16O
Elastic
2.99
E
-
02
1.19
E
-
02
4.18
E
-
02
16O
n,n′
1.37
E
-
05
-
1.95
E
-
08
1.37
E
-
05
16O
n,2n
4.32
E
-
12
-
2.24
E
-
17
4.32
E
-
12
16O
n,γ
-
8.21
E
-
07
1.00
E
-
08
-
8.11
E
-
07
16O
n,p
-
1.90
E
-
07
-
2.02
E
-
11
-
1.90
E
-
07
16O
n,d
-
4.02
E
-
08
-
1.07
E
-
12
-
4.02
E
-
08
16O
n,t
-
5.75
E
-
13
-
5.14
E
-
18
-
5.75
E
-
13
16O
n,alpha
1.12
E
-
04
2.69
E
-
08
1.12
E
-
04
234U
Total
-
4.27
E
-
05
2.05
E
-
06
-
4.07
E
-
05
234U
Scattering
-
9.13
E
-
08
1.47
E
-
06
1.37
E
-
06
234U
Elastic
-
9.43
E
-
07
1.47
E
-
06
5.23
E
-
07
234U
n,n′
8.50
E
-
07
3.27
E
-
10
8.50
E
-
07
234U
n,2n
1.99
E
-
09
-
1.12
E
-
13
1.99
E
-
09
234U
Fission
4.15
E
-
08
4.34
E
-
09
4.58
E
-
08
234U
n,γ
-
4.27
E
-
05
5.83
E
-
07
-
4.21
E
-
05
234U
Nubar
4.36
E
-
08
—
4.36
E
-
08
234U
Chi
-
3.13
E
-
15
—
-
3.13
E
-
15
235U
Total
-
1.01
E
-
01
3.36
E
-
03
-
9.75
E
-
02
235U
Scattering
1.08
E
-
03
8.97
E
-
04
1.97
E
-
03
235U
Elastic
-
9.45
E
-
05
8.96
E
-
04
8.01
E
-
04
235U
n,n′
1.16
E
-
03
8.01
E
-
07
1.16
E
-
03
235U
n,2n
1.38
E
-
05
-
1.14
E
-
09
1.38
E
-
05
235U
Fission
-
7.56
E
-
02
1.39
E
-
03
-
7.42
E
-
02
235U
n,γ
-
2.63
E
-
02
1.07
E
-
03
-
2.53
E
-
02
235U
Nubar
-
2.75
E
-
04
—
-
2.75
E
-
04
235U
Chi
1.22
E
-
11
—
1.22
E
-
11
238U
Total
7.85
E
-
01
-
8.43
E
-
02
7.01
E
-
01
238U
Scattering
3.25
E
-
02
-
5.69
E
-
02
-
2.44
E
-
02
238U
Elastic
4.68
E
-
03
-
5.69
E
-
02
-
5.22
E
-
02
238U
n,n′
2.75
E
-
02
8.91
E
-
06
2.75
E
-
02
238U
n,2n
3.38
E
-
04
8.26
E
-
08
3.38
E
-
04
238U
Fission
1.04
E
-
01
-
3.16
E
-
07
1.04
E
-
01
238U
n,γ
6.48
E
-
01
-
2.74
E
-
02
6.21
E
-
01
238U
Nubar
1.38
E
-
04
—
1.38
E
-
04
238U
Chi
-
1.12
E
-
12
—
-
1.12
E
-
12Figure 4
The five reaction cross-sections with the highest contribution to thekeff uncertainty for the test cases within Exercise I-1.Also, the most relevant sensitivity profiles for the BWR PB-2 and GEN-III unit-cell cases are shown in Figures5 and 6, respectively. One of the benchmark requirements was the evaluation of the uncertainty associated to the calculation of the one-group absorption and fission microscopic cross-sections for 235U and 238U within the test cases of Exercise I-1. The results are given in Table 4. The uncertainty of the microscopic cross-section values is around one order of magnitude higher than the one of the keff cases, ranging in between ~1% and ~4%. The highest uncertainty value was systematically found for the fission cross-section of 238U. As far as the test cases of Exercise I-2 are concerned, as required from the benchmark specification, some selected homogenized macroscopic cross sections with the associated uncertainties have been evaluated in a two-group structure with a cut-off energy of 0.625 eV. Results are provided in Table 5. The first energy group (E≤0.625 eV) was systematically found to be the one with the lower associated uncertainties. Also, as a general trend the uncertainties have been evaluated to be very consistent within all the test cases, and higher values (in the order of 1.35%) were computed for the homogenized absorption cross-section.Table 4
Exercises I-1: one-group microscopic cross-sections.
Exercise I-1test cases
Microscopic cross-section (barn)
235U absorption
238U absorption
235U fission
238U fission
(uncertainty %)
(uncertainty %)
(uncertainty %)
(uncertainty %)
BWR
HZP
2.01
E
+
01
(
1.04
E
+
00
)
3.10
E
-
01
(
9.80
E
-
01
)
1.65
E
+
01
(
1.06
E
+
00
)
3.01
E
-
02
(
4.07
E
+
00
)
HFP
4.15
E
+
01
(
1.22
E
+
00
)
8.80
E
-
01
(
9.70
E
-
01
)
3.34
E
+
01
(
1.22
E
+
00
)
8.57
E
-
02
(
4.79
E
+
00
)
PWR
HZP
4.41
E
+
01
(
1.09
E
+
00
)
9.42
E
-
01
(
9.65
E
-
01
)
3.57
E
+
01
(
1.10
E
+
00
)
9.91
E
-
02
(
3.92
E
+
00
)
HFP
4.30
E
-
01
(
1.10
E
+
00
)
9.64
E
-
01
(
9.71
E
-
01
)
3.47
E
+
01
(
1.11
E
+
00
)
9.91
E
-
02
(
3.94
E
+
00
)
VVER
HZP
5.98
E
+
01
(
1.02
E
+
00
)
1.02
E
+
00
(
9.88
E
+
01
)
4.90
E
+
01
(
1.04
E
+
00
)
9.26
E
-
02
(
3.86
E
+
00
)
HFP
5.84
E
+
01
(
1.03
E
+
00
)
1.04
E
+
00
(
9.90
E
-
01
)
4.78
E
+
01
(
1.05
E
+
00
)
9.27
E
-
02
(
3.88
E
+
00
)
KRITZ-2 : 1
Cold
8.23
E
+
01
(
1.07
E
+
00
)
1.01
E
+
00
(
1.02
E
+
00
)
6.83
E
+
01
(
1.09
E
+
00
)
8.90
E
-
02
(
4.20
E
+
00
)
Hot
2.71
E
+
01
(
1.13
E
+
00
)
3.84
E
-
01
(
1.01
E
+
00
)
2.24
E
-
01
(
1.15
E
+
00
)
3.38
E
-
02
(
4.51
E
+
00
)
GEN-III
HFP
1.53
E
+
01
(
1.38
E
+
00
)
9.18
E
-
01
(
1.01
E
+
00
)
1.10
E
+
01
(
1.22
E
+
00
)
1.17
E
-
01
(
3.83
E
+
00
)Table 5
Exercise I-2 results: macroscopic cross-sections.
Cross-section
Energy group
Value (cm−1) (Uncertainty %)
PWR
BWR
GEN-IIIType 1
GEN-IIIType 2
GEN-III Type 3
GEN-IIIType 4
Total
1
1.44
E
+
00
1.58
E
+
00
1.31
E
+
00
1.32
E
+
00
1.33
E
+
00
1.50
E
+
00
(
1.38
E
-
01
)
(
1.29
E
-
01
)
(
1.41
E
-
01
)
(
1.39
E
-
01
)
(
1.39
E
-
01
)
(
1.39
E
-
01
)
2
5.69
E
-
01
5.79
E
-
01
5.33
E
-
01
5.34
E
-
01
5.34
E
-
01
5.24
E
-
01
(
8.78
E
-
01
)
(
8.40
E
-
01
)
(
9.04
E
-
01
)
(
9.03
E
-
01
)
(
9.01
E
-
01
)
(
9.73
E
-
01
)
Absorption
1
1.11
E
-
01
5.72
E
-
02
1.07
E
-
01
1.17
E
-
01
3.45
E
-
1
1.24
E
-
01
(
8.77
E
-
01
)
(
6.06
E
-
01
)
(
7.00
E
-
01
)
(
5.61
E
-
01
)
(
9.79
E
-
01
)
(
5.11
E
-
01
)
2
1.06
E
-
02
7.32
E
-
03
1.04
E
-
02
1.06
E
-
02
5.09
E
-
01
1.07
E
-
02
(
1.33
E
+
00
)
(
1.38
E
+
00
)
(
1.35
E
+
00
)
(
1.34
E
+
00
)
(
1.47
E
+
00
)
(
1.34
E
+
00
)
Fission
1
7.95
E
-
01
2.94
E
-
02
6.86
E
-
02
6.60
E
-
02
1.90
E
-
01
6.42
E
-
02
(
3.17
E
-
03
)
(
3.23
E
-
01
)
(
3.23
E
-
01
)
(
3.24
E
-
01
)
(
6.26
E
-
01
)
(
3.24
E
-
01
)
2
3.59
E
-
03
1.95
E
-
03
3.17
E
-
03
3.10
E
-
03
4.97
E
-
01
3.05
E
-
03
(
3.55
E
-
01
)
(
6.81
E
-
01
)
(
3.71
E
-
01
)
(
3.75
E
-
01
)
(
4.45
E
-
01
)
(
3.80
E
-
01
)
Nufission
1
1.94
E
-
01
7.02
E
-
02
1.67
E
-
1
1.61
E
-
01
5.45
E
-
01
1.56
E
-
01
(
4.44
E
-
01
)
(
4.49
E
-
01
)
(
4.48
E
-
01
)
(
4.49
E
-
01
)
(
1.09
E
+
00
)
(
4.49
E
-
01
)
2
9.08
E
-
03
4.69
E
-
03
8.02
E
-
03
7.86
E
-
03
1.44
E
-
02
7.73
E
-
3
(
5.12
E
-
01
)
(
1.01
E
+
00
)
(
5.71
E
-
01
)
(
5.82
E
-
01
)
(
7.75
E
-
01
)
(
5.92
E
-
01
)Figure 5
k
eff sensitivity profiles for the BWR PB-2 unit cell case (HFP conditions).Figure 6
k
eff sensitivity profiles for the GEN-III unit cell case (HFP conditions).


## 5. The Reference Solutions

As part of the activities within the UAM benchmark. KIT and PSU are also working on the development of reference solutions for the test cases by means of Monte Carlo methodologies. As it is well known, the Monte Carlo method can be considered as a “numerical experiment” that represents a high quality reference solution for the validation of deterministic codes. For this purpose we are mainly using the SERPENT code [22], a Monte Carlo code designed for lattice physics calculations [23]. A comparison between absorption and fission microscopic cross sections values as computed with SERPENT and SCALE and relative to three test cases within Exercise I-1 is summarized in Table 6. A very good agreement between the two approaches can be observed. keff values computed by means of the SERPENT code (v. 1.1.16) with different sets of NDLs are summarized in Table 7. As far as the modeling, a good statistics in the criticality calculations has been achieved by simulating 5.0E+6 neutron histories over 1000 active cycles (5000 neutron source × 1000 cycles). Continuous cross-section libraries in SERPENT are given in steps of 300 K (300 K, 600 K, etc.), and the ‘‘actual’’ problem temperatures have been simulated by defining a mix of two materials with proportions equal to the weighted averages in between the two closest temperatures sets as available in SERPENT package. One can observe quite a good agreement in between the results with the different NDLs. Also, a systematic slight increase in the prediction of the keff when switching from JEFF3.1 to JEFF3.1.1 up to ENDF/B-VII, respectively, can be noted. keff computed with JEFF3.1.1 are larger by ~100 pcm with respect to those obtained from JEFF3.1, and this is mainly due to the improvement of the 16O(n,α) cross section which was decreased by 30% in the JEFF3.1.1 evaluation [24]. The comparison of these values with the correspondent ones obtained from the SCALE6.1 code (see Table 2) shows some discrepancies. One of the reasons for these differences was found to be due to the different Boundary Conditions (BCs) used in the SERPENT modeling; the “white” BC is not yet implemented in SERPENT, and therefore the “reflective” BC has been adopted for our simulations. In order to assess the impact of the different BCs on the computed reactivities, three test cases at HZP conditions have been modeled with the MCNP5 code (v. 1.40) [25] with “white” BCs. The corresponding results, which are given in Table 7 and have been obtained by simulating 1500 generation of 2000 neutrons each, show how the discrepancies are then reduced.Table 6
Phase I-1 Exercises.
Micro-XS
SCALE 6.1
SERPENT
Uncertainty due to nuclear data (%)
Unit cell
235U absorption
41.48
40.41 ± 0.0086
1.22
BWR
238U absorption
0.88
0.80 ± 0.0011
0.97
235U fission
33.43
32.56 ± 0.00069
1.22
238U fission
0.086
0.089 ± 0.00097
4.79
235U absorption
42.95
42.18 ± 0.00088
1.09
PWR
238U absorption
0.96
0.93 ± 0.0011
0.97
235U fission
34.72
34.10 ± 0.00064
1.11
238U fission
0.099
0.10 ± 0.00096
3.94
235U absorption
58.13
57.26 ± 0.00085
1.03
VVER
238U absorption
1.042
1.005 ± 0.0012
0.99
235U fission
47.84
47.76 ± 0.00063
1.05
238U fission
0.093
0.095 ± 0.00100
3.88Table 7
Phase I-1 Exercise:keff evaluation with SERPENT and MCNP*.
Test case
k
eff
JEFF3.1
JEFF 3.1.1
ENDFB-7
VVER
HZP
1.34764 ± 0.00028
1.34937 ± 0.00026
1.34986 ± 0.00027
1.34542 ± 0.00029*
HFP
1.33152 ± 0.00028
1.33356 ± 0.00029
1.33435 ± 0.00029
PWR
HZP
1.42785 ± 0.00027
1.42888 ± 0.00025
1.42923 ± 0.00027
1.42018 ± 0.00030*
HFP
1.41136 ± 0.00026
1.41315 ± 0.00028
1.41401 ± 0.00026
BWR
HZP
1.34541 ± 0.00027
1.34673 ± 0.00025
1.34691 ± 0.00026
1.34191 ± 0.00034*
HFP
1.23046 ± 0.00032
1.23080 ± 0.00032
1.23295 ± 0.00032
KRITZ-2 : 1
Cold
1.23762 ± 0.00028
1.23846 ± 0.00027
1.23984 ± 0.00027
Hot
1.22632 ± 0.00028
1.22864 ± 0.00026
1.22863 ± 0.00027
GEN-III
HFP
1.01485 ± 0.00039
1.01602 ± 0.00039
1.01805 ± 0.00037


## 6. Conclusions

The neutronics modeling of some selected test cases within the “Neutronics Phase” of the OECD UAM Benchmark has been presented. A S/U analysis on the impact of the uncertainty in the basic nuclear data on the calculation of the multiplication factor and microscopic and macroscopic cross-sections have been performed using the perturbation methodologies implemented in the TSUNAMI code. Uncertainties were found to be~0.5% on keff and higher (up to ~4%) for the cross-sections. The 238U capture cross-section and the 239Pu nubar were found to be the highest contributors to the total uncertainty for the UOX and MOX LWR’s representative test cases. The deterministic solutions were also compared with the corresponding reference solutions obtained using Monte Carlo methods, and a good agreement has been found for several cases, particularly in the evaluation of the microscopic cross-sections.



---

*Source: 573697-2013-02-07.xml*