High Burnup Effects Programme (HBEP) A1/8-4, A3/6-4 and H8/36-4

Overview

The purpose of the High Burnup Effects Programme (HBEP) Task 3 was to provide high burnup effects evaluations, fission gas sampling, and parameter effects (HBEP, 2002). Rods A1/8-4, A3/6-4 and H8/36-4 were selected for the assessment of BISON code; they were irradiated to 48.5 MWd/kgU, 47.8 MWd/kgU and 46.6 MWd/kgU respectively, in the TVO-1 boiling water reactor (BWR) (HBEP, 2002).

Test Description

The three rods in this series, manufactured by ASEA, have similar design characteristics but slightly different pellet and cladding diameter. The fuel stack consisted of solid pellets with an active length of 3.68 meters. This portion of the HBEP experiments designated as Task 3 was to provide well-characterised data on the evaluation of high burnup effects on different fuel characteristics with an emphasis on fission gas release.

Table 1: HBEP Test Rod Specifications (A1/8-4, A3/6-4, and H8/36-4)

Fuel RodMeasurementUnit
Fuel stack height3.68m
Nominal plenum height
A1/8-4246mm
A3/6-4246mm
H8/36-4244mm
Number of pellets per rod337
Pellet Height
A1/8-410.9mm
A3/6-410.9mm
H8/36-410.9mm
Fill gas compositionHe
Fill gas pressure.37MPa
FuelMeasurementUnit
MaterialUO
Enrichment
A1/8-42.33%
A3/6-43.08%
H8/36-41.36%
Density95.6%
Outer diameter
A1/8-49.939mm
A3/6-410.439mm
H8/36-49.939mm
Grain diameter
All Rods6.6m
Fuel densification
All Rods.43%
Roughness
All Rods1.1m
CladdingMeasurementUnit
MaterialZr-2
Outer diameter
A1/8-411.73mm
A3/6-412.25mm
H8/36-411.73mm
Inner diameter
A1/8-410.14mm
A3/6-410.65mm
H8/36-410.14mm
Roughness
All Rods.35m

Operating Conditions and Irradiation History

Rods A1/8-4, A3/6-4 and H8/36-4 were irradiated in the TVO-1 commercial BWR Cycles 2 through 7 (HBEP, 2002). Reactor operating coditions were, coolant pressure of 7 MPa and coolant inlet temperature of 552 K (HBEP, 2002). Peak-in-life LHGR values occurred during the first and third operating cycles, with maximum rod-average LHGR values being approximately 18 to 27 KW/m. Rod-average burnups ranged from 44 to 50 MWd/kgM with rod-average FGR values ranging from 0.3 to 17.1%. Clad temperatures and local power histories were taken at 12 axial locations and obtained from the IFPE database. The plot below, Figure 1, shows the power histories for the A1/8-4, A3/6-4, and H8/36-4 rods.

Figure 1: The linear heat rates for A1/8-4, A3/6-4 and H8/36-4.

Geometry and Mesh

All three fuel rods were meshed using 2-D axisymmetric models with eight-node quadratic elements. For simplicity, the pellet stack was modeled as a single continuous fuel column (smeared mesh). The meshes were generated using BISON's automatic mesh generator. The fuel pellets had 200 axial elements and 12 radial elements, and the cladding consisted of 200 axial elements and 4 radial elements.

Figure 2: A BISON generated mesh showing the fuel and cladding mesh, axially minimized by 100 for visual purposes.

Material and Behavioral Models

The following material and behavioral models were used for the fuel:

  • UO2Sifgrs: fission gas release model used with the gaseous swelling model UO2VolumetricSwellingEigenstrain

For the clad material, a constant thermal conductivity of 16 W/m-K was used and both thermal (primary and secondary) and irradiation creep were considered using the Limback creep model (Limbäck and Andersson, 1996). The following material and thermal behavior models were used for the cladding:

Details and references for all of these models listed above can be found on the linked BISON documentation pages.

Input files

The BISON input and all supporting files (power histories, axial power profile) are provided with the code distribution at bison/assessment/LWR/validation/HBEP/analysis.

Results

HBEP rods A1/8-4, A3/6-4, and H8/36-4 were designed and built to utilize well-characterized data on the effects of fuel temperature, burnup, power history, and different fuel characteristics with an emphasis on FGR. To achieve this the rods were built to the exact specification as described in the report. These rods were irradiated in the TVO-1 reactor for 6 cycles, and were then taken to a hot cell for fission gas characterization and diametral measurements.

Figure 3 and Figure 4 compare BISON calculations for both fission gas release and cladding diametral changes following irradiation. BISON's fission gas calculations slightly over predict for two of the cases and dramatically under predict for the H8/36-4 case. Cladding diametal comparisons vary from 12-55 microns.

Figure 3: Comparison of measured fission gas release and the BISON calculations

Figure 4: Comparison of measured cladding displacements and the BISON calculation; measured diameterchange was averaged over the region 500 to3500 mm above the bottom of the fuel rod

Table 2 makes comparisons between BISON and measured fission gas, plenum pressure, clad elongation, free gas volume, and burnup. Table 3 and Table 4 compare BISON pellet-cladding gap calculations to the corisponding axilly measured pellet-clad gap values. Furthermore, Table 4 compares pellet OD's between experimental measurements and BISON calculations to better determine the differences between Table 3 pellet-clad gap comparisons.

Table 2: BISON Comparisons to Experimental Results

A1/8-4ExperimentBISONUnits
Fission Gas Release.30.54%
Plenum Pressure.48.7MPa
Clad Elongation3.9663.942m
Free Gas Volume27.921.1cm
Burnup48.545.2MWd/kgU
A3/6-4ExperimentBISONUnits
Fission Gas Release11.68%
Plenum Pressure.64.92MPa
Clad Elongation3.9673.942m
Free Gas Volume28.522.8cm
Burnup47.844.5MWd/kgU
H8/36-4ExperimentBISONUnits
Fission Gas Release17.3.27%
Plenum Pressure2.86.630MPa
Clad Elongation3.9643.938m
Free Gas Volume29.221.64cm
Burnup46.643.2MWd/kgU

Table 3: BISON Comparisons to Gap Experimental Measurements

H8/36-4
Axial Location (cm)Avg Exp Gap (microns)BISON (microns)
220.52414.3
1702019.2
48.24846.2
234855.3
A3/6-4
Axial Location (cm)Avg. Exp Gap (microns)BISON (microns)
1572513

Table 4: BISON Comparisons Pellet OD to Experimental Measurements

H8/36-4
Axial Location (cm)Exp (mm)Exp (mm)BISON (mm)
220.510.06710.11610.058
17010.10510.12110.049
48.210.01410.06610.002
2310.06810.079.988
A3/6-4
Axial Location (cm)Exp (mm)Exp (mm)BISON (mm)
157.510.53610.55410.577

Figure 5 shows the strong axial power dependence at times corresponding to axial control blade movements.

Figure 5: Axial Power Profile for H8/36-4 at selected times during the control blade insertion

Discussion

Modeling these three rods was useful to the BISON development to validate the capability of modeling high burnup fuel rods. The HBEP report (HBEP, 2002) stated, fission gas release was related to the rod location to the control blades. Rods closer to the control blades had significantly higher fission gas release than rods further away. Figure 5 shows, there is significant power changes caused by the blade movement, which could result in high fission gas release. BISON code appears to under calculate FGR under such conditions. Rods A1/8-4 and A3/6-4 were not located near a control blade.

During the PIE the cladding diameter was measured, and this value is the average value from 500-3500 mm. Furthermore, the cladding diameter does include the oxide thickness in the measurement. BISON's calculations are with in .07 to .02 mm of the measured value. There is no oxidation model in BISON to model the oxide growth of BWR fuel cladding, and this could also have implication on the calculations.

BISON calculated a lower free gas volume for all cases, and calulated a higher plenum pressure for A1/8-4, and A3/6-4 and a lower plenum pressure for H8/36-4. The difference between the calculated and measurement data on plenum pressure could be contributed to the difference in fission gas release. As for cladding elongation, BISON calculation is lower on all accounts. BISON currently does not consider the plenum spring or friction between the fuel and cladding. Both areg driving mechanisms for cladding elongation.

References

  1. HBEP. Summary of the High Burn-up Effects Programme as abstracted from the Programme Final Report. Technical Report, High Burn-up Effects Programme, December 2002.[BibTeX]
  2. 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]