skip navigation links 
 
 Search Options 
Index | Site Map | FAQ | Facility Info | Reading Rm | New | Help | Glossary | Contact Us blue spacer  
secondary page banner Return to NRC Home Page

Regulatory Guide 3.54 - Spent Fuel Heat Generation in an Independent Spent Fuel Storage Installation


(Draft was issued as DG-3010)

Revision 1
January 1999

Availability Notice


A. Introduction

In 10 CFR Part 72, "Licensing Requirements for the Independent Storage of Spent Nuclear Fuel and High-Level Radioactive Waste," paragraph (h)(1) of Section 72.122, "Overall Requirements," requires that spent fuel cladding be protected during storage against degradation that leads to gross ruptures or the fuel must be otherwise confined such that degradation of the fuel during storage will not pose operational safety problems with respect to its removal from storage. It has been shown that, under certain environmental conditions, high storage temperatures can cause degradation and gross rupture of the fuel rods to occur very rapidly. It is necessary to know what storage temperatures are anticipated during the life of the storage installation and that these temperatures will not significantly degrade the cladding to a point that causes gross ruptures. The temperature in an independent spent fuel storage installation is a function of the heat generated by the stored fuel assemblies. The spent fuel storage system is required by 10 CFR 72.128(a)(4) to be designed with a heat removal capability consistent with its importance to safety.

This regulatory guide presents a method acceptable to the Nuclear Regulatory Commission (NRC) staff for calculating heat generation rates for use as design input for an independent spent fuel storage installation. The original guide, issued in September 1984, was based on validated analyses performed for pressurized-water reactors (PWRs), and boiling-water reactors (BWRs) were considered only as a simple conservative extension of the PWR data base. In this revision, the procedure for determining heat generation rates for both PWRs and BWRs is based on analyses of each reactor type using calculational methods that have been validated against measured heat generation data from PWR and BWR assemblies. This revision presents a methodology that is simpler and is therefore expected to be more useful to applicants and reviewers.

This regulatory guide contains no information collection requirements and therefore is not subject to the requirements of the Paperwork Reduction Act of 1980 (44 U.S.C. 3501 et seq.).

B. Discussion

The methodology of NUREG/CR-5625(1) is appropriate for computing the heat generation rates of fuel assemblies from light-water-cooled power reactors as a function of burnup, specific power, and decay time. The computed heat generation results are used in the next section in a procedure for determining heat generation rates for PWR and BWR assemblies.

Calculations of decay heat have been verified by comparison with the existing data base of experimentally measured decay heat rates for PWR and BWR spent fuel. The range of parameter values in the procedure is considered to lie in the mainstream of typical burnup, specific power, enrichment, and cooling time. A detailed example is shown in Appendix A.

The following terms and units have been used in this guide.

Terms and Units Used in Guide

Term/Unit Definition
Be - burnup in last cycle, MWd/kgU
Be-1 - burnup in next-to-last cycle, MWd/kgU
Bi - fuel burnup increase for cycle i, MWd/kgU
Btot - total burnup of discharged fuel, MWd/kgU
Es - initial fuel enrichment, wt-% 235U
P - specific power of fuel as in Equations 2 and 3, kW/kgU
Pave - average cumulative specific power during 80% uptime, kW/kgU
Pave,e-1 - average cumulative specific power (at 80%) through cycle e-1, the next-to-last cycle
Pe - fuel-specific power during the last cycle e
Pe-1 - fuel-specific power during cycle e-1, the next-to-last cycle
PL, PH - lower and higher values of specific power that bracket the specific power value of Pave
Ptab - heat generation rate that is obtained from the table by interpolation between the lower and higher bracketing values
Pfinal - final heat generation rate determined by applying all adjustment factors, followed by the safety factor to the value Ptab
TL, TH - lower and higher time values in a table that bracket the cooling time of interest, TC
S - percentage safety factor applied to decay heat rates, ptab
Tc - cooling time of an assembly, in years
Te - cycle time of last cycle before discharge, in days
Te-1 - cycle time of next-to-last cycle, in days
Ti - cycle time of ith reactor operating cycle, including downtime for all but last cycle of assembly history, in days
Tres - reactor residence time of assembly, from first loading to shutdown for discharge, in days
f7 - last-cycle short cooling time modification factor
f'7 - next-to-last cycle short cooling time factor
fe - 235U initial enrichment modification factor
fp - excess power adjustment factor
p - heat generation rate of spent fuel assembly, W/kgU

C. Regulatory Position

The following method for determining heat generation rates of reactor spent fuel assemblies is acceptable to the NRC staff. There may be fuel assemblies with characteristics that are sufficiently outside the mainstream of typical operations that they need a separate computation of the heat generation rate. A discussion of the characteristics of assumed typical reactor operations is given in Appendix B.

The first part of this section contains the definitions and derivations, as used in this guide, of parameters needed in the determination of the heat generation rate of a fuel assembly. The second part contains the procedure used in deriving the final heat rate of an assembly. Although allowance has been made to use simple adjustment factors for cases that are somewhat atypical, many cases will probably not require any adjustment of the table heat rate other than the safety factor.

Heat generation rate tables for actinides, fission products, and light elements are given in Appendix C for informational purposes only. They are not used directly in this guide's method for determining heat rates.

1. Definitions and Derivations of Parameters

The following definitions and derivations of parameters of the spent fuel assembly are used in the procedure in this guide.

1.1 Heat Generation Rate (p)

The heat generation rate of the spent fuel assembly is the recoverable thermal energy (from radioactive decay) of the assembly per unit time per unit fuel mass. The units for heat generation rate used in this guide are watts per kilogram U (W/kgU), where U is the initial uranium loaded. Heat generation rate has also been referred to as decay heat rate, afterheat, or afterheat power.

1.2 Cycle and Cycle Times (Ti)

A cycle of the operating history for a fuel assembly is, with one exception, the duration between the time criticality is obtained for the initially loaded or reloaded reactor to the time at which the next reloaded core becomes critical. The exception is for the last cycle, in which the cycle ends with the last reactor shutdown before discharge of the assembly. Ti denotes the elapsed time during cycle i for the assembly. Specifically, the first and last cycles are denoted by i  = s (for start) and i = e (for end), respectively. Tres, the total residence time of the assembly, is the sum of all Ti for i = s through e, inclusive. Except for the last cycle for an assembly, the cycle times include the downtimes during reload. Cycle times, in this guide, are in days.

1.3 Fuel Burnup of the Assembly (Bi and Btot)

The fuel burnup of cycle i, Bi, is the recoverable thermal energy per unit fuel mass during the cycle in units of megawatt days per metric ton (tonne) initial uranium (MWd/tU), or in the SI units(2) of mass used in this guide, megawatt days per kilogram U (MWd/ kgU). Bi is the best maximum estimate of the fuel assembly burnup during cycle i. Btot is the total operating history burnup:

Equation 1 - B sub tot shown as the sum over i of B sub i

1.4 Specific Power of the Fuel (Pi, Pe, Pave, and Pave,e-1)

Specific power has a unique meaning in this guide. The reason for developing this definition is to take into account the differences between the actual operating history of the assembly and that used in the computation of the tabulated heat generation rates. The calculational model applied an uptime (time at power) of 80% of the cycle time in all except the last cycle (of the discharged fuel assembly), which had no downtime. The definition of specific power, used here, has two basic characteristics. First, when the actual uptime experienced by the assembly exceeds the 80% applied in the SAS2H/ORIGEN-S calculations, the heat rate derived by the guide procedure maintains equivalent accuracies within 1%. Second, when the actual uptime experienced is lower than the 80% applied in the calculations, the heat rate is reduced. The technical basis for these characteristics is presented in NUREG/CR-5625.

The specific power of cycle i, or e (last cycle), in kW/kgU, using burnup in MWd/kgU, is determined by

Equation 2 - showing specific power of cycle for i less than e and for i equal to e

The average specific power over the entire operating history of a fuel assembly, using the same units as in Equation 2, is determined by:

Equation 3 - average specific power P sub ave

The average specific power through the next-to-last cycle is used in applying the adjustment factor for short cooling time (see Regulatory Position 2.2). This parameter is determined by:

Equation 4 - expression for P sub ave sub (e-1)

Note that Btot and Pave, as derived here, are used in determining the heat generation rate with this guide. Also, for cooling times 7 years, Pe is used in an adjustment formula. The method applied here accommodates storage of a fuel assembly outside the reactor during one or two cycles and returning it to the reactor. Then, Bi = 0 may be set for all intermediate storage cycles. If the cooling time is short (i.e., <10 years), the results derived here may be excessively high for cases in which the fuel was temporarily discharged. Other evaluation methods that include the incorporation of storage cycles in the power history may be preferable.

1.5 Assembly Cooling Time (Tc)

The cooling time, Tc, of an assembly is the time elapsed from the last downtime of the reactor prior to its discharge (at end of Te) to the time at which the heat generation rate is desired. Cooling times, in this guide, are in years.

1.6 Assembly Initial Fuel Enrichment (Es)

The initial enrichment, Es, of the fuel assembly is considered to be the average weight percent 235U in the uranium when it is first loaded into the reactor. Heat generation rates vary with initial enrichment for fuel having the same burnup and specific power; the heat rate increases with lower enrichment. If the enrichment is different from that used in the calculations at a given burnup and specific power, a correction factor is applied.

2. DETERMINATION OF HEAT GENERATION RATES

Directions for determining the heat generation rates of light-water-reactor (LWR) fuel assemblies from Tables 1 through 8 are given in this section. First, a heat rate, ptab, is found by interpolation from Tables 1 through 3 or Tables 5 through 7. Next, a safety factor and all the necessary adjustment factors are applied to determine the final heat generation rate, pfinal . There are three adjustment factors (see Regulatory Positions 2.2 to 2.4) plus a safety factor (see Regulatory Position 2.5) that are applied in computing the final heat generation rate, pfinal, from ptab. In many cases, the adjustment factors are unity and thus are not needed. An alternative to these directions is the use of the light-water-reactor afterheat rate calculation (LWRARC) code on a personal computer; the code is referred to in Regulatory Position 2.7. This code evaluates ptab and pfinal using the data and procedures established in this guide.

Table 1 : BWR Spent Fuel Heat Generation Rates, Watts Per Kilogram U, for Specific Power = 12 kW/kgU

Cooling Time, Years Fuel Burnup, MWd/kgU
20 25 30 35 40 45
1.0 4.147 4.676 5.121 5.609 6.064 6.531
1.4 3.132 3.574 3.955 4.370 4.760 5.163
2.0 2.249 2.610 2.933 3.281 3.616 3.960
2.8 1.592 1.893 2.174 2.472 2.764 3.065
4.0 1.111 1.363 1.608 1.865 2.121 2.384
5.0 0.919 1.146 1.371 1.606 1.844 2.087
7.0 0.745 0.943 1.142 1.349 1.562 1.778
10.0 0.645 0.819 0.996 1.180 1.369 1.561
15.0 0.569 0.721 0.876 1.037 1.202 1.370
20.0 0.518 0.656 0.795 0.940 1.088 1.240
25.0 0.477 0.603 0.729 0.861 0.995 1.132
30.0 0.441 0.556 0.672 0.792 0.914 1.039
40.0 0.380 0.478 0.576 0.678 0.781 0.886
50.0 0.331 0.416 0.499 0.587 0.674 0.764
60.0 0.292 0.365 0.438 0.513 0.589 0.666
70.0 0.259 0.324 0.387 0.454 0.520 0.587
80.0 0.233 0.291 0.347 0.405 0.464 0.523
90.0 0.212 0.263 0.313 0.365 0.418 0.470
100.0 0.194 0.241 0.286 0.333 0.380 0.427
110.0 0.179 0.222 0.263 0.306 0.348 0.391

Table 2 : BWR Spent Fuel Heat Generation Rates, Watts Per Kilogram U, for Specific Power = 20 kW/kgU

Cooling Time, Years Fuel Burnup, MWd/kgU
20 25 30 35 40 45
1.0 5.548 6.266 6.841 7.455 8.000 8.571
1.4 4.097 4.687 5.173 5.690 6.159 6.647
2.0 2.853 3.316 3.718 4.142 4.540 4.950
2.8 1.929 2.296 2.631 2.982 3.324 3.673
4.0 1.262 1.549 1.827 2.117 2.410 2.705
5.0 1.001 1.251 1.501 1.760 2.024 2.292
7.0 0.776 0.985 1.199 1.420 1.650 1.882
10.0 0.658 0.838 1.023 1.215 1.413 1.616
15.0 0.576 0.731 0.890 1.056 1.227 1.403
20.0 0.523 0.663 0.805 0.954 1.107 1.263
25.0 0.480 0.608 0.737 0.871 1.009 1.150
30.0 0.444 0.560 0.678 0.800 0.925 1.053
40.0 0.382 0.481 0.579 0.682 0.786 0.893
50.0 0.332 0.417 0.501 0.588 0.677 0.767
60.0 0.292 0.365 0.438 0.513 0.589 0.666
70.0 0.259 0.324 0.386 0.452 0.518 0.585
80.0 0.233 0.290 0.345 0.403 0.460 0.519
90.0 0.211 0.262 0.311 0.362 0.413 0.465
100.0 0.193 0.239 0.283 0.329 0.375 0.421
110.0 0.178 0.220 0.260 0.302 0.343 0.385

Table 3 : BWR Spent Fuel Heat Generation Rates, Watts Per Kilogram U, for Specific Power = 30 kW/kgU

Cooling Time, Years Fuel Burnup, MWd/kgU
20 25 30 35 40 45
1.0 6.809 7.786 8.551 9.337 10.010 10.706
1.4 4.939 5.721 6.357 7.006 7.579 8.169
2.0 3.368 3.958 4.463 4.979 5.453 5.938
2.8 2.211 2.651 3.050 3.460 3.855 4.256
4.0 1.381 1.705 2.016 2.339 2.663 2.991
5.0 1.063 1.335 1.605 1.885 2.172 2.462
7.0 0.797 1.015 1.239 1.471 1.713 1.958
10.0 0.666 0.850 1.039 1.237 1.443 1.653
15.0 0.579 0.737 0.898 1.067 1.242 1.422
20.0 0.525 0.667 0.811 0.962 1.117 1.276
25.0 0.482 0.611 0.741 0.877 1.017 1.160
30.0 0.445 0.563 0.681 0.805 0.931 1.061
40.0 0.382 0.482 0.581 0.685 0.790 0.898
50.0 0.332 0.418 0.502 0.589 0.678 0.769
60.0 0.292 0.366 0.438 0.513 0.589 0.666
70.0 0.259 0.323 0.386 0.451 0.517 0.584
80.0 0.232 0.289 0.344 0.401 0.459 0.517
90.0 0.210 0.261 0.310 0.361 0.411 0.463
100.0 0.192 0.238 0.282 0.327 0.372 0.418
110.0 0.177 0.219 0.259 0.300 0.340 0.382

Table 4 : BWR Enrichments for Burnups in Tables

Fuel Burnup, MWd/kgU Average Initial Enrichment, wt-% U-235
20 1.9
25 2.3
30 2.7
35 3.1
40 3.4
45 3.8

Table 5 : PWR Spent Fuel Heat Generation Rates, Watts Per Kilogram U, for Specific Power = 18 kW/kgU

Cooling Time, Years Fuel Burnup, MWd/kgU
25 30 35 40 45 50
1.0 5.946 6.574 7.086 7.662 8.176 8.773
1.4 4.485 5.009 5.448 5.938 6.382 6.894
2.0 3.208 3.632 4.004 4.411 4.793 5.223
2.8 2.253 2.601 2.921 3.263 3.595 3.962
4.0 1.551 1.835 2.108 2.398 2.685 2.997
5.0 1.268 1.520 1.769 2.030 2.294 2.576
7.0 1.008 1.223 1.439 1.666 1.897 2.143
10.0 0.858 1.044 1.232 1.430 1.633 1.847
15.0 0.744 0.905 1.068 1.239 1.414 1.599
20.0 0.672 0.816 0.963 1.116 1.272 1.437
25.0 0.615 0.746 0.879 1.018 1.159 1.308
30.0 0.566 0.686 0.808 0.934 1.063 1.197
40.0 0.487 0.588 0.690 0.797 0.904 1.017
50.0 0.423 0.510 0.597 0.688 0.780 0.875
60.0 0.372 0.447 0.522 0.601 0.680 0.762
70.0 0.330 0.396 0.462 0.530 0.599 0.670
80.0 0.296 0.355 0.413 0.473 0.534 0.596
90.0 0.268 0.321 0.372 0.426 0.480 0.536
100.0 0.245 0.293 0.339 0.387 0.436 0.486
110.0 0.226 0.270 0.312 0.356 0.399 0.445

Table 6 : PWR Spent Fuel Heat Generation Rates, Watts Per Kilogram U, for Specific Power = 28 kW/kgU

Cooling Time, Years Fuel Burnup, MWd/kgU
25 30 35 40 45 50
1.0 7.559 8.390 9.055 9.776 10.400 11.120
1.4 5.593 6.273 6.836 7.441 7.978 8.593
2.0 3.900 4.432 4.894 5.385 5.838 6.346
2.8 2.641 3.054 3.435 3.835 4.220 4.642
4.0 1.724 2.043 2.352 2.675 2.999 3.346
5.0 1.363 1.637 1.911 2.195 2.486 2.793
7.0 1.045 1.271 1.500 1.740 1.987 2.248
10.0 0.873 1.064 1.261 1.465 1.677 1.900
15.0 0.752 0.915 1.083 1.257 1.438 1.627
20.0 0.677 0.823 0.973 1.128 1.289 1.457
25.0 0.619 0.751 0.886 1.027 1.171 1.322
30.0 0.569 0.690 0.813 0.941 1.072 1.208
40.0 0.488 0.590 0.693 0.800 0.909 1.023
50.0 0.424 0.511 0.599 0.689 0.782 0.877
60.0 0.372 0.447 0.523 0.601 0.680 0.762
70.0 0.330 0.396 0.461 0.529 0.598 0.668
80.0 0.295 0.354 0.411 0.471 0.531 0.593
90.0 0.267 0.319 0.371 0.424 0.477 0.531
100.0 0.244 0.291 0.337 0.385 0.432 0.481
110.0 0.225 0.268 0.310 0.352 0.396 0.440

Table 7 : PWR Spent Fuel Heat Generation Rates, Watts Per Kilogram U, for Specific Power = 40 kW/kgU

Cooling Time, Years Fuel Burnup, MWd/kgU
25 30 35 40 45 50
1.0 8.946 10.050 10.900 11.820 12.580 13.466
1.4 6.514 7.400 8.111 8.863 9.514 10.254
2.0 4.462 5.129 5.692 6.284 6.821 7.418
2.8 2.947 3.441 3.884 4.346 4.787 5.267
4.0 1.853 2.212 2.554 2.910 3.265 3.647
5.0 1.429 1.728 2.021 2.327 2.639 2.970
7.0 1.067 1.304 1.543 1.793 2.052 2.325
10.0 0.881 1.078 1.278 1.488 1.705 1.936
15.0 0.754 0.921 1.091 1.268 1.452 1.645
20.0 0.678 0.827 0.978 1.136 1.298 1.469
25.0 0.619 0.754 0.890 1.032 1.178 1.331
30.0 0.570 0.693 0.816 0.945 1.077 1.215
40.0 0.488 0.592 0.695 0.803 0.912 1.026
50.0 0.423 0.512 0.599 0.691 0.783 0.879
60.0 0.371 0.448 0.522 0.601 0.680 0.762
70.0 0.329 0.396 0.461 0.529 0.597 0.668
80.0 0.294 0.353 0.410 0.470 0.530 0.592
90.0 0.266 0.319 0.369 0.422 0.475 0.530
100.0 0.243 0.290 0.336 0.383 0.430 0.479
110.0 0.224 0.267 0.308 0.351 0.393 0.437

Table 8 : PWR Enrichments for Burnups in Tables

Fuel Burnup, MWd/kgU Average Initial Enrichment, wt-% U-235
25 2.4
30 2.8
35 3.2
40 3.6
45 3.9
50 4.2

2.1 Computing Heat Rate Provided by Tables

Tables 1 through 3 are for BWR fuel, and Tables 5 through 7 are for PWR fuel. The heat rates in each table pertain to a single average specific power and are listed as a function of total burnup and cooling time. After determining Pave, Btot, and Tc as above, select the next lower (L-index) and next higher (H-index) heat rate values from the tables so that:

PL Pave PH

BL Btot BH

and

TL Tc TH

Compute ptab, the heat generation rate, at Pave, Btot, and Tc, by proper interpolation between the tabulated values of heat rates at the lower and higher parameter limits. A linear interpolation should be used between heat rates for either burnup or specific power interpolations. In computing the heat rate at Tc, the interpolation should be logarithmic in heat rate and linear in cooling time. Specifically, the interpolation formulas for interpolating in specific power, burnup, and cooling time are, respectively,

Equation 5 - p for interpolating in specific power

Equation 6 - p for interpolating in burnup

Equation 7 - p for interpolating in cooling time

where pL and pH represent the tabulated or interpolated heat rates at the appropriate parameter limits corresponding to the L and H index. If applied in the sequence given above, Equation 5 would need to be used four times to obtain p values that correspond to BL and BH at values of TL and TH. A mini-table of four p values at Pave is now available to interpolate burnup and cooling time. Equation 6 would then be applied to obtain two values of p at TL and TH. One final interpolation of these two p values (at Pave and Btot) using Equation 7 is needed to calculate the final ptab value corresponding to Pave, Btot, and Tc. The optional Lagrangian interpolation scheme offered by the LWRARC code is also considered an acceptable method for interpolating the decay heat data.

If Pave or Btot falls below the minimum table value range, the minimum table-specific power or burnup, respectively, may be used conservatively. If Pave exceeds the maximum table value, the table with the maximum specific power (Table 3 for BWR fuel and Table 7 for PWR fuel) may be used in addition to the adjustment factor, fp, described in Regulatory Position 2.3.

The tables should not be applied if Btot exceeds the maximum burnup in the tables, or if Tc is less than the minimum (1 year). If Tc exceeds the maximum (110 years) cooling time of the tables, the 110-year value is acceptable, although it may be too conservative.

2.2 Short Cooling Time Factors f7 and f7

The heat rates presented in Tables 1 through 3 and Tables 5 through 7 were computed from operating histories in which a constant specific power and an uptime of 80% of the cycle time were applied. Expected variations from these assumptions cause only minor changes (<1%) in decay heat rates beyond approximately 7 years of cooling. However, if the specific power near the end of the operating history is significantly different from the average specific power, Pave, ptab needs to be adjusted if Tc 7. The ratios Pe/Pave and Pe-1/Pave,e-1 are, respectively, used to determine the adjustment factors f7 and f7. The factors reduce the heat rate ptab if the corresponding ratio is less than 1 and increase the heat rate ptab if the corresponding ratio is greater than 1. The formulas for the factors are below.

f7 = 1 when Tc>7 years or e = s (i.e., 1 cycle only) (Equation 8)
f7 = 1 + 0.35 R/ square root of T sub c when 0 < R < 0.3
f7 = 1 + 0.25 R/Tc when -0.3 < R < 0
f7 = 1 - 0.075/Tc when R < -0.3

where

Equation 9 - expression for R

f7 = 1 when Tc>7 years or e < 3 (Equation 10)
f7 = 1 + 0.10 R'/square root of T sub c when 0 < R' < 0.6
f7 = 1 + 0.08 R'/Tc when -0.5 < R' < 0
f7 = 1 - 0.04/Tc when R'<-0.5

where

Equation 11 - expression for R'

It can be observed that there are upper limits to R and R in Equations 8 and 10. It is recommended not to use the decay heat values of this guide if any of the following conditions occur:

if Tc 10 years and Pe /Pave > 1.3,
if 10 years < Tc 15 years and Pe /Pave > 1.7,
if Tc 10 years and Pe-1 /Pave,e-1 > 1.6.

Although it is safe to use the procedures in this guide, the heat rate values for pfinal may be excessively high when

Tc 7 years and Pe /Pave < 0.6,
Tc 7 years and Pe-1 /Pave,e-1 < 0.4.

2.3 The Excess Power Adjustment Factor fp

The maximum specific power, Pmax, used to generate the data in Tables 1 through 3 and Tables 5 through 7 is 40 kW/kgU for a PWR and 30 kW/kgU for a BWR. If Pave, the average cumulative specific power, is more than 35% higher than Pmax (i.e., 54 kW/kgU for PWR fuel and 40.5 kW/kgU for BWR fuel), the guide should not be used. When 1 < Pave/Pmax < 1.35, the guide can still be used, but an excess power adjustment factor, fp, must be applied. The excess power adjustment factor is

Equation 12 - f sub p equal to the square root of the quatity P sub ave divided by P sub max (Equation 12)

For Pave Pmax, fp = 1

2.4 The Enrichment Factor fe

The decay heat rates of Tables 1 through 3 and Tables 5 through 7 were calculated using initial enrichments of Tables 4 and 8. The enrichment factor fe is used to adjust the value ptab for the actual initial enrichment of the assembly Es. To calculate fe, the data in Tables 4 (BWR) or 8 (PWR) should be interpolated linearly to obtain the enrichment value Etab that corresponds to the assembly burnup, Btot. If Es/Etab < 0.6, the NRC staff recommends not using this guide. When Es/Etab 0.6, set the enrichment factor as follows:

fe = 1 + 0.01[a + b(Tc - d)][1 - Es/Etab ]

when Es/Etab 1.5,

fe = 1 - 0.005 [a + b(Tc - d)]

when Es/Etab > 1.5,

(Equation 13)

where the parameters a, b, and d vary with reactor type, Es, Etab, and Tc. These variables are defined in Tables 9 and 10.

Table 9 : Enrichment Factor Parameter Values for BWR Assemblies

Parameter in Equation 13 Parameter Value
Es /Etab < 1 Es /Etab > 1
1 Tc 40 Tc > 40 1 Tc 15 Tc > 15
a 5.7 5.7 0.6 0.6
b -0.525 0.184 -0.72 0.06
d 40 40 15 15

Table 10 : Enrichment Factor Parameter Values for PWR Assemblies

Parameter in Equation 13 Parameter value
Es /Etab < 1 Es /Etab > 1
1 Tc 40 Tc > 40 1 Tc 20 Tc > 20
a 4.8 4.8 1.8 1.8
b -0.6 0.133 -0.51 0.033
d 40 40 20 20

2.5 Safety Factor S

Before obtaining the final heat rate pfinal, an appropriate estimate of a percentage safety factor S should be determined. Evaluations of uncertainties performed as part of this project indicate that the safety factor should vary with burnup and cooling time.

For BWR assemblies:

S = 6.4 + 0.15 (Btot - 20) + 0.044 (Tc - 1) (Equation 14)

For PWR assemblies:

S = 6.2 + 0.06 (Btot - 25) + 0.050 (Tc - 1) (Equation 15)

The purpose of deriving spent fuel heat generation rates is usually to apply the heat rates in the computation of the temperatures for storage systems. A preferred engineering practice may be to calculate the temperatures prior to application of a final safety factor. This practice is acceptable if S is accounted for in the more comprehensive safety factors applied to the calculated temperatures.

2.6 Final Heat Generation Rate Evaluation

The equation for converting ptab, determined in Regulatory Position 2.1, to the final heat generation rate of the assembly, is

pfinal = (1 + 0.01S) f7 f7 fp fe ptab (Equation 16)

where f7, f7, fp, fe, and S are determined by the procedures given in Regulatory Positions 2.2 through 2.5.

2.7 Heat Rate Evaluation by LWRARC Code

The LWRARC (light-water-reactor afterheat rate calculation) code is an MS-DOS PC program that performs the calculations in this guide. The only input for cases in which the cooling time exceeds 15 years are Btot, Tres, Es, and Tc. Additionally, the short cooling time factors require Bi and Ti of the last and next-to-last cycles. The code features a pull-down menu system with data entry screens containing context-sensitive help messages and verification dialog boxes. The menus may be used with either a keyboard or a mouse. The code printout (one page per case) contains the input data, the computed safety and adjustment factors, and the interpolated and final computed decay heat rates. The output file may be printed, observed on a monitor, or saved. Input cases may be saved, retrieved, duplicated, or stacked in the input file.

The LWRARC code may be requested from the Radiation Safety Information Computational Center (RSICC).

Radiation Safety Information Computational Center
Oak Ridge National Laboratory
P.O. Box 2008
Oak Ridge, TN 37831-6362
Telephone: (423)574-6176
FAX: (423)574-6182
Electronic Mail: PDC@ornl.gov


Appendix A : Sample Case Using Heat Generation Rate Tables

A BWR fuel assembly with an average fuel enrichment of 2.6 wt-% 235U was in the reactor for four cycles. Determine its final heat generation rate with safety factors, using the method in this guide, at 4.2 years cooling time. Adequate details of the operating history associated with the fuel assembly are shown in Table A.1.

Table A.1 : Sample Case Operating History

Relative
Fuel
Cycle
Time from Startup of Fuel, Days Accumulated Burnup
(Best Maximum Estimate),
MWd/kgU
Cycle Startup Cycle Shutdown
1 0 300 8.1
2 340 590 14.7
3 630 910 20.9
4 940 1240 26.3

Note that the output of the LWRARC code for this case is shown in the first case of Appendix B of NUREG/CR-5625.(3)

Using Regulatory Position 1:

The following were given in the sample case (see Regulatory Position 1 for definitions):

Tres = 1240 d

Btot = 26.30 MWd/kgU

Tc = 4.2 y

Es = 2.6 wt-% 235U

Compute Te, Be, Pe, Te-1, Pe-1, Pave,e-1, and Pave from Regulatory Position 1 and Equations 2 through 4.

Te = 1240 - 940 = 300 d

Be = 26,300 - 20,900 = 5,400 kWd/kgU

Pe = (26,300 - 20,900)/300 = 18.00 kW/kgU

Te-1 = 940 - 630 = 310 d

Be-1 = 20,900 - 14,700 = 6,200 kWd/kgU

Pe-1 = 6,200/[0.8(310)] = 25.0 kW/kgU

Pave,e-1 = 20,900/[0.8(940)] = 27.793 kW/kgU

Pave = 26,300/[300 + 0.8 (940)] = 25.00 kW/kgU

Using Regulatory Position 2

Ptab should be determined from Pave, Btot, and Tc, as described in Regulatory Position 2.1. First, select the nearest heat rate values in Tables 2 and 3 for the following limits:

PL = 20 Pave PH = 30

BL = 25 Btot BH = 30

TL = 4 Tc TH = 5

Next, use the prescribed interpolation procedure for computing ptab from the tabular data. Although the order is optional, the example here interpolates between specific powers, burnups, and then cooling times. Denote the heat rate, p, as a function of specific power, burnup, and cooling time by p(P,B,T). The table values at PL and PH for BL and TL are

p(PL, BL, TL) = p(20,25,4) = 1.549

p(PH, BL, TL) = p(30,25,4) = 1.705

First, interpolate the above heat rates to Pave using

p(Pave,25,4) = p(20,25,4) + Fp [p(30,25,4) - p(20,25,4)]

where

Fp = (Pave - PL )/(PH -PL ) = 0.5

The result at p(Pave,25,4) is

p(Pave,25,4) = 1.549 + 0.5 (1.705 - 1.549)
= 1.627

The other three values at Pave are computed with a similar method:

p(Pave,30,4) = 1.827 + 0.5 (2.016 - 1.827)
= 1.9215
p(Pave,25,5) = 1.293
p(Pave,30,5) = 1.553

These are heat rates at the burnup and time limits.

Second, interpolate each of the above pairs of heat rates to Btot from the values at BL and BH :

FB = (Btot - BL )/(BH - BL ) = 0.26
p(Pave,Btot,4) = 1.627 + 0.26 (1.9215 - 1.627)
= 1.7036
p(Pave,Btot,5) = 1.3606

Third, compute the heat rate at Tc from the above values at TL and TH by an interpolation that is logarithmic in heat rate and linear in time:

FT = (Tc -TL )/(TH -TL ) = 0.2
log[p(Pave,Btot,Tc)] = log 1.7036 + 0.2 (log 1.3606 - log 1.7036)
= 0.2118
ptab = p(Pave,Btot,Tc) = 10 0.2118 = 1.629 W/kgU

With the value for ptab, the formulas of Regulatory Positions 2.2 through 2.6 can be used to determine pfinal. Since Tc < 7 y, use Equations 8 through 11 to calculate the short cooling time factors:

R = Pe /Pave - 1 = (18/25) - 1 = - 0.28

f7 = 1 + [0.25(-0.28)]/4.2 = 0.983

R = Pe-1 /Pave,e-1 - 1 = -0.1005

f7 = 1 + [0.08(-0.1005)]/4.2 = 0.998

Since Pave < PH = Pmax, the excess power factor, fp, is unity. Interpolating Table 4 enrichments to obtain the enrichment associated with the burnup yields

Etab = 2.3 + (2.7 - 2.3)(26.3 - 25)/(30 - 25)
= 2.404

The enrichment factor, fe, is then calculated using Equation 13:

fe = 1 + 0.01 (8.376)(1 - 2.6/2.404) = 0.993

because Es > Etab

The safety factor, S, for a BWR is given in Equation 14:

S = 6.4 + 0.15 (26.3 - 20) + 0.044 (4.2 - 1)
= 7.49%

Then, using Equation 16,

pfinal = (1 + 0.01 S) f7 f7 fp fe ptab

with the above adjustment factors and ptab yields

pfinal = 1.0749 � 0.983 � 0.998 � 1 � 0.993 � 1.629
= 1.0749 � 1.587 = 1.706 W/kgU

Thus, the final heat generation rate, including the safety factor, of the given fuel assembly is 1.706 W/kgU.

Appendix B : Acceptability and Limits of the Guide

Inherent difficulties arise in attempting to prepare a heat rate guide that has appropriate safety factors, is not excessively conservative, is easy to use, and applies to all commercial reactor spent fuel assemblies. In the endeavor to increase the value of the guide to licensees, the NRC staff made an effort to ensure that safe but not overly conservative heat rates were computed. The procedures and data recommended in the guide should be appropriate for most power reactor operations with only minor limitations in applicability.

In general, the guide should not be applied outside the parameters of Tables 1 through 8. These restrictions, in addition to certain limits on adjustment factors, are given in the text. The major table limits are summarized in Table B.1.

Table B.1 : Parameter Range for Applicability of the Regulatory Guide

Parameter BWR PWR
Tc(year) 1-110 1-110
Btot (MWd/kgU) 20-45 25-50
Pave (kW/kgU) 12-30 18-40

In using the guide, the lower limit on cooling time, Tc, and the upper limit on burnup, Btot, should never be extended. An adjustment factor, fp, can be applied if the specific power, Pave, does not exceed the maximum value of the tables by more than 35%. Thus, if Pave is greater than 54 kW/gU for PWR fuel or 40.5 kW/kgU for BWR fuel, the guide should not be applied. The minimum table value of specific power or burnup can be used for values below the table range; however, if the real value is considerably less than the table minimum, the heat rate derived can be excessively conservative. Also, the upper cooling time limit is conservative for longer cooling times.

In preparing generic depletion/decay analyses for specific applications, the most difficult condition to model is the power operating history of the assembly. Although a power history variation (other than the most extreme) does not significantly change the decay heat rate after a cooling time of approximately 7 years, it can have significant influence on the results in the first few years. Cooling time adjustment factors, f7 and f7, are applied to correct for variations in power history that differ from those used in the generation of the tables. For example, the heat rate at 1 year is increased substantially if the power in the last cycle is twice the average power of the assembly. The limits on the conditions in Regulatory Position 2.2 on ratios of cycle to average specific power are needed; first, to derive cooling time adjustment factors that are valid, and second, to exclude cases that are extremely atypical. Although these limits were determined so that the factors are safe, a reasonable degree of discretion should be used in the considerations of atypical assemblies-- particularly with regard to their power histories.

Another variable that requires attention is the 59Co content of the clad and structural materials. Cobalt-59 is partly transformed to 60Co in the reactor and subsequently contributes to the decay heat rate. The 59Co content used in deriving the tables here should apply only to assemblies containing Zircaloy-clad fuel pins. The 60Co contribution can become excessive for 59Co contents found in stainless-steel-clad fuel pins. Thus, the use of the guide for stainless-steel-clad assemblies should be limited to cooling times that exceed 20 years. Because 60Co has a 5.27-year half-life, the heat rate contribution from 60Co is reduced by the factor of 13.9 in 20 years.

In addition to the parameters used here, decay heat rates are a function of other variables to a lesser degree. Variations in moderator density (coolant pressure, temperature) can change decay heat rates, although calculations indicated that the expected differences (approximately 0.2% heat rate change per 1% change in water density, during any of the first 30-year decay times) are not sufficient to require additional corrections. The PWR decay heat rates in the tables were calculated for fuel assemblies containing water holes. Computed decay heat rates for assemblies containing burnable poison rods (BPRs) did not change significantly (<1% during the first 30-year decay) from fuel assemblies containing water holes.

Several conditions were considered in deriving the safety factors (Equations 14 and 15) that were developed for use in the guide. Partial uncertainties in the heat generation rates were computed for selected cases by applying the known standard deviations of half-lives, Q-values, and fission yields of all the fission product nuclides that make a significant contribution to decay heat rates. This calculation did not account for uncertainties in contributions produced by the neutron absorption in nuclides in the reactor flux, or from variations in other parameters. In addition to the standard deviations in neutron cross sections, much of the uncertainty from neutron absorption arose from approximations in the model used in the depletion analysis. In developing the safety factors, these more indirect uncertainties were determined from comparisons of the calculated total or individual nuclide decay heat rates with those determined by independent computational methods, as well as comparisons of heat rate measurements obtained for a variety of reactor spent fuel assemblies. Note from the equations that the safety factors increase with both burnup and cooling time. This increase in the safety factor is a result of the increased importance of the actinides to the decay heat with increased burnup and cooling time together with the larger uncertainty in actinide predictions caused by model approximations and limited experimental data.

Whenever the design or operating conditions for a spent fuel assembly exceed the parameter ranges accepted in this guide, another well-qualified method of analysis that accounts for the exceptions should be used. A well-qualified method would be one that has a technical basis that is validated against measured heat-rate data and has been demonstrated to provide conservative heat-rate estimates (i.e., per justified safety factors consistent with the measured data) for the extended design or operating conditions.

Appendix C : Contributions to Decay Heat Rates by Actinides, Fission Products, and Light Elements

The decay heat rates determined by the methods recommended in this guide are totals resulting from all sources of radioactive decay. In the tables of this Appendix C, the contributions to these totals from actinides, fission products, and light elements are listed separately. These values were used to construct the totals given in Tables 1-7. The values in this Appendix C represent some of the many results available from the codes described in NUREG/CR-5625.

Table C.1 : BWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 12 kW/kgU, Set 1

Table C.2 : BWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 12 kW/kgU, Set 2

Table C.3 : BWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 12 kW/kgU, Set 3

Table C.4 : BWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 20 kW/kgU, Set 1

Table C.5 : BWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 20 kW/kgU, Set 2

Table C.6: BWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 20 kW/kgU, Set 3

Table C.7 : BWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 30 kW/kgU, Set 1

Table C.8 : BWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 30 kW/kgU, Set 2

Table C.9 : BWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 30 kW/kgU, Set 3

Table C.10 : PWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 18 kW/kgU, Set 1

Table C.11 : PWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 18 kW/kgU, Set 2

Table C.12 : PWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 18 kW/kgU, Set 3

Table C.13 : PWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 28 kW/kgU, Set 1

Table C.14 : PWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 28 kW/kgU, Set 2

Table C.15 : PWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 28 kW/kgU, Set 3

Table C.16 : PWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 40 kW/kgU, Set 1

Table C.17 : PWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 40 kW/kgU, Set 2

Table C.18 : PWR Decay Heat Rates (W/kgU) of Light Elements, Actinides, and Fission Products, for Specific Power = 40 kW/kgU, Set 3

Value/Impact Statement

A Value/Impact Statement was published with Regulatory Guide 3.54 when it was issued in September 1994. No changes are necessary, so a separate value/impact statement for this proposed Revision 1 has not been prepared. A copy of the value/impact statement is available for inspection or copying for a fee in the Commission's Public Document Room at 2120 L Street NW., Washington, DC, under Regulatory Guide 3.54. The PDR's mailing address is Mail Stop LL-6, Washington, DC 20555; telephone (202)634-3273; fax (202)634-3343.

1. Technical Support for a Proposed Decay Heat Guide Using SAS2H/ORIGEN-S Data,  NUREG/CR-5625 (ORNL-6698), September 1994. Copies are available for  inspection or copying for a fee from the NRC Public Document Room at 2120 L  Street NW., Washington, DC; the PDR's mailing address is Mail Stop LL-6,  Washington, DC 20555; telephone (202)634-3273; fax (202)634-3343. Copies may  be purchased at current rates from the U.S. Government Printing Office, P.O.  Box 37082, Washington, DC 20402-9328 [telephone (202)512-1800]; or from the  National Technical Information Service by writing NTIS at 5285 Port Royal  Road, Springfield, VA 22161.

2. The International System of units.

3. Technical Support for a Proposed Decay Heat Guide Using SAS2H/ORIGEN-S Data, NUREG/CR-5625 (ORNL-6698), September 1994. Copies are available for inspection or copying for a fee from the NRC Public Document Room at 2120 L Street NW., Washington, DC; the PDR's mailing address is Mail Stop LL-6, Washington, DC 20555; telephone (202)634-3273; fax (202)634-3343. Copies may be purchased at current rates from the U.S. Government Printing Office, P.O. Box 37082, Washington, DC 20402-9328 [telephone (202)512-1800]; or from the National Technical Information Service by writing NTIS at 5285 Port Royal Road, Springfield, VA 22161.



Privacy Policy | Site Disclaimer
Thursday, February 15, 2007