Research Article  Open Access
O. Cabellos, "Presentation and Discussion of the UAM/Exercise I1b: “PinCell BurnUp Benchmark” with the Hybrid Method", Science and Technology of Nuclear Installations, vol. 2013, Article ID 790206, 12 pages, 2013. https://doi.org/10.1155/2013/790206
Presentation and Discussion of the UAM/Exercise I1b: “PinCell BurnUp Benchmark” with the Hybrid Method
Abstract
The aim of this work is to present the Exercise I1b “pincell burnup benchmark” proposed in the framework of OECD LWR UAM. Its objective is to address the uncertainty due to the basic nuclear data as well as the impact of processing the nuclear and covariance data in a pincell depletion calculation. Four different sensitivity/uncertainty propagation methodologies participate in this benchmark (GRS, NRG, UPM, and SNU&KAERI). The paper describes the main features of the UPM model (hybrid method) compared with other methodologies. The requested output provided by UPM is presented, and it is discussed regarding the results of other methodologies.
1. Introduction to the UAM/Exercise I1b “PinCell BurnUp Benchmark”
The general frame of the OECD LWR UAM benchmark consists of three phases with different exercises for each phase [1]. In the Phase I (“Neutronics Phase”), the Exercise 1 (I1) “Cell Physics” is focused on the derivation of the multigroup microscopic crosssection libraries. Since the OECD LWR UAM benchmark establishes a framework for propagating crosssection uncertainties in LWR design and safety calculations, the objective of the extension of this Exercise I1 to I1b (cell burnup physics) is to address the uncertainties in the depletion calculation due to the basic nuclear data as well as the impact of processing of nuclear and covariance data. The SCALE6.0/1 covariance library [2] is the recommended source of crosssection data uncertainty. However, covariance data coming from other source of uncertainty together with evaluated nuclear data files can be used without any inconvenience.
To address this problem different sensitivity/uncertainty (S/U) tools can be used to propagate nuclear data (e.g., crosssection) uncertainties. The requested output of Exercise I1b is criticality value, reactions rates, collapsed crosssections and nuclide concentrations as well as their uncertainties for depletion in a PWR pincell model.
1.1. Specifications of the “PinCell BurnUp Benchmark”
The specification of this pincell benchmark is given in Tables 1 and 2 (geometry and material specifications), showing a typical configuration of a TMI1 PWR unit cell.


The linear fuel density (gU/cm) calculated according to values taken from Tables 1 and 2 is 6.2784 gU/cm. The average power density (W/gU) can be assumed to be equal to 33.58 W/gU. The fuel sample is burned for a unique complete cycle, the length of the burn time, and subsequent cooling time is given in Table 3. The specific power and the final cumulative burnup are also given, 61.28 GWd/MTU.

Concerning boundary conditions, the following type of boundary conditions can be used: (a) for a “cylindrical pincell” model, reflective boundary conditions are utilized at the centerline boundary while white boundary conditions are applicable at the peripheries of the cell model; (b) for a “square pincell” model, reflective boundary conditions on all surfaces are applied. For depletion, it can be considered an infinite burnup spectrum mode.
1.2. Requested Output of the “PinCell BurnUp Benchmark”
Results and associated uncertainties are provided at eight burnup steps: 0, 10, 20, 30, 40, 50, 60, and shutdown (61.28) GWd/MTU. And, six additional decay steps are required at 1, 3, 5, 10, 50, and 100 years of cooling time. The requested output can be summarized in the following three sets of information:(i)criticality values: Kinf and nuclide reactions that contribute the most to the uncertainty in kinf;(ii)reaction rates and collapsed macroscopic crosssections:(a)Reaction rates (capture and fission) and uncertainties for major isotopes: ^{235,238}U and ^{239,240,241}Pu;(b)Twogroup macroscopic crosssections, fast and thermal, and associated uncertainties for the homogenized pin cell: absorption, fission, nufission, and diffusion coefficient. The thermal energy cutoff is 0.625 eV. (iii)Number densities:(a)actinides (15): ^{233,234,235,236,238}U; ^{237}Np; ^{238,239,240,241,242}Pu; ^{241,243}Am; ^{244,246}Cm;(b)fission products (36): ^{95}Mo; ^{99}Tc; ^{101,106}Ru; ^{103}Rh; ^{109}Ag; ^{133,134,135,137} Cs; ^{139}La; ^{140,142,144}Ce; ^{142,143,145,146,148,150}Nd; ^{147,148,149,150,151,152,154}Sm; ^{151,153,154,155}Eu; ^{154,155,156,158,160}Gd.
2. Summary of Propagation Uncertainty Methodologies in BurnUp Calculations
The first phase of participation in this exercise was completed in April 2012 with a total of 4 participants: GRS, NRG, UPM, and SNU&KAERI. Table 4 summarizes the main calculation methodologies and nuclear data libraries and their uncertainties. The results were presented at the Sixth Workshop (UAM6) of OECD Benchmark for Uncertainty Analysis in BestEstimate Modelling (UAM).

On one hand, depletion calculations are performed by GRS and UPM with SCALE6 code system [3], while NRG uses SERPENT code [4] and SNU&KAERI participates in the benchmark with its own McCARD code [5], both Monte Carlo codes. On the other side, for uncertainty calculations, GRS and NRG use Monte Carlo techniques, GRS with a sampling methodology (XSUSA [6]) of multigroup crosssection libraries provided in SCALE6 format and NRG using the technique of Total Monte Carlo [7] with TENDL2011. UPM applies a hybrid method [8] based on determining the sensitivity coefficients with TSUNAMI code [9] and performing a Monte Carlo sampling to determine the uncertainty of the number densities; these uncertainties are computed with ACAB code [10]. McCARD code makes use of the technique of Adjoint Weighted Perturbation (AWP) method to predict the sensitivity coefficients.
Regarding crosssection covariance data, GRS, SNU&KAERI, and UPM use SCALE6/COVA44 groups. In addition, SNU&KAERI provides results with uncertainties coming from JENDL3.3 and ENDF/BVII.0. Figure 1 shows an example of crosssection covariance data taken from SCALE6.1/COVA44G. In this figure, the original ^{235}U COVERX/SCALE6.1 file is processed with ANGELO, LAMBDA, and NJOY codes to visualize the correlation matrix. NRG uses TENDL2011 and their uncertainty for crosssection data libraries. In addition, NRG and UPM have carried out some calculations with the uncertainty provided in Fission Yields (TENDL2011, JEFF3.1.1) and Decay Data (JEFF3.1.1) libraries.
Next, the main characteristics of the uncertainty propagation methodologies used in this Benchmark are summarized, and the uncertainty propagation in number density is used as an example in the following Figures 2, 3, and 6. Figure 2 shows the calculation scheme of the Monte Carlo methodologies. NRG uses for each sampling a different nuclear data library TENDL2011; the generation of this library is done using the TASMAN code [7]. TASMAN is a computer code for the production of covariance data using results of the nuclear model code TALYS, and for automatic optimization of the TALYS results with respect to experimental data. It is assumed that each nuclear model (i.e., TALYS input) parameter has its own uncertainty; running TALYS many times, it provides a sampling of ENDF files or a single file with full covariance information. GRS will generate a set of multigroup libraries in SCALE6 format; this sampling is done with the SCALE6.1/44groups covariance library using XSUSA code. The sensitivity/uncertainty procedure is based on a first order Taylor series approach. So, the number density can be written as where .
We can define the sensitivity coefficients as , and is the error in the 1group effective crosssections. This 1group error depends explicitly on the uncertainty of crosssections, and implicitly on the neutronflux uncertainty, Here, is the error due to nuclear data and is the error due to neutronflux. The variance in the number density can be obtained using the sandwich formula: The first term propagates the multigroup crosssection uncertainty with no uncertainty in the neutron flux. And, the second term propagates the effect of this uncertainty with the uncertainty in the neutron flux.
If the uncertainty in the neutron flux can be considered negligible, a simple scheme of S/U can be illustrated in Figure 3. In this case, TRITON code [3] is run to determine the number densities at different burnup steps, as a reference or nominal calculation without uncertainties. And, the number densities calculated in the nominal case are used to generate TSUNAMI [9] inputs at each burnup step. With TSUNAMI code, S/U analysis can be provided for criticality , twogroup crosssections and reaction rates . However, number density sensitivities are not calculated with TSUNAMI code.
Once, the sensitivity coefficients are calculated by TSUNAMI code, the criticality uncertainty analysis based on “nuclear data uncertainties” can be formulated as follows: it is explicitly dependent on the nuclear data (e.g., crosssections, nubar, …) and implicitly dependent on the number density which characterizes the system: is the sensitivity coefficient explicitly of crosssections and is the sensitivity coefficient of number density, ; both are calculated by TSUNAMI code. Figures 4 and 5 show the eff integrated sensitivity coefficients for crosssection and number density at each burnup step. In Figure 4, the evolution of shows the importance of ^{239}Pu at high burnups, mainly for nubar nuclear reaction. For ^{238}U, and reactions are the most important for all burnup. For ^{235}U, sensitivity decreases with burnup, being nubar with the highest value. Evolution of ^{135}Xe is also shown. Some “fissiongamma” crosscorrelations for ^{239}Pu and ^{235}U are also illustrated. Figure 5 shows the integrated sensitivities, , for the most important isotopes related with criticality: ^{239,240, 241}Pu, ^{235,238}U. Also, some important fission products are shown: ^{135}Xe and ^{103}Rh.
is the covariance crosssection data taken from SCALE6.1/COVA, and is the covariance number densities predicted by ACAB code. It can be calculated with the uncertainty due to crosssection, fission yield and/or decay data. Our ACAB code is used to propagate nuclear data uncertainty (crosssection, fission yield, and decay data) in the prediction of number density uncertainty: ACAB accounts for the impact of nuclear data uncertainty as follows (see Figure 6). (i) In a first step, a coupled neutrondepletion calculation (without uncertainties) is carried out only once, taken the bestestimated values for all the parameters involved in the problem. (ii) In a second step, ACAB performs a simultaneous random sampling of the probability density functions (PDF) of all these variables: crosssection, fission yield, and decay data. Then, ACAB computes the isotopic concentrations at the end of each burn step, taking the fluxes halfway through each burn step determined in the bestestimated calculation. Then, only the depletion calculations are repeated or run many times. A statistical analysis of the results allows assessing the uncertainty in the calculated number density and determining . Table 5 shows an example of this type of information.

3. Results with the Hybrid Method
In Table 6, and their associated uncertainty for PWR unitcell are summarized at four different burnups. The five most important nuclide reactions that contribute to uncertainty are identified: (i) for fresh fuel, , and , , , and (ii) for high burnup: and , , and . In addition, the contribution of number density uncertainty, var, is evaluated, being the crosssections and fission yields the most important contributions, and it can be concluded that the contribution of decay data uncertainty is negligible.

Table 7 shows the uncertainty of twogroup crosssections: , and (subscript 1 refers to fast group and subscript 2 to the thermal group). The low contribution of the uncertainty due to number density uncertainty except for thermal groups can be seen. The total uncertainty is about 1%, and the contribution due to the uncertainty in fission yields is negligible.

As an example of integrated sensitivities of macroscopic twogroup crosssections, Figures 7 and 8 show these values for . ^{238}U is the most important contributor with the and reactions.
Table 8 shows the uncertainty for the following capture and fission reaction rates: ^{235,238}U and ^{239,240,241}Pu. The total uncertainty is in the range of 1%–3%. In general, the uncertainty contribution due to the uncertainty in the number density is below the contribution due to crosssection , except for ^{240}Pu and ^{241}Pu reaction rates where this contribution is larger.

In Table 9, it can be seen that the number density uncertainty for some major and minor actinides due to crosssection data remains below 3%. Larger uncertainties are predicted for minor actinides (e.g., ^{246}Cm) and the uncertainty throughout irradiation period rises. And, it can be concluded that the uncertainty due to decay data uncertainty is negligible.

In Table 10, the uncertainty in the number of fission products due to crosssections, decay, and fission yields data has been predicted. Some isotopes, ^{155}Gd, ^{154,155}Eu, and ^{149}Sm show a relative error above 10%, being the high uncertainty in crosssection data, the reason of this large uncertainty. In general, the uncertainty due to fission yields remain below 3%, except for ^{95}Mo with 4.5% (with high sensitivity to ^{95}Zr fission yield) and ^{149}Sm with 4.7% (with high sensitivity to ^{149}Pm fission yield) [11]. For decay data uncertainties, the isotope ^{151}Eu reaches a maximum uncertainty of 3.2% as a consequence of the 6.7% relative error in the halflife of ^{151}Sm.
