# A Detailed Analysis of the H.B. Robinson-2 Reactor Pressure Vessel Dosimetry Benchmark

Romain Vuiart, Mariya Brovchenko, Julien Taforeau, Eric Dumonteil

► To cite this version:

Romain Vuiart, Mariya Brovchenko, Julien Taforeau, Eric Dumonteil. A Detailed Analysis of the H.B. Robinson-2 Reactor Pressure Vessel Dosimetry Benchmark. Energies, 2022, 15 (14), pp.5098. 10.3390/en15145098. irsn-04122073

HAL Id: irsn-04122073

https://irsn.hal.science/irsn-04122073

Submitted on 8 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# A Detailed Analysis of the H.B. Robinson-2 Reactor Pressure Vessel Dosimetry Benchmark

Romain Vuiart<fnref n="1" /> ID, Mariya Brovchenko<fnref n="1" />*, Julien Taforeau<fnref n="1" /> and Eric Dumonteil<fnref n="2" />

**Abstract:** The operation of many nuclear pressurized water reactors is being extended beyond their design lifetime limit. From the perspective of possible further lifetime extension, safety requirements are a priority. Therefore, the quantification of the neutron irradiation embrittlement of the reactor pressure vessel (RPV) is an important issue, as this is a guiding parameter that influences the reactor lifetime. In this context, the Institut de Radioprotection et de Sûreté Nucléaire developed a calculation scheme for the analysis of RPV aging under neutron irradiation, named VACS (vessel aging calculation scheme). VACS couples a deterministic approach (CASMO5 and SIMULATE5) to evaluate the full-core fission neutron source term and a Monte Carlo modeling (MCNP6) approach to model the neutron attenuation from the core to sites of interest (RPV, surveillance capsules, etc.). To ensure the reliability of aging predictions, this paper describes a detailed analysis of the neutron H.B. Robinson-2 reactor pressure vessel dosimetry benchmark. The results indicate that VACS shows satisfactory accuracy when the ENDF-B/VII.1 or JEFF-3.3 nuclear data libraries are used in the attenuation calculation. However, the use of ENDF-B/VIII.0 leads to significantly worse results.

**Keywords:** nuclear safety; vessel aging; vessel embrittlement; validation; dosimetry benchmark; H.B. Robinson-2; neutronics

## 1. Introduction

The lifespan of many nuclear pressurized water reactors (PWRs) is being extended beyond their design lifetime. With the prospect of further lifetime extensions, meeting safety requirements is a priority. In particular, the characterization of the structural integrity of the reactor pressure vessel (RPV) is an important issue. Indeed, under the operation of a PWR, the RPV is bombarded by neutrons from fission reactions occurring in the reactor core. These incident neutrons (especially fast neutrons, i.e., those with energies above 1 MeV [1]) induce numerous atomic displacements in the RPV material, causing the mechanical properties of the RPV to gradually deteriorate. To monitor this embrittlement, surveillance programs have been implemented for the majority of RPVs around the world [2–5]. These surveillance programs include the so-called capsules, which contain vessel material specimens as well as dosimeters. The capsules are usually attached to the core barrel or the thermal neutron shield, in order to increase their exposure to neutron fluxes with regard to the RPV. Thus, the material specimens age faster than the RPV inner wall. Mechanical tests are performed on the specimens, which are withdrawn regularly, and the acquired data, in conjunction with the data provided by the dosimeters, allow an empirical correlation to be established between the accumulated fast neutron fluence and the deterioration of the mechanical properties of the material. The neutron fluence in the vessel is based on a simulation, and determining the precision of this simulation is the aim of this study.

In this context, efforts have been made at the Institut de Radioprotection et de Sûreté Nucléaire (IRSN) to develop a versatile, best-estimate calculation scheme for the analysis

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Vuiart, R.; Brovchenko, M.; Taforeau, J.; Dumonteil, E. A Detailed Analysis of the H.B. Robinson-2 Reactor Pressure Vessel Dosimetry Benchmark. *Energies* **2022**, *15*, 5098.
https://doi.org/10.3390/en15145098

Academic Editor: Hiroshi Sekimoto

Received: 28 May 2022

Accepted: 27 June 2022

Published: 12 July 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

Copyright: © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSN-RES/SNC/LN, F-92260 Fontenay-aux-Roses, France; romain.vuiart@gmail.com (R.V.); julien.taforeau@irsn.fr (J.T.)</footnote>

<footnote n="2">² Commissariat à l'Énergie Atomique et aux Énergies Alternatives (CEA), Institut de Recherche sur les Lois Fondamentales de l'Univers, Université Paris Saclay, F-91191 Gif-sur-Yvette, France; eric.dumonteil@cea.fr</footnote>

<footnote>* Correspondence: mariya.brovchenko@irsn.fr</footnote>
of RPV embrittlement under neutron irradiation [6], named VACS (vessel aging calculation scheme). The methodology relies on a two-step approach, coupling deterministic and Monte Carlo calculations. The first step estimates the full-core fission neutron source term and the second step considers the transport of neutrons from the core to the RPV, using the fission neutron distribution determined in the previous step. Throughout the construction of the methodology, several studies were conducted to minimize the inherent methodological biases [6–9] and to pave the way towards the quantification of the global uncertainty in its predictions [9,10]. In order to validate the numerical predictions using experimental dosimeter data from ex-core locations, a detailed analysis of the H.B. Robinson-2 (HBR-2) RPV dosimetry benchmark [11] was conducted. This benchmark was chosen because it includes measurements performed for an operating PWR using in-vessel (capsule location) and ex-vessel (reactor cavity) dosimeters. This feature is of major interest because it allows the accuracy of aging predictions to be characterized, not only for in-vessel applications but also, to some extent, for ex-vessel ones. In addition, the HBR-2 benchmark provides a fission neutron source term which is intended to be directly used to validate the accuracy of the neutron attenuation calculation only (second step). However, the data provided were used to analyze the accuracy of both steps of the calculation scheme (fission neutron source term and attenuation) and thus to provide a complete validation.

The H.B. Robinson-2 Reactor and the data provided in the benchmark are presented in Section 2. Section 3 highlights the modeling and the methods used to calculate the fission neutron source term and model the attenuation of these neutrons at the sites of interest. Section 4 presents and discusses the results obtained at each step of the analysis. Finally, conclusions are drawn in Section 5.

## 2. The H.B. Robinson-2 Reactor Pressure Vessel Dosimetry Benchmark

The H.B. Robinson-2 (HBR-2) neutron dosimetry benchmark [11] from the SINBAD database [12] is widely used in the scientific community for the validation of calculation schemes designed to predict vessel aging, as shown in [13–16]. It features in-vessel (surveillance capsule) and ex-vessel (reactor cavity) dosimetry data acquired during the 9th cycle of the H.B. Robinson-2 reactor.

### 2.1. *Description of the H.B. Robinson-2 Reactor*

The H.B. Robinson-2 reactor is a 2300 MWth PWR designed by Westinghouse and owned by the Carolina Power and Light Company. The reactor is located near Hartsville in South Carolina, United States. Its operation began in March 1971.

The reactor core is loaded with 157 UOX assemblies, initially enriched to 2.9% of $^{235}$U. The assemblies contain 204 fuel rods plus 21 guide tubes that serve as inserts for control rods. Fuel-rod and guide-tube claddings are made of Zircaloy. A visualization of the geometry of a fuel assembly is given in Figure 1. The reactor core is surrounded by the baffle, the moderator bypass region, the core barrel, a first moderator downcomer region, the thermal neutron shield, a second moderator downcomer region, the RPV, the reactor cavity and the biological shield. These elements can be found in Figure 2, which is taken from [14] and presents a simplified view of an eighth of the reactor.

On the one hand, the dimensions of the majority of the components at or above assembly scale (assembly, baffle, barrel, thermal neutron shield, RPV, etc.) are provided in the benchmark data [11]. On the other hand, the dimensions of some components below the assembly scale (fuel pins, guide tubes, etc.) are not specified in [11] but can be found in [17]. Descriptions of the materials constituting the main elements of the reactor (fuel, baffle, cladding, thermal neutron shield, vessel, etc.) are also provided in the benchmark data [11].
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Radial cut of a UOX fuel assembly from the H.B. Robinson-2 reactor. Fuel rods are shown in green.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Schematic radial cut of an octant of the H.B. Robinson-2 (taken from [14]).</figcaption></figure>

The HBR-2 reactor has four groups of black control banks and one bank of shutdown rods. The insertion layout of these banks can be found in [17]. However, no information could be found in the literature on the insertion sequence of these control rods during core power variations.

The reactor is equipped with several sets of dosimeters, whose locations are marked by crosses in Figure 2. A first set of dosimeters is placed in a surveillance program capsule attached to the thermal neutron shield (azimuth 20°; 191.15 cm from the reactor center). These dosimeters are inserted in holes drilled in the center of the capsule. A second set of dosimeters is placed in a cavity behind the vessel (azimuth 0°; 238.02 cm from the reactor center). Cavity dosimeters are held by an aluminum 6061 holder, which is supported by stainless steel gradient wires mounted vertically in the gap between the insulation and the biological shield (see [11]). Both sets of dosimeters were inserted into the reactor at the beginning of its 9th irradiation cycle. Specific activities were measured at the end of this cycle, for the two dosimeter locations and for the six types of dosimeters inserted: $^{237}$Np(n,f)$^{137}$Cs, $^{238}$U(n,f)$^{137}$Cs, $^{58}$Ni(n,p)$^{58}$Co, $^{54}$Fe(n,p)$^{54}$Mn, $^{46}$Ti(n, p)$^{46}$Sc and $^{63}$Cu(n, α)$^{60}$Co). Only specific activities in the equatorial plane of the reactor (0 cm axial position) are detailed in [11]. These experimental data (M) constitute the reference values to which the calculated values (C) must be compared in order to validate the numerical predictions of a calculation scheme. The provided specific activities are detailed in Table 1 [11]. However, it is worth mentioning that the authors in [11] do not provide uncertainty values for these experimental data.
<table><caption>**Table 1.** Specific activities measured in the equatorial plane, at the end of cycle no. 9, for the dosimeters placed in the capsule of the surveillance program and in the cavity of the H.B. Robinson-2 reactor [11].</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Dosimeter Type</th><th colspan="2">Dosimeter Location</th></tr><tr><th>Capsule</th><th>Reactor Cavity</th></tr></thead><tbody><tr><td rowspan="6">Specific activities<br>(Bq/mg)</td><td><sup>237</sup>Np(n,f)<sup>137</sup>Cs</td><td>3.671·10<sup>+2</sup></td><td>2.236·10<sup>+1</sup></td></tr><tr><td><sup>238</sup>U(n,f)<sup>137</sup>Cs</td><td>5.345·10<sup>+1</sup></td><td>8.513·10<sup>-1</sup></td></tr><tr><td><sup>58</sup>Ni(n,p)<sup>58</sup>Co</td><td>1.786·10<sup>+4</sup></td><td>1.959·10<sup>+2</sup></td></tr><tr><td><sup>54</sup>Fe(n,p)<sup>54</sup>Mn</td><td>9.342·10<sup>+2</sup></td><td>8.711·10<sup>+0</sup></td></tr><tr><td><sup>46</sup>Ti(n,p)<sup>46</sup>Sc</td><td>3.500·10<sup>+2</sup></td><td>3.310·10<sup>+0</sup></td></tr><tr><td><sup>63</sup>Cu(n,α)<sup>60</sup>Co</td><td>2.646·10<sup>+1</sup></td><td>2.645·10<sup>-1</sup></td></tr></tbody></table>

To help the reader better understand the following results, Figure 3 shows the activation cross sections for the six types of dosimeters considered, processed at 600K using the NJOY2016 code [18] with ENDF files taken from the JEFF-3.3 data library [19]. Figure 3 shows that all dosimeters are sensitive to fast neutrons (i.e., those with energies above 1 MeV), making them highly relevant for assessing the reliability of vessel aging assessments (given that fast neutrons are the primary contributors to radiation damage [6]). Finally, it is interesting to note that <sup>58</sup>Ni(n,p), <sup>54</sup>Fe(n,p), <sup>46</sup>Ti(n,p) and <sup>63</sup>Cu(n,α) represent reactions with a strong energy threshold.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Activation cross sections taken from the JEFF-3.3 library and processed at 600 K.</figcaption></figure>

## 2.2. *Description of Cycle No. 9 of the H.B. Robinson-2 Reactor*

The data provided in [11] allow a partial reconstruction of the course of the 9th irradiation cycle of the HBR-2 reactor to be performed. The initial burn-up distribution of assemblies is not perfectly symmetrical. However, such slight asymmetries are expected in a real reactor, since fuel assemblies are loaded in different areas of the core depending on the number of irradiation cycles they have undergone. Nevertheless, not all assemblies with the same number of irradiation cycles have experienced exactly the same history, which induces asymmetries in the initial assembly burn-up distribution.

After its loading, the reactor was in operation between 21 August 1982 and 26 January 1984 [15]. The power history of the HBR-2 reactor over this time period is provided in [11]. From this data, it is possible to know that the HBR-2 reactor did not operate at nominal power during cycle no. 9, but operated for most of the time at 80% of its nominal power
(NP) and experienced numerous power drops. The average power of the reactor over the cycle was equal to 75.2% NP, which corresponds to roughly 1729.6 MWth. At such power levels, it is very likely that the control rods were partially inserted during the cycle. However, in the absence of information on the insertion sequence of the control-rod banks, it is not possible to reconstruct the average insertion position of each rod bank over cycle no. 9.

Other data included in the benchmark [11] give the average power distribution of the reactor over cycle no. 9. This power distribution is provided at the assembly level for the entire core and at the fuel-pin level for the upper right quarter. Axially, a single power distribution with 12 divisions distributed over the active height of the reactor is provided. The latter is equal to the average of the axial power distributions of the two outermost rows of assemblies. The authors of [11] specify that the assembly-level cycle-averaged power distribution was calculated using the assembly burn-up distributions at the beginning and end of cycle no. 9, taking into account the uranium content of each assembly. However, no further details are provided on how these burn-up distributions were obtained. In the absence of additional information, it can be assumed here that the burn-up distributions were evaluated by combining experimental data and neutronics calculations.

In addition, [11] provides the cycle-average fuel-pin power distribution, which was obtained by performing a weighted average of the pin powers given at eight steps during cycle no. 9 (147, 417, 1632, 3363, 5257, 7595, 9293 and 10,379 MWd/t). In the same way as before, no further details are provided on the evaluation of these radial pin power distributions, and it is reasonable to assume here that these distributions were evaluated by means of neutronics calculations. The analysis of the provided data lead to notice that the power values of several fuel pins in assemblies J15, K15, L15, N13, J12, R11, R10, M09 and R09 are very low compared to the power values of neighboring fuel pins (the assembly positions are presented in Figure 7 of Section 4.1). It was found in Version 21 of the Final Safety Analysis Report for the HBR-2 reactor [17] that, at some point, integral burnable absorber rods were thought to be used in the reactor. These rods were composed of gadolinia ($Gd_2O_3$) in various concentrations in $UO_2$, and the goal was to control the peaking and to reduce the critical boron at the beginning of the cycle. Although there is no mention of the presence of any absorber fuel rods in the benchmark data [11], it seems plausible that the absorber rods were present in the core during the experimental campaign, and this could explain the observed lower power values of the above-mentioned fuel pins. Nevertheless, due to the lack of a more detailed description, no gadolinia-bearing fuel pins were modeled using VACS (Section 3.1.3).

It should also be mentioned that, apart from the reactor dimensions, no uncertainty estimates are provided in [11] for the specific activity measurements or the reactor power distributions. Therefore, the provided power distributions will not be considered here as an absolute reference but rather as a point of comparison, to put into perspective the fission neutron source term calculated using VACS.

### 3. Modeling and Methods

In order to validate VACS, specific activity estimations for dosimeters placed at the monitoring capsule and in the cavity were performed. In addition, a second methodology that performed the same analysis using the assumptions and data provided in [11] was employed for comparative purposes. Thus, in this section, these two methods are described, as well as the models developed to perform the study.

#### 3.1. *Evaluation of the Fission Neutron Source Term*

As outlined previously, the first step of the decoupled calculations aims at evaluating the fission neutron source term, which describes the spatial and energy distributions of fission neutrons in the reactor. The fission neutron source term can be evaluated by means of numerical calculations or from experimental power distributions measured in the reactor.

In this analysis, two approaches were used to assess the fission neutron source term:
* The α approach: the fission neutron source term was evaluated using only the data provided in the benchmark (no additional calculations were made). This first approach aims to follow as closely as possible the assumptions recommended in [11]. It will be used to put into perspective the results obtained using VACS.

* The β approach: the fission neutron source term was evaluated using VACS, which couples the CASMO5 [20] and SIMULATE5 [21] codes to perform the fission neutron source term assessment.

The following sections provide a description of the α and β approaches. Prior to that, a reminder of the typical approximations made in decoupled calculation schemes is given.

### 3.1.1. Typical Approximations in Decoupled Approaches

As well as the angular dependency, the number of fission neutrons emitted per second at a time t, per unit volume around a point $\vec{r}$ in space and per unit energy around an energy $E'$, can be defined as:

$$S(\vec{r}, E', t) = \sum_{i} \int_{E} R_{f,i}(\vec{r}, E, t) \cdot v_{i}(E) \cdot \chi_{i}(E \rightarrow E') \cdot dE \quad (1)$$

where, for an isotope i, $R_{f,i}$ is the fission rate (cm$^{-3} \cdot s^{-1} \cdot eV^{-1}$), $v_i$ is the number of neutrons emitted per fission (n), $\chi_i$ is the fission spectrum (eV$^{-1}$), E is the energy of the neutron inducing the fission (eV) and $E'$ is the energy of the emitted neutron (eV).

By identifying the space-energy-dependent fission neutron emission rate of isotope i, $F_i(E) = R_{f,i}(E)v_i(E)$ in Equation (1), it can be rewritten according to Equation (2):

$$S(\vec{r}, E', t) = \sum_{i} \int_{E} F_{i}(E) \cdot \chi_{i}(E \rightarrow E') \cdot dE \quad (2)$$

where $F_i$ is the fission neutron emission rate (n · cm$^{-3} \cdot s^{-1} \cdot eV^{-1}$), $\chi_i$ is the fission spectrum (eV$^{-1}$), E is the energy of the neutron inducing the fission (eV) and $E'$ is the energy of the emitted neutron (eV). F represents the spatial distribution of fission neutrons, whereas $\chi$ describes the energetic distribution of these neutrons.

However, two-step calculation schemes such as those used in this article usually perform spatial, energy and isotopic averaging approximations [7], leading to fission neutron source terms of the form of Equation (3).

$$S_{M_1}(E', t) = F_{M_1}(t) \cdot \chi_{M_2}(E', t) \quad (3)$$

where, for a medium $M_1$ and a medium $M_2$ (which may or may not be the same), $F_{M_1}$ is the space-energy-integrated fission neutron emission rate (n · s$^{-1}$), $\chi_{M_2}$ is the average fission spectrum and $E'$ is the energy of the emitted neutron (eV). For example, $M_1$ and $M_2$ can represent a fuel pin, a fuel assembly, the entire reactor core, etc.

Equations (2) and (3) imply that F is integrated over all incident neutron energies and every fissile isotope, as well as spatially integrated over the medium $M_1$. In the same manner, $\chi$ is averaged over the medium $M_2$ and is also considered to be integrated over all incident neutron energies. The two media $M_1$ and $M_2$ are introduced here because the spatial dependence of F is more pronounced than that of $\chi$. Therefore, it is common to use a spatial discretization with a larger mesh for $\chi$ than for F.

Since approaches α and β are performed in two steps, their fission neutron source terms follow the formalism of Equation (3). However, both methodologies make additional approximations when evaluating the spatial distribution of fission neutrons (F) and their energy distribution ($\chi$). These particularities are described in the following sections.
### 3.1.2. The $\alpha$ Approach: Use of the Data Provided in the Benchmark

In this first approach, the fission neutron source term is defined to be as close as possible to the data and recommendations given in [11], which leads to the following description of the cycle-averaged fission neutron source term:

$$ \overline{S_p(E')} = \overline{F_p} \cdot \overline{\chi_{core}(E')} \approx \overline{P_{th,p}} \cdot C_{f,bench} \cdot \chi_{bench}(E') \quad (4) $$

where, for a fuel pin $p$, $\overline{F_p}$ is the cycle-averaged, space-energy-integrated, fission neutron emission rate (n · s⁻¹), $\overline{\chi_{core}}$ is the cycle-averaged fission spectrum of the reactor core, $\overline{P_{th,p}}$ is the cycle-averaged pin power (MWth), $C_{f,bench}$ is the power-to-emitted-neutrons conversion factor (n · s⁻¹ · MW⁻¹), $\chi_{bench}(E')$ is the fission spectrum (eV⁻¹) and $E'$ is the energy of the emitted neutron (eV).

Equation (4) implies that the spatial distribution of fission neutrons is approximated by assuming that fission neutrons are distributed in the reactor in a power-like manner. This means that this approach makes the assumption that the same number of neutrons are emitted per fission in each fuel pin of the reactor. The impact of this hypothesis is therefore of increasing importance as heterogeneities between fuel-pin compositions increase (since ²³⁹Pu fission reactions emit on average more neutrons than ²³⁵U fission reactions). For cycle no. 9 of the HBR-2 reactor, this assumption is likely to underestimate the number of neutrons emitted in depleted assemblies (such as S08, S09, S10, etc.), which contain on average more ²³⁹Pu than other loaded assemblies. In contrast, this hypothesis implies that the fission neutron emission rates of fresh assemblies (such as R06, R07, R08, R09, etc.) are likely to be overestimated (for the opposite reason to the case outlined in the previous sentence). Additionally, the $\alpha$ approach assumes that the recoverable energy per fission is equal for every fission reaction occurring in the reactor. The impact of this hypothesis on the fission neutron distribution is expected to be less important than the previous one, since both ²³⁵U and ²³⁹Pu fission reactions emit energies of about 200 MeV (on average).

As explained in Section 2.2, [11] provides the cycle-averaged pin power distribution of the HBR-2 reactor. Axially, a single distribution was defined for the whole reactor, including 12 divisions distributed over the active height. This distribution is equal to the average of the axial distributions of the assemblies of the two outermost rows. Then, the three-dimensional pin powers were converted into fission neutron emission rates by multiplying them with the power-to-emitted-neutrons conversion factor $C_{f,bench}$. The value of this parameter considered in approach $\alpha$ was taken directly from [11], which provides a value of 8.175 · 10¹⁶ n · s⁻¹ · MW⁻¹. This value was evaluated by I. Remec and F.B.K Kam in 1997, using the DOTSOR code [22] and considering a burn-up equal to 28.596 GWd/t. This burn-up corresponds to the average exposure over cycle no. 9 of the assemblies labeled S08, S09 and S10.

Moreover, the $\alpha$ approach considers a single fission spectrum for the entire reactor core. This fission spectrum is equal to the average of the fission spectra of ²³⁵U and ²³⁹Pu, as recommended in [11]. The isotopic fission spectra are approximated here by Watt spectra, using the $\alpha$ and $\beta$ coefficients indicated in Table 2 (see Section 3.1.3).

<table><caption>**Table 2.** Data for fission neutron spectrum modeling.</caption><thead><tr><th>Nuclide</th><th>Incident<br>Neutron Energy</th><th>Fission<br>Spectrum Type</th><th>a (MeV)</th><th>b (MeV<sup>-1</sup>)</th></tr></thead><tbody><tr><td><sup>235</sup>U</td><td>Thermal</td><td>Watt</td><td>0.988</td><td>2.249</td></tr><tr><td><sup>238</sup>U</td><td>2.6 MeV</td><td>Watt</td><td>0.920</td><td>3.121</td></tr><tr><td><sup>239</sup>Pu</td><td>Thermal</td><td>Watt</td><td>0.966</td><td>2.842</td></tr><tr><td><sup>241</sup>Pu</td><td>Thermal</td><td>Maxwell</td><td>1.3597</td><td>-</td></tr></tbody></table>

Finally, the fission neutron source term was converted into a spatial and energetic source description for an MCNP6 [23] model in which the neutron transport was simulated from the core to the vessel (see Section 3.2). In practice, the spatial distribution of fission
neutrons given to MCNP6 is normalized, so that the sum of fission neutron emissions in the reactor core is equal to 1. Therefore, the MCNP6 tally outputs are given per source neutron and are post-processed following the methodology described in Section 3.3, to obtain the dosimeter specific activities.

### 3.1.3. The β Approach: Use of VACS

In this second approach, the fission neutron source term averaged over the 9th cycle of the HBR-2 reactor was computed using VACS. A detailed description of this methodology can be found in [6], and in this section we only recall the main principles of the calculation scheme. In addition, this section details the few approximations not present in the methodology as described in [6] that had to be made in order to carry out the analysis of the HBR-2 benchmark.

#### Main Principles of VACS

The data flow diagram in Figure 4 represents the different steps of VACS and shows the codes and data used. The calculation steps are sequential and are not automatic. Dedicated scripts are used to prepare parts of the MCNP6 [23] input file, based on the outputs of the SIMULATE5 [21] calculation.

In VACS, the cycle-averaged fission neutron source term is described by Equation (5).

$$ \overline{S_p(E')} = \overline{F_p} \cdot \overline{\chi_a(E')} \approx \sum_i \overline{R_{f,i,a}} \cdot \overline{v_i} \cdot \overline{X_p} \cdot \overline{\chi_a(E')} \quad (5) $$

where, for an isotope i and a fuel pin p within a fuel assembly a, $\overline{F_p}$ is the cycle-averaged, space-energy-integrated fission neutron emission rate (n · s⁻¹), $\overline{R_{f,i,a}}$ is the cycle-averaged, space-energy-integrated fission rate over assembly a (s⁻¹), $\overline{v_i}$ is the average number of neutrons emitted per fission (n), $\overline{\chi_a}$ is the average fission spectrum of fuel assembly a (eV⁻¹), $\overline{E'}$ is the energy of the emitted neutron (eV) and $\overline{X_p}$ is the pin power reconstruction factor.

Here, the spatial dependency of the fission neutron source term ($R_f \cdot v$ product) is evaluated by means of core-following calculations performed using a deterministic approach, coupling the CASMO5 [20] and SIMULATE5 codes. The CASMO5 code is used to generate the fuel-assembly homogenized 8-group macroscopic cross-section data. Then, the SIMULATE5 code uses these data as a multi-parametric database (since the energy-group cross sections depend on moderator density, boron concentration, etc.) to compute the full-core fission neutron source term.

On the CASMO5 side, the 2D stationary transport equation (the axial dimension being considered as infinite) is solved for assemblies in an infinite lattice, using a 586-group data library based on the ENDF/B-VII.1 nuclear data [24]. Firstly, the collision probability ($P_{ij}$) method [25] is used with 586 energy groups, on a coarse spatial mesh. Secondly, the method of characteristics (MOC) [25] is used on a refined spatial mesh with 19 energy groups for UOX assemblies and 35 energy groups for MOX assemblies. On the SIMULATE5 side, the diffusion equation is solved for 8 energy groups on 4 radial nodes and 20 axial divisions in each assembly. SIMULATE5 evaluates the total fission neutron emission rate distribution ($F_a(t) = \sum_i R_{f,i,a}(t) \cdot v_i$) in the reactor core, at an assembly scale, at multiple burn-up steps throughout the reactor cycle. Then, fuel-pin-level fission neutron emission rate distributions are reconstructed from the assembly level using the 3D form factors of the pin power reconstruction of the SIMULATE5 $X_p(t)$ [26]. These are multiplicative coefficients that specify how the power of an assembly is distributed among the different pins. In this study, 26 burn-up steps scattered along the cycle were modeled.
<figure><img src="image_7.png" /><figcaption>Figure 4. Flowchart of the different codes and steps used in VACS.</figcaption></figure>

Then, the pin-scale fission neutron emission distributions calculated by following the above methodology are used to compute the cycle-averaged fission neutron distribution within the reactor, according to Equation (6).

$$ \bar{X} = \frac{\sum_{n=0}^{N-1} X_n \cdot (bu_{n+1} - bu_n)}{bu_N} \qquad (6) $$

where, for a burn-up step *n* and a total number of burn-up steps *N*, *X* is the quantity to be averaged and *bu* is the reactor core exposure. Equation (6) means that the quantity *X* computed at the final burn-up step (*N*) is not used in the cycle-averaging process. In this study, *N* was equal to 26.

In the same manner, Equation (6) is used to compute the cycle-averaged, assembly-scale fission neutron emission rates $\bar{R}_{f,i,a}$ from the SIMULATE5 data. These are used to compute the cycle-averaged fission neutron spectrum of each fuel assembly $\bar{\chi}_{a}$, according to Equation (7).
$$ \overline{\chi_a(E')} = \sum_i \frac{\overline{R_{f,i,a}} \cdot v_i \cdot \chi_i(E')}{\sum_i \overline{R_{f,i,a}} \cdot v_i} \qquad (7) $$

where, for an isotope i and a fuel assembly a, $\overline{R_{f,i,a}}$ is the cycle-averaged, space-energy-integrated, fission rate (s⁻¹), $v_i$ is the average number of neutrons emitted per fission (n), $\chi_i$ is the fission spectrum (eV⁻¹) and $E'$ is the energy of the emitted neutron (eV).

In Equation (7), the fission neutron spectrum of each isotope is weighted by the probability that a neutron originates from the fission of a particular isotope. Here, only $^{235}$U, $^{238}$U, $^{239}$Pu and $^{241}$Pu are accounted for in the fission spectrum reconstruction. The fission spectra of these four isotopes are pre-evaluated on the basis of either a Watt spectrum or a Maxwell spectrum, using the a and b spectrum coefficients indicated in Table 2. The fission spectrum was compared against the Monte Carlo eigenvalue calculations in [7], and no significant discrepancies were found in the fast neutron fluence assessments at the RPV. It should be mentioned that the use of the Maxwell or Watt spectrum to model isotopic fission spectra is an approximation to the use of experimental spectra. This modeling point may be improved in future versions of VACS.

Finally, the fission neutron source term is automatically converted by a linking tool into a spatial and energetic source description for an MCNP6 model where the neutron transport is performed from the core to the vessel. The spatial representation is treated by using 20 axial divisions in each fuel pin (instead of the 200 used in the original VACS [6]). This approximation was made because it reduces the calculation time of the attenuation calculation (Section 3.2) by a factor of 10, from approximately 1,050,000 min.cpu to 105,000 min.cpu. The authors believe that this simplification of the modeling does not induce a significant bias in the assessments of the dosimeters' specific activities, since a mesh size of 20 divisions can already be considered as relatively fine. Regarding the energy, the fission neutron spectra of each assembly are described over 2000 energy bins distributed logarithmically between 10⁻⁴ MeV and 20 MeV.

In the same way as for approach a, the spatial distribution of fission neutrons provided to MCNP6 is, in practice, normalized so that the sum of fission neutron emissions in the reactor core is equal to 1. Therefore, the MCNP6 tally outputs are given *per source neutron* and are post-processed following the methodology described in Section 3.3, to obtain the dosimeter specific activities.

## Study-Specific Modeling Approximation

One asset of VACS is its ability to take into account the effects related to the history (power variations, control-rod insertions, etc.) of each reactor [6]. This capability is important, because it was shown in [8] that not taking into account the history of the reactor (i.e., modeling a reactor operating continuously at nominal power) can induce a significant deviation in the estimates of RPV aging. However, in order to accurately model the reactor core power variations, it is mandatory to have a detailed description of the control-rod composition, the insertion map of the various rod banks in the reactor and their insertion sequences.

In the case of the HBR-2 reactor, the authors could not find any public data about the insertion sequence of the control-rod banks as a function of the reactor power level. Therefore, the reactor was modeled as operating at full power with all control-rod banks fully extracted. Considering the study detailed in [8], this approximation could bias the dosimeter activity estimates by up to about 8%. However, this value must be considered with caution, as it was obtained with control-rod management that was probably very different from that used in the HBR-2 reactor. Indeed, the control-rod type and the insertion map of the different banks given in [8] (which correspond to those used in several French 900 MW$_e$ PWRs) significantly differ from the data found for the HBR-2 reactor in [17]. It is therefore only possible to say here that this assumption could lead to a non-negligible bias in the evaluations of the specific activities of the dosimeters.
Based on this assumption, the end of the cycle was automatically determined by SIMULATE5 by seeking for the point in time at which the critical boron concentration equaled 0 ppm. This process resulted in an estimated cycle length equal to 10.258 GWd/t. This value is 3.5% lower than the value given in the benchmark data [11], which seems reasonable given the approximations made regarding the reactor operation.

## 3.2. Modeling Fission Neutron Attenuation from the Core to Sites of Interest

The second step of the calculation deals with the transport of fission neutrons from the core to sites of interest (i.e., the dosimeters located at the monitoring program capsules or in the cavity), using the fission neutron distribution determined in the previous step. To perform this step, a MCNP6 fixed source calculation [27] was performed using a model developed on the basis of the "HBRHE81" MCNP5 input file, available on the SINBAD database, which was initially provided by P. Ortego (Shielding Engineering and Analysis, Spain).

The following sections describe the reactor model as provided in the "HBRHE81" input file and highlight the modifications made by the authors. The model developed in this process was used to perform fixed source calculations using the fission neutron source terms generated using the $\alpha$ and $\beta$ approaches. The MCNP6 fixed source calculation allows the evaluation of the per source neutron isotopic reaction rates in the dosimeters, which are then post-processed to obtain the specific activities.

### 3.2.1. Reactor Geometry

The reactor geometry described in the "HBRHE81" input file was used without modification. It describes one eighth of the HBR-2 reactor with reflection boundary conditions on the surfaces acting as symmetry cuts. Figure 5 displays a radial cut of the model in the equatorial plane (i.e., the axial position 0 cm). Similarly, Figure 6 presents an axial cut of the model in the 20° azimuthal plane (i.e., the plane in which the capsules are centered).

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Radial cut of the MCNP6 modeling of the H.B. Robinson-2 reactor.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>Figure 6. Axial cut of the MCNP6 modeling of the H.B. Robinson-2 reactor.</figcaption></figure>

The reactor core is described at a pin level, with assemblies described according to the geometry presented in Figure 1. (In Figures 5 and 6, the reactor core appears in black due to graphical artifacts resulting from the very large number of surfaces to be displayed in this area.) The inside of the guide tubes is filled with borated water, in order to model the reactor as operating with all control rods fully extracted. The boron concentration in the moderator here was adjusted to match 500 ppm (which is the recommended value in the benchmark data [11]). Moreover, in the "HBRHE81" input file, all the fuel rods are composed of fresh uranium oxide (UO₂) enriched to 2.9% by mass of ²³⁵U. Although the impact of this approximation on the dosimeter activity estimates is expected to be small (see [10]), it was decided to adjust this fuel composition in order to take into account the burn-up of the fuel in the isotopic vector. The fuel composition was taken here from [28], in which an ORIGEN-6.1 [29] evolution calculation was performed to compute the isotopic vector of the spent fuel for a burn-up equal to 27.4 GWd/t. This fuel exposure corresponds to the average of the cycle-averaged burn-up of assemblies S09 and S10. These assemblies were selected because they are located in the core periphery and are therefore strong contributors to the dosimeters' specific activities (especially in the cavity).

In the "HBRHE81" input file, the moderator density is homogeneously equal to 0.766 g·cm⁻³ in the reactor core. This density corresponds to a temperature of about 280 °C under the pressure conditions of cycle no. 9 (approximately 155 bar). In the same manner, the moderator density is homogeneously equal to 0.776 g·cm⁻³ (~273 °C) in the moderator bypass region (between the baffle and the barrel, see Figure 2) and to 0.787 g·cm⁻³ (~267 °C) in the downcomer regions (see Figure 2). These values correspond to those recommended in the benchmark [11] and have not been modified.

In the surveillance capsule, A533B steel plates (equivalent to the RPV steel) are placed at the -28 cm, 0 cm and +28 cm axial positions (three plates per axial location). These plates are coated with a thin layer of gadolinium and mounted on a stainless steel (SS-304) support (see Figure 5). Each steel plate is 15 cm high and approximately 1 cm thick. In the "HBRHE81" input file, the total width of the surveillance capsule (including its support) is approximately equal to 5°. This value is not consistent with the data provided in [11], where a total width of 2.13° is specified. Nevertheless, it seems reasonable to consider here that this modeling discrepancy has only a small impact on the assessment of the dosimeter specific activities compared to the overall uncertainty of such assessments. As an
illustration, the ±1σ uncertainty in the estimates of the HBR-2 dosimeter specific activities was quantified as equal to approximately 15% in [13]. Therefore, the monitoring capsule modeling was not modified when performing the validation study. For the geometry cell associated with the cavity dosimeter, the modeling provided in the "HBRHE81" input file was not modified. The cavity is located just behind the RPV insulation and contains air (see Figure 5). The cell measures 10 cm in width and 76 cm in height, with a thickness of 1.1 cm. The air volume is surrounded by a thin layer of cadmium, which is also surrounded by a thin layer of aluminum (each about 0.05 cm thick).

Finally, all the nuclear data used to model neutron transport in the "HBRHE81" input file were considered at 300 K. In the present study, the MCNP input file was modified so that the cross sections used for the transport calculation were considered:

* At 900 K in the fuel;

* At 600 K in the moderator;

* At 600 K in all regions between the reactor core and the vessel (including the RPV);

* At 300 K in all regions outside the vessel.

For the α approach (Section 3.1.2), only the ENDF/B-VII.1 library [24] was used to model neutron attenuation. For the β approach (Section 3.1.2) the following nuclear data libraries were used: ENDF/B-VII.1, ENDF/B-VIII.0 [30] and JEFF-3.3 [19]. All cross sections used were processed using the NJOY2016 code [18].

### 3.2.2. Modeling of the Dosimeters

The dosimeters located in the surveillance capsule and in the cavity are not "physically" modeled in the MCNP6 model. Therefore, the isotopic reaction rates (which are used to compute the specific activities) were evaluated in the central steel plates of the surveillance capsules as well as in the geometry cell corresponding to the cavity (see Figure 5). The two considered cells were axially located in the equatorial plane (Z = 0 cm), since the experimental activities (see Table 1) are provided only for this axial plane. The isotopic reaction rates are evaluated from the average neutron fluxes in these cells. For this purpose, the "F4" MCNP6 cell flux tallies [27] are combined with the "FM" cards [27], which are used to multiply the evaluated neutron flux by the cross sections of the reactions of interest. The dosimeter composition is not explicitly given in the SINBAD database. However, a fictitious material (the term "fictitious material" refers here to a material which is not used in the definition of the geometry) representative of the dosimeter composition was defined in the "HBRHE81" input. This material can therefore be associated with the "FM" cards in order to calculate the isotopic reaction rates. The composition of this material is shown in Table 3. However, it is important to note that this dosimeter composition could not be found elsewhere in the benchmark data. Regarding the dosimetry cross sections, the following libraries were tested to calculate the reaction rates:

* The nuclear data library used to model the neutron transport (ENDF/B-VII.1, ENDF/B-VIII.0 or JEFF-3.3);

* The IRDFF-II (International Reactor Dosimetry and Fusion File) library [31].

<table><caption>**Table 3.** Dosimeter composition as defined in the "HBRHE81" input file [11].</caption><thead><tr><th>Isotope</th><th>Atomic Density (Atom · Barn<sup>-1</sup> · cm<sup>-1</sup>)</th></tr></thead><tbody><tr><td><sup>54</sup>Fe</td><td>5.13·10<sup>-3</sup></td></tr><tr><td><sup>63</sup>Cu</td><td>5.88·10<sup>-2</sup></td></tr><tr><td><sup>58</sup>Ni</td><td>6.26·10<sup>-2</sup></td></tr><tr><td><sup>46</sup>Ti</td><td>4.707·10<sup>-3</sup></td></tr><tr><td><sup>237</sup>Np</td><td>5.081·10<sup>-2</sup></td></tr><tr><td><sup>238</sup>U</td><td>4.75·10<sup>-2</sup></td></tr></tbody></table>
### 3.2.3. Variance Reduction and Simulation Parameters

To improve the calculation efficiency of the MCNP6 fixed source calculation, the weight-window variance reduction method [27] was used. The weight-window parameters were optimized for 47 neutron energy bins covering the range [10⁻⁵ eV-17.3 MeV]. This evaluation was performed using the ADVANTG-3.2 tool [32], which uses the deterministic code DENOVO [33]. A set of optimized weight windows for both dosimeter locations (capsule and cavity) and for all reactions was generated using the forward-weighted CADIS method [34]. Finally, an MCNP6 attenuation calculation was performed for each of the nuclear data library combinations in the studied approach (e.g., the β approach with the ENDF/B-VII.1 library). In each fixed source calculation, 8·10⁹ neutron histories were simulated, leading to a statistical uncertainty σ of about 0.2-2% in the reaction rate estimates.

### 3.3. Post-Processing of MCNP6 Tally Outputs

As explained previously, the fission neutron source term provided and used in the MCNP6 attenuation calculation was normalized so that the sum of fission neutron emissions in the reactor core was equal to 1. Therefore, the MCNP6 calculation provides an estimate of the dosimeter reaction rates $R_{c,k}$ in units of reaction·source neutron⁻¹·atom⁻¹. In order to obtain $R_{c,k}$ in units of reaction·s⁻¹·atom⁻¹, the results for the MCNP6 tallies must be post-processed according to Equation (8).

$$R_{c,k} = \frac{\overline{P_{th,core}} \cdot C_f \cdot T_{c,k}}{8} \qquad (8)$$

where, for a reaction c and a dosimeter location k, $T_{c,k}$ is the MCNP6 tally output (reaction·neutron⁻¹·atom⁻¹), $\overline{P_{th,core}}$ is the average thermal power of the reactor core in cycle no. 9 (which is approximately equal to 1730 MWth) and $C_f$ is the power-to-emitted-neutrons conversion factor (neutrons·s⁻¹·MW⁻¹). The factor 1/8 is added to account for the fact that only one eighth of the reactor is modeled in the MCNP6 attenuation calculation.

On the one hand, for the α approach (Section 3.1.2), the $C_f$ coefficient value was taken from the study performed by I. Remec and F.B.K Kam [11], in which a value of 8.175·10¹⁶ neutrons·s⁻¹·MW⁻¹ is specified. As explained before, this value was evaluated using the DOTSOR code [22] for a burn-up equal to 28.596 GWj/t, which corresponds to the average exposure of the S08, S09 and S10 assemblies in cycle no. 9. On the other hand, for the β approach (Section 3.1.3), the $C_f$ coefficient was evaluated according to Equation (9), using the cycle-averaged data calculated by SIMULATE5.

$$C_f = \frac{\sum_a \sum_i \overline{R_{f,i,a}} \cdot v_i}{\overline{P_{th,S5}}} \qquad (9)$$

where, for an isotope i and an assembly a, $\overline{R_{f,i,a}} \cdot v_i$ is the cycle-averaged fission neutron production (neutrons·s⁻¹) and $\overline{P_{th,S5}}$ is the average thermal power of the reactor over the cycle modeled with SIMULATE5 (2300 MWth = reactor nominal power).

Equation (9) led to an estimate of $C_f$ equal to 7.93302·10¹⁶ neutrons·s⁻¹·MW⁻¹. This value is slightly lower than the one given by I. Remec and F.B.K Kam. This observation is consistent with the observations made in [13], suggesting that the estimate of $C_f$ evaluated by I. Remec and F.B.K Kam is an overestimate of the actual $C_f$ value. Finally, the $C_f$ coefficient evaluated using Equation (9) was used to post-process the results of the MCNP6 attenuation calculation for the β approach.

Once the reaction rates $R_{c,k}$ are evaluated, the dosimeter specific activities $A_k$ (in units of Bq·mg⁻¹) are deduced from $R_{c,k}$, according to Equation (10).
$$A_{c,k} = f_{c,k} \cdot R_{c,k} \qquad (10)$$

where, for a reaction c and a dosimeter location k, $f_{c,k}$ is the reaction rate for a specific activity conversion factor (atom·mg⁻¹) and $R_{c,k}$ is the reaction rate (reaction·s⁻¹·atom⁻¹). The values of $f_{c,k}$ used in this study were taken from [15] and are summarized in Table 4.

<table><caption>**Table 4.** Reaction rates for a specific activity conversion factor, taken from [15].</caption><thead><tr><th colspan="4">Dosimeter Location</th></tr><tr><th></th><th>Reaction</th><th>Surveillance Capsule</th><th>Reactor Cavity</th></tr></thead><tbody><tr><td rowspan="6">f<sub>c,k</sub><br>(atom·mg<sup>-1</sup>)</td><td><sup>237</sup>Np(n,f)<sup>137</sup>Cs</td><td>3.082·10<sup>15</sup></td><td>3.085·10<sup>15</sup></td></tr><tr><td><sup>238</sup>U(n,f)<sup>137</sup>Cs</td><td>2.0950·10<sup>15</sup></td><td>2.952·10<sup>15</sup></td></tr><tr><td><sup>58</sup>Ni(n,p)<sup>58</sup>Co</td><td>3.576·10<sup>18</sup></td><td>3.844·10<sup>18</sup></td></tr><tr><td><sup>54</sup>Fe(n,p)<sup>54</sup>Mn</td><td>2.494·10<sup>17</sup></td><td>2.558·10<sup>17</sup></td></tr><tr><td><sup>46</sup>Ti(n,p)<sup>46</sup>Sc</td><td>5.364·10<sup>17</sup></td><td>5.737·10<sup>17</sup></td></tr><tr><td><sup>63</sup>Cu(n,a)<sup>60</sup>Co</td><td>6.697·10<sup>17</sup></td><td>6.726·10<sup>17</sup></td></tr></tbody></table>

Finally, the calculated specific activities were directly compared to the experimental data using the C/E ratio, defined by Equation (11).

$$C/E = \frac{A_{calc}}{A_{exp}} \qquad (11)$$

where $A_{calc}$ is the activity calculated with the α or β approaches (Bq·mg⁻¹) and $A_{exp}$ is the experimental specific activity provided in the benchmark data (Bq·mg⁻¹, see Table 1).

The C/E ratios are the quantities that measure the accuracy of the calculation scheme. C/E ratios close to 1 indicate accurate modeling. The further away from 1 the C/E ratios are, the more the modeling quality can be questioned.

# 4. Results

## 4.1. Spatial Distribution of Fission Neutrons

First, the spatial distribution of fission neutrons obtained with the β method was compared to that reconstructed from the data provided in the benchmark (α approach). Figure 7 shows the relative differences between the two approaches.

At the assembly scale, the observed differences are mostly below 3% (in absolute value) in the center of the core and tend to increase up to 8.76% for peripheral assemblies. A significant asymmetry in the relative differences can be observed between the left and right sides of the core. This behavior originates from the cycle-average powers provided by [11], which are asymmetrical between the left and right sides of the core. This observation suggests that there is a power tilt in the reactor core, which cannot therefore be reproduced with the β approach. At the pin scale, the differences are relatively homogeneous within the assemblies located in the core center, with some local increases (in absolute value) at the junctions between two assemblies with significantly different burn-up values (thus with differing <sup>239</sup>Pu contents). In contrast, a strong increase in the relative differences is observed for the outermost rows of pins, for which the β approach systematically overestimates the fission neutron production rates compared to the α approach. These overestimates reach 200% for some fuel pins (not shown in Figure 7 for viewing convenience) and can be explained by the fact that the pin powers provided in [11] for these fuel pins are very low.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Relative differences (%) between the cycle-averaged fission neutron distributions evaluated with the α and β approaches (β/α). Differences are presented at the assembly scale (top) and at the fuel-pin scale (bottom). For pin-scale differences, only values less than or equal to 20% are displayed, for viewing convenience.</figcaption></figure>
It is also interesting to note that the $\beta$ approach overestimates the fission neutron production rates in the S08, R08, R09, S09, R10 and S10 assemblies, which are the main contributors to the neutron flux in the cavity. In contrast, the fission neutron production rates in the other assemblies contributing to the neutron flux in the capsules (e.g., R12, P13, N14, etc.) tend to be underestimated by the $\beta$ approach compared to the $\alpha$ approach. Considering these results, it is likely that the $\beta$ approach will lead to higher specific activity predictions in the cavity than the $\alpha$ approach. Moreover, only a small variation in the predictions of the capsule specific activities is expected when switching from the $\alpha$ to the $\beta$ approach, due to compensation effects (e.g., the fission neutron production rate of assembly S10 is overestimated while that of assembly R12 is underestimated).

Overall, it is likely that most of the differences explained above result from model differences between the $\alpha$ approach and the $\beta$ approach, as follows:

* The methodology employed to evaluate the fission neutron source term in the $\beta$ approach uses the SIMULATE5 code which is based on the diffusion approximation. In comparison, the cycle-averaged power distribution provided in [11] was evaluated from burn-up distributions that were probably obtained by combining experimental data and neutronics calculations.

* The fission neutrons are distributed in the same way as the thermal power in the $\alpha$ approach. In contrast, the $\beta$ approach computes the spatial distribution of source neutrons using the fission neutron production rates (product $R_f \cdot v$), which is a more accurate approach in the context of a reactor containing assemblies with various exposures.

* In the $\beta$ approach, the reactor core was modeled as operating at nominal power. However, the power history of the reactor was probably taken into account in the $\alpha$ approach.

* Assemblies R09, R10 and R11 are among those that were likely to contain gadolinium-bearing fuel rods, which were not modeled in the $\beta$ approach. Nevertheless, the trends observed here are consistent with the fact that these assemblies could contain gadolinia-bearing fuel rods.

* The nuclear data library used in the $\beta$ approach to evaluate the fission neutron source term (based on ENDF/B-VII.1) probably differs from the one used in the $\alpha$ approach.

At this point in the paper, it is not possible to discuss the quality of the two fission neutron source terms. Such an analysis will be performed in Section 4.3, which discusses the C/E obtained with both approaches.

## 4.2. Dosimeter Reaction Rates

Reaction rates were calculated using the average fission neutron distribution in cycle no. 9, evaluated using the $\alpha$ or $\beta$ approaches. Post-processing of the MCNP6 tally outputs was performed as described in Section 3.3. The obtained results are presented in Table 5 for the capsules and in Table 6 for the cavity. In order to have a point of comparison, these tables also summarize the reaction rates calculated in 1997 by I. Remec and F.B.K Kam [11], using the two-dimensional transport code DORT [35] and the nuclear data library BUGLE-96 [36]. Putting Tables 5 and 6 side by side, it can be seen that the standard deviation of the MCNP6 tallies are systematically lower in the cavity than in the capsule. This behavior is due to the parameterization of the weight windows generated by the ADVANTG-3.2 code, which was more successful for neutron histories heading toward the cavity than for those heading toward the surveillance capsule. In general, the obtained standard deviations for the reaction rate estimates were inferior or equal to 1%. However, they were slightly more pronounced for the $^{46}$Ti(n,p) and $^{63}$Cu(n,$\alpha$) reactions, for which they reached 1.65% and 2.21%, respectively (in the cavity). This behavior is induced by the fact that these reactions have an energy threshold of the order of 4 MeV (see Figure 3). Indeed, the number of high-energy neutrons in the surveillance capsule and in the cavity is rather low compared to the population of low-energy neutrons.
<table><caption>**Table 5.** Reaction rates obtained in the surveillance capsule. Values in parentheses indicate the standard deviations of MCNP6 tally responses (in units of %).</caption><thead><tr><th rowspan="2">Methodology</th><th rowspan="2">Nuclear Data—<br>Attenuation</th><th rowspan="2">Nuclear Data—<br>Dosimetry</th><th colspan="6">Surveillance Capsule Reaction Rates (Reactions ·s<sup>-1</sup>·atom<sup>-1</sup>)</th></tr><tr><th><sup>237</sup>Np(n,f)</th><th><sup>238</sup>U(n,f)</th><th><sup>58</sup>Ni(n,p)</th><th><sup>54</sup>Fe(n,p)</th><th><sup>46</sup>Ti(n,p)</th><th><sup>63</sup>Cu(n,a)</th></tr></thead><tbody><tr><td>Remec and Kam-DORT [11]</td><td>BUGLE-96</td><td></td><td>1.06·10<sup>-13</sup></td><td>1.54·10<sup>-14</sup></td><td>4.74·10<sup>-15</sup></td><td>3.51·10<sup>-15</sup></td><td>5.62·10<sup>-16</sup></td><td>3.57·10<sup>-17</sup></td></tr><tr><td rowspan="2">α approach and MCNP6</td><td rowspan="2">ENDF/B-VII.1</td><td>ENDF/B-VII.1</td><td>1.07·10<sup>-13</sup><br>(0.55)</td><td>1.37·10<sup>-14</sup><br>(0.68)</td><td>4.12·10<sup>-15</sup><br>(0.87)</td><td>3.04·10<sup>-15</sup><br>(0.94)</td><td>5.04·10<sup>-16</sup><br>(1.41)</td><td>2.96·10<sup>-17</sup><br>(1.85)</td></tr><tr><td>IRDFF-II</td><td>1.07·10<sup>-13</sup><br>(0.55)</td><td>1.39·10<sup>-14</sup><br>(0.67)</td><td>4.15·10<sup>-15</sup><br>(0.90)</td><td>3.00·10<sup>-15</sup><br>(0.97)</td><td>5.47·10<sup>-16</sup><br>(1.40)</td><td>3.02·10<sup>-17</sup><br>(1.84)</td></tr><tr><td rowspan="6">β approach and MCNP6</td><td rowspan="3">ENDF/B-VII.1</td><td>ENDF/B-VII.1</td><td>1.04·10<sup>-13</sup><br>(0.58)</td><td>1.41·10<sup>-14</sup><br>(0.72)</td><td>4.17·10<sup>-15</sup><br>(0.93)</td><td>3.06·10<sup>-15</sup><br>(0.99)</td><td>5.07·10<sup>-16</sup><br>(1.54)</td><td>3.06·10<sup>-17</sup><br>(2.05)</td></tr><tr><td>IRDFF-II</td><td>1.11·10<sup>-13</sup><br>(0.58)</td><td>1.43·10<sup>-14</sup><br>(0.72)</td><td>4.19·10<sup>-15</sup><br>(0.95)</td><td>3.01·10<sup>-15</sup><br>(1.02)</td><td>5.49·10<sup>-16</sup><br>(1.53)</td><td>3.19·10<sup>-17</sup><br>(2.05)</td></tr><tr><td>ENDF/B-VIII.0</td><td>9.97·10<sup>-14</sup><br>(0.61)</td><td>1.22·10<sup>-14</sup><br>(0.76)</td><td>3.45·10<sup>-15</sup><br>(1.03)</td><td>2.51·10<sup>-15</sup><br>(1.09)</td><td>4.27·10<sup>-16</sup><br>(1.65)</td><td>3.84·10<sup>-17</sup><br>(2.21)</td></tr><tr><td rowspan="2">ENDF/B-VIII.0</td><td>IRDFF-II</td><td>1.00·10<sup>-13</sup><br>(0.62)</td><td>1.22·10<sup>-14</sup><br>(0.76)</td><td>3.50·10<sup>-15</sup><br>(1.02)</td><td>2.51·10<sup>-15</sup><br>(1.09)</td><td>4.63·10<sup>-16</sup><br>(1.64)</td><td>2.65·10<sup>-17</sup><br>(2.18)</td></tr><tr><td>JEFF-3.3</td><td>1.13·10<sup>-13</sup><br>(0.57)</td><td>1.44·10<sup>-14</sup><br>(0.71)</td><td>4.26·10<sup>-15</sup><br>(0.92)</td><td>2.74·10<sup>-15</sup><br>(1.02)</td><td>5.29·10<sup>-16</sup><br>(1.53)</td><td>2.89·10<sup>-17</sup><br>(2.06)</td></tr><tr><td rowspan="2">JEFF-3.3</td><td>IRDFF-II</td><td>1.13·10<sup>-13</sup><br>(0.58)</td><td>1.47·10<sup>-14</sup><br>(0.71)</td><td>4.19·10<sup>-15</sup><br>(0.94)</td><td>3.00·10<sup>-15</sup><br>(1.02)</td><td>5.30·10<sup>-16</sup><br>(1.53)</td><td>2.96·10<sup>-17</sup><br>(2.03)</td></tr></tbody></table>

<table><caption>**Table 6.** Reaction rates obtained in the reactor cavity. Values in parentheses indicate the standard deviations of MCNP6 tally responses (in units of %).</caption><thead><tr><th rowspan="2">Methodology</th><th rowspan="2">Nuclear Data—<br>Attenuation</th><th rowspan="2">Nuclear Data—<br>Dosimetry</th><th colspan="6">Cavity reaction rates (reactions ·s<sup>-1</sup>·atom<sup>-1</sup>)</th></tr><tr><th><sup>237</sup>Np(n,f)</th><th><sup>238</sup>U(n,f)</th><th><sup>58</sup>Ni(n,p)</th><th><sup>54</sup>Fe(n,p)</th><th><sup>46</sup>Ti(n,p)</th><th><sup>63</sup>Cu(n,a)</th></tr></thead><tbody><tr><td>Remec and Kam-DORT [11]</td><td>BUGLE-96</td><td></td><td>4.13·10<sup>-15</sup></td><td>2.11·10<sup>-16</sup></td><td>4.79·10<sup>-17</sup></td><td>3.23·10<sup>-17</sup></td><td>5.16·10<sup>-18</sup></td><td>3.63·10<sup>-19</sup></td></tr><tr><td rowspan="2">α approach and MCNP6</td><td rowspan="2">ENDF/B-VII.1</td><td>ENDF/B-VII.1</td><td>5.28·10<sup>-15</sup><br>(0.21)</td><td>2.28·10<sup>-16</sup><br>(0.39)</td><td>5.00·10<sup>-17</sup><br>(0.55)</td><td>3.32·10<sup>-17</sup><br>(0.64)</td><td>5.560·10<sup>-18</sup><br>(1.04)</td><td>3.63·10<sup>-19</sup><br>(1.34)</td></tr><tr><td>IRDFF-II</td><td>5.32·10<sup>-15</sup><br>(0.21)</td><td>2.34·10<sup>-16</sup><br>(0.39)</td><td>4.80·10<sup>-17</sup><br>(0.59)</td><td>3.25·10<sup>-17</sup><br>(0.67)</td><td>6.00·10<sup>-18</sup><br>(1.04)</td><td>3.69·10<sup>-19</sup><br>(1.33)</td></tr><tr><td rowspan="6">β approach and MCNP6</td><td rowspan="3">ENDF/B-VII.1</td><td>ENDF/B-VII.1</td><td>5.65·10<sup>-15</sup><br>(0.22)</td><td>2.41·10<sup>-16</sup><br>(0.41)</td><td>5.19·10<sup>-17</sup><br>(0.58)</td><td>3.43·10<sup>-17</sup><br>(0.68)</td><td>5.63·10<sup>-18</sup><br>(1.10)</td><td>3.67·10<sup>-19</sup><br>(1.40)</td></tr><tr><td>IRDFF-II</td><td>5.70·10<sup>-15</sup><br>(0.22)</td><td>2.48·10<sup>-16</sup><br>(0.21)</td><td>4.97·10<sup>-17</sup><br>(0.62)</td><td>3.34·10<sup>-17</sup><br>(0.71)</td><td>6.08·10<sup>-18</sup><br>(1.09)</td><td>3.74·10<sup>-19</sup><br>(1.40)</td></tr><tr><td>ENDF/B-VIII.0</td><td>4.75·10<sup>-15</sup><br>(0.23)</td><td>1.92·10<sup>-16</sup><br>(0.46)</td><td>3.61·10<sup>-17</sup><br>(0.72)</td><td>2.48·10<sup>-17</sup><br>(0.81)</td><td>4.43·10<sup>-18</sup><br>(1.27)</td><td>4.45·10<sup>-19</sup><br>(1.59)</td></tr><tr><td rowspan="2">ENDF/B-VIII.0</td><td>IRDFF-II</td><td>4.78·10<sup>-15</sup><br>(0.23)</td><td>1.92·10<sup>-16</sup><br>(0.46)</td><td>3.71·10<sup>-17</sup><br>(0.70)</td><td>2.48·10<sup>-17</sup><br>(0.81)</td><td>4.76·10<sup>-18</sup><br>(1.27)</td><td>3.11·10<sup>-19</sup><br>(1.58)</td></tr><tr><td>JEFF-3.3</td><td>5.85·10<sup>-15</sup><br>(0.22)</td><td>2.63·10<sup>-16</sup><br>(0.40)</td><td>5.57·10<sup>-17</sup><br>(0.55)</td><td>3.22·10<sup>-17</sup><br>(0.67)</td><td>6.00·10<sup>-18</sup><br>(1.11)</td><td>3.64·10<sup>-19</sup><br>(1.45)</td></tr><tr><td rowspan="2">JEFF-3.3</td><td>IRDFF-II</td><td>5.89·10<sup>-15</sup><br>(0.22)</td><td>2.69·10<sup>-16</sup><br>(0.40)</td><td>5.25·10<sup>-17</sup><br>(0.59)</td><td>3.49·10<sup>-17</sup><br>(0.68)</td><td>6.01·10<sup>-18</sup><br>(1.11)</td><td>3.71·10<sup>-19</sup><br>(1.43)</td></tr></tbody></table>

Table 6 shows that the reaction rates obtained in the monitoring capsules with the β approach are globally similar to those obtained with the α approach (with equivalent nuclear data libraries). In contrast, Table 6 highlights that the reaction rates obtained with the β approach in the cavity are systematically higher than those obtained with the α approach. These results are consistent with what was expected from analyzing the differences between the fission neutron source terms of the two approaches (see Section 4.1). Moreover, the results obtained with the ENDF/B-VII.1 and JEFF-3.3 libraries are globally consistent. In contrast, the reaction rates obtained with ENDF/B-VIII.0 are largely lower
than those obtained with the other libraries. This trend is not improved by using the IRDFF-II library to evaluate reaction rates. These observations therefore indicate that the modeling of neutron transport from the core to the sites of interest is significantly different with the ENDF/B-VIII.0 library compared to the other libraries. These tendencies were also observed in a more extended study of ex-core dosimeter validations [37]. Finally, the reaction rates obtained in the surveillance capsule with the ENDF/B-VII.1 and JEFF-3.3 libraries are, on average, slightly lower than those obtained by I. Remec and F.B.K Kam. Conversely, the reaction rates evaluated with the ENDF/B-VII.1 and JEFF-3.3 libraries are slightly higher than those calculated by I. Remec and F.B.K Kam. Given that the specific activity predictions of I. Remec and F.B.K Kam tend to be inferior to the experimental values [11], an improvement in the C/E in the cavity and a slight degradation in the C/E in the surveillance capsule can be expected.

### 4.3. Dosimeter Specific Activities

The specific activities of dosimeters and the associated C/E ratios were calculated from the reaction rates, as described in Section 3.3. Figure 8 presents the results obtained using the α and β approaches with the ENDF/B-VII.1 nuclear data library. In the same way as before, the results obtained by I. Remec and F.B.K Kam are shown here in order to provide a point of comparison.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** C/E ratios obtained with the α and β approaches using the ENDF/B-VII.1 nuclear data library. The results obtained by I. Remec and F.B.K. Kam [11] are shown for comparison. The $^{237}$Np(n,f) specific activity is not taken into account in the calculation of the average C/E value in the cavity.</figcaption></figure>

In Figure 8 (and in those following), we chose not to represent the statistical uncertainties of the C/E values, as these were not significant. Only the standard deviation of the average C/E is shown, which is significantly greater than the statistical uncertainty of the Monte Carlo calculations. Nevertheless, it is still worth mentioning that the uncertainty propagation work done by R. Pittarello et al. in [13] showed that, in the case of the HBR-2 benchmark, the ±1σ uncertainty in the specific activity assessments can be of the order of 15%.
For the surveillance capsule, the C/E ratios obtained with the α and β approaches are relatively similar and remain above 0.8 overall. The average C/E value is, in both cases, approximately equal to 0.84. The left side of Figure 8 also shows that the C/E values obtained for the surveillance capsule are globally worse than those obtained by I. Remec and F.B.K Kam. For the reactor cavity, the C/E values obtained with the α and β approaches are mostly higher than 0.9, except for the ²³⁷Np(n,f) reaction, where the C/E coefficient is significantly lower. Similar observations were made in [13–15] regarding the ²³⁷Np(n,f) reaction in the cavity. Therefore, it is reasonable to consider here that there is a measurement error in this experimental value. Thus, the specific activity of ²³⁷Np(n,f) was not taken into account in the calculation of the average C/E value in the cavity (in Figure 8 and in those following). Considering these results, it seems that the α and β approaches improve the specific activity predictions in the cavity compared to those in the study of I. Remec and F.B.K Kam. The average C/E ratios obtained in the cavity with the α and β approaches are about 0.94 and 0.98, respectively. Therefore, it also seems that the β model is slightly more accurate than the α model. Nevertheless, it remains important to point out that the discussion on the quality of the various models must be tempered, given the uncertainties inherent in such C/E evaluations (±1σ = 15% according to [13]) and the standard deviation of the calculated average C/E values.

Overall, the results obtained with VACS (the β approach) are consistent with those of [13], where it was demonstrated that a fine analysis of the HBR-2 benchmark using the MCNP code induced, compared to the work of I. Remec and F.B.K Kam, a degradation of the activation predictions in the capsules and an improvement of those in the cavity. However, as mentioned in [13], although the C/E values obtained by I. Remec and F.B.K Kam are closer to the experimental values in the cavity [11], the quality of these results is at least partially fortuitous, i.e., it is not guaranteed that such good results would be obtained for another configuration.

In order to put into perspective the C/E ratios obtained with the β approach and the ENDF/B-VII.1 library, Figure 9 presents these results together with those obtained with the β approach and the ENDF/B-VIII.0 and JEFF-3.3 libraries. This figure shows a good consistency between the results obtained with the JEFF-3.3 library and those obtained with the ENDF/B-VII.1 library. In contrast, the C/E ratios evaluated with the ENDF/B-VIII.0 library are largely degraded in comparison with the other tested nuclear data libraries. It is also possible to note that the C/E ratio obtained with ENDF/B-VIII.0 for the ⁶³Cu(n,α) specific activity is very different from those evaluated for the other specific activities. This observation suggests that there is also an anomaly in the ENDF/B-VIII.0 ⁶³Cu(n,α) cross section. In order to confirm this last statement, Figure 10 presents the obtained C/E values using the ENDF/B-VIII.0 library to model neutron transport and the IRDFF-II library for dosimetry.

This figure shows that using the IRDFF-II dosimeter cross sections does not signifi-
cantly improve the C/E ratio obtained, when using the ENDF/B-VIII.0 library to model
neutron attenuation from the core to the sites of interest. The only significant difference
observed here is for the ⁶³Cu(n,α) reaction, whose C/E value regains consistency with the
other reactions using the IRDFF-II library (which was not the case with ENDF/B-VIII.0).
Therefore, these results suggest that the ENDF/B-VIII.0 library does not allow accurate
modeling of the neutron transport from the core to the sites of interest (the capsule and the
cavity) to be performed and that an anomaly exists in the ⁶³Cu(n,α) cross section in this
library. In contrast, the C/E ratios obtained with the ENDF/B-VII.1 and JEFF-3.3 libraries
are consistent with those obtained using IRDFF-II for the dosimetry step, as shown in
Figure 11.
<figure><img src="image_12.png" /><figcaption>**Figure 9.** C/E ratios obtained with the *β* approach using the ENDF/B-VII.1, ENDF/B-VIII.0 and JEFF-3.3 nuclear data libraries. The ²³⁷Np(n,f) specific activity is not taken into account in the calculation of the average C/E value in the cavity.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 10.** C/E ratios obtained with the *β* approach using the ENDF/B-VIII.0 library for transport and the ENDF/B-VIII.0 or IRDFF-II library for dosimetry. The ²³⁷Np(n,f) specific activity is not taken into account in the calculation of the average C/E value in the cavity.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Figure 11.** The C/E ratios obtained with the *β* approach and various combinations of nuclear data libraries. The <sup>237</sup>Np(n,f) specific activity is not taken into account in the calculation of the average C/E value in the cavity.</figcaption></figure>

In summary, the predictions of the specific activities of the dosimeters obtained using VACS (β approach) can be considered satisfactory, considering the uncertainties in the C/E assessments (about 15% according to [13]). Moreover, the consistency between the results obtained here and those of [13–15] also reinforces the confidence in the vessel aging predictions using VACS. Nevertheless, the correction of some assumptions made in this study could result in an improvement in the C/E predictions:

* In the *β* approach, the fission neutron source term was evaluated considering a reactor operating at nominal power. However, it was shown in [8] that predictions of vessel aging are sensitive to the history of the reactor. Therefore, consideration of the reactor power history (and associated control-rod insertions) in the evaluation of the fission neutron source term could improve the C/E predictions.

* As explained in Section 2.2, it is likely that gadolinia-bearing fuel pins were used during cycle no. 9, and these were not modeled in the *β* approach. Modeling these gadolinia-bearing fuel pins could result in an improvement in the C/E predictions.

* In the MCNP6 model, the moderator density gradients in the reactor core and in the moderator bypass zones were not modeled. However, in [6], the authors indicate that such assumptions can lead to a significant bias in vessel aging predictions. Therefore, modeling the moderator density gradients could improve the accuracy of C/E estimates.

* The surveillance capsule modeling could be reworked to better match the data provided in [11], which could also lead to a very slight improvement in the obtained results.

It is worth mentioning that these approximations are not usually made in VACS [6,9]. In particular, the reactor history is usually taken into account in the fission neutron source term assessment. Moreover, a fine description of the moderator density gradients in each assembly, as calculated by SIMULATE5, is usually added in the MCNP6 fixed source calculation, and the moderator density in the bypass is usually fixed as the average of those
of the outermost row of assemblies. This description of the moderator density gradient was
added after the publication of [6], and the impact of this change on aging assessments can
be found in [9].

**5. Conclusions**

Reactor lifetime extension motivates the development of precise and accurate simulation tools for vessel aging assessments and in-depth analysis of the parameters that may influence this aging (fuel management variability, core power history, etc.). Recently, a calculation scheme dedicated to the analysis of PWR vessel aging was developed, named VACS (vessel aging calculation scheme). VACS couples a deterministic approach (CASMO5 and SIMULATE5 codes) to evaluate the full-core fission neutron source term with Monte Carlo modeling (MCNP6 code) to perform the neutron attenuation from the core to sites of interest. This article featured an in-depth analysis of the H.B. Robinson-2 reactor pressure vessel dosimetry benchmark, constituting a first step towards a complete validation of VACS. The relative qualities of several nuclear data libraries (ENDF-B/VII.1, ENDF-B/VIII.0 and JEFF-3.3) for modeling the neutron attenuation from the core to the capsules and cavity were also studied.

The H.B. Robinson-2 dosimetry benchmark was intended to be used to validate the attenuation of the neutron flux only, which corresponds to the second step of the VACS approach. Nevertheless, information was gathered to use this benchmark to validate both steps, including the fission distribution inside the core. Thus, this work provides some recommendations regarding information that is still missing for performing a precise validation of the core calculation. The authors expect that with some additional information on the core management, the worth of this benchmark will be extended.

The results obtained in the frame of the current study indicated a satisfactory accuracy for the specific activity predictions obtained with VACS and the ENDF-B/VII.1 or JEFF-3.3 nuclear data libraries in the attenuation calculation. Indeed, in these cases, the final calculation/experiment (C/E) ratios remained mostly within 20% of unity in the capsules and were close to unity in the cavity. These results are satisfactory, given the global uncertainty remaining in such evaluations (±1σ ≃ 15% [13]). In contrast, the results showed that using the ENDF/B-VIII.0 nuclear data library to model fission neutron attenuation from the core to the capsule or the cavity resulted in a systematic underestimation of the dosimeter specific activities of about 30%. This has also been reported in other studies [37]. Moreover, an anomaly was observed in the ⁶³Cu(n,α) cross section in the ENDF-B/VIII.0 library.

Although the validation work described in this paper indicates a satisfactory perfor-
mance for VACS, it is not sufficient to make a definitive statement on the general accuracy
of the methodology.

**Author Contributions:** Conceptualization, R.V., M.B., J.T. and E.D.; methodology, R.V. and M.B.; software, R.V. and J.T.; formal analysis, R.V.; validation, M.B., J.T. and E.D.; writing—original draft preparation, R.V.; writing—review and editing, M.B., J.T. and E.D.; supervision, M.B., J.T. and E.D.; project administration, M.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** The authors would like to thank L. Leal from the Institut de Radioprotection et de Sûreté Nucléaire (IRSN) for fruitful discussions on the sensitivity of the results to nuclear data libraries.

**Conflicts of Interest:** The authors declare no conflict of interest.
## Abbreviations

The following abbreviations are used in this manuscript:

HBR-2 H.B. Robinson-2

IRSN Institut de Radioprotection et de Sûreté Nucléaire

PWR Pressurized Water Reactor

RPV Reactor Pressure Vessel

VACS Vessel Aging Calculation Scheme developed at IRSN

## References

1. Pichon, C.; Brillaud, C.; Deydie, D.; Alberman, A.; Soulat, P.; Chapter Neutron Spectrum Effect and Damage Analysis on Pressure Vessel Steel Irradiation Behaviour; ASTM International: West Conshohocken, PA, USA, 2000; pp. 87–97.

2. Wagner, J.C.; Haghighat, A.; Petrovic, B.G. Monte Carlo Transport Calculations and Analysis for Reactor Pressure Vessel Neutron Fluence. Nucl. Technol. **1996**, *114*, 373–398. [CrossRef]

3. Vasiliev, A.; Ferroukhi, H.; Zimmermann, M.A.; Chawla, R. Development of a CASMO-4/SIMULATE-3/MCNPX calculation scheme for PWR fast neutron fluence analysis and validation against RPV scraping test data. Ann. Nucl. Energy **2007**, *34*, 615–627. [CrossRef]

4. Matsushita, K.; Kurosawa, M. Neutron fluence analyses around the reactor pressure vessel of BWR using MCNP with a heterogeneous and homogeneous mixed core model. Prog. Nucl. Sci. Technol. **2014**, *4*, 463–466. [CrossRef]

5. Bourganel, S.; Petit, O.; Diop, C. Three-Dimensional Particle Transport Using Green's Functions in TRIPOLI-4 Monte Carlo Code: Application to PWR Neutron Fluence and Ex-Core Response Studies. Nucl. Technol. **2013**, *184*, 29. [CrossRef]

6. Vuiart, R.; Brovchenko, M.; Taforeau, J.; Jaiswal, V.; Dumonteil, E. A Versatile Methodology for Reactor Pressure Vessel Aging Assessments. Nucl. Sci. Eng. **2022**, *196*, 455–477. [CrossRef]

7. Vuiart, R.; Console Camprini, P.; Burn, K.; Brovchenko, M.; Taforeau, J.; Dumonteil, E. Impact of Fission Spectrum Modeling on Aging Estimations of Pressurized Water Reactor Vessels. Trans. Am. Nucl. Soc. **2021**, *125*, 1106–1109. [CrossRef]

8. Vuiart, R.; Brovchenko, M.; Taforeau, J. Impact of core power variations on the fast neutron flux incident on pressurized water reactor vessels. In Proceedings of the 2019 International Congress on Advances in Nuclear Power Plants (ICAPP 2019), Juan-les-Pins, France 12–15 May 2019.

9. Vuiart, R. Analyse du Vieillissement Sous Irradiation Neutronique des Cuves des Réacteurs Nucléaires à eau sous Pression. Ph.D Thesis, Université Paris-Saclay, Fontenay-aux-Roses, France, 2021.

10. Benedet, B.; Vuiart, R.; Taforeau, J.; Brovchenko, M. Sensitivity of Pressurized Water Reactors Vessels Aging to Fuel Composition. Trans. Am. Nucl. Soc. **2021**, *125*, 1110–1113. [CrossRef]

11. Remec, I.; Kam, F. H.B. Robinson-2 Pressure Vessel Benchmark; NUREF/CR-6453, ORNL/TM-13204, Oak Ridge National Laboratory: Oak Ridge, TN, USA, 1997.

12. Kodeli, I.; Sartori, E. SINBAD-Radiation shielding benchmark experiments. Ann. Nucl. Energy **2021**, *159*, 108254. [CrossRef]

13. Pittarello, R.; Vasiliev, A.; Ferroukhi, H.; Chawla, R. Refined Monte Carlo analysis of the H.B. Robinson-2 reactor pressure vessel dosimetry benchmark. Ann. Nucl. Energy **2011**, *38*, 1842–1851. [CrossRef]

14. Vasiliev, A.; Ferroukhi, H.; Kolbe, E. Performance of a Monte-Carlo solution for the H.B. Robinson-2 Pressure Vessel Dosimetry Benchmark. Ann. Nucl. Energy **2010**, *37*, 1404–1410. [CrossRef]

15. Orsi, R. H.B. Robinson-2 pressure vessel dosimetry benchmark: Deterministic three-dimensional analysis with the TORT transport code. Nucl. Eng. Technol. **2019**, *52*, 448–455. [CrossRef]

16. Bourganel, S.; Raskiniteb, I.; Soldevilac, M. Analysis of the H.B. ROBINSON-2 reactor pressure vessel dosimetry benchmark using TRIPOLI-4 Monte Carlo code. Prog. Nucl. Sci. Technol. **2014**, *4*, 312–316. [CrossRef]

17. HBR-2-Updated Final Safety Analysis Report-CHAPTER 4, Carolina Power & Light Document, Hosted by NRC, Washington D.C.; Revision No. 21. Available online: https://www.nrc.gov/docs/ML1729/ML17298A851.pdf.

18. Macfarlane, R.; Muir, D.; Boicourt, R.M.; Kahler, A.C., III; Conlin, J.L. The NJOY Nuclear Data Processing System, Version 2016; Technical Report for Los Alamos National Lab. (LANL): Los Alamos, NM, USA, 2017. [CrossRef]

19. Plompen, A.J.; Cabellos, O.; De Saint, J.C.; Fleming, M.; Algora, A.; Angelone, M.; Archier, P.; Bauge, E.; Bersillon, O.; Blokhin, A.; et al. The joint evaluated fission and fusion nuclear data library, JEFF-3.3. Eur. Phys. J. **2020**, *56*. [CrossRef]

20. Rhodes, J.; Hykes, J.; Ferrer, R. CASMO5 A Fuel Assembly Burnup Program User's Manual; Studvik Scandpower Inc.: Idaho Falls, ID, USA, 2017; Rev. 12.

21. Bahadir, T.; Lindahl, S. Studsvik's next generation nodal code SIMULATE-5. In Proceedings of the Advances in Nuclear Fuel Management IV (ANFM 2009), Hilton Head Island, SC, USA, 12–15 April 2009.

22. Williams, M.L. DOTSOR: A Module in the LEPRICON Computer Code System for Representing the Neutron Source Distribution in LWR Cores; EPRI Research Project: 1399-1 Interim Report, RSIC Peripheral Shielding Routine Collection PSR-277; EPRI: Washington, DC, USA, 1985.

23. Goorley, T.; James, M.; Booth, T.; Brown, F.; Bull, J.; Cox, L.; Durkee, J.; Elson, J.; Fensin, M.; Forster, R.; et al. Initial MCNP6 Release Overview - MCNP6 Version 1.0. Nucl. Technol. **2013**, *164*. [CrossRef]
24. Chadwick, M.; Herman, M.; Obložinský, P.; Dunn, M.E.; Danon, Y.; Kahler, A.C.; Smith, D.L.; Pritychenko, B.; Arbanas, G.; Arcilla, R.; et al. ENDF/B-VII.1 Nuclear Data for Science and Technology: Cross Sections, Covariances, Fission Product Yields and Decay Data. *Nucl. Data Sheets* **2011**, *112*, 2887–2996. [CrossRef]

25. Hébert, H. Applied Reactor Physics. In *Presses Internationales Polytechnique*; Worldcat: Montréal, QC, Canada, 2009; ISBN 978-2-553-01436-9.

26. Bahadir, T.; Lindahl, S.; Scandpower, S. SIMULATE-4 Pin Power Calculation. In Proceedings of the PHYSOR-2006 American Nuclear Society's Topical Meeting on Reactor Physics: Organized and Hosted by the Canadian Nuclear Society, Vancouver, BC, Canada, 10–14 September 2006.

27. Initial MCNP6 Release Overview. *Nucl. Technol.* **2012**, *180*, 298–315. [CrossRef]

28. Fonte, L. Modélisation et Analyse du Benchmark SINBAD du Réacteur H.B. Robinson-2; AERGON, Internal IRSN Report: Fontenay-aux-Roses, France, 2021.

29. Skutnik, S. ORIGEN-Based Nuclear Fuel Inventory Module for Fuel Cycle- Assessment Final Project Report; Department of Nuclear Engineering, University of Tennessee-Knoxville: Knoxville, TN, USA, 2017.

30. Brown, D.; Chadwick, M.; Capote, R.; Kahler, A.; Trkov, A.; Herman, M.; Sonzogni, A.; Danon, Y.; Carlson, A.; Dunn, M.; et al. ENDF/B-VIII.0: The 8th Major Release of the Nuclear Reaction Data Library with CIELO-project Cross Sections, New Standards and Thermal Scattering Data. *Nucl. Data Sheets* **2018**, *148*, 1–142. [CrossRef]

31. Trkov, A.; Griffin, P.; Simakov, S.; Greenwood, L.; Zolotarev, K.; Capote, R.; Aldama, D.; Chechev, V.; Destouches, C.; Kahler, A.; et al. IRDFF-II: A New Neutron Metrology Library. *Nucl. Data Sheets* **2020**, *163*, 1–108. [CrossRef]

32. Mosher, S.; Johnson, S.; Bevill, A.; Ibrahim, A.; Daily, C.; Evans, T.; Wagner, J.; Johnson, J.; Grove, R. ADVANTG—*An Automated Variance Reduction Parameter Generator*, *Rev. 1*; Oak Ridge National Lab. (ORNL): Oak Ridge, TN, USA, 2015.

33. Evans, T.; Stafford, A.; Slaybaugh, R.; Clarno, K. DENOVO: A New Three-Dimensional Parallel Discrete Ordinates Code in SCALE. *Nucl. Technol.* **2010**, *171*, 171–200. [CrossRef]

34. Wagner, J.; Haghighat, A. Automated Variance Reduction of Monte Carlo Shielding Calculations Using the Discrete Ordinates Adjoint Function. *Nucl. Sci. Eng.* **1998**, *128*, 186–208. [CrossRef]

35. Rhoades, W.; Childs, R. The DORT Two-Dimensional Discrete Ordinates Transport Code. *Nucl. Sci. Eng.* **1988**, *99*, 88–89. [CrossRef]

36. White, J.; Ingersoll, D.; Skater, C.; Roussin, R. BUGLE-96: A Revised Multigroup cross Section Library for LWR Applications based On ENDF/B-VI Release 3. In Proceedings of the ANS Radiation Protection & Shielding Topical Meeting, Falmouth, MA, USA, 21–25 April 1996.

37. Fischer, A. Benchmarking of ENDF/B-VIII.0 Cross Sections for LWR Ex-Core Transport Applications. In Proceedings of the PHYSOR 2022, Pittsburgh, PA, USA, 15–20 May 2022.