Modelling of as-fabricated porosity in UO2 fuel by MFPR code

For consistent modelling of behaviour of as-fabricated porosity in UO2 fuel irradiated under various conditions of in-pile and out-of-pile tests as well as under normal and abnormal conditions of nuclear reactor operation, the additional analysis of experimental observations and critical assessment of available models are presented. On this base, the mechanistic MFPR code, including physically-grounded models for the fuel porosity evolution in UO2 fuel under various irradiation and thermal regimes, is refined. These modifications complete the consistent description of the fuel porosity evolution in theMFPR code and result in a notable improvement of the code predictions.


Introduction
The in-pile dimensional behaviour of oxide fuels in nuclear reactors is a well-known phenomenon of great technological interest.It is generally established that at the beginning of irradiation the fuel densifies due to shrinking of the asfabricated pores remaining from the fuel sintering process with a wide distribution of their sizes [1,2].The densification is most pronounced in low density fuel, especially in the case of fine-dispersed porosity with pores typically less than one micron diameter.
Re-sintering in the furnace can be generally understood and described analytically by thermal diffusion processes, but not so in-pile densification: it was additionally assumed by Stehle and Assmann [3] that in-pile densification is a mixed athermal/thermal process, including the thermal evaporation of vacancies from pores (which dominates at relatively high temperatures above ≈1200 °C), and the athermal atomization of pores into lattice vacancies by fission spikes.
For consistent modelling of porosity behaviour in UO 2 fuel irradiated under various conditions of in-pile and outof-pile tests as well as under normal and abnormal conditions of nuclear reactor operation, the critical assessment of available models, their modification and development of more advanced models for implementation in the mechanistic codes, become rather an important task.The code MFPR (Module for Fission Products Release) was developed for analysis of fission products (FP) release from irradiated UO 2 fuel in collaboration between IBRAE and IRSN (Cadarache, France) [4,5].The mechanistic approach applied in this code allows the realistic consideration of fuel porosity evolution, self-consistently with analysis of FP release, based on physically-grounded parameters and mechanisms.
Some important modifications of the existing models of MFPR and development of new models for the fuel porosity evolution in UO 2 fuel under various irradiation and thermal regimes are presented in this paper.

Initial fuel porosity
Optical microscopy reveals that the majority of the internal cavities are located on grain boundaries [6]; the pores are generally non-spherical in shape.In the current analysis, the pores are considered as intergranular lenticular voids with the dihedral angle u = 50°.Their volume, V, and surface area, S, are [7]: where the pore curvature radius, R, relates to the experimentally measured median radius, r (which below will be simply referenced to as 'radius') as: where j % 1:29 (see Appendix A of Ref. [8]).
In the typical fresh fuel, the pores sizes are distributed within the interval 0.1-10 mm, their density distribution function can be satisfactory approximated as [8]: where C 0 is the total pore concentration and r is the mean radius.The maximum contribution to the total porosity makes pores with r ¼ 3r.
Figure 1 illustrates approximations of the experimental distributions observed in reference [9] for normal grain (mean grain diameter d gr = 8 mm) and large grain (23 mm) samples.The normal grain data are approximated by equation (3) with the mean radius of 0.35 mm and total porosity of 4.7%.Note that pores with r from 0.5 to 2 mm contribute near 75% to the total porosity.The large grain data are approximated by a superposition of two exponents corresponding to two pore populations, P1 and P2, with the mean radii of 0.45 and 3.0 mm, the partial porosities being 1.2 and 4.3% respectively.

Mechanisms of pore size relaxation
If the grain boundary self-diffusion is the rate controlling mechanism of the thermal pore relaxation, the pore volume change is described by the equation [10,11]: where D gb is the grain boundary diffusivity, w ≈ 0.5 nm is the thickness of the grain boundary layer, V = 4.09 Â 10 -29 m 3 is the atomic volume, k is the Boltzmann constant, T is the temperature.The pressure difference, dP, is given by equation: where N p is the number of gas atoms in the pore, P h is the external pressure, B is the van der Waals constant, ' is the fractional coverage of the grain boundary by pores.
The factor F SB was derived in reference [10] for the case of small identical pores uniformly distributed over an infinitely large grain boundary, if the vacancy diffusion in the grain boundary is rate controlling: This function is of order of 1 for typical ' values of 10-20%, however it has logarithmic singularity at '→0 and cubic singularity at '→1.Moreover, applicability of equation ( 6) is unclear in the case of ensemble of different pores as well as in the case of large pores, which size is comparable with inter-pore distance or with grain face size.Therefore, for simplification it was assumed in this paper that F SB = 1.
As for the grain boundary diffusivity, considerable uncertainty still exists in the literature.It was shown in reference [8] that the best fit to the re-sintering data of references [9] and [12] is provided by the Arrhenius correlation for the diffusivity with parameters of Reynolds and Burton [13].For instance, simulations of the re-sintering conditions (24 h at 1700 °C) in the Harada and Doi test resulted in the density change of 1.15% for the normal grain and 0.175% for the large grain samples (including reduction by 0.173% for population P1 and 0.002% for P2), which should be compared with the experimental values of 1.08% and 0.19% [9].Dollins and Nichols [10], following Stehle and Assmann [3], concluded that the thermal vacancy emission alone is not sufficient to explain the healing of pores under irradiation, especially at low temperatures.This is illustrated in Figure 2 where the results are presented of simulation of porosity evolution with equation ( 4) under steady irradiation conditions in the Harada and Doi test (the line denoted as 'thermal').In these calculations, the mean pellet temperature was supposed to be equal to 1100 K, in accordance with reference [9].To overcome this difficulty, the irradiation-induced vacancy knock-out mechanism was introduced in references [3,10]: where h is the number of vacancies that escape the pore per each hit, l ∼ 1 mm is the "viable" track length of the fission fragment, and G is the fission rate.As for the key parameter, h, the authors referenced the value of 600 deduced from the fuel sputtering experiments [14]; however, they considered this value as the upper limit and set h = 100.Note for comparison that in the subsequent sputtering experiments [15] the value of h ∼ 20 was observed at typical stopping power of 20 keV/nm.The above equation ( 7) has a trivial solution: With h = 100 and typical fission rate of 10 19 m -3 s -1 , this parameter equals to ≈10 -13 m/s, so that the pores with r < 10 mm would disappear during standard LWR campaign.With the initial exponential pore size distribution, equation (3), the total fuel porosity decays exponentially: where x ¼ v rad t=r.In particular, this equation predicts decrease of the initial porosity in the typical LWR fuel by an order of magnitude at burnup of 1 GWd/t, which is considerably faster than the experimental observations.Even if to decrease the parameter h down to 20, the kinetics of fuel densification remains strongly overestimated (the curves 'Dollins & Nichols' in Fig. 2).
In particular, equation (7) does not predict the saturation of the densification process, which can be explained considering that pores with the size greater than some threshold value do not shrink, the threshold being associated with the grain size [16].Therefore, to take into account this threshold effect, the cut-off of the irradiation term was suggested in reference [8], which can be implemented in equation (7) in the smoothed form: where r pr is the pore projection radius and L edge ≈ 0.69R gr is the typical length of the grain edge 1 .With the choice h = 50, this allows reasonably reproducing not only the experimental correlation [9] for the densification kinetics, Figure 2 (curve 'modified'), but also the kinetics of the total fuel density due to both pores and inter-and intragranular fission gas bubbles, measured in references [17,18] and presented in reference [9] (Fig. 3).

Fission gas capture by pores
The initial number of gas atoms in pores per one grain can be evaluated as: where p 0 is the initial porosity, T sint and P sint are the temperature and pressure during fuel sintering.The total number of gas atoms, N rel , released from one grain during reactor campaign is equal to kbf g V gr , where k ≈ 0.
, and different variants of the irradiation term, equation (7); markers correspond to the experimental correlation [9], 'modified' corresponds to equation (10).fission gas yield per one fission, f g is the fractional release of the gas atoms to the grain boundaries and b is the burnup (number of fissions per unit volume).Therefore, one evaluates that: For the typical values T sint = 2000 K, P sint = 10 5 Pa, b = 10 27 m -3 and f g = 0.1, one evaluates this ratio as 0.6% (whereas the ratio of N 0 to the total generated gas is an order of magnitude less than this estimate).
The pores can capture the fission gas escaping from fuel grains.The capture rate is estimated multiplying the pore area, equation (1), by the gas flux density, F: where the gas flux density can be found as the time derivative of the number of gas atoms released from the grain per unit area of grain surface: In the case of constant fission rate, b (t) = Gt and thus: , whereas in the case of high burnup f g t ð Þ % f g t ð Þ.On the other hand, the gas atoms can be knocked out from pores by passing fission fragments (irradiationinduced resolution).Following Nelson's model [19], the resolution rate for intergranular pores is estimated in MFPR as [20]: where b 0 is the resolution constant, l is the average distance the ejected atom travels in pore, d is the width of the resolution layer [19].As explained in reference [21], the original Nelson model is used for intergranular porosity without modifications, suggested in reference [22] for intragranular bubbles (in order to avoid duplication of the backward flux of atoms, struck from pores, to the grain boundary).

Qualitative analysis
At the beginning of irradiation, the pores are generally under-compressed so that they tend to shrink due to both thermal and irradiation mechanisms, equations ( 4) and (10).As a result the internal gas pressure in pores increases and eventually the pressure difference dP, equation ( 5), approaches zero.Neglecting the van der Waals correction (required for small bubbles with R < 5 mm), one derives the relationship between the number of gas atoms in pore N eq and its curvature radius at equilibrium R eq : In the limiting cases one estimates the resolution term, equation ( 16), as: Therefore, one gets estimates for the rhs of the complete equation for N p in the equilibrium: where R 0 , and f 0 are the constants depending on the model parameters and external conditions: For the typical parameter values R gr = 5 mm and T = 1100 K, one evaluates that R 0 ≈ 7.6 mm, f 0 ≈ 3.2.It follows from these estimates for small pores (which quickly equilibrate so that f g << 1) that dN p =dt < 0, therefore after equilibration small pores definitely lose the gas and hence continue shrinking.The same conclusion can be drawn for large equilibrated pores; however, with less reliability in view of uncertainties of the resolution model and the relevant parameters.Therefore, one expects that pores of all sizes lose their gas after equilibration.The opposite trend cannot be excluded under some extreme conditions (high burnups at high temperature of the large grain fuels).Note that these conclusions were drawn for the equilibrated pores whereas the gas content in nonequilibrium pores can be either reducing or growing.

Quantitative analysis
The above qualitative considerations are illustrated in Figure 4 by MFPR numerical simulations of evolution of pores with initial radii of 0.1, 1 and 10 mm (curves labelled in the graph as 1, 2 and 3, respectively) under irradiation conditions of the Harada and Doi test [9]: normal grain fuel (d gr = 8 mm) with the initial porosity of 4.7% under temperature of 1100 K, pressure of 3 MPa and fission rate of 10 19 m -3 s -1 ; the initial pore distribution was approximated by equation ( 3) with r ¼ 0:35 mm.
It is seen that the relatively small pores quickly equilibrate, monotonically shrink and eventually (equilibration times are ∼6 Â 10 5 , 5 Â 10 6 and 5 Â 10 7 s for pores with the initial radii of 0.1, 0.3 and 1 mm).The large pores (≥ 3 mm) are practically unchanged in their sizes.The final fuel porosity averaged over the pore ensemble turned out to be 2.8%.
As for the gas content, it monotonically decreases in small pores (r 0.5 mm in these calculations) because the resolution mechanism dominates in the initial stage of irradiation when the small pores effectively shrink, e.g.see curve 1 in the right panel.For the larger pores, N p (t) can be non-monotonic function, however it finally decreases (at least after equilibration, see previous Sect.5), e.g.curve 2. The gas content in pores with r > 1.2 mm increases to the end of irradiation in comparison with the initial value; however these pores remain under-pressurized (curve 3).The maximum relative gas increase (∼37%) is attained in pores with the initial radius of 2 mm.As for the overall gas content in the pore ensemble, first it rapidly decreased by a factor of ∼15 at burnup of ≈3 GWd/t and then slowly the gas content was near 80% of the initial value to the end of the campaign (54 GWd/t); this change can be estimated using equation (12) as ∼0.1% of the gas released from the fuel grains.
To clarify the role of gas capture/resolution effects, the calculations were repeated with fixed N p (dashed curves in Fig. 4).It is seen on the left panel that the solid and dashed curves are very close to each other except of the stage of quick shrinking.However, at this stage the volumes of the pores are much less in comparison with the initial values and hence do not essentially contribute to the fuel porosity.As for the large pores, their sizes are practically constant, so the two approaches are close to each other too.In addition, the contribution of the largest pores to the total porosity is exponentially small (it can be evaluated by Eq. ( 9) with x ¼ r=r).These qualitative considerations were confirmed by our calculations which showed that the gas capture/ resolution effect influenced the final porosity by ≈0.01%only.
These calculations demonstrated that the pores cannot be considered as effective traps of the fission gas released from the nuclear fuel.In addition, it was justified that the neglect of gas content variation in pores is a good approximation in numerical calculations, as qualitatively discussed in reference [8].To check these conclusions, the additional calculations have been performed with the same initial pore distribution but varying one of the external parameters: grain size, fission rate or irradiation temperature.In all these cases, the pore kinetics were found to be qualitatively similar to that presented in Figure 4.
Simulations with the increased fission rate (5 Â 10 19 m -3 s -1 ) have shown that an increase of the overall gas content (following a fast initial decrease by a factor of 8) resulted in full compensation of the gas content in pores at burnup of 75 GWd/t.These variations were within 0.1% of the gas amount released from the grains, which is comparable with the above considered cases.The gas capture/resolution mechanisms were found to contribute to the final fuel porosity (which was equal to 2.33%) very similarly to the above examples.The greatest differences were found in calculations with the increased temperature.At T = 1500 K, the initial decrease of the gas content was followed (at burnup of ≈1 GWd/t) by a slight increase, which in turn followed again (at ≈4 GWd/t) by decrease up to the end of the campaign; the final content was found to be of 16% of the initial value.In addition, the simulations of the fuel volume evolution (including both pores and fission gas bubbles) under irradiation conditions of the Harada and Doi test [9] were performed (Fig. 5).The following external conditions were chosen: the mean irradiation temperature 1100 K, fission rate 10 19 m -3 s -1 , and external pressure 3 MPa.In the case of the large grain fuel, the realistic bimodal initial pore size distribution was simulated as superposition of P1 and P2 populations, see Section 2.
It was found that the pore populations P1 and P2 lost 84% and 1.4% of their initial volume respectively so that the total pore densification turned out to be of 1.07%.The final fuel density change (at the burnup of 21 GWd/t) was estimated as -1.5% and -0.9% for the normal and large grain fuels respectively.Comparison with experiment [9] demonstrates good agreement for the very initial stages of irradiation (up to ∼1 and 3 GWd/t for the large and normal grain fuels respectively).However, in later stages there is a qualitative difference in the fuel volume variation kinetics.Moreover, the similar disagreement takes place between the experimental data in reference [9] and references [17,18], also cf. Figure 3.The authors of reference [9] supposed that the discrepancy could be caused by "fuel fragment relocation at the early stage of the irradiation".This effect was not simulated in our calculations, but is foreseen in the forthcoming version of the code.Besides, for large grain samples a more detailed experimental information is required on pore size distribution in grains of different sizes (rather than available averaged data presented in Fig. 1), since the real grain size distribution is rather flat (see Ref. [9]) and this may strongly influence the threshold sizes for shrinking pores in equation (10).

Conclusions
The MFPR model for intergranular pore evolution was updated and verified against experimental data [9,17,18], and then applied to analysis of the sintering porosity behaviour under various conditions of in-pile irradiation.The performed analysis demonstrated that generally the resolution of gas atoms from pores prevails over capture of the fission gas released from grains in early stages of irradiation, which somewhat accelerates fuel densification.In later stages, the gas content of the survived large pores can increase, but not significantly.As a result, pores lose their gas content during typical reactor campaign; however the effect being vanishingly small.This implies that pores can be hardly considered as effective traps for the fission gas.
The model predicts a comparatively rapid fuel densification due to shrinkage of small pores with projection radius less or comparable with the grain face size, whereas the coarse pores remain unchanged, in agreement with numerous observations, e.g.references [6,16].In particular, this implies that the second population of relatively large pores in fuel with large grains, fabricated with a pore former, provides a rather small densification (mainly due to the first population of small pores) and thus can be hardly used for accommodation of the fuel swelling.At high burnups this can result in significant pellet-cladding mechanical interaction caused from the swelling due to retained gases and solid FPs.However, this preliminary conclusion should be thoroughly verified against additional experimental data.exp [9], large grain exp [9], normal grain exp [17,18], normal grain

Fig. 1 .
Fig.1.Initial pore size distributions in fuel samples of the Harada and Doi test[9] and their approximations.

Fig. 2 .
Fig.2.Simulation of fuel porosity under irradiation in the Harada and Doi test[9] with the thermal relaxation term, equation (4), and different variants of the irradiation term, equation(7); markers correspond to the experimental correlation[9], 'modified' corresponds to equation(10).

Fig. 3 .
Fig. 3. Kinetics of the fuel density calculated by MFPR for normal grain samples.

Fig. 4 .
Fig.4.Time dependence of pore radii and gas atom content in pores under irradiation conditions of the test[9]; dashed lines correspond to calculations with N p = const.

Gh
Nomenclature b burnup (number of fissions per unit volume) B van der Waals volume of the gas atom C (r) pore size distribution function normalized to the pore concentration fission rate (number of fissions in unit volume per unit time) mean number of vacancies knocked out from pore per one hit

Fig. 5 .'
Fig. 5. Kinetics of fuel volume change under irradiation for the large and normal grain samples.