Riso AN8
Overview
The Riso AN8 experiment was conducted at the Riso DR3 water-cooled HP1 rig and utilized an unopened rod from the Biblis A pressurized water reactor (PWR). The original fuel rod, named CB10, was irradiated over four reactor cycles up to about 44 GWd t, then underwent ramp testing in unopened conditions, hence this rod did not experience any re-fabrication process.
BISON time domain results on this page have not been updated to reflect the latest modeling and simulation formulations.
Test Description
Rod Design Specifications
The CB10 rod geometry is tabulated in Table 1, in which details regarding fuel characteristics are also reported.
Table 1: Riso AN8 Test Rod Specifications
| Fuel Rod | Measurement | Unit |
|---|---|---|
| Overall length | 0.5418 | m |
| Fuel stack height | 0.514 | m |
| Nominal plenum height | 61.0 | mm |
| Fill gas composition | He | |
| Fill gas pressure | 2.31 | MPa |
| Fuel | Measurement | Unit |
| Material | UO | |
| Enrichment | 2.95 | |
| Density | 93.74 | |
| Outer diameter | 9.053 | mm |
| Grain diameter | 6.0 | m |
| Cladding | Measurement | Unit |
| Material | Zr-4 | |
| Outer diameter | 10.81 | mm |
| Inner diameter | 9.261 | mm |
| Wall thickness | 0.772 | mm |
Operating Conditions and Irradiation History
The power history for the base irradiation carried out at the Biblis A PWR is shown in Figure 1. The experiment power history carried out at the Riso DR3 facility is shown in Figure 2. The rod axial power profile, clad surface temperature, coolant pressure and fast flux during experiments have been prescribed coherently with the data available in the International Fuel Performance Experiments database (Sartori et al., ) and in the Bump Test AN8 final report (R, 1990). Details on reactor operation parameters are tabulated in Table 2.
Table 2: Operational input parameters.
| Base Irradiation | ||
|---|---|---|
| Coolant inlet temperature | C | 284.7 |
| Coolant pressure | MPa | 15.52 |
| Fast neutron flux | n cm s per (kW m) | 3.4 |
| Power Ramps | ||
| Coolant pressure | MPa | 15.3 |
| Fast neutron flux | n cm s per (kW m) | 4.0 |

Figure 1: Base irradiation history for fuel segment CB10, carried out at Biblis A PWR.

Figure 2: Riso DR3 irradiation period for test AN8 (CB10).
Model Description
Geometry and Mesh
The fuel rod geometry was modeled for the entire simulation with a smeared column of flat ended pellets.
A 2-dimensional axisymmetric quadratic (Quad8 elements) mesh was used to model the geometry of the rod used in the AN8 experiment. The mesh script files used to generate the mesh of the considered rod are included in the rod analysis subfolder.
Input files
The BISON input and all supporting files (power histories, axial power profile, fast neutron flux history, etc.) for this case are provided with the code distribution at bison/assessment/LWR/validation/ and bison/assessment/LWR/validation/Riso_AN8/analysis.
To avoid code duplication, the input files are structured in this format: A first base input file contains characteristics common to all the Riso cases: Riso_Base.i. A second base input file contains characteristics common to all the Riso case that do not use the action: Riso_Base_sub.i. Input files containing information specific to the fuel rod and the type of problem solving use the !include function to build a complete input file with the base files.
To run a specific assessment, such as the Riso AN8, run: Riso_AN8.i.
Material and Behavioral Models
The following material and behavioral models for UO fuel were used:
UO2Thermal - NFIR: NFIR model for temperature and burnup dependent thermal properties
ComputeFiniteStrainElasticStress and UO2ElasticityTensor: elastic mechanical behavior
UO2RelocationEigenstrain: relocation strains, relocation activation threshold power set to 5 kW/m
ComputeThermalExpansionEigenstrain: thermal expansion with a constant instantaneous thermal expansion coefficient
UO2VolumetricSwellingEigenstrain : volumetric expansion due to solid and gaseous swelling
UO2Sifgrs: fission gas release model used with the gaseous swelling model
UO2VolumetricSwellingEigenstrain(Pastore et al., 2015)
For the cladding material, a constant thermal conductivity of 16 W/m-K was used and both thermal and irradiation creep were considered using the Limback model (Limbäck and Andersson, 1996). The following material and thermal behavior models were used for the Zircaloy-4 cladding:
HeatConductionMaterial: Thermophysical material properties
ZryCreepLimbackHoppeUpdate and ComputeIsotropicElasticityTensor: mechanical creep and elastic deformation behavior for Zircaloy-4, with a Young's modulus of 1.010 and a Poisson's ratio value of 0.3
ZryIrradiationGrowthEigenstrain: ESCORE model for volumetric swelling due to irradiation exposure
ZryThermalExpansionMATPROEigenstrain: thermal expansion with an instantaneous thermal expansion coefficient that varies according to the MATPRO cthexp model
ZryOxidation: cladding oxidation thickness growth model for metal-water reactions
Details and references for all of these models listed above can be found on the linked BISON documentation pages.
Boundary and Operating Conditions
The power history considered comprehends the base irradiation, shown in Figure 1, and the subsequent ramp test, reported in Figure 2. The temperature boundary conditions as a function of time was prescribed according to the files available in the IFPE database (Sartori et al., ).
Results Comparison
The Riso AN8 experiment is used to assess the code's capability to predict integral fission gas release (FGR) at EOL, as well as to predict cladding deformation after the ramp test.

Figure 3: BISON total fission gas release comparison to Riso experimental data.

Figure 4: BISON cladding diameter results compared to measured data.
Fission Gas Release
The calculated integral fuel rod FGR is compared to the final, measured data in Figure 3. BISON results show a good agreement with the experimental data, in line with the intrinsic modeling uncertainties (Pastore et al., 2015).
EOL Cladding Diameter
The BISON results of cladding outer diameter at the end of irradiation as a function of axial position are presented in Figure 4 and compared to the experimental measurements. The figure shows an acceptable agreement between simulated and measured cladding diameter at EOL. The deviation of BISON calculations appear to be in line with state-of-the-art fuel performance codes accuracy (IAEA, 2008-2012). The systematic over-estimation of the cladding diameter may be reduced if fuel creep is considered in the simulations. This issue will be further investigated in future work.
References
- IAEA.
Improvement of Computer Codes Used for Fuel Behaviour Simulation (FUMEX-III): Report of a Coordinated Research Project 2008-2012.
Technical Report IAEA-TECDOC-1697, International Atomic Energy Agency, 2008-2012.[BibTeX]
- M. Limbäck and T. Andersson.
A model for analysis of the effect of final annealing on the in- and out-of-reactor creep behavior of zircaloy cladding.
In Zirconium in the Nuclear Industry: Eleventh International Symposium, ASTM STP 1295, 448–468. 1996.[BibTeX]
- G. Pastore, L.P. Swiler, J.D. Hales, S.R. Novascone, D.M. Perez, B.W. Spencer, L. Luzzi, P. Van Uffelen, and R.L. Williamson.
Uncertainty and sensitivity analysis of fission gas behavior in engineering-scale fuel modeling.
Journal of Nuclear Materials, 465:398–408, 2015.[BibTeX]
- Giovanni Pastore, L.P. Swiler, J.D. Hales, S.R. Novascone, D.M. Perez, B.W. Spencer, L. Luzzi, P. Van Uffelen, and R.L. Williamson.
Uncertainty and sensitivity analysis of fission gas behavior in engineering-scale fuel modeling.
Journal of Nuclear Materials, 456:398–408, 2015.
doi:10.1016/j.jnucmat.2014.09.077.[BibTeX]
- R.
The Third RisøFission Gas Project: Bump Test AN8 (CB10).
Technical Report Risø-FGP3-AN8, Risø, september 1990.[BibTeX]
- E. Sartori, J. Killeen, and J. A. Turnbull.
International Fuel Performance Experiments (IFPE) Database.
OECD-NEA, 2010, available at http://www.oecd-nea.org/science/fuel/ifpelst.html.[BibTeX]