https://doi.org/10.1051/epjn/2017001
Regular Article
Neutronics characterization of an erbia fully poisoned PWR assembly by means of the APOLLO2 code
ENEA,
Via Martiri di Monte Sole 4,
40129
Bologna, Italy
^{⁎} email: roberto.pergreffi@enea.it
Received:
19
October
2016
Received in final form:
19
December
2016
Accepted:
24
January
2017
Published online: 21 March 2017
Recently, increasing demands on the reduction of fuel cycle costs have led to higher burnup fuel designs. According to the erbiacredit super high burnup fuel concept, developed by mixing low content of erbia to UO_{2} powder directly after reconversion process so that all fuel pins in a given fuel assembly are homogeneously doped, the present study aims to characterize, from a neutronic point of view, a 17 × 17 pressurized water reactor assembly enriched to 10.27 wt.% in ^{235}U with an erbia content of 1 at.% (i.e. 0.7 wt.%) by means of the deterministic neutronic code APOLLO2. For this purpose, a simplified thermalhydraulic analysis was performed in order to evaluate the effects on fuel thermal conductivity of adding erbia to uranium oxide. The results obtained allow to conclude that an Erdoped assembly enriched to >5 wt.% in ^{235}U represents an advantageous solution for very long fuel cycles, and it is so suited for very high burnups.
© R. Pergreffi et al., published by EDP Sciences, 2017
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The idea of using neutron poison materials was originally developed in order to increase the allowable initial core fuel enrichment. In fact, the high neutron absorption cross sections of such materials permit to compensate, during the early stages of core life, the excess reactivity due to higher initial enrichment. Moreover, they burn out somewhat faster than fuel so that their contribution in core life in terms of negative reactivity is negligible [1].
Erbia's (Er_{2}O_{3}) role as a poison in light water reactor (LWRs) was first highlighted in 1970, but only late in the 1980s that it was recognized as an alternative absorber to gadolinia (Gd_{2}O_{3}) and, like that, mixed with fissile material in a small number of fuel pins of an assembly [2]. From a safety point of view, erbium was found to be very effective at minimizing radial power peaking thanks to relatively low thermal absorption effective cross section (162 ± 8 barns vs. 49,000 ± 1000 barns of gadolinium) and at controlling transients thanks to the higher resonance integral (740 ± 10 vs. 390 ± 10 barns). Moreover, erbium has a rather short evolution chain and, differently from gadolinium, its efficiency as a function of the content and number of poisoned rods, is nearly linear [3,4]. Over the years, erbia has been used in a certain number of pressurized water reactors (PWRs) [5].
Recently, increasing demands on the reduction of fuel cycle costs has led to higher burnup fuel designs. In UOx–LWRs, extended burnups are achieved by higher initial uranium enrichment and consequently higher amounts of gadolinia. For this purpose, one of the main issues before very high assembly average burnups (>70 GWd/MTU) can be achieved is represented by the current enrichment limit for commercialtype LWR fuel, that is 5 wt.%, due to criticality safety requirements related to the design of fabrication plants.
The erbiacredit super high burnup (ErSHB) fuel concept developed a few years ago goes beyond the so called “5 wt.% barrier” without requiring significant modifications and relicensing of fuel cycle facilities. The innovative idea is to mix low amount of erbia to UO_{2} powder directly after the reconversion process so that all fuel pins, in a given fuel assembly, are homogeneously doped [6,7]. In this way, notwithstanding the ^{235}U enrichment exceeds the current limit, the initial reactivity is equivalent to ≤5 wt.% fuel. It is worth noting that analyzing the multiplication factor of a fuel assembly enriched to 6 wt.% in ^{235}U with a poison content set to 0.2 wt.%, it was found that gadolinia does not represent an alternative solution to erbia for such high burnup fuel concept because of a too large suppression of reactivity at beginning of life (BoL) due to its very large absorption cross section [7]. The ErSHB fuel concept was not only studied from a neutronic viewpoint, in fact the effects of erbia addition on the thermal and mechanical properties of (U_{1−x}Er_{x})O_{2} (with 0 ≤ x ≤ 0.1) were measured [8].
The present work aims to characterize a 17 × 17 PWR assembly poisoned with erbia from a neutronic point of view comparing its neutronic performance to that of conventional UO_{2} fuel assembly enriched to 5 wt.% in ^{235}U. According to ErSHB fuel concept, erbium oxide is homogeneously dispersed into all fuel pins. More in detail, some crucial parameters of both fuel configurations such as multiplication factor, residual reactivity penalty, neutron spectrum, spectral index, boron reactivity worth and temperature reactivity coefficients were calculated at steadystate conditions. For this purpose, it was decided to use the APOLLO2 code that is recognized worldwide as a standard for deterministic cell calculations for LWRs. Moreover, a simplified thermalhydraulic analysis was performed in order to evaluate the effects on fuel thermal conductivity of adding erbia to uranium oxide.
2 Materials and methods
2.1 Geometry and material data
The geometry adopted in this study corresponds to an EPR™ like assembly containing 265 fuel rods, each with an active length of 420 cm, and 24 guide tubes. The fuel assembly is assumed to be composed of fuel cladding and moderator material only. Because no control rods have been taken into account, guide tubes have been filled with water (measurement and evaluation of the control rod worth is postponed to a future study). Helium into the gap between fuel and cladding has been neglected except for thermalhydraulic analysis (see next paragraph). The water blade around the assembly is set to a value of 0.084 cm. Dimensions of fuel, clad and guide tubes at room temperature are provided in Table 1.
The fuel is assumed to be pure UO_{2}, only containing ^{235}U and ^{238}U, mixed with small quantities of Er_{2}O_{3}; no ^{234}U was included, being considered negligible at first approximation for making comparisons between systems. Except for isocriticality curves that were drawn at varying ^{235}U enrichment and erbium content, for the rest of the analysis these values were fixed at 10.27 wt.% and 0.7038 wt.% (= 1 at.%), respectively. If it is not differently specified, the boron concentration in the moderator is set to 1000 ppm by mass (value for the BoC of the EPR™ first cycle). Cladding and guide tube are taken to be M5^{®}, a ternary alloy licensed by AREVA as fuel cladding material up to burnup values of 80 GWd/MTU. In this study, M5^{®} was modeled as zirconium (98.875 wt.%), niobium (1 wt.%) and oxygen (0.125 wt.%). The M5^{®} chemical composition at BoL at operating temperature is described in Table 2.
The presence of any structural materials (e.g., spacer grids) has been ignored.
The expanded dimensions of fuel, clad and guide tubes at operating temperatures have been obtained using linear thermal expansion coefficients. Value, reference temperature and source for each coefficient are detailed in Table 3.
Material densities at operating temperatures are provided in Table 4. In particular, pellet density depends on the content of erbia. As reported in [8], increasing the erbium content from 0 at.% to 5 at.% the percentage of the theoretical density decreases going from 95% to 90%. As almost all simulations refer to 1 at.% Erdoped UO_{2} pellet, a percentage of the theoretical density of 94% was considered. In this case, the fuel density is 9.6955 g/cm^{3}. To calculate it, the fuel mass in one rod at room temperature has been divided by the internal volume of the pin at operating temperature. The effect of boric acid on the moderator density has been neglected.
Table 5 shows the value of average linear power density required for the normalization of the neutron flux. This value has been computed multiplying the total core thermal power (4500 MWth) by the power percentage generated in the fuel (97.4%) and dividing by the number of assemblies (241).
Fuel assembly geometry data.
Clad and guide tube composition.
Linear thermal expansion coefficients.
Densities of fuel, clad and moderator.
Average linear power density.
2.2 Description of neutronic code
All simulations were performed by means of the deterministic transport code APOLLO2 version 2.83.E [10]. APOLLO2 is a modular cell code for 2dimensional multigroup transport calculations. Microscopic cross sections refer to 281 groups' master library with a SHEM group structure CEA2005V4.1.2.patch, based on JEFF 3.1.1 evaluations. Calculations were done assuming reflective boundary conditions and taking advantage of all symmetries of the system, so that only 1/8 of assembly was considered. A schematic of 1/8 of modeled assembly is given in Figure 1. Elementary fuel cell was divided into four concentric zones corresponding, from center to periphery, to 50%, 30%, 15%, and 5% of the total volume, respectively, one clad zone and one moderator zone (the search for an optimized number of fuel regions to rigorously take into account the rim effect for erbia is postponed for a future study). The elementary guide tube cell was modeled in two moderator zones, divided by the guide tube. Collision Probability method was used in APOLLO2 to resolve the transport equation to compute the neutron flux. Selfshielding calculations were done bearing in mind the following order of isotopes: ^{238}U, ^{235}U, ^{167}Er, ^{239}Pu, ^{240}Pu, Zr_nat. TR (all resonances) approximation was used over the entire energy range. The effect of leakage on the neutron spectrum was taken into account by means of homogenous B1 model that computes the buckling to assure criticality.
Fig. 1 Model of 1/8 of the fuel assembly. 
3 Thermalhydraulic analysis
The purpose of the analysis is to calculate the impact on PWR fuel pin temperature distribution due to erbia poisoning. Therefore, a pin containing 1 at.% of erbia and enriched to 10.27% in ^{235}U was compared with a reference case with enrichment of 5% and no poisoning. From a preliminary qualitative analysis, two main antagonistic effects are foreseen which are as follows:

reduction of fuel conductivity in the poisoned pin is expected to cause a rise in centerline fuel temperature;

increase in neutron absorption in the peripheral area of the pellet due to the presence of erbia and to the higher fissile concentrations is expected to cause an increase in centre pellet flux depression with consequent decrease of centerline fuel temperature.
The thermalhydraulic analysis was based on the data from an EPR™ preconstruction safety report [11]. Some thermalhydraulic parameters are described in Table 6.
The analysis was carried out at BoL. The power to be transferred from the fuel rods to the coolant was calculated as described in previous paragraph. No hottest pin analysis was necessary because the neutronic calculations were performed, as usual, using reflective boundary conditions; therefore, no difference of flux and, consequently, of power is possible with only one type of fuel pins in the assembly.
Heat transfer from the pins to the coolant is a typical case of convection involving heat transfer between a surface and an adjacent fluid at different temperatures. The dominant contribution to such heat transfer is the bulk motion of fluid, with a minor contribution by conduction. The convective heat flux is modeled by Newton's law of cooling, with all the variations due to fluid properties, surface geometry and flow conditions included in the convection coefficient h_{w}: (1) where T_{we} is the external clad temperature, T_{B} the bulk temperature (averaged over the core), q^{′} the average linear power, and D the subchannel hydraulic diameter.
The convection coefficient depends upon many variables; the approach to its determination consists in identifying universal functions of dimensionless groups with physical meaning for convective flow. The Nusselt number Nu = hL/k, where L and k represent the characteristic length of the surface and the thermal conductivity, respectively, physically represents the dimensionless temperature gradient at the surface. In forced convection Nu can be correlated with the Reynolds number Re = ρvL/μ, physically representing the ratio of inertia and viscous forces, and the Prandtl number Pr = c_{p}μ/k, representing the ratio of the momentum and thermal diffusivity. In water cooled reactors for purely single phase flow, the Nusselt number can be evaluated by employing the DittusBoelter correlation: (2)
All the physical properties and correlations for water were derived by interpolating the data from NIST standard reference data. In the case considered Nu has a value of 844 and consequently h_{w} = 37,157 W/(m^{2} K). The external cladding temperature results T_{we} = 325 °C.
Heat transfer through the cladding material was modeled by the steady state heat conduction equation reduced to a onedimensional equation in the radial direction in absence of heat generation. After integration the equation becomes: (3) where t_{c} is the cladding material thickness (M5^{®}). The conductivity of M5^{®} was taken from [12] and set to the value of k_{c} = 18 W/(m K). The internal cladding temperature was found to be T_{wi} = 344 °C.
Heat transfer through the gap was modeled through the standard general equation: (4)
According to the RossStaute gap conductance model, reported in [13], the gap conductance h_{G} was calculated as: (5) where k_{G} is the thermal conductivity of the gap gas (pure helium at the beginning of cycle) and G is the gap thickness at the beginning of cycle. At the beginning of cycle the effects of burnup on the pellet are not detectable yet, but the pellet expansion can't be calculated as simple thermal expansion because, as soon as the reactor power increases, pellets crack due to thermal stress induced by radial temperature gradients [13]. This phenomenon causes an increase of the apparent pellet diameter and consequently a significant gap reduction. This gap thickness at the beginning of cycle during irradiation was calculated by application of the model proposed by Oguma [13] as (6)where G is the gap of the beginning of cycle full power pin, G_{0} is the initial asfabricated gap, D is the pellet diameter all in microns and q^{′} is the linear power in W/cm. In such model, the experimental constants were derived from in pile gap analyses of instrumented fuel rods irradiated in a boiling water reactor (BWR). Therefore, to apply this model to a PWR, a correction was introduced to take into account the difference in coolant pressure and its effect on elastic deformations of the cladding material. The correction was calculated using, for the calculation of the Young modulus of Zircaloy, the MATPRO formula [12]: (7)by neglecting the modifications due to the effect of oxidation (K_{1} = 0) and the effect of fast neutron fluence (K_{3} = 1) and calculating the modification to account for the effect of cold work (the fractional reduction in crosssection area due to processing) as: (8)with the cold work assumed as C = 0.2 (default value of FRAPCON code) [14]. The pellet surface temperature determined this way was T_{ps} = 453 °C.
Heat transfer through the fuel was modeled by the steady state heat conduction equation reduced to a onedimensional equation in the radial direction: (9) where q^{‴} is the power density generated by the nuclear reactions within the fuel. For the solution, the q^{‴} radial distribution needs to be known. This distribution can be obtained by parabolic interpolation of power densities calculated by APOLLO2 in ten different concentric fuel regions, but a value of the effective fuel temperature T_{eff} is required by APOLLO2 for effective cross sections determination. To provide a first attempt value of T_{eff}, the maximum centerline temperature T_{pc} was assumed to have the value resulting by solving the heat equation in the assumption of a radially constant q^{‴}. This first attempt value of T_{eff} was calculated as follows [15]: (10)
With this value of T_{eff} a first attempt neutronic simulation was run to obtain a first attempt power distribution. These results were interpolated by a curve in the form: (11)
Such expression was introduced in (9). After integration and separation of variables and considering that the minimum fuel temperature is reached on the pellet surface T_{ps} and the maximum in the centerline T_{pc}_{,} the heat equation becomes: (12) where T_{pc} is the unknown value to be found, and R is the fuel pin radius whose thermal expansion was calculated by the following correlation [16]: (13)
Every time a new value of T_{eff} was calculated, R(T_{eff}) was updated. According to [8], the temperature dependence of the thermal conductivity of erbia enriched uranium dioxide (U_{1−x}Er_{x})O_{2} as a function of Er content, x is: (14) after substituting and integrating over T, the centerline temperature T_{pc} could be determined in the cases of x = 0 and x = 0.01 (i.e. 1 at.% of erbia). With this temperature a second attempt value of T_{eff} was calculated as before and a second more accurate neutronic simulation was carried out. The procedure was iterated until convergence was reached; starting with the value of T_{eff} calculated assuming a radially constant q^{‴}, convergence to within ±0.1°C was reached with three iterations.
In the poisoned pin the power distribution showed a slightly higher variation between the peripheral area and the pellet center with respect to the reference case. The power calculated for the outermost of the ten regions considered in the simulation resulted 4% higher than the reference case, while in the inner region it was 3% lower. Consequently, the last iteration values of the parabolic interpolation coefficients of the power distribution were a = 2.39 × 10^{12} W/cm^{5} and c = 2.85 × 10^{8} W/cm^{3} for the reference case and a = 3.74 × 10^{12} W/cm^{5} and c = 2.74 × 10^{8} W/cm^{3} for the erbia poisoned pin. The two values of maximum centerline temperatures were as follows:

T_{pc} = 740 °C for the reference case with x = 0,

T_{pc} = 744 °C for the poisoned pin with x = 0.01.
The increase in temperature is very limited given the small decrease (−3%) of the conductivity of the poisoned fuel at the higher enrichment. This effect is further alleviated by the more favorable power distribution, characterized by a lower power generation in the inner regions where the heat produced is more difficult to dispose of. The radial temperature profile of the erbia enriched fuel pellet is plotted in Figure 2.
All relevant temperatures used for the neutronics analysis described in the next paragraph are summarized in Table 7.
Thermalhydraulic data.
Fig. 2 Radial temperature profile of the erbia enriched fuel pellet. 
Temperatures of pellet, clad and moderator.
4 Neutronic analysis
The main results of the neutronics characterization of a (U_{,} Er)O_{2} fuel assembly are presented below.
First of all, three isocriticality curves at BoL related to three different boron concentrations are plotted in Figure 3. Each curve correlates the erbium content with the uranium enrichment so that the multiplication factor of the (U_{1−x}Er_{x})O_{2} fuel is equivalent (difference less than ±5 ppm) to that of the corresponding UO_{2} fuel (i.e. enriched to 5 wt.% in ^{235}U). It is worth noting that, increasing the boron concentration, the erbium content required to compensate the initial excess reactivity increases. This means that the negative reactivity worth of erbium is partially reduced by that of boron.
In order to make the interpretation of the following results easier, the three considered fuel configurations have been summarized in Table 8. The choice of a UO_{2} fuel enriched to 5 wt.% in ^{235}U as Reference case depends on the fact that this enrichment represents a superior limit for the commercialtype power reactors.
The increase of initial ^{235}U enrichment beyond the “5 wt.% barrier” requires an accurate assessment of criticality implications. This includes verifying that the positive extrareactivity due to higher enrichment is controlled by the negative reactivity due to erbia at every burnup step. Figure 4 shows the infinite neutron multiplication factor as a function of burnup of the three different configurations. As it can be noted, Case2 has the same criticality at BoL of Reference case (5 wt.% fuel) but, differently from this, it achieves k = 1 when the burnup value is about 76 GWd/MTU. The fact that, notwithstanding the initial enrichment of Case2 is to 10.27 wt.%, its initial criticality does not exceed the maximum criticality value of the Reference case (1.32400 vs. 1.32402), is a remarkable aspect from a safety point of view. It is also noteworthy to observe the fact that the reduction of initial reactivity of Case2 with respect to Case1 (1.32400 vs. 1.45179) is completely due to the erbia content. Furthermore, the burnup for which k = 1 is quite similar to the value indicated in some studies as optimal point in terms of generation cost for a conventional PWR with >5 wt.% uranium enrichment level [17,18].
The effects of erbia on the assembly reactivity can also be taken into account in terms of residual reactivity penalty. The reactivity penalty is defined as the difference of reactivity between Case1 and Case2, i.e. between two different fuels both enriched to 10.27 wt.% in ^{235}U with and without erbium oxide. As plotted in Figure 5, the residual reactivity penalty decreases with increasing burnup and it is of about 1000 pcm when the burnup exceeds 76 GWd/MTU. As suggested in [6], the residual penalty could be removed by modifying erbium isotopically, i.e. eliminating ^{166}Er. In fact the residual reactivity penalty at end of life (EoL) is mainly due to the presence of ^{167}Er isotope which is built up from ^{166}Er. This fact is confirmed by ^{167}Er and ^{166}Er concentrations as a function of burnup plotted in Figure 6. In addition, if the ^{167}Er concentration significantly lessens in the first part of the curve as a result of its absorption macroscopic cross section, it is substantially constant over time starting from a burnup value of about 40 GWd/MTU.
Another relevant neutronic aspect is related to the production of ^{239}Pu. The concentration of ^{239}Pu as a function of burnup is plotted in Figure 7 with reference to the two fuel configurations: Reference case (curve in blue) and Case2 (curve in red).
This concentration expresses the balance between production and destruction of ^{239}Pu instant by instant. It is well known that:

^{239}Pu production is roughly proportional to fast flux and ^{238}U effective macroscopic cross section;

^{239}Pu destruction is roughly proportional to thermal flux and ^{239}Pu effective macroscopic cross section.
Therefore, denoting the production and destruction terms as: (15) (16) and analyzing the two curves in Figure 7, we can observe that these two terms in the two configurations are very similar up to a burnup value of about 20 GWd/MTU. In other words the following relation can be written as: (17)
It is however important to underline that the equivalence of the differences in the two cases does not imply the equivalence term to term. On the contrary, analyzing the data in Tables 9 and 10 it can clearly be seen that ^{238}U concentration as well as the fast and thermal fluxes at BoL are very different in the two configurations. In addition, easy physical considerations permit to extend the same conclusion also to the effective microscopic cross sections of uranium and plutonium. Between 20 GWd/MTU and 76 GWd/MTU, ^{239}Pu concentration of Case2 builds up differently from the Reference case. This fact means that in the assembly with lower enrichment in ^{235}U, ^{239}Pu reaches equilibrium at a much lower burnup value because of its more important contribution to the total power produced.
In Figure 8, the normalized fluxes per unit lethargy at 0 GWd/MTU of both configurations are plotted. As it can clearly be seen, the flux of the assembly poisoned with erbia is less thermalized than the other. In fact, the increased thermal neutron absorption due to erbium causes a hardening of the flux. This effect is confirmed by analyzing the spectral index values at different burnup steps in Table 10. The spectral index, which is defined here as the ratio between fast and thermal flux, is always higher in the Erdoped assembly, owing to the increased thermal neutron absorption or, similarly, to the larger amount of absorbers at thermal energy in that assembly. But, if each single term is analyzed, an opposite trend can be observed: at each burnup step, both the fast flux and the thermal flux are larger in the reference assembly than in the assembly poisoned with erbia. After all, even if the neutron flux is less hardened at low enrichment, the total flux has to be bigger in order to produce the same power.
In Figure 9, the evolution of the flux per unit lethargy of Case2 configuration is plotted with reference to four different burnup values. As a rule, in order to produce the same power, the total flux has to increase because of the reduction of fissile material. It is interesting to observe that, while going from low to high burnup values, the thermal part of the flux is progressively more relevant. In addition, the dip at thermal energies (∼1 eV) due to ^{239}Pu resonances can clearly be seen.
The reactivity worth of 1 ppm of boron according to the following definition: (18) where w_{B} is the boron worth, k is the multiplication factor, and c_{B} is the boron concentration, is plotted in Figure 10 as a function of burnup. This graph refers to the assembly enriched to 10.27 wt.% in ^{235}U and poisoned with erbia. The boron reactivity worth has a moderate variation over time up to a burnup value of about 65 GWd/MTU. In fact after a rapid reduction due to the xenon and samarium poisoning, the boron worth goes from a value of about 5 pcm/ppm to a value of 4.5 pcm/ppm as a result of the spectrum hardening due to the plutonium buildup.
In terms of boron controlled reactivity, no differences arise with respect to a standard PWR, given the fact that k_{∞} at BoL is the same for the two systems. However, the Erdoped fuel assembly needs almost twice boric acid to achieve the same reactivity value. In fact the boron reactivity worth of Case2 is about half that of Reference case, even if initial concentration values are the same (1000 ppm). This fact may have some drawbacks as far as the primary coolant chemistry is concerned; the analysis of this aspect is however out of the scope of the present paper.
One of the most important aspects of the fission process from the reactor control viewpoint, is the presence of effective delayed neutrons. The effective delayed neutron fraction − β_{eff} − is calculated weighting delayed neutrons on the adjoint neutron flux. For this purpose, adjoint flux calculations were performed together with direct flux calculations at each burnup step. In Figure 11, three examples of normalized adjoint neutron flux at three different burnup values have been plotted. The effects of ^{239}Pu resonances are very well shown in the two curves corresponding to burnup values of 30 and 70 GWd/MTU.
β_{eff} as a function of burnup with reference to two fuel configurations is shown in Figure 12. All calculations were performed with APOLLO2 considering production and decay of eight groups of delayed neutron precursors. In both curves, increasing the burnup value, the effective beta decreases, going from about 800 to about 400 pcm as a result of the ^{239}Pu buildup and ^{235}U depletion. But the effective beta of the assembly poisoned with erbia decreases more slightly than the Reference case.
Temperature reactivity coefficients are crucial parameters in transients of LWRs. As the temperature does not change uniformly throughout the assembly, two different temperature reactivity coefficients related fuel and moderator were calculated:

fuel temperature coefficient (FTC), denoted as α_{Dop} and defined as the fractional change in k per unit change in effective fuel temperature;

moderator temperature coefficient (MTC), denoted as α_{Mod} and defined as the fractional change in k per unit change in moderator temperature considering a change in water density [19].
All temperature reactivity coefficients have been computed according to the following relation [20]: (19) where ρ is the reactivity, k_{i} and k_{f} are the initial and final multiplication factors, and T is the temperature.
The FTC for Erdoped assembly, in going from initial effective fuel temperature (540 °C) to final temperatures, produces a reactivity change very similar to other fuels and in particular to the Reference case, as shown in Table 11. Considering that all calculations were performed at zero burnup, the largest contribution to the FTC is due to the increase in resonant capture by ^{238}U. Nevertheless, the ^{167}Er and ^{166}Er contribution in Erdoped assembly is not negligible. In particular, the positive contribution of the erbium isotopes to the FTC can be appreciated by comparing the two configurations enriched to 10.27 wt.% in ^{235}U.
The same reasoning can be extended to the MTC (Tab. 12). In this case too, a change of temperature causes a change of the opposite sign in the reactivity. In addition, in going from initial effective moderator temperature (314 °C) to final temperatures, the reactivity of the Erdoped assembly changes much more than those of the other fuels also including the Reference case. As the temperature effect on reactivity is due to moderator density change, Table 13 summarizes water density values used for MTC calculations. It is worth noting that the water density depends only on temperature (i.e., the effect of boric acid concentration on water density was neglected) [21].
It can be concluded that the performance of Erdoped assembly from the safety viewpoint is at least as good as that of the Reference case, i.e., of UO_{2} fuel enriched to 5 wt.% in ^{235}U. The effect of high Erdoping on the worth of control rods remains to be evaluated in a future study.
Fig. 3 Isocriticality curves at BoL. 
Fuel configurations.
Fig. 4 Multiplication factor vs. burnup. 
Fig. 5 Residual reactivity penalty. 
Fig. 6 ^{167}Er and ^{166}Er concentrations vs. burnup. 
Fig. 7 Burnup varying ^{239}Pu concentration. 
^{235}U and ^{238}U concentrations at BoL.
Spectral index.
Fig. 8 Normalized flux per unit lethargy. 
Fig. 9 Flux per unit lethargy. 
Fig. 10 Boron worth as a function of burnup. 
Fig. 11 Normalized adjoint neutron fluxes. 
Fig. 12 Burnup varying effective beta for two fuel configurations. 
Fuel temperature coefficient.
Moderator temperature coefficient.
Water density.
5 Conclusions
The results of the comparison between a 17 × 17 PWR assembly enriched to 10.27 wt.% in ^{235}U with an erbia content of 1 at.% (i.e. 0.7 wt.%) and a conventional UO_{2} fuel assembly (enriched to 5 wt.% in ^{235}U) in terms of neutronics parameters such as multiplication factor, ^{239}Pu concentration, neutron spectrum, spectral index, beta effective and temperature reactivity coefficients are summarized below:

Erdoped assembly can reach a burnup of 76 GWd/MTU, more than twice the 36 GWd/MTU burnup of conventional fuel assembly but without exceeding the limit represented by the maximum criticality value.

EoL ^{239}Pu concentration in Erdoped assembly is twice less than that in the conventional fuel assembly.

Effective beta, as a function of burnup, decreases more slowly in Erdoped assembly than in the conventional fuel assembly.

Fuel and moderator temperature reactivity coefficients of Erdoped assembly are at least as good as the conventional ones.
These results allow us to infer that, from a neutronic point of view, the performance of the erbia poisoned fuel assembly is at least as good as that of a conventional fuel assembly. In addition, the numerical analysis performed to assess the effects of erbia in a PWR fuel pin showed that, from a thermalhydraulic point of view, the addition of 1 at.% erbia to uranium oxide produces a very limited increase in maximum centerline temperature.
It can therefore be concluded that this kind of fuel assembly may represent an advantageous solution to achieve very high burnups, provided that mechanical interaction between fuel and cladding does not constitute a strong limiting constraint. However, this aspect is out of the scope of this study.
Nomenclature
at.%: percentage of number of atoms
β_{eff}: effective delayed neutron fraction
D: subchannel hydraulic diameter
FTC: fuel temperature coefficient
GWd/MTU: gigawatt days per metric ton of uranium
h: convective heat transfer coefficient of the flow
k_{c}: thermal conductivity of the cladding material
k_{f}: thermal conductivity of the fuel
k_{G}: thermal conductivity of the gap gas
M5^{®}: lowcorrosion zirconium alloy
MTC: moderator temperature coefficient
pcm: per cent mille − unit of reactivity corresponding to 10^{−5} Δk/k
ppm: partspermillion, 10^{−6}
PWR: Pressurized water reactor
q^{′}: average linear power density
T_{B}: average bulk temperature
t_{c}: cladding material thickness
T_{eff}: effective fuel temperature
T_{pc}: pellet centerline temperature
T_{ps}: pellet surface temperature
T_{we}: external clad temperature
T_{wi}: internal clad temperature
µ: dynamic viscosity of the fluid
Acknowledgments
The APOLLO2 code is developed by CEA and coowned by CEA, EDF and AREVA NP.
References
 J. Duderstadt, L. Hamilton, Nuclear reactor analysis (John Wiley & Sons, New York, 1976) (In the text)
 J. Porta, M. Asou, Erbium: alternative poison? Stabilisation additive? What future? Prog. Nucl. Energy 38, 347 (2001) [CrossRef] (In the text)
 M. Asou, J. Porta, Prospects for poisoning reactor cores of the future, Nucl. Eng. Des. 168, 261 (1997) [CrossRef] (In the text)
 J. Porta et al., Qualification of the neutronic efficiency of erbium at zero burnup, Prog. Nucl. Energy 38, 355 (2001) [CrossRef] (In the text)
 OECD, Nuclear energy agency report no. 6224, very high burnups in light water reactors (OECD, Paris, 2006) (In the text)
 M. Yamasaki, The study on erbia credit superhighburnup fuel with isotopically modified erbia, in ANS 2010 Winter Meeting, Las Vegas, November 7–11, 2010 (2010) (In the text)
 M. Yamasaki et al., Development of erbiacredit super high burnup fuel: experiments and numerical analyses, Nucl. Technol. 177, 63 (2012) [CrossRef] (In the text)
 S. Yamanaka et al., Thermal and mechanical properties of (U, Er)O_{2}, J. Nucl. Mater. 389, 115 (2009) [CrossRef] (In the text)
 IAEATHPH, Thermophysical properties of materials for nuclear engineering: a tutorial and collection data (IAEATHPH, Vienna, 2008) (In the text)
 R. Sanchez et al., APOLLO2 Year 2010, Nucl. Eng. Technol. 42, 474 (2010) [CrossRef] (In the text)
 UK EPR, The preconstruction safety report, subchapters 4.2 (fuel system design) and 4.4 (thermal and hydraulic design) (UK EPR, 2012), available at: http://www.eprreactor.co.uk/scripts/ssmod/publigen/content/templates/show.asp?P=290&L=EN (In the text)
 W.G. Luscher, K.J. Geelhood, Material Property Correlations: Comparisons between FRAPCON3.4, FRAPTRAN 1.4 and MATPRO, USNRC NUREG/CR7024, August 2010 (2010) (In the text)
 M. Oguma, Cracking and relocation behaviour of nuclear fuel pellets during rise to power, Nucl. Eng. Des. 76, 35 (1983) [CrossRef] (In the text)
 K.J. Geelhood, W.G. Luscher, C.E. Beyer, FRAPCON3.4: A Computer Code for the Calculation of SteadyState ThermalMechanical Behavior of Oxide Fuel Rods for High Burnup, NUREG/CR7022, PNNL19418, Washington D.C. (2011), Vol. 1 (In the text)
 T. Kozlowski, T. Downar, Pressurized water reactor MPX/UO_{2} core transient benchmark final report, NEA/NSC/DOC(2006)20, Technical Report (OECD/NEA and US NRC, Paris, 2006) (In the text)
 U.S.NRC, Material Property Correlations: Comparisons between FRAPCON3.4, FRAPTRAN 1.4, and MATPRO, March 2011 (2011) (In the text)
 R. Gregg, A. Worrall, Effect of highly enriched/highly burnt UO_{2} fuels on nuclear design parameters and economics, in Advances in Nuclear Fuel Management III (ANFM 2003), Hilton Head Island, South Carolina, October 5–8, 2003 (2003) (In the text)
 J.R. Secker et al., Optimum Discharge Burnup and Cycle Lenght for PWRs, in Advances in Nuclear Fuel Management III (ANFM 2003), Hilton Head Island, South Carolina, October 5–8, 2003 (2003) (In the text)
 J.R. Lamarsh, A.J. Baratta, Introduction to nuclear engineering (Prentice Hall, Upper Saddle River, New Jersey, 2001) (In the text)
 R.D. Mosteller, The Dopplerdefect benchmark: overview and summary of results, in Joint International Topical Meeting on Mathematics & Computation and Supercomputing in Nuclear Applications, Monterey, April 15–19, 2007 (2007) (In the text)
 S.C. McCutcheon, J.L. Martin, T.O. Barnwell Jr., Water quality, in Handbook of hydrology, edited by D.R. Maidment (McGrawHill, New York, 1993) (In the text)
Cite this article as: Roberto Pergreffi, Davide Mattioli, Federico Rocchi, Neutronics characterization of an erbia fully poisoned PWR assembly by means of the APOLLO2 code, EPJ Nuclear Sci. Technol. 3, 8 (2017)
All Tables
All Figures
Fig. 1 Model of 1/8 of the fuel assembly. 

In the text 
Fig. 2 Radial temperature profile of the erbia enriched fuel pellet. 

In the text 
Fig. 3 Isocriticality curves at BoL. 

In the text 
Fig. 4 Multiplication factor vs. burnup. 

In the text 
Fig. 5 Residual reactivity penalty. 

In the text 
Fig. 6 ^{167}Er and ^{166}Er concentrations vs. burnup. 

In the text 
Fig. 7 Burnup varying ^{239}Pu concentration. 

In the text 
Fig. 8 Normalized flux per unit lethargy. 

In the text 
Fig. 9 Flux per unit lethargy. 

In the text 
Fig. 10 Boron worth as a function of burnup. 

In the text 
Fig. 11 Normalized adjoint neutron fluxes. 

In the text 
Fig. 12 Burnup varying effective beta for two fuel configurations. 

In the text 