ML20154Q283
ML20154Q283 | |
Person / Time | |
---|---|
Site: | Byron |
Issue date: | 05/31/1988 |
From: | COMMONWEALTH EDISON CO. |
To: | |
Shared Package | |
ML20154Q280 | List: |
References | |
NUDOCS 8806060236 | |
Download: ML20154Q283 (183) | |
Text
{{#Wiki_filter:- LI ENSING REPORT 07; ON HIGH DENSITY SPENT FUEL RACKS FOR BYRON UNITS I AND 2 NRC DOCKET NO. 50 - 454 50 - 455 i e COMMONWEALTH EDISON COMPANY l CHICAGO, ILLINOIS 60603 l MAY 1988 (Revision 2) O_- np r 1888e 888881s4 DCD
O
\~ / TABLE OF CONTENTS Section .P, ale,
1.0 INTRODUCTION
1- 1 2.0 GENERAL ARRANGEMENT 2-1 3.0 RACK CONSTRUCTION 3-1 3.1 Fabrication Details 3-1 3.1.1 Region 1 3-1 3.1.2 Region 2 3-3 3.2 Codes, Standards, and Practices for the Spent 6 Fuel Pool Modification 4.0 CRITICALITY SAFETY ANALYSES 4-1 4.1 Design Bases 4-1 4.2 Summary of Criticality Analyses 4-2 4.2.1 Normal Operating Conditions 4-2 4.2.2 Abnormal and Accident Conditions 4-3 4.2.3 New Fuel Storage 4-4 4.3 Reference Fuel Storage Cell 4.3.1 Reference Fuel Assembly 4-5 4.3.2 Region 1 Storage Cells 4-5 4.3.3 Region 2 Storage Cells 4-6 4.4 Analytical Methodology 4-6 4.4.1.1 Reference Analytical Methods and Bias 4-6 4.4.1.2 Analytical Methods and Bias for Modified Cases 4-8 4.4.2 Fuel Burnup Calculations 4-9 4.4.3 Ef fect of Axial Burnup Distribution 4-11 4.4.4 Long-Term Decay 4-12 4.5 Region 1 Criticality Analysis and Tolerance Variations 4-12 4.5.1 Nominal Design Case 4-12 4.5.2 Boron Loading Variation 4-13 O' 4.5.3 Storage Cell Lattice Pitch Variation 4-13 4-14 4.5.4 Stainless Steel Thickness Tolerances 4.5.5 Fuel Enrichment and Density Variation 4-14 i
4 k) TABLE OF CONTENTS (continued) Section Page 4.5.6 Boraflex Width Tolerance Variation 4-14 4.5.7 Axial Cutback of Boraflex 4-15 4.6 RLgion 2 Criticality Analysis and Tolerance Variations 4-15 4.6.1 Nominal Design Case 4-15 4.6.2 Boron Loading Variation 4-16 4.6.3 Storage Cell Lattice Pitch Variations 4-17 4.6.4 Stainless Steel Thickness Tolerance 4-17 4.6.5 Fuel Enrichment and Density Variation 4- 17 4.6.6 Boraflex Width Tolerance 4-17 4.7 Abnormal and Accident Conditions 4- 17 4.7.1 Eccentric Positioning of Fuel Assembly in Storage Rack 4-17 4.7.2 Temperature and Water Density Ef fects 4-18 4.7.3 Dropped Fuel Assembly Accident 4-18 4.7.4 Abnormal Location of a Fuel Assembly 4-19 ( }) 4.7.5 Lateral Rack Movement 4-19 , 4.8 New Fuel Storage 4-20 4.8.1 Storage in Region 1, Dry 4-20 4.8.2 Storage in Region 2 Flooded 4-20 4.8.3 Storage in Region 2, Dry 4-21 REFERENCES 4-22 5.0 THERMAL-HYDRAULIC ANALYSES 5-1 5.1 Decay Heat Calculations for the Spent Fuel 5-1 5.1.1 Basis 5-2 5.1.2 Model Description 5-5 5.1.3 Decay Hea Calculation Results 5-7 5.2 Thermal-Hydraulic Analyses for Spent Fuel Cooling 5-8 5.2.1 Basis 5-9 5.2.2 Model Description 5-9 + 5.2.3 Results 5-11 Q'- REFERENCES 5-13 ii _ _ _ _ . _ _ - _ _ . _ ~,__ _ _ ____, _ __
TABLE OF CONTENTS (V3 (continued) Section Page 6.0 STRUCTURAL ANALYSIS 6-1 6.1 Analysis Outline 6-1 6.2 Fuel Rack - Fuel Assembly Model 6-3 6.2.1 Outline of Model 6-3 6.2.2 Model Description 6-5 6.2.3 Fluid Coupling 6-5 6.2.4 Damping 6-6 6.2.5 Impact 6-7 6.3 Assembly of the Dynamic Model 5-7 6.4 Time Integration of the Equations of Motion 6-9 6.4.1 Time-History Analysis Using Multi- 6-9 Degree of Freedom Rack Model 6.4.2 Evaluation of Potential for Inter- 6-11 Rack Impact 6.5 Structural Acceptance Criteria 6-11 i 6.6 Material Properties 6-13 6.7 Stress Limits for Various Conditions 6-13 6.7.1 Normal and Upset Conditions (Level A or Level B) 6-13 6.7.2 Level D Service Limits 6-15 6.8 Results 6-15 6.9 Impact Analyses 6-17 6.9.1 Impact Loading Between Fuel Assembly and Cell Wall 6-17 6.9.2 Impacts Between Adjacent Racks 6-17 6.10 Weld Stresses 6-17 6.10.1 Baseplate to Rack Welds and Cell-to-Cell Welds 6-18 6.10.2 Heating of an Isolated Cell 6-18 O 6.11 Su-ery of Mecnanicei Anaisses 6-19
; 6.12 Definition of Terms Used in Section 6 6-20 iii I
'l C) i TABLE OF CONTENTS (continued)
Page Sectig REFERENCES 6-21 7.0 ENVIRONMENTAL EVALUATION 7-1 7.1 Summary 7-1 7.2 -Characteristics of Stored Fuel 7-1 7.3 Related Industry Experience 7-2 7.4 Byron Nuclear Power Station Experience 7-3 7.5 Spent Fuel Cooling and Cleanup System 7-3 7.6 Fuel Pool Radiation Shielding 7-4 7.6.1 Source Terms 7-4 7.6.2 Radiation Shielding 7-5 Q 7-6 7.7 Radiological Consequences 7.8 Reracking Operation 7-6 7.9 Conclusions 7-7 REFERENCES 7-9 8.0 IN-SERVICE SURYEILLANCE PROGRAM FOR BORAFLEX NEUTRON ABSORBING MATERIAL 8-1 8.1 Program Intent 8-1 8.2 Description of Specimens 8-1 8.3 Specimen Evaluation 8-1 9-1 9.0 COST / BENEFIT ASSESSMENT 9.1 Specific Needs for Spent Fuel Storage 9-1 9.2 Cost of Spent Fuel Storage 9-1 9.3 Alternatives to Spent Fuel Storage 9-2 9-3 9.4 Resource Comitments iv ,
O TA8u 0F CONTEuS
.(continued) l l
Section g REFERENCES 9-4 l 10.0 QUALITY ASSURANCE PROGRAM 10-1 10.1 Introduction 10-1 10.2 General 10-1 10.3 System Highlights 10-1 10.4 Sumary 10-2 i APPENDIX A BENCHMARK CALCULATIONS A-1 APPENDIX 8 NRC QUESTIONS FROM DECEMBER 22, 1986 : MEETING 8-1 l O l l r i r a ? a 4 1 I i L t i O i y ! l
(]
1.0 INTRODUCTION
This report describes the design, fabrication, and safety analysis of high density spent fuel storage racks manuf actured by Joseph Oat Corporation (Oat) for the Byron Station Unit 1 and Unit 2. The plant, which is located two miles east of the Rock River and approximately three miles southwest of Byron in Ogle County, is owned and operated by Comonwealth Edison Company (CECO). Byron is a two-unit pressurized water reactor (PWR) with a net design capacity of 1120 megawatts electric for each unit. Each of the two reactor cores contains 193 fuel assemblies and is rated to produce 3411 thermal megawatts (MWt). At present, there is one (normal core offload) batch spent fuel assemblies stored in the spent fuel pool. Unit 1 went into comercial operation in September of 1985. Unit 2 went into comercial operation in August, 1987. The two units share one common spent fuel storage pool which is currently
- "]
r licensed for the storage of 1060 spent fuel assemblies. As shown in Table 1.1, the storage pool would lose full core discharge capability in 1994. The proposed reracking will increase the number of pool storage locations to 2370 (includes six failed fuel locations). Table 1.1 indicates that the new racks will provide adequate storage with full core discharge capability well into the next century (circa 2009). Table 1.1 is based on an estimated 18-month fuel cycle. Current trends toward longer cycles, extended burnup, and higher enrichment would further extend the time span of onsite storage. The proLosed racks are free-standing and self-supporting. The principal construction materials are ASTM A-240, Type 304L stainless steel for the structural members and shapes, with "Boraflex", a paterted product of 815C0 (a division of Brand, Inc.), and Boral, a patented product of AAR Brooks & Perkins for neutron attenuation. Whereas the fixed height support legs employ 304L series austenitic stainless plate and pipe material, the adjustable support legs are constructed from SA351-CF3 and SA217-CA15 casting stock, and for certain support legs 400 series stainless (SA479-410) steel material . ,/~} 1-1
O Tne specificatioes for design, construction, and auai4ty assurance for t e high density spent fuel storage racks were prepared by Sargent & Lundy l Engineers (S&L) of Chicago, Illinois. The mechanical design and fabrication of the hardware was done by Oat. Seismic / structural analysis, thermal-hydraulic analysis, and other related calculations were performed by Holtec International of Mount Laurel, New Jersey and S&L. S&L provided the seismic
.'e evaluation. ~
response spectra and performed the spent fuel pool st S&L performed the radiation shielding analysis. Southern Scieace, a division of Black and Yeatch, served as a consultant to Oat in the area of criticality analysis. Further criticality analysis was performed by Westinghouse ' Electric. The nathematical analysis responsibility for this contract was taken over by Holtec International of Mount Laurel, New Jersey, who prepared the December,1986 amendment to this licensing report on behalf of Oat. The analyses performed by Holtec in conjunction with Black and Yeatch and S&L demonstrate that acceptable margins of safety exist with respect to appropriate NRC and ASME acceptance criteria. A cost-benefit comparison of several potential spent fuel disposition alternatives indicates that reracking of the Byron pool is the lowest risk and most cost-ef fective alternative, and that neither the reracking operation nor the increased onsite storage of irradiated material pose an undue hazard to the plant staff or the public. The following sections provide a synopsis of the design, fabr':ation, nuclear criticality analysis, thermal / hydraulic analysis, structural analysis, accident analysis, environmental analysis, and cost-benefit appraisal of the high density spent fuel racks. In particular, the integrity of the rack structure under the specified combinations of inertial, seismic, and mechanical loads and thermal gradient per NUREG-0800 is demonstrated. Also included are descriptions of the rack In-Service Surveillance Program and the Oat Quality Assurance Program. This Quality Assurance Program has been reviewed and found ace.eptable for engineered fabrication of ASME Section III, Class 1, 2 and 3 and MC Components by both ASME and the NRC. O 1- 2
t die t t _O BYR0R UNIT 1 & UNIT 2 FUEL ASSEMBLY DISCHARGE (TENTATIVE SCHEDULE) Remaining Capacity Total Number of Discharged Without With*** Assemblies Batch Proposed Proposed Uni t Date Discharged Site Expansion Ex.pansion
#1 Feb. 1987 88 88 972 2782 .
l
#1 Aug. 1988 88 176 884 2694 #2 Dec. 1988 88 264 796 2006 #1 Feb. 1990 88 352 708 2518 #2 June 1990 88 440 620 2430 #1 Aug. 1991 88 528 532 2342 q
L
#2 Dec. 1991 88 616 444 2254 #1 rib. 1993 84 700 360 2170 f2 June 1993 88 788 272 2082 l #1 Aug. 1994 84 872 188** i .* iB #2 Dec. 1994 84 956 104 1914 #1 Feb. 1996 84 1040 20* 1830 #2 June 1996 84 1124 1746 #1 Aug. 1997 84 1208 1662 I #2 Dec. 1997 84 1292 1578
- Partial core discharge capability lost - 84 assemblies
** Full core discharge capability lost - 193 assemblies *** Includes 6 failed fuel storage positiors O
1-3
. [#] Table 1.1 (Continued) BYRON UNIT 1 & UNIT 2 FUEL ASSEMBLY DISCHARGE (TENTATIVE SCHEDULE) Remaining Capacity Total ' Number of Discharged Without With*** Assemblies Batch Proposed Proposed Unit Date Discharged Size Expansion Expansion
#1 Feb. 1999 84 1376 1494 #2 June 1999 84 1460 1410 #1 Aug. 2000 84 1544 1326 #2 Dec. 2000 84 1628 1242 #1 Feb. 2002 84 1712 1158 #2 June 2002 84 1796 1074 #1 Aug. 2003 84 1880 990 #2 Dec. 2003 84 1964 906 #1 Feb. 2005 84 2048 822 #2 June 2003 84 2132 738 #1 Aug. 2006 84 2216 654 #2 Dec. 2006 84 2300 570 #1 Feb. 2008 84 2384 486 #2 June 2008 84 2468 402 Aug. 2009 84 2552 318 #1 Dec. 2009 84 2636 234 #2 #1 Feb. 2011 84 2720 150** #2 June 2011 84 2804 66* *** Includes 6 failed fuel storage positions.
1-4
1
,~
h, 2.0 GENERAL ARRANGEMENT The high density spent fuel racks consist of individual cells with 8.85-inch (nominal) square cross-section, each of which accommodates a single Westinghouse PWR fuel issembly or equivalent. A total of 2864 cells and six f ailed fuel storage cells are arranged in 23 distinct modules of varying sizes in two regions. Region 1 is designed for storage of new fuel assemblies with enrichments up to 4.2 weight percent U-235. Region 1 is also designed to store fuel assemblies with enrichments up to 4.2 weight percent U-235 that have not achieved adequate burnup for Region 2. The Region 2 cells are capable of accommodating fuel assemblies with various initial enrichments which have accumulated minimum burnups within an acceptable bound as depicted in Figure 4.1. Figure 2.1 shows the arrangement of the rack modules in the spent fuel pool. The high density racks are engineered to achieve the dual objective of maxirm protection against structural loadings (arising from ground motion, thermal O) L stresses, etc.) and the maximization of available storage locations. In general, a greater width-to-height aspect ratio provides greater margin against rigid body tipping. Hence, the modules are made as large as possible within the constraints of transportation and site handling capabilities. As shown in Figure 2.1, there are 23 discrete modules arranged in the fuel pool. Each rack module is equipped (see Figures 2.2a and 2.2b),with girdle bars, one-inch-thick by 3-1/2 inches high. The nominal gap between cell walls of adjacent modules is two inches. The modules make surface contact between their contiguous sides at the girdle bar locations and thus maintain a specified gap between them. Table 2.1 summarizes the typical physical data for each Region 1 and Region 2 rack. Table 2.2 summarizes other pertinent information on each rack module. 1 0 2-1
l O Table 2.1 DESIGN DATA (Cell Pitch) Flux Trap Region Nominal Min. Boral Min. B-10 Gap (nominal) i in. Loading Loading in.
~
1 Nominal 10.32 N&S N/A .020 gm/cm 2 1.16* , 410.42 E&W N/A 1.26 Boral 10.32 N&S .020gm/cm 2 .020gm/cm 2 1.01 Insert a10.42 E&W 1.11 1 Boral 10.32 N&S .025gm/cm 2 N/A 1.18 Only &10.42 E&W 1.28 2 9.03** N/A .010 gm/cm 2 0.0 4 The minimum flux trap gap can be .040" less than the nominal gap. The minimum pitch can be .040" less than the nominal pitch. O 2-2
... - - - . - - . , . . . - . . , . - - - - . ~ . - . . - . . . . _ . . . - - - . - - . . . - . . . . . . - . . - - . . . - - - - .
Table 2,2 MODULE DATA Number Cells Approximate Module of per Module Weight Region Type Modules Module Size (ib/ module) I Al 1 104 13x8 20,800 I B1-3 3 96 12x8 19,200 11 Cl-6 6 168 14x12 26,900
!! 01&3 2 126 14x9 20,150 .II 04 1 113 14x9 18.360 -(2x2+3x3) 14x9-(4x3) 18,250
{} II D2 El 1 1 114 112 14x8 17,909 II II F1 1 165 11x15 26,600 II G1 1 90 10x9 14,700 t' 7x8 8,950 II H1 1 56 II J1 1 35+6 7x5 10,150 failed fuel containers II K1 1 117 13x9 19,000 11 L1-L2 2 156 13x12 25,200 l 11 M1 1 98 14x7 16,000 l O 2-3 i
. . _ . _ . . . . _ , . _ _ . - . . . . . . ~ . . . . _ . _ , . . _ . . . . . . . . . . _ . , . _ . - . - . . . . - . . . _ . - - - _ . . .-
1 3 4 O 'y RE '"\""['"~j""Y"' 7
.N e' -
Di Cs Cl Ai g i 14X12 14 X ia 13X 8 Y R 14x 9 d~ l REGION i ( I N D4 Cs C2 Bi % bD 14X12 14 X 12 LEX 8 9 I4X 9 -13 REGION h z 1 H _- CL . 8 'i4 D, 14X 9-12 C4 C3 14XI2 B2 12X 8 s 1 5 14X 1?. REGION % i
. R t " - I * < 1 - ~
l 0 1 I : K, L i. La B3 % 13 Xia 12X 8 l $ $ /3X9 13 X12 RE610N s t t
' I
! F ( Ji Hi Gi ez~ *~' l Q~~ 7X5 7x8 Fi ~ ~~ 9 j il X l.: y s ; 4 io x .9
',/
I i= / l /
- 0 M' I D3 1 k SHIPPING //
/CASK AREA 14X7 14X 8 14X9 g ; s /, a/ / REF7 - ,u_ w ,w
- PLUS (6) /42" 43.5I ! 72 f'~ e/ 55' - .f.# 2*
O FAILEO ASSEMBLY; J 0 CELLS. FIG U R E 2.1 POOL L AYOUT 2-4
GIRDLE 3AR C)
.f ~ ~ , ' i e-p- fi 7 ~2 ? -{j u '- -; l - r . . i y
- --j ;-u i l i l i i . i f
. i ,
i i -
, i e
- i. ' I !
I I I t l , , e , [ ' I i l , q
- t l i
l l
. I i ' I I i e i I l . 3 i , - ,
i i , l , i i , ,
. , , i i , , , i i i i !
I I i , , i i i .
, I i , I I I e i ,
i i i
, ' r I I , l i i i i i i i , I I '
I l 1 e . , ;
, I, i , , , , i , ,
l l l l l l l
' l l La , , , __i i ..t i ..; i . _-
i i_ _2. m y..- lD i ! l ; i ,i O l i l i l i i i i i
' I l
i
, i i
l i i
! , i i I ( ' ' I i 1 ' I e i l l ' i l i i i , i ; ', i i ,
i i i i i i i i I i , i i I i i i ; i i i l I i i l g l ' I I
' ( i i I 3 , l i l l i i l i e l l , ' i 1 ' ! I l i i I i l l i i I i i i .i i .
l 1 I i i . , i
- I l i i ! I i g i l ! . g I I I i i i i I I i i i ' I i l ,'
l I 1 I
; i l i ,
t i i i i ! g
'i I ;
IO !O ! OI l ' j O! ;, Ol ; i l 0
- i;
'i i
- i 1
i, i i 1 f
'\ bf g
_+-~ j I _ i J i O FIGURE 2.2 a TYPICAL RACK ELEVATION-REGION 1 2-5
GIRDLE BAR
= . . -
i 1 i i i i i.
/ n $ Rg0 ( . ' ~ ' l ' ~ ' ' - ! l !l p R % l t .
bII * ' I t ,' i ' i e s ,' l(
%) - h ,'i j i l ,, , - i i'i . i i n i, i ..., '
I iWI 18 l l 8 s, e
.l ,
A - e
' I t i i r t : l g d
l I ' l I I l I I I l l l
\
i I 1 I l l l i l l i i e I
' I I ! l l 1
3
; l l
l e I l l l I l ' I ' I I l s l ; l l l ! ' l I l l l 1 l l l i I i i i I i l' I . i i I l - n I 1 I l l l i l
.. l i .
l , i , i . . i 2m
! ! '__l _
j l l ,!. _ l _ .l l
~T ~1 i T i ir- ; w I - 'j i l i
i ; l i l i l l l
- l .
l l j e 1 l 8 l 3 8 g l ' 8 8 l e l i 8 I l , I p l I e l l 8
, l h l 1
I l i l i i
- l l 1 8 l 1
[ l l l l i i i l l 1
' ' 8 I I I ;
l ,
, i l ! l I
i l I I l l , I I I l I l ', ' I , I I g l I i I I
! ' l l , I I ! l t '
l - ' I I I i l l I I ' I i , I . i 1 . 8 I , l l l l I
' I l l l
I I l I l I l r l e i ' l r I ' i l i l I i i I i mi i i *g 4 ,1 j ,j 8
-( l- >i "" 7" l- .)[ {* 4.j j
4
~~r. 4 ]7,1 1 f I t i Ik _ ,i f l I I
f
, p l_
i e J
}~ ~ i-I I 1
l TYPICAL RACK ELEV ATIOrl-REGIOff 2 O FIGURE 2.2 b 2-6
g 3.0 RACK CONSTRUCTION 3.1 FABRICATION DETAILS 3.1.1 Region 1 The rack module is fabricated from ASTM A-240-304L austenitic stainless steel sheet and plate material, ASME SA351-CF3 and SA217-CA15 casting material and SA 479-410 material. = The weld filler material utilized in body welds is ASME SFA-5.9, Type 308L and 308LSI. Boraflex and Boral serve as the neutron l absorber material. The detailed neutronic properties of Boraflex may be found in Section 4. The Boraflex experience list is given in Table 3.1. The Boral experience is given to Table 3.1.a. A typical module contains storage cells which have an 8.85-inch nominal square cross-sectional opening. This dimension ensures that fuel assemblies with maximum expected axial bow can be inserted and removed from the storage cells without any damage to the fuel assemblies or the rack modules. Figure 3.1 shows a horizontal cross-section of a 3x3 array. The cells provide a smooth and continuous surf ace for lateral contact with the fuel assembly. The anatomy of the rack modules is best explained by describing the components of the design, namely: 0 Internal squure box 0 Neutron absorber material (Boraflex Boral) j l 0 Poison sheathing 0 Gap element 0 Baseplate 0 Support assembly 0 Top Lead-in
- a. Internal Square Box This element provides the lateral bearing surface to the fuel assembly, it is f Tbricated by joining two formed cnannels (Figure 3.2) using a controlled seam welding operation. This element is an 8.85-inch square (nominal) cross-section by 168-7/8 inches long.
O 3-1 l 1
1 1 l ( b. Neutron Absorber Material (Boraflex, Boral) l Boraflex is placed on all four sides of a square tube over a length of 139-1/2 inches, which covers the active fuel length except the top and bottom 2.25 inches, for Region 1 Nominal Case. In addition to the Boraflex, Boral sheets are added to the flux traps of two Region 1 racks (B2 and 83) which have already been built. Boral sheets shall be used in place of Boraflex to cover the entire active fuel on Region 1 racks ( Al and B1) which are not yet built. For Region 2 the Boraflex length is 144 inches and it covers the entire length of 1.he active fuel,
- c. Poison Sheathing The poison sheathing (cover plate) shown in Figure 3.4, serves to position and retain the poison material in its designated space.
This is accomplished by spot welding the cover sheet to the square tube along the former's edges at numerous (at least 20) locations. This manner of attichment ensures that the poison material will not sag or laterally displace during fabrication processes and under any subsequent loading condition.
- d. Gap Element Gap elements, illustrated in Figure 3.3, position two inner boxes at a predetermined distance to maintain the minimum flux trap gap
(>) s required between two boxes. The gap element is welded to the inner box by fillet welds. An array of composite box assemblies welded as indicated in Figure 3.1 form the honeycomb gridwork of cells which harnesses the structural strength of all sheet and plate type members in an efficient manner. The array of composite boxes has overall bending, torsional, and axial rigidities which are an order of magnitude greater than configurations utilizing grid bar type of construction,
- e. Baseplate The baseplate is a 5/8-inch thick plate type member which has 6-inch diameter holes concentrically located with respect to the internal square tube, except. at support leg locations, where the hole size is 5 inches in diameter. These holes provide the primary path for coolant flow. Secondary flow paths are available between adjacent cells via tne lateral flow holes (1 inch in diameter) near the root of the honeycomb (Figures 3.5a and 3.5b) which preclude flow
- blockages. The honeycomb is welded to the baseplate with 3/32-inch fillet welds,
- f. Support Assembly Each module has at least four support legs. Most supports are g adjustable in length to enable levelling of the rack. The variable
{J height support assembly consists of a flat-footed spindle which rides into an internally-threaded cylindrical member. The cylindrical member is attached to the underside of the baseplate 3-2
l O tnrousa ri11et e"o pertie' 9eaetretioa weias. Tae bese or tne riet-footed spindle sits on the pool floor. Levelling of the rack 1 I I modules is accomplished by turning the square sprocket i ' the spindle using a long arm (approximately 46 feet long) scJare head wrench. Figure 3.6A shows a vertical cross-section of ohe adjustable support assembly. The supports elevate the module baseplate approximately 7-1/2 inches above the pool floor, thus creating the water plenum for coolant flow. The lateral holes in the cylindrical member provide the coolant entry path leading into the bottom of the storage locations. Short adjustable legs are used on racks El, L2, C2 (for detail see Fig. 3.6B). In addition, one small fixed support leg is used on rack D' (for detail see Fig. 3.6C).
- g. Top Lead-in Lead-ins are provided on each cell to facilitate fuel assembly insertion. The lead-ins of contiguous walls of adjacent cells are structurally connected at the lead-in. These lead-in joints aid in reducing the lateral deflection of the inner square tube due to the impact of fuel assemblies during the ground motion (postulated seismic motion specified in the FSAR). This type of construction leads to natural venting locations for the inter-cell space where the neutron absorber material is located.
3.1.2 Region 2 The rack modules in Region 2 are fabricated from the same material as that used for Region 1 modules, i.e., ASTM A-240-304L austenitic stainless steel. A typical Region 2 module storage cell also has an 8.85-inch nominal square cross-sectional opening. Figure 3.5b shows a partial vertical section of Region 2 module. Figure 3.7 shows a horizontal cross-section of a 3 x 3 array. The rack construction varies from that for Region 1 f nasmuch as the sta .11ess steel cover plates, gap elements, and top lead-ins are eliminated. Hence, the basic components of this design are as follows: 0 Internal square box 0 Neutron absorber material l 0 Side strips l 0 Baseplate 0 Support assembly 3-3
o V In this construction, two channel elements form a box with an 8.85-inch nominal square cross-sectional opening. The poison material is placed between two boxes as shown in Figure 3.7. Stainless steel side strips are inserted in both sides of the poison material to firmly locate it in the lateral direction. The bottom strip positions the poison material in the vertical direction to envelope the entire active fuel length of a fuel assembly (Figure 3.5b). Two adjacent boxe; and the side strip between boxes are welded togetner as shown in Figure 3.7, to form the honeycomb rack module. The bast.olate and adjustable support assemblies are incorporated in exactly the same nanner as described for Region 1 in the preceding section. 3.2 CODEC. STANDARDS, AND PRACTICES FOR THE SPENT FUEL P0OL MODIFlTATION The fabrication of the rack modules is performed under a strict quality assurance system suitale for ASME Section III, Class 1, 2 and 3 manufacturing O wnicn nes been in piece at Oet for over 10 veers. The following codes, standards and practices were used as applicable for the design, construction, and assembly of the spent fuel storage rackc and analysis of the pool structure. Additional specific references related to detailed analyses are given in each section,
- a. Design Codes (1) AISC Manual of Steel Construction, 8th Edition, 1980.
(2) ANSI N210-1976, "Design Objectives for Light Water Reactor Spent Fuel Storage Facilities at Nuclear Power Stations." (3) American Society of Mechanical Engineers ( ASME), Boiler and Pressure Vessel Code, Section III, 1983 Edition up to and including Summer 1983 Addenda (Subsection NF). (4) ASNT-TC-1A June,1980 American Society for Nondestructive Testing (Recomended Practice for Personnel Qualifications). ,.s
- b. Material Codes (1) American Society for Testing and Materials (ASTM) Standards -
A-240. 3-4
l'b \_/ (2) American Society of Mechanical Engineers ( ASME), Boiler and Pressure Vessel Code, Section II - Parts A and C,1983 Edition, up to and including Summer 1983 Addenda,
- c. Welding Codes ASME Boiler and Pressure Vessel Code, Section IX - Welding and Brazing Qualifications,1983 Edition up to and including Summer, 1983 Addenda,
- d. Quality Assurance, Cleanliness, Packaging, Shipping, Receiving, Storage, and Handling Requirements (1) ANSI N45.2.2 - Packaging, Shipping, Receiving, Storage and Handling of Items for Nuclear Power Plants.
(2) ANSI 45.2.1 - Cleaning of Fluid Systems and Associated Components during Construction Phase of Nuclear Power Plants. (3) AS:iE Boiler and Pressure Vessel, Section V, Nondestructive Examination,1983 Edition, including Summer and Winter 1983 Addenda. (4) ANSI - N16.1-75 Nuclear Criticality Safety Operations with Fissiunable Materials Outside Reactors. (]) (5) ANSI - N16.9-75 Validation of Calculation Methods for Nuclear Criticality Safety. (6) ANSI - N45.2.11,1974 Quality Assurance Requirements for the Design of Nuclear Power Plants.
- e. Other References (1) NRC Regulatory Guides, Division 1, Regulatory Guides 1.13, Rev.
2 (proposed); 1.29, Rev. 3; 1.31, Rev. 3; 1.61, Rev. 0; 1.71, Rev. 0; 1.85, Rev. 22; 1.92, Rev.1; 1.124, Rev.1; and 3.41, Rev. 1. . (2) General Design Criteria for Nuclear Power Plants, Code of Federal Regulations, Title 10, Part 50, Appendix A (GDC Nos.1, 2, 61, 62, and 63). (3) NUREG-0800, Standard Review Plan (1981). (4) "0T Position for Review and Acceptance of Spent Fuel Storage and Handling Applications," dated April 14, 1978, and the modifications to this document of January 18, 1979. 3-5
.) '# Table 3.1 BORAFLEX EXPERIENCE FOR HIGH DENSITY RACKS Unit Plant NRC Licensing Site No. Type Docket No. Status Point Beach 1&2 PWR 50-226 & 301 Licensed Nine Mile Point 1 BWR 50-220 Licensed Oconee 1&2 PWR 50-269 & 270 Licensed Prairie Island 1&2 PWR 50-282 & 306 Licensed Calvert Cliffs 2 PWR 50-318 Licensed Quad Cities
- 1&2 BWR 50-254 & 265 Licensed Watts Bar 1&2 PWR 50-390 & 391 Pending Waterford 3 PWR 50-382 Pending fh Fermi, 50-341 Licensed V 2 BWR H.B. Robinson 2 PWR 50-261 Licensed River Bend 1 BWR 50-458 Licensed Rancho Seco 1 PWR 50-312 Licensed Nine Mile Point 2 BWR 50-410 To Be applied for Shearon Harris 1 PWR 50-400 To be applied for Millstone 3 PWR 50-423 To be applied for Grand Gulf
- 1 BWR 50-416 Pending Oyster Creek BWR 50-219 Licensed V.C. Sumer PWR 50-395 Licensed Diablo Canyon
- 1&2 PWR 50-275 & 323 Licensed O -
aosePn Oet CorPoretion febriceted recks 3-6
i- , (,/, - Table 3.1.a BORAL EXPERIENCE FOR HIGH DENSITY RACKS Unit Plant NRC Licensing Site No. Type Docket No. Status Bellefonte 1&2 PWR 50-438 D. C. Cook 1&2 PWR 50-316 Indian Point 3 PWR 50-286 Maire Yankee 50-309 Salem 1&2 PWR 50-272 Seabrook PWR 50-443 Sequoyah 1&2 PWR 50-327 Yankee Rowe PWR 50-029 Zion 1&2 PWR 50-295 Browns Ferry 1,283 BWR 50-259,260, 296 Brunswick 1&2 BWR 50-325 l (%)') Clinton BWR 50-461 Cooper BWR 59-298 Dresden 283 BWR 50-237, 249 Duane Arnold BWR 50-331 f J. A. Fitzpatrick BWR 50-333 E.I. Hatch 1&2 BWR 50-231, 366 Hope Creek BWR 50-354 Humbolt Bay 50-133 La Crosse BWR 50-409 Limerick 1&2 BWR 50-352 l Monticello BWR 50-263 Peachbottom 2&3 BWR 50-277 Perry 182 BWR 50-440, 441 Pilgrim BWR 50-293 Shoreham BWR 50-322 Susquehanna 132 BWR 50-387 () Vermont Yankee BWR 50-271 3-7
O #'# # NEUTRON . M. 3 2 ABSORBER INTERNAL SQUARE BOX (BORAL SHEET) PITCH x.m um ._x mxo m x x _- - _ y p a g-i a '
't M e ' < a 1 * . "- F' 5
a "l y-
, , l i , , < m . ! l %
a o i . l M i i w p ul ! I
.m ,
m , w,mm J 1., J Li
~mg .i 3 g.
o N,- N9ipiFilik .n r
- m ,- -m, t r , BUT97ERIPilMM _
z 4 m'
'M's y '
s l
'l ; " -
l h l .
! e . t.n ,L
- , y,' ,
, nu a
l l
! ' [':[ ' II, y,2 l ,
n, l Q er ! ,
,a 'l__
i.,_,.m _ _ . . _m,, Jf n .- __ n
'-c "i
- - * - + -
-1, - ~ ~ ~ = ~ ,
n.~__, i o n
. t -
1m 1 ri ep
- e. c .
a s 1 e , e s ,
. =
- I __
i u 3 !. .l =
; 4 :
1, g b f kh.M [bl a o it w a i wu , w _ _; c a 1.
.y.m name= -
CONNECTING
. COVER GAP ELEMENT SHEET NEUTBON ABSORBER SHEET (BORAFLEX)
O REGION 1 F(GURE 3,3 3 x 3 TYP. ARRAY 3-8
-- - - - - - . - . . - - _ _ __ ~"'"-Newww,y--p_ _ ^-~*h~-w,_,,_
O .: . c s ' NEUTRON / C. 3 2 ,, ABSORBER
~
INTERNAL SQUARC BOX PITCH (BORAFLEX)\ . = _=i . ,- . . . . t_---
-.-, lq j . ,
1 -. .
' a e
l s D
._.T7 - . T-, . 'I hN
(
-. x , r ; 2
_F , g
- q _ f.he i . . ' ' 1 .
i O _? ]J
't. _
3 i A j.C ._ {
'm '~' -
l , , 1 ' i . I
- l e : ; ,
, 4
[~ '1 a. l I
". ~. _ ..
i ~ f .;f -
* ~.
I CONNECTING COVER GAP ELEMENT SHEET NEUTRON ASSORSER (BORAL SHEET) O REGION 1 FI G U R E 3 ,J a 3 x 3 TYP. ARRAY 3-9
O '
+.016 4.435 .
r
- ,4 -
- en - * - m
- O
+ i - - .06 o e
1 - S. l - m lO !
- .0 6 R" 1 '1 ks i .
l FIGURE 5.2 ' CHANNEL ELEMENT - REGION 1 & 2 (2 FOR SQU - - \
. . m -
s W3 C w w a C o h C
; s ec :
- c e ,' m c m,
2
-~g- ' 0 7 5,,
w -;- ,o75
- e ; s o l-o--r 6
l O ricuRE2.3 are EtsuEnrr REGION 1 3-10
--.,,,<--,,-..r,-,,,,-.,,,.----,-,_w-,,,,,,-,-,-,,,-_.,n- ,, e ,,,.,,--,,,,,v.-, - ,,-,-,_w-,n-- , , , , _,, .. - ,--w-,,,.----
O SQUARE 80X (FROM 2 FORMED CHANNELS )
. q 'R I > .. i'O 3 2"so, ~
8.'850 _
~ ~ ~
O 9 x COVER SHEET x 0 2 di,0 0 3" f> x
- w. .
f (BORAFLEX ONLY) tu:
' e x m xl vu N EUTRON ABSORBER (BORAL SHEET OR BORAFLEX)
I COMPOSITE BOX A S S E M B LY - R E G I O N FI G. 3.4 O 3-11
- . __ _b 365V 3, '; j j .m; ,
l-WELO 7 /4 LG. g [, [T ' ' '; ' y( (4) 51 DE S b b A_ l -[3-O ,,# An g y . y fl-"r{- y -
% l- , , ,
sus 2 . . .
-~~_u;_, .. _ . . . _ _...
B ASE PLATE 8.850"r ^.,, ~
= S o.
( T Y P. )
.z-
_z FUEL ASSEM BLY l o
-\tbREF.
O 11I.I_D TLi M V
/yyp' ,I !, I I
i I!') o a
, ,l L' --! i ;- - BORAFLEX '; I. .
- '- -1 O- l l l--- i ,
ni-BORAL SHEET a . L _i -, __, .
)
- f m POISON /
, . I - H i, e COVER SHEET $ ,e l
b r : _al : i . u m : l l
' I .
il
', 7C ~- 9 l l ' 'l M '
i .
- 1 i_ _ __~l ,. s. --l 12 : l F- - !<': 25 .. , . I l DIA.
a
; l l l '! l l_. _ -
I '
- FLOW HOLE l !: , - -l : l' lj , 1 .
[,- l
't[
f ' l ' l i .l B ASE PL ATE
,e , , .. ,y:- , ,
U a , M' ry , , o 'j --
}
N " DI 6 A H O L E ( T Y P. ) Q FIG. 3.5a TYPICAL CELL ELEV ATION -REGION 1 . 3-12
TYP) / TYP. BORAFL YP) , 3
.U . . . = ,,-, 1 , ,; - W E L D 7 /4 LG.
(4) SIDES l(T of j j. (T . []-- - O ,,# g< v.4 v , 7- (f.o B ASE PLATE
'"' '~ ' ^
_ l' . . _ , N 65 = T, s o.0" ^ . ,,. g ,,z,, ( T Y P. ) - FUEL ASSEMBLY
- I I. /2 REF, u " 7,1I V 'TYP. \ $_1_I_D -
i l l L _ '- 4 3 i r- -7 L- --d .O I , ;
-j - BORAL SH'EET 'g ~
b _t , I I e POISON l
@ su hr r ._1 i -> ' I e I M i f-7-5 _ - l
_ i
'F2 ~l , "Zs 7Z' l . ..
I I DIA. l L_ . ._] FLOW HOLE r _, f l l B ASE PL ATE
- i. --,
,e _,
3_ y=~,' U e , o ( } E
' S" D I A H O L E ( T Y P. )
O FIG. 3.5a,I TYPIC AL .CEL L ELEV ATION -REGION 1 .. 3 13
.. J/ g '\(T Y P.) *p d WELD 7.7 5' LG.
e 3 l . .., . - . . -, .-- h _- , (4) SIDES
~ '
i i l 8.85* ' ' ' ' 2.===2 SQ. BASE PLATE (T Y P.) FUEL ASSEMBLY 1.5' R E F. iji qi
./TT i s. _
L!..L Lj I l
. i ! L ,-.-d ', POISON l F- - ~-I l O i i,, - -d , h' 2M I - l e '
l 8 * - e v- l t . , l e .
, C' -.-Q \ \ ..
e E l E l
l I
r_-_q l 5 I '
'e l l l>
s . a l ' e i e l f-q p--- l; l1. " D I A . F L O W H O L E -
' ' I .i l_ i.- l g B ASE PL ATE l
j
> ! F~ 2 ~-f T . / \ .
t O Ne Di A. sOte (TYe.) FIG _3. 5 b T_Y P I C A _L__C E L L E L E V A T I O N - R E G I O N 2 3-14
1 0 - r I 5'/A 0fA. f S ASE PL ATE l, i l' E - 5"O 1 A.W . i t l l 2" i 1 s y 3 *i
^ s kN ,
i v N / 4
~
(81-1OIA. - 6 ae stow soteS
- O .
u , .. 2 : i o ' ' 1 (8)3o/4 G U S S E T S E , l,
=
p 4"- 4 U N C
*s ' 7 -
l t
- ,! ' CLASS lA ->i < ! :N < ~2 ) i +b .. o ,
y 1 W\ _(i 'I. V
,6, i K I DI A. H O LE o ,
I s /2 O f A . = l l_
=
l 8" OfA. l SUPPORT FI G. 3. S A A DJ U ST A B L E Q 3-15
I Q . l 5"D I A , _ BASE PL ATE l
' l- . ) 1 ' I 5 .A j/ / ). l\ l, i
o 2 4~ 1
//
f// n: . . , s>s.
/
ee ,e a e/ !A A, f _q u so. IV l , me i =$ M ',' 8-
@ -) j'% \
N O ' r ' ' ' 72 4unc CLASS lA (4) I DI A. HOLES tl
- 10 '/2 0,D. =
FIG 3,6 B SHORT ADJUSTABLE SUPP T, 1 e O mw-, 3-16
1 l O 8ASEk Sf . a g ,,, v BOS p/PE
- & . hs ).
7 ,z
'/e 7 -
I o \(4) / '$ Notts p, o .. O* A l# @ 90
- aper U y
/l6 l 4 %'mx,7(b.o.
( So 4 f.)
/
SmLL preo somr 1 FIGllRE 3.6C O 3-17
1 l \ O d. 04
- 9. c 3 INTERNAL SQ. BOX !
7 PITCH L I 1 I I l K?W W D 1 PC C t h M l.'- ~ [ C0CCDCCd I ' f , f5 h 'l s f> ' L I % i . ' - 5 -. o - r l L . _ I BCrMXXQ l INI I BcOOXXNS I - l l t i la ' ouNo4R3 0 - i met - i xxcwws i . - r-i - p _ I I E COVER
) SHEET . ~ > e n [ mm'(A 1 I N'YM 1 I WN I I l i i i I I i i $ l l
l ! NEUTRON ASSORBER CONNECTING STRIPS
- o 4 *og c8 rxx SHEET (BORAFLEX) 2 FIGURE 3 7 3X3 TYPICAL ARRAY REGION O
3-18 NOTE: Connecting strips vary in thickness from .040" to .080", resulting in average wall to wall separation of 06"
- .n as-welded cond :icn.
() 4.0 NUCLEAR CRITICALITY ANALYSIS 4.1 DESIGN BASES The high density spent fuel storage racks for the 3yron Nuclear Power Station are designed to assure that the neutron multiplication factor (keff) is equal to or less than 0.95 with the racks fully loaded with fuel of the highest anticipated reactivity in each of two regions, and flooded with unb", rated , water at temperature corresponding to the highest reactivity. The maximum calculated reactivity includes a margin for uncertainty in reactivity calculations and in mechanical tolerances, statistically combined, such that the true keff will be equal to or less than 0.95 with a 95% probability at a 95% confidence level. Applicable codes, standards, and regulations, or pertinent sections thereof, including the following: o 0 General Design criterion 62, Prevention of Criticality in Fuel V Storage and Handling. 0 USNRC Standard Review Plant, NUREG-0800, Section 9.1.1, New Fuel Storage, and Section 9.1.2, Spent Fuel Storage. O USNRC letter of April 14, 1978, to all Power Reactor Licensees - OT Position for Review and Acceptance of Spent Fuel Storage and Handling Applications, including modificaticn letter dated January 18, 1979. 0 USNRC Requiatory Guide 1.13, Spent Fuel Storage Facility Design Basis, Rev. 2 (proposed), December 1981. 0 USNRC Regulatory Guide 3.41, Yalidation of Calculational Methods for Nuclear Criticality Safety (and related ANSI N16.9-1975). O ANS I/ AN S-57. 2- 1983, Design Requirements for Light Water Reactor Spent Fuel Storage Facilities at Nuclear Powar Plants. O ANSI N210-1976, Design Objectives for Light Water Reactor Spcnt Fuel Storage Facilities at Nuclear Power Plants. : 0 ANS I N18.2-1973, Nuclear Safety Criteria for the Design of Stationary Pressurized Water Reactor Plants. to essure tne true reectivity ~411 ei evs be 1ess tnen tne ceicuieted O reactivity, the following conservative assumptions were made: 4-1
O O Moderator is pure, unborated water at a temperature corresponding to the highest reactivity. O Lattice of storage racks is assumed infinite in all directions, i.e., no credit is taken for axial or radial neutron leakage (except in the assessment of certain abnormal / accident conditions). O Neutron absorption in minor structural me2ers is neglected, i.e., spacer girds are replaced by water. The design basis fuel assembly is a 17x17 Westingnouse optimized fuel assembly containing U02 at a maximum initial enrichment of 4.2% U-235 by weight, corresponding to 48.6 grams U-235 per axial centimeter of fuel assembly. Two separate storage regions are provided in the spent fuel storage pool, with separate criteria defining the highest anticipated reactivity in each of the two regions as follows: 0 Region 1 is designed to accomodate new fuel with a maximum enrichment of 4.2 wt% U-235, or spent fuel regardless of the discharge fuel burnup. -Q 0 Region 2 is designed to accommodate fuel of various initial enrichments which have accumulated minimum burnups within an acceptable bound as depicted in Figure 4.1. 4.2
SUMMARY
OF CRITICALITY ANAL)SES 4.2.1 Normal Operating Conditions The criticality analyses of each of the two separate regions of the spent fuel storage pool previously described are summarized in Table 4.1 for the anticipated normal storage conditions. The calculated maximum reactivity in Region ? includes a burnup-dependent allowable for uncertainty in depletion calculations and, further:nore, provides an additional margin of more than 25 W k below the limiting ef fective multiplication factor (keff) of 0.95. As cooling time increases in long-term storage, decay of Pu-241 results in a significant decrease in reactivity, which will provide an increasing suberiticality margin and tends to further compensate for any uncertainty in depletion calculations. Spacing between the two dif ferent rack modulus is suf ficient to preclude adverse nuclear interaction between modules. 4-2
O aesion 2 cea eccommodete fuei of verious iaitiei enrichmeats ead discner9e fuel bJrnups, provideo the combination f alls within the acceptable domain illustrated in Figure 4.1. For convenient reference, the minimum burnup values in Figure 4.1 have been fitted by linear tangents at various values and the results are tabulated below. Initial Mini .am Initial Minimum Enrichment, % Burnup, MWD /MTV Enrichment, 5 Burnup, MWD /MTU 1.52 0 3.00 22,490 1.80 5,230 3.20 25,000 2.00 8.750 3.40 27,510 2.20 11,340 3.60 30,020 2.40 14,390 3.80 32,540 2.60 17,050 4.00 34,960 2.80 19,700 4.20 37,370 Linear interpolation between the tabulated values will always yield values on or conservatively above the curve of limiting burnups. g These data will be implemented in appropriate administrative procedures to d assure verified burnup as specified in draf t Regulatory Guide 1.13, Revision
- 2. Administrative procedures will also be employed to confirm and assure the presence of soluble poison in the pool water during fuel handling operations, as a further margin of safety and as a precaution in the event of fuel misplacement during fuel handling operations as discuss in Section 4.2.2.
4.2.2 Abnormal and Accident Conditions Although credit for the soluble poison normally present in the spent fuel pool I water is permitted under abnormal or accident conditions,* most abnormal or l accident conditions will not result in exceeding the Limiting reactivity (keff of 0.95) even in the absence of soluble poison. Tne effects on reactivity ot j credible abnormal and accident conditions are summarized in Table 4.2. Of these abnormal / accident conditions, only one has the potential for a more than negligible positive reactivity ef fect, l l i O
- Double contingency ' principle of ANSI N16.1-1975, as specified in the April 14, 1978, NRC letter (Section 1.2) and implied in the proposed revision (draf t) to Regulatory Guide 1.13 (Section 1.4, Appendix A).
4-3 l t
r ()\ The inadvertent misplacement of a new fuel assembly (either into a Region 2 storage cell or outside and adjacent to a rack module) has the potential for exceeding the limiting reactivity should there be a concurrent and independent accident condition resulting in the loss of all soluble poison, Administrative procedures to assure the presence of soluble poison during fuel handling operations will preclude the possibility of the simultaneous occurrence of these two independent accident conditions. The largest reactivity increase occurs from accidentally placing a new fuel assembly into
- Region 2 storage cell with all other cells fully loaded. Under this
..cnd' :;on, the presence of only 300 ppm soluble boron assures that the infinite multiplication factor would not exceed the design basis reactivity for Region 2. With the nominal concentration of soluble poison present (2000 ppm baron), the maximum reactivity, Bk , is less than 0.95 even if Region 2 were to be fully loaded with fresh fuel of 4.2% enrichment.
4.2.3 New Fuel Storage
/^\
U Region 1 of the storage racks is designed to safely accommodate new unirradiated fuel of 4.2% enrichment, when fully flooded with clean unborated water. Under certain circumstances, it may be desirable to store new fuel in the dry condition in Region 1 or to utilize Region 2 for the temporary storage of new fuel, either dry or fully flooded. These conditions were analyzed to assure the acceptability of Region 1 in the dry condition and to determine an arrangement in Region 2 that would assure criticality safety in conformance with the requirements of SRP 9.1.1, "New Fuel Storage." rriticality analyses confirmed that Region 1 does not exhibit a peak in reactivity at low moderator densities (e.g., for or foam mderation) and that the optimum nederation (highest keff) occurs for the fully flooded condition. This condition is the design basis for Region 1 where the maximum k8 , including all uncertainties, is less than 0.947. In Region 2, it was determined that a checkerboard pattern (fuel assemblies 0 d aligned diagonally) provided an acceptable kB , in either the fully flooded or l the dry (low density moderation) condition for new fuel assemblies of 4.2% 4-4
1 b enri chment. These calculations indicated a nominal kB of 0.813 + 0.014 (1s) when fully flooded with clean unborated water--a value substantially less than the limiting keff of 0.95, even with an additional allowance for uncertainties (maximum Bk , of 0.86 at 95%/95% tolerance limits). Calculations, using Monte Carlo techniques, did not reveal a peak in reactivity at low moderator densities, and fully flooded condition corresponds to the highest reactivity (optimum moderation). Thus, the checkerboard pattern of new 4.2% enriched fuel in Region 2 represents a safe configuration in conformance with SRP 9.1.1. and 9.1.2. 4.3 REFERENCE FUEL STORAGE CELL 4.3.1 Reference Fuel Asserrbly The design basis fuel assembly, illustrated ir Figure 4.2, is a 17x17 array of fuel rods with 25 rods replaced by 24 c9ntrol rod guide tubes and 1 instrument O b thi rrbl e. Table 4.3 summarizes the Westinghouse optimized fuel assembly (0FA) design specifications and the expected range of significant variations. 4.3.2 Region 1 Storage Cells The nominal spent fuel storage cells used for the original critically analyses of Region 1 storage cells is shown in Figure 4.2.1. The rack is composed of Boraflex absorber material sandwiched between a 0.060-inch inner stainless steel box and a 0.020-inch outer stainless steel (SS) coverplate (0.125-inch coverplate for module periphery cell walls). The fuel assemblies are ( central!y located in each storage cell on a nominal lattice spacing of 10.320 l
+ 0.050 inches in one direction and 10.420 + 0.050 inches in the other direction. Stainless steel gap channels connect one storage cell box to another in a rigid structure and define an outer water space between boxes.
This outer water space constitutes a flux-trap between the two Boraflex absorber sheets that are essentially opaque (black) to thermal neutrons. The Boraflex absorber has a thickness of 0.075 + 0.007 inch and a nominal B-10 areal density of 0.0238 gram per cm2, 4-5
C) Figure 4.2.2 shows the nominal cell modified with Boral inserts. The Boral 2 has a minimum B-10 loading of .020gm/cm and thickness of .075 1 004 inch. The minimum boron loading is combined With the maximum thickness as a "worst case" conservatism. The Boral plates are to be flush with the Boraflex cover plate. However a manufacturing tolerance allowing .125 inch clearance is assumed for the analysis. Figure 4.2.3 shows the cell with Boral replacing Boraflex. The cover plate has been eliminated as it serves no useful purpose. The Boral has a minimum 2 B-10 loading of 0.025 gm/cm and a thickness of .085 *.004 inch. These are again used simultaneously for a "worst case" calculation. 4.3.3 Region 2 Storage Cells Region 2 storage cells were initially designed for fuel of 3.2 wt% U-235 initial enrichment burned to 25,000 MWD /MTU and extended to encompass fuel of m Q 4.25 initial enrichment burned to 37,370 MWD /MTV. In this region, the storage cells are composed of a single Boraflex absorber sandwiched between the 0.060-inch stainless steel walls of adjacent storage cells. These cells shown in Figure 4.3, are located on a lattice spacing of 9.03 1 0.04 inches. 4.4 ANALYTICAL METHODOLOGY 4.4.1.1 Reference Analytical Methods and Bias for Nominal Case The CASMO-2E computer code (References 1, 2 and 3), a two-dimensional multigroup transport theory code for fuel assemblies, has been benchmarked (see Appendix A) and is used both as a primary method of analysis and as a means of evaluning small reactivity increments associated with manufacturing tolerances. CASM0-2E benchmarking resulted in a calculational bias of 0.0013 1 0.0018 (95%/95%). In fuel rack analyses, for independent verification, criticality analyses of O tne nign density spent fuei storege recks were ei - performed itn tne ^Mex-KENO computer package (References 4 and 5), using the 27-group SCALE
- cross-4-6
section library (Reference 6) with the NITAWL subroutine for U-238 resonance shielding ef fects (Nordheim integral treatment). Details of the benchmark calculations with the 21-group SCALE cross-section library are also presented in Appendix A. These benchmark calculations resulted in a bias of 0.0106 _+ _ 0.0048 (95%/955). In the geometric model used in KENO, each fuel rod and its cladding were described explicitly. For two-dimensional X-Y analysis, a zero current (white albedo) boundary condition was applied in the axial direction and, for Region 1, at the centerline through the outer water space (flux-trap) on all four sides of the cell, effectively crating an infinite array of storage cells. In Region 2, the zero current boundary condition was applied at the center of the Boraflex absorber sheets between storage cells. The AMPX-KEN 0 Monte Carlo calculations inherently include a statistical u.. certainty due to the random nature of neutron tracking. To minimize the statistical uncertainty of the KENO-calculated reactivity, a total of 50,000 neutron histories is normally accumulated for each calculation, in 100 generations of 500 neutrons each. CASM0-2E is also used for burnup calculations, with independent verification by EPRI-CELL and NULIF calculations. In tracking long-term (30 year) reactivity ef fect of spent fuel stored in Region 2 of the fuel storage rack, EPRI-CELL calculations indicate a continuous reduction in reactivity with time (af ter Xe decay) due primarily to Pu-241 decay and Am-241 growth. A third independent method of criticality analysis, utilizing diffusion / blackness theory, was also used for additional confidence in results of the primary calculational methods, although no reliance for criticality safety is placed on the reactivity value from the diffusion /blac< ness theory technique. This technique, however, is used for auxiliary calculations of small incremental reactivity ef fects (e.g., axial cutback or mechanical tolerances) that would otherwise be lost in normal KEN 0 statistical variations, or would be inconsistent with CASMO-2E geometry limitations. lO
- SCALE is an acronym for Standard Computer _A,nalysis for Licensing Evaluation.
4-7
Il v Cross sections for the dif fusion / blackness theory calculations were derived frnm CaSM0-2E or calculated by the NULII computer code (Reference 7), supplemented by a blackness theory routine that effectively imposes a transport theory boundary condition at the surface of the Boraflex neutron absorber. Two different spatial diffusion theory codes, PDQ07 (Reference 8 in two dimensions and SNEID* in one dimension, were used to calculate reactivities. The two-dimensional PDQ07 code was used to describe the actual storage cell geometry, with NULIF cell-homogenized constdnts representing each fuel rod and its associated water moderator. SNEID is one-dimensional model, in cylindrical or slab geometry, used for the calculation of axial cutback reactivity ef fects and in the assessment of abnormal occurrences. 4.4.1.2 Analytical Methods and Bias for Modified Cases The design method used to demonstrate the criticility safety of fuel in the modified fuel storage racks uses the AMPX system of computer codes [4] for
) cross-section generation (NITAWL-S and XSDRNPM-S), and the KENO-IV (CRC) code
[5] for reactivity detennination. Tolerance variations of the nominal case were then applied. An ORNL 227 energy group ENDF/8-V cross-section library [15,16] was utilized for the criticality analysis. The NITAWL program computes resonance cross-sections and prepares a 227 neutron group P3 cross-section data file for use in the XSDRNPM one-dimensionalnS neuron transport code. The Nordheim Integral Treatment [17] used. Energy and spatial weighting of cross-sections to reduce the 227 group data to 27 group (SCALE) date is performed by the XSDRNPM program. Zone weighted cross-sections are in used in the fuel regions that are modeled as discrete fuel rods in KENO. Discrete modeling of structural regions such as the storage cell structures also use zone weight Cross-sections.
*SNEID is one-dimensional diffusion theory routine developed by Black & Veatch and verified by comparison with PDQ07 one-dimensional calculations.
4-8
O KENO-IV is used in the evaluation of criticality of the Byron storage racks. KENO-IV is a three-dimensional Monte Carlo theory program designed for reactivity calculations. The calculational method and cross-section values which were used in the criticality analysis of the fuel storage racks have been verified by comparison with critical experiment data for assemblies similar to those for which the racks are designed. This benchmarking data is sufficiently diverse to establish that thte method bias and uncertainty will apply to fuel storage rack conditions which include strong neutron absorbers and low moderator densities. A set of 33 critical experiments has been analyzed to demonstrate its applicability to criticality analysis and to establish the method bias and variability. The experiments range from water moderated oxide fuel array separated by various materials (BORAL, steel and water) that simulate LWR fuel O shipping and storage conditions [18], to dry, harder spectrum uranium metal cylinder arrays with various interspersed materials [19] (Plexiglas, paraffin and air) that demonstrate the wide range of applicability of the method. The results and some descriptive factors about each of the 33 benchmark critical experiments are give in Table 4.9. The average Keff of the benchmarks is 0.9917 which demonstrates that there is a 0.0083 bias associated with the method. The standard deviation of the Keff values is 0.00082 celta-
- k. The 95/95 one-side tolerance limit factor for 33 values is 2.19.
There is a 95 percent probability with 95 percent confidence level that the uncertair.ty in reactivity, due to the method, is not greater than 0.0018 delta-k. 4.4.2 Fuel Burnup Calculatio,7s Fuel burnup calculations in the hot operating condition were performed primarily with the CASM0-2E code. However, to enhance the credibility of the burnup calculations (in lieu of critical experiments), the CASM0-2E results were independently checked by calculations with the NULIF code (Reference 7) and with EPRI-CELL (Reference 9). Figure 4.4 compares results of these 4-9
m b independent methods of burnup analysis under hot reactor operating conditions. The results agree within 0.003 Wk in the hot operating condition. In addition to depletion calculations under not operating conditions, reactivity comparisons under conditions more representative of fuel to be stored in the racks (cold, xenon-free) are also significant in storage rack criticality analyses. Table 4.4 compares the clod, xenon-free reactivities calculated by CASM0-2E, NULIF/PDQ07, and EPRI-CELL. In the cold condition, the CASMO-2E calculations gave a slightly higher reactivity value for the Region 2 fuel storage cell and the good agreement generally observed lends credibility to the calculations, particularly in view of the know bias and uncertainty in CASMO-2E calculations (Appendix A). No definitive methods exists for determining the uncertainty in burnup-dependent reactivity calculations. All of the codes discussed above have been used to accurately follow reactivity loss rates in operating reactors. CASMO-2E has been extensively benchmarked ( Appendix A; References 1, 2, 3 and 'd 10) against cold, clean, critical experiments (including plutonium-bearing fuel), Monte Carlo calculations, reactor operations, and heavy-element concentrations in irradiated fuel. In particular, the analyses (Reference 10) of 11 critical experiments with plutonium-bearing fuel gave an average keff of 1.002_+ 0.011 (95%/95%), showing adequate treatmant of the plutonium nuclides. In addition, Johansson (Reference 11) has obtained very good agreement in calculations of close-packed, high-plutonium-content, experimental configurations. Since critical-experiment data with spent fuel is not available, it is necessary to assign an uncertainty in reactivity based on other consideratior.s, supported by the close agreement between different calculational methods and the general industry experience in predicting reactivity loss rates in operating plants. Over a considerable portion of the burnup, the reactivity loss rate in PWRs is approximately 0.01 Wk for each 1,000 MWD /MTU, becoming somewhat smaller at the higher burnups. By conservatively assuming an uncertainty in reactivity
- of 0.5 x 10-6 times the burnup in MWD /MTV, a burnup-dependent uncertainty is defined that increases 4-10
( with increasing fuel burnup, as would be reasonably expected. This assumption provides an estimate of the burnup uncertainty that is more conservative and bounds estimates frequently employed in other fuel rack licensing applications (i.e., 5% of the total reactivity decrement). Table 4.5 summarizes results of the burnup analyses and estimated uncertainties. These uncertainties are appreciably larger, in general, than would be suggested by the industry experience in predicting reactivity loss rates and boron let-down curves over many cycles in operating plants. The increasing level of conservatism at the higher fuel burnups provides an adequate margin in the uncertainty estimate to accomnodate the possible existence of a small positive reactivity increment from the axial distribution in burnup (see Section 4.4.3). In addition, although the burnup uncertainty may be either positive or negative, it is treated as an additive term rather than being cont >ined statistically with other uncertaincies. Thus, the allowance for uncertainty in burnup calculations is believed to be a conservative estimate, p;rticularly in view of the substantial reactivity decrease with aged fuel as discussed in Section 4.4.4. 4.4.3 Ef fect of Axial Burnup Distribution Initially, fuel loaded into the reactor will burn with a slightly skewed cosine power distribution. As burnup progresses, the burnup distribution will tend to flatten, becoming more highly burned in the central regions than in the upper and lower ends. This ef fect may be clearly seen in the curves compiled in Reference 12. At high burnup, the more reactive fuel near the ends of the fuel assembly (less than average burned) occurs in regions of lower reactivity worth due to neutron leakage. Consequently, it is expected that distributed-burnup fuel assemblies would exhibit a slightly lower reactivity than that calculated for the average burnup. As burnup progresses, the distribution, to some extent, tends to be self-regulating as controlled by the axial power distribution, precluding the existence of large regions of significantly reduced burnup. A nunt)er of one-dimensional dif fusion theory analyses have been made based upon calculated and measured axial burnup distributions, these analyses confirm the minor and generally negative reactivity ef fect of axially 4-11
' distributed burnup. The trends observed, however, suggest the possibility of a small positive reactivity ef fect at the high burnup values, and the uncertainty in kB due to burnup, assigned at the higher burnups (Section 4.4.2), is adequately conservative to encompass the potential for a small positive reactivity effect of postulated axial burnup distributions.
Furthermore, reactivity decreases with time in storage (Section 4.4.4), and, in addition, there is a large margin in reactivity (>0.02 Wk) below the limiting keff value (0.95) which can accommodate any reasonable reactivity effects that might be larger than expected. 4.4.4 Long-Term Decay Since the fuel racks in Region 2 are intended to contain spent fuel for long periods of time, calculations nr* made using EPRI-CELL (which incorporates the CINDER code) to follow the kor,g-term changes in reactivity of spent fuel over a 30-year period. CINDER tracks the decay and burnup dependence of some 179 fission products. Early in the decay period, xenon grows in (reducing reactivity) and subsequently decays, with the reactivity reaching a maximum at 100-200 hours. The decay of Pu-241 (13-year half-life) and growth of Am-241 substantially reduce reactivity during long term storage, as indicated in Table 4.6. The reference design criticality calculations do not take credit for this long-term reduction in reactivity, other than to indicate an increasing subcriticality margin in Region 2 of the spent fuel storage pool. 4.5 REGION 1 CRITICALITY ANALYSIS AND TOLERANCE VARIATIONS 4.5.1 Nominal Design Case Under normal conditions, with nomir.21 dimensions, the kB values calculated by the three methods of analysis are as follows: Maximum k Analytical Method Bias-corrected kg (95%/95%)B CASM0-2E 0.9387 1 0.0018 0.9405 g 0.9362 Q AMPX-XENO Diffusion blackness theory 0.9301 1 0.0061 0.9393 0.9393 4-12
O Tne ^MeFx-xEw0 ceicuietioes iaciede e one-sided toler>nce fector (aeference
- 13) of 1.799 corresponding to 957, probability at a 9M confidence limit. For the nominal design case, the CASMO-2E calculation yields the highest (most conservative) reactivity, and, therefore, the independent verification calculations substantiate CASM0-2E as the primary calculational method.
4.5.2 Boron Loading Verification The Boraflex absorber sheets used in Region 1 storage cells are nominally 0.075-inch thick, with a B-10 areal density of 0.0238 g/cm2 . Independent 2 manufacturing tolerance limits are +0.007 inch in thickness and +0.0017 g/cm in B-10 content. This assures that at any point where the minimum boron concentration (0.0221 grsm B-10/cm2 ) and minimum Baraflex thickness (0.068 inch) may coincide, the boron-10 areal density will not be less than 0.020 gram /cm22 Differential CASMO-2E calculations indicate that these tolerance limits result in an incremental reactivity uncertainty of +0.0021 Wk for boron content and +0.0047 for Boraflex thickness variations. (q> ~ 4.5.3 Storage Cell Lattice Pitch Variation The design storage cell lattice spacing between fuel assenblies in Region 1 is 10.32 inches in one direction and 10.42 inches in the other direction. A decrease in storage cell lattice spacing may or may not increase reactivity depending upon other dimensional changes that may be associated with the decrease in lattice spacing. Increasing the water thickness between the fuel and the inner stainless steel box results in a small increase in reactivity. The reactivity ef fect of the flux-trap water thickness, however, is more significant, and decreasing the flux-trap v.ater thickness increases rea cti vi ty. Both of these effects have been evaluated for independent design tolerances. The inner stainless steel box dimension, 8.850 _+ 0.032 inches, defines the inner water thickness between the fuel and the inside box. For the tolerance limit, the uncertainty in reactivity is +0.00'.8 Wk as determined by Q _, dif ferential CASM0-2E calculations, with kB increasing as the inner stainless steel box dimension (and derivative lattice specing) increases. 4-13 l --
O. V The design flux-trap water thicknesses are 1.160 + 0.040 inches and 1.260 + 0.040 inches, which result in an uncertainty of +0.0033 Wk due to the tolerance in flux-trap water thickness, assuming the water thickness is simultaneously reduced on all four sides. Since the manufacturing tolerances on each of the four sides are statistically independent, the actual reactivity uncertainties would be less than +0.0038, although the nore conservative value has been used in the criticality evaluation. 4.5.4 Stainless steel Thickness Tolerances The nominal stainless steel thickness in Region 1 is 0.60 inch for the inner stainless steel box and 0.020 inch for the Boraflex coverplate (0.125 inch on module boundary). The maximum positive reactivity effect of the expected stainless steel thickness tolerance variations, statistically combined, was calculated (CASMO-2E) to be +0.0025 Wk. 4.5.5 Fuel Enrichment and Density Variation The design maximum enrichment is 4.20 + 0.05 wt% U-235. Calculations of the sensitivity to small enrichment variations by CASM0-2E yielded a coefficient of 0.0047 Wk per 0.01 wt% U-235 at the design enrichment. For a tolerance on U-235 at the design enrichment. For a tolerance on U-235 enrichment of +0.05 in wt%, the uncertainty on 'rg is +0.0024 Wk. fuel density increased to maximum value of Calculations were made with the U02 97% theoretical density (TD). For the midrange value (95% TD) used for the reference design calculations, the uncertainty in reactivity is +0.0026 Wk over the range of U02 densities expected. 4.5.6 Boraflex Width Tolerance Variation The reference storage cell design for Region 1 (Figure 4.2) uses a Boraflex blade width of 7.75 + 0.0625 inches. A positive increment in reactivity L occurs for a decrease in Boraflex absorber width. For reduction in width of the maximum tolerance 0.0625 inch, the calculated positive reactivity increment is +0.0007 Wk. 4-14
(3 V 4.5.7 Axial Cutback of Boraflex The axial length of the Boraflex poison material is less than the active fuel length by three inches at the top and at the bottom of the Revision 1 storage rack modules. To account for the reactivity effect of this axial cutback, one-dimensional (slab) diffusion theory calculations were made using flux-weighted homogenized diffusion theory constants edited from CASC-2E calculations of the array of storage cells, with and without Boraflex present. In the one-dimensional calculations, an infinite (30-cm) water reflector was used above and below the fuel assembly, with the lengths of the unpoisoned "cutback" regions, top and bottom, varied in a series of parametric calculations. Results of these calculations showed that the keff remains less than the kB of the referer.ce central storage cell region, until the axial cutback exceeds four inches top and bottom. Thus, the actual axial neutron leakage nore than compensates for the three-inch design cutback, and the reference infinite multiplication factor, k , Bremains a conservative O overest4=ete or tae tr e react <v4tv 4.6 REGION 2 CRITICAllTY ANALYSIS AND TOLERANCE VARIATIONS 4.6.1 Nominal Design Case The principal method of analysis in Region 2 was CASMO-2E code, using the restart option in CASMO to transfer fuel of a specified burnup into the O storage rack configuration at a reference temperature of O C. Calculations were made for fuel of several dif ferent initial enrichments and, for each enrichment, a limiting kg value established which included an additional f actor for uncertainty in the burnup analysis and for the axial burnup distribution. The restart CASM0-2E calculations (cold, clean, rack geometry) were then interpolated to define the burnup value yielding the limiting kB value for each enrichment, as indicated in Table 4.7. These converged burnup values define the boundary of the acceptable domain shown in Figure 4.1. O 4-15
O At e buraun of 37.000 swo/stu. the seasitivity to burn n is ceicuietea to be
-0.0079 Wk per 1000 MWD /MTu, During long-term storage, the kB values of the Region 2 fuel rack will decrease continuously from decay of Pu-241 as indicated in Section 4.4.4.
Two independent calculational methods were used to provide additional confidence in the reference Region 2 criticality analyses. Fuel of 1.5% initial enrichment (approximately equivalent to the reference rack design for burned fuel) was analyzed by AMPX-KENO (27-group SCALE cross-section library) and by the CASM0-2E m) del used for the Region 2 rack analysis. For this case CASM0-2E kB (0.9014) was within the statistical uncertainty of the bias-corrected value (0.904310.0030 (1s) obtained In the AMPX-KEN 0 calculations. This agreement confirms the validity of the primary CASM0-2E calculations. The second independent method of analysis used the NullF code for burnup analysis, and for generating diffusion theory constants (cold, clean) for the NulIF-calculated composition at 25,000 Mwo/MTu with fuel of 3.2s initial O enri chment. These constants, together with blackness theory constants for the Boraflex absorber, where then used in a two-dimensional PDQ07 calculation for the storage rack configuration. Results of this calculation (kB of 0.9017) compared favorably with the CASM0-2E calculation for the same conditions (kB of 0.9061) and thus tend to confirm the validity of the Primary calculational method. 4.6.2 Boron Loading Variation The Boraflex absorber sheets used in the Region 2 storage cells are nominally 0.041 inch thick with B-10 areal density of 0.0130 g/cm2 . Independent 2 manufacturing limits are +0.007 inch in thickness and +0.0009 g/cm in B-10 content. This assures that at any point where the minimum boron concentration (0.01206 B-10 g/cm2 ) and the minimum Boraflex thickness (0.034 inch) may 2 coincide, the boron-10 areal density will not be less than 0.010 g/cm , Differential CASM0-2E calculations indicate that these tolerance limits result O in en incrementei reectivity uncerte4nty of 10.0028 wk for boren content end 10.0078 Wk for Boraflex thickness. 4- 16
G V 4.6.3 Storage Cell Lattice Pitch Variations The value used for the storage sell lattice spacing between fuel asseslies in Region 2 is 9.0310.040 inches, corresponding to an uncertainty in reactivity of 0.0011 Wk. 4.6.4 Stainless Steel Tnickness Tolerance The nominal thickness of the stainless steel box wall is 0.060 inch with a tolerance limit of 10.005 inch, resulting in an uncertainty in reactivity of 10.0001 Wk. 4.6.5 Fuel Enrichment and Density Variation Uncertainties in reactivity due to tolerances on fuel enrichment and UO2 density in Region 2 are assumed to be the same as those determined for Region 1. 4.6.6 Boraflex Width Tolerance The reference storage cell design for Region 2 (Figure 4.3) uses a Boraflex absorber width of 7.25 1 0.625 inches. For a reduction in width of the maximum tolerance, the calculated positive reactivity increment is 0.0009 Wk. 4.7 ABNORMAL AND ACCIDENT CONDITIONS 4.7.1 Eccentric Positioning of Fuel Asse21y in Storage Rack The fuel asse21y is normally located in the center of the storage rack cell with bottom fittings and spacers that mechanically limit lateral mvement of the fuel asse 211es. Nevertheless, calculations were made with the fuel assemblies moved into the corner of the storage rack cell (four-assembly cluster at closest approach). These calculations indicated that the reactivity decreases very slightly in both regions, as determined by PDQ07 calculations with diffusion coefficients
- generated by NULIF and a blackness 4- 17 t
theory routine. The highest reactivity therefore corresponds to the reference l] design with the fuel assemblies positioned in the center of the storage cells. 4.7.2 Temperature and Water Density Effects The moderator temperature coef ficient of reactivity in both regions is negative; a moderator temperature of 00 0, with a water density of 1.0 g/cm3 , was assumed for the reference designs, which assures that the true reactivity will always be lower, regardless of temperature. Tenperature effects on reactivity have been calculated and the results are shown in Table 4.8. Introducing voids in the water internal to the storage cell (to simulate boiling) decreased reactivity, as shown in the table. Voids due to boiling will not occur in the outer (flux-trap) water region of Region 1. With soluble poison present, the temperature coefficients of reactivity would O be exaected to differ fro = tnose '"'errea fro = tne aet ia tebie 4 8-However, the reactivities would also be substantially lower at all temperatures with soluble boron present, and the data in Table 4.8 is pertinent to the higher-reactivity unborated case. 4.7.3 Droped Fuel Assembly Accident To investigate the possible reactivity effect of a postulated fuel assent >1y drop accident, calculations were made for unpoisoned assemblies separated only by clean unborated water. Figure 4.5 shows the results of these calculations. From these data, the reactivity (kB ) will be less than 0.95 for any water gap spacing greater than 6 to 7 inches in the absence of any absorber nuterial, other than water, between assemblies. For a drop on top of the rack, the fuel assent)1y will come to rest horizontally on top of the rack with a minimum separation distance of more than 12 inches. Maximum expected O inis ceicuietionei engroecn wes aecessers since tne reectivit, effects ere too small to be calculated by KENO, and CASMO-2E geometry is not readily amendable to eccentric positioning of a fuel assent >ly. 4- 18
O deformation uader seismic or eccideat coaditions wiii aot reduce tne miaimom spacing between fuel assemblies to less than 12 inches. Consequently, fuel asse21y drop accident will not result in an increase in reactivity above that. calculated for the infinite nominal design storage rack. Furthermore, soluble boron in the pool water would substantially reduce the reactivity and assure that the true reactivity is always less than the limiting value for any conceivable fuel handling accident. 4.7.4 Abnormal Location of a Fuel Assembly The abnormal location of a fresh unirradiated fuel assedly of 4.25 enrichment could, in the absence of soluble poison, result in exceeding the design reactivity limitation (kB f 0.95). This could occur if the assembly were to be either positioned outside and adjacent to a storage rack module or loaded into a Region 2 storage cell, with the latter condition producing the larger positive reactivity increment. Soluble poison, however, is normally present in the spent fuel pool water (for which credit is permitted under these O coaditions) aad wouid mei#taia tne reectivitv substaatialiv iess tnaa tae design limitation. The largest reactivity increase occurs for accidentally placing a new fuel assembly into a Region 2 storage cell with all other cells fully loaded. under this condition, the presence of 300 ppm soluble boron assures that the infinite multiplication factor would not exceed the design basis reactivity. With the nominal concentration of soluble poison present (2000 ppm boron), the maximum reactivity, kB , is less than 0.95 even if Region 2 were to be fully loaded with fresh fuel of 4.2% enrichment. Administrative procedures will be f used to confirm and assure the continued presence of soluble poison in the ! spent fuel pool water during fuel handling operations. l 4.7.5 Lateral Rack Movement Lateral notion of the rack modules under seismic onditions could potentially j alter the spacing between rack modules. However, lirdle bars on the modules prevent closing the spacing to less than 2.0 inches, which is greater than the normal flux-trap water gap in the Region I reference design. Region 2 storage 4- 19
O cells do not use flux-trap and the reactivity is insensitive to the spacing between modules. Furthermore, soluble poison would assure that a reactivity less than the design limitation is maintained under all conditions. 4.8 NEW FUEL STORAGE 4.8.1 Storage in Region 1, Dry Region 1 is normally designed to acconnodate new unirradiated fuel asse211es in the fully flooded condition. For storage in the dry condition, the racks must also conform to the requirements of SRP 9.1.1 which specify a limiting keff value of 0.98 under optimum low density moderation. Calculations were made, using AMPX-KEN 0, for several hypothetical low-moderator densities down to 0.05 g/cc simulating fog or foam moderation. These calculations showed a continuously decreasing kB as the moderator density decreased, yielding a kB of 0.546 + 0.008 (1s) at 105 moderator density. Axial leakage was neglected in these alculations, but would substantially reduce the already low kB values. These results are consistent with the general observation that a low-density optimum-moderation peak in reactivity does not exist in poisoned racks (Reference 14). 4.8.2 Storage in Region 2, Flooded In a succession of trial-and-error calculations, it was found that a checkerboard storage pattern in Region 2 would allow new fuel assemblies of 4.25 enrichment to be safely accormodated without exceeding the limiting 0.95 In this checkerboard loading pattern, the fuel assemblies are keff value. located on a diagonal array, as illustrated below, with alternate storage cells empty of any fuel. 4-20
/N.
U
//// '///n //// //// ///, //// '///i //// //// '///<
r///, //// //// r///. ///s r///. //// //// r/// ////
//// '///i //// //// ///< //// '///i //// //// '///i ;/ USAUUE / . '/// // / / '/ / '///. //// //// '///; //// //// '///, //// //// '///4 '/// ///s //// r/// //// //// NON-r/// : ///s //// r///. '///> M USABLE //// '///a //// //// ' / / /. //// '///, //// //// '///,
(empty)
'///. ///s //// r/// '////
r/// //// //// f/// '///s
//// '///. //// //// '///< //// '///, //// //// '///.
r///, //// //// r///, ///s r///. ///J //// r///. '//// Monte Carlo calculations (AMPX-XENO) resulted in a kg of 0.8133 + 0.0138. O witn e one sided x-fector (aefereace 13) for 9ss probebiitty et e est confidence level and a Wk of 0.009 for uncertainties (Table 4.1 for Region 2), the maximum kg is 0.863, which is substantially less than the 0.95 limiting value. Thus, Region 2 may be safely used for the temorary storage of new fuel assemblies provided the storage configuration is restricted to te checkerboard pattern indicated above. 4,8.3 Storage in Region 2, Dry As indicated in Section 4.8.1 above, a peak in reactivity (keff) at low nederator densities is not expected for poisoned rack designs. AMPX-KEN 0 j calculations confirmed the absence of a low-moderator-density peak in Region 2 with 4.2% enriched fuel arranged in the checkerboard pattern. At 105 l l moderator density, the calculated kB was 0.552, which would be substantially reduced if axial leakage were to be included. Thus, Region 2 conforms to the requirements of SRP 9.1.1 (kB < 0.98 at optimum moderation) for the safe storage of 4.21 enriched fuel, dry, in the checkerboard loading pattern. 4-21 l
REFERENCES
- 1. A. Ahlin, M. Edenius, H. Haggblom, "CASMO - Fuel Assembly Burnup Program," AE-RF-76-4153, Studsvik report (proprietary).
- 2. A. Ahlin and M. Edenius, "CASMO - A Fast Transport Theory Depletion Code for LWR Analysis," ANS Transactions, Vol. 26, p. 604, 1977.
- 3. M. Edenius et al ., "CASMO Benchmark Report," Studsvik/RF-78-6293, Aktiebolaget Atomenergi, March 1978.
4 Green, Lucious, Petrie, Ford, White, Wright, "PSR-63/AMPX-1 (code package) AMPX Modular Code System for Generating Coupled Multigroup Neutron-Gamma Libraries from ENDF/B," ORNL-TM-36-706, Oak Ridge National Laboratory, March 1976.
- 5. L. M. Petrie and N. F. Cross, "KENO-IV, An Improved Monte Carlo Criticality Program," 0RNL-4933, Oak Ridge National Laboratory, Novetter 1975.
- 6. R. M. Westf all et al., "SCALE: A Modular Code System for Performing Standardized Computer Analyses for Licensing Evaluation," NUREG/CR-0200, 1979.
- 7. W. A. Wittkopf, "NULIF -Neutron Spectrum Generator, Few-Group Constant C' Generator and Fuel Depletion Code," BAW-426, The Babcock & Wilcox Company, August 1976.
- 8. W. R. Cadwell, PDQ07 Reference Manual, WAFD-TM-678, Bettis Atomic Power Laboratory, January 1967.
- 9. W. J. Eich, "Advanced Recycle Methodology Program, CEM-3," Electric Power Research Institute,1976.
- 10. E. E. Pilat, "Methods for the Analysis of Boiling Water Reactors (Lattice Physics)," YAEC-1232, Yankee Atomic Electric Co., December 1980.
- 11. E. Johansson, "Reactor Physics Calculations on Close-Packed Pressurized Water Reactor Lattices," Nuclear Technology, Vol. 68, pp. 263-268, February 1985.
- 12. H. Richings, Some Notes on PWR (W) Power Distribution Probabilities for LOCA Probabilistic Analyses, NRC Memorandum to P. S. Check, dated July 5, 1977.
- 13. M. G. Natrella, Experimental Statistics National Bureau of Standards, Handbook 91, August 1963.
- 14. J. M. Cano et al., "Supercriticality Through Optimum Moderation in Nuclear Fuel Storage," Nuclear Technology, Vol. 48, pp. 251-260, May
(] \ 1980. 4-22
- () 15. W. E. Ford !!!, et al, "CRSL-V
- Processed ENDF/B-V 227 Neutron Group and Pointwise Cross Section Libraries for Criticality Safety Reactor and Shielding Studies," ORNL/CSO/TM-160, June,1982.
- 16. R. Kinsey, Ed. "ENDF/B Summary Documentation," BNL-NCS-17541 (ENDF-201),
3rd Edition, Brookhaven National Laboratory,1979.
- 17. T. W. Nordneim, "The Theory of Resonance Absorption," Proceeding of Symposia in Applied Mathematics, Vol. XI, p. 58, Garrett Birkhoff and Eugene Wigner, Eds., AM. Math. Soc.,1961.
- 18. M. N. Baldwin et al, "Critical Experiment Supporting Close Proximity i Water Storage of Power Reactor Fuel," BAE-1484-7, July,1979. ..
- 19. J. T. Thomas, "Critical Three-Dimensional Arrays of U(93.2) Metal Cylinders "Nuclear Science and Engineering, Volume 52, pages 350-359, 1973.
i 5 ^ l I P
~
4-23 p
- v. .. .__ -. -
4 1 O rabie 4.1
SUMMARY
OF- CRITICALITY SAFETY ANALYSES Region 1 Region 1 Nominal Boral Inserts Minimum acceptable burnup 0 0 MWD /MTU G 4.2% initial enrichment 0 0 Temperature assumed for analysis 0C 0C Reference k8 (n minal) 0.9374 0.9215 Calculational bias 0.0013 0.0083 Uncertainties 81a s +0.0018 +0.0018 KENO T/A T0.0041 O 8-10 coaceateatioa (8erafiex) Boraflex thickness
+0.0021 70.0047 +0.0021 70.0047 BorafIex width T0.0007 70.0007 Inner box dimension 70.0018 70.0018 Water gap thickness 70.0038 70.0038 SS thickness 70.0025 70.0025 Fuel enrichment 70.0024 70.0024 i Fuel density 70.0026 70.0026 Eccentric assembly position negative Tegative Statistical combination (1) +0.0082 +0.0091 Allowance for burnup uncertainty N/A N/A Total 0.9387 + 0.0082 0.9298 + 0.0091
- Maximum reactivity 0.9469 0.9389 T
4 O IIISquare root of sum of squares. 4-24
4 Table 4.1 Cont'd)
SUMMARY
OF CRITICALITY SAFETY ANALYSES Region 1 Region 2 Boral l i Minimum acceptable burnup- 0 37,370 MWD /MTV 9 4.25 initial enrichment 0 0
. Temperature assumed for analysis 0C 0C Reference kg (nominal) 0.9282 0.8999 -Calculational bias 0.0083 0.0013 Uncertainties Bias +0.0018 +0.0018 KENO I.0041 Y/A B-10 concentration (Boraflex) N/A +0.0028 ]' Boraflex thickness N/A 70.0078 Boraflex width N/A 70.0009 Inner box dimension +0.0018 70.0011 Water gap thickness 70.0038 T/A SS thickness 70.0025 +0.0001 Fuel enrichment 70.0024 70.0024 Fuel density 70.0026 70.0026 Eccentric assembly position negative negative III +0.0073 +0.0093 Stati stical combination Allowance for burnup uncertainty N/A N/A Total 0.9365 + 0.0073 0.9199 + 0.0093 Maximum reactivity 0.9438 0.9292
- o IIISquare root of sum of squares.
4-25 I
i
..() Table 4.2 REACTIVITY EFFECTS OF ABNORMAL AND ACCIDENT CONDITIONS Accident / Abnormal Conditions Reactivity Effect Temperature increase Negative in both regions Void (boiling) Negative in both regions Assently dropped on top of rack Negligible Lateral rack module movement Negli gible r
Misplacement of a fuel assembly Posi tive O l l l O 4-26 l
~.
O Tebie 4 3 FUEL ASSEMBLY DESIGN SPECIFICATIONS Fue1~ Rod Data Outside diameter, in. 0.360 Cladding thickness, in. 0.0225 Cladding material 2frcaloy-4 Pellet diameter, in. 0.3083 UO2 pellet density, % TO 95 1 2 3 U02 stack density, g/cm 10.288 1 0.217 Enrichment, wt1 U-235 4.2 1 0.05 Fuel Assembly Data Numer of fuel rods 264 (17 x 17 array) ] Fuel rod pitch, in. 0.496 Control rod guide tube Nut er 24 Outside diameter, in. 0.474 Thickness, in. 0.016 Material 2frcaloy-4 Instrument thiele 1 Numer Outside diameter, in. 0.474 Thickness, in. 0.016 Material Zircaloy-4 U-23S loading g/ axial em of asse21y 48.6 1 1.0 0 4-27
s O- Tabie 4.4 COMPARIS0N OF COLD, CLEAN REACTIVITIES CALCULATED AT 25,000 MWD /MTU BURNUP AND 3.2% ENRICHMENT kB Xe-free 0 0C Calculational Method Fuel Assertly In Region 2 Cell CASM0-2E 1.1206 0.9061 NULIF/PDQ07 1.1294 0.9017 EPRI-CELL 1.1201(II - f (1)EPRI-CELL kB at maximum value during long-term (30-year) storage. 1 O Tabie 4.5 ESTIMATED UNCERTAINTIES IN REACTI'i!TY 00E TO FUEL DEPLETION EFFECTS i , Design 0.5 x 10-6 Burnup Times Reactivity Loss, Wktl) Initial Enrichment MWD /MTU Burnup, Wk l 1.8 5,230 0.0026 0.0475 2.5 15,720 0.0079 0.1575 3.2 25,000 0.0125 0.2337 j 3.7 31,280 0.0156 0.2757 4.2 37,370 0.0187 0.3107 (1) Total reactivity decrease, calculated for the cold, Xe-free condition in t- the fuel storage rack, from the beginning-of-life to the design burnup, i O l 4-28 l [.
L Tasie 4.6
-O LONG-TERM CHANGES IN REACTIVITY IN STORAGE RACK (XENON-FREE) 4 Wk from Shutdown (Xenon-free) ,
Storage 3.2%E 4.25E Time, years 025,000 MWD /MTV 037,000 MWD /MTU
. 0.5 -0.0046 -0.0057 1.0 -0.0080 -0.0103 10.0 -0.0406 -0.0529 20.0 -0.0588 -0.0756 30.0 -0.0692 -0.0886 Table 4.7 ;
FUEL BURNUP VALUES FOR REQUIRED REACTIVITIES (kB ) WITH FUEL 0F VARIOUS INITI AL ENRICHMENTS Fuel Initial Reference Uncertainty III Design Burnup, Enrichment kg In Burnup, W Limit kB MWD /MTU l - 0.9186 0 0.9186 0 1.58 1.8 0.9186 0.0026 0.9160 5.230 l 1 0.0079 0.9107 15,720 2.5 0.9186 , j 0.0125 0.9061 25,000 ; 3.2 0.9186 0.0156 0.9030 31,280 3.7 0.9186 0.9186 0.0187 0.8999 37,370 r 4.2 l l (1)See Section 4.4.2. l 4-29 f r
. _ - . . . _ . . - - . _ . _ . - _ , _ , _ , _ _ _ . _ _ , . _ _ _ _ , _ . . , _ _ _ _ _ _ _ , _ _ _ _ _ . . , _ ~ . _ , _ , . . , _ , _ . _ . _ _ _ _ _ , , , _ _ . _ . . . _ _ . . , _ . -
O Tabie 4 8 : EFFECT OF TEMPERATURE AND V0!0 ON CALCULATED REACTIVITY OF STORAGE RACK Incremental Reactivity Change, Wk Case Region 1 Region 2 0 Reference 0C Reference 200C -0.0022 -0.0047 50 C -0.0084 -0.0081 800C -0.0165 -0.0121 1200C -0.0298 -0.0178 120 0C + 207, void -0.0953 -0.0520 O O 4-30
\ \ 'J TABLE 4.9 BENCHHARK EXPERIHENTS Soluble General U-235 Separating Boron Description W/0 Reflector Material Iggi Keff
- 1. U02 rod lattice 2.46 Water Water 0 0.9857 +/ .0028
- 2. U02 rod lattice 2.46 Water Water 1037 0.9906 +/ .0018
- 3. UO2 rod lattice 2.46 Water Water 764 0.9895 +/ .0015 4 U02 rod lettice 2.46 Water B4C pins 0 0.9914 +/ .0025
- 5. UO2 rod lattice 2.46 Water B4C pins 0 0.9891 +/ .0026
- 6. U02 rod lattice 2.46 Water B4C pins 0 0.9955 +/ .0020
- 7. U02 rod lattice 2.46 Water B4C pins 0 0.9889 +/ .0026
- 8. UO2 rod lattice 2.46 Water B4C pins 0 0.9983 +/ .0025
- 9. U02 rod lattice 2.46 Water Water 0 0.9931 +/ .0028
- 10. U02 rod lattice 2.46 Water Water 143 0.9928 +/ .0025
- 11. U02 rod lattice 2.46 Water SS 514 0.9967 +/ .0020
- 12. U02 rod lattice 2.46 Water SS 217 0.9943 +/ .0019
- 13. U02 rod lattice 2.46 Water Borated A1 15 0.9892 +/ .0023 14 002 rod lattice 2.46 Water Borated A1 92 0.9884 +/ .0023
- 15. UO2 rod lattice 2.46 Water Borated Al 395 0.9832 +/ .0021
- 16. 002 rod lattice 2.46 Water Borated Al 121 0.9848 +/ .0024 2.46 Water Borated A1 487 0.9895 +/ .0020 gs 17 . UO2 rod lattice 2.46 Water Borated A1 197 0.9885 +/ .0022
(,) 18. U02 rod lattice Borat:d Al 634 0.9921 +/ .0019
- 19. UO2 rod lattice. 2.46 Water
- 20. UO2 rod lattice 2.46 Water Borated Al 320 0.9920 +/ .0020 2.46 Water Borated A1 72 0.9939 +/ .0020
- 21. U02 rod lattice 0.9905 +/ .0020
- 22. U metal cylinders 93.2 Bare Air 0 93.2 Bare Air 0 0.9976 +/ .0020
- 23. U metal cylinders
- 24. U metal cylinders 93.2 Bare Air 0 0.9947 +/ .0025
- 25. U metal cylinders 93.2 Bare Air 0 0.9928 +/ .0019
- 26. U metal cylinders 93.2 Bare Air 0 0.9922 +/ .0026 27 . U metal cylinders 93.2 Bare Air 0 0.9950 +/ .0027
- 28. U metal cylinders 93.2 Bare Plexi glas 0 0.9941 +/ .0030 29 . U metal cylinders 93.2 Paraffin Plexiglas 0 0.9928 +/ .0041
- 30. U metal cylinders 93.2 Bare Plexiglas 0 0.9968 +/ .0018 U metal cylinders 93.2 Paraffin Plexiglas 0 1.0042 +/ .0019
- 31. 0.9963 +/ .0030
- 32. U metal cylinders 93.2 Paraffin Plexiglas 0
! 33. O metal cylinders 93.2 Paraffin Plexiglas 0 0.9919 +/ .0032 4
/~% ,
4-31
. . . . i. . . . . -
35 - 2 - F - AC OEP T Al3LE - E - BURNUP DC' MAIN N 30 / _ o
\ , 25 -
20 - o - z - a: - o - m ) %
\ - * - UNAC CEPTA IBLE - ~ .J BURNtJP DOF AIN m
x
]
g 10 , - f - I e ~
~
l 6 - \ - 1 - 0 3.0 3.5 4.0 4.5 1.5 2.0 2.5 O INITIAL ENRICHMENT,WT% U-235 Fig. 4,1 ACCEPTABLE BURNUP DOMAIN IN REGION 2 OF THE BYRON STATION SPENT FUEL STORAGE RACKS. 4-32
i 8.05" 1 0.032"
, ~
j' BOX I.D.
,!f f '= 90RAFLEX -.i_i- --
7 3/4" i 1/16" (
. ._. j _ _ __ . .
j 0.075" 1 0.007" Thick i 0.060" 1 0.005 0.0238 1 0.0017 g bio / cat Inner ss Box I 00000000000000000 00000000000000000 i
! 8882RRRRR8RRR2888 ! 0o0555555c556000o !
o i l@EEEEl0 D0g;,1 :.*>- OElE@l@@0000000000 0000000 00000000000000000 i-
! 00000000000000000 !
i OOeOOeOOeOOeOOeOO i j 00000000000000000 1.160110.040" water cap " 000e0000000009000 =
! OOOOO900eOOeOOOOO : . 00000000000000000 I 00000000000000000 j i j
l i i o,ngo. Thick 1.26' i 0.040" Water Cap y Cap Channel O Fig . 4.2.1 REGION 1 STORAGE CELL GEOMETRY. 4-33
i l 1 l O ! l 8.85" 2 0.032" - 00X I.D. C BORAFLEX j f
!-= ~ ~"
O.075 00 hick o, g . t 0.005" gg 0.023810.0017 g 810/ cat ] l ' I 00000000000000000 00000000 0000 3 g BORAL SHEET 7 .5',( - 0",4 0.125") j sanganaa 88ag ./! I e flillill.lill l= 00000000000000000 00000000000000000 ! i i OOeOOeOOeOOeOOeOO O.75,2"min. j 00000000000000000 Wate cae
- COOeOOOOOOOOOeOOO = " ! OOOOOeOOeOOeOOOOO .
00000000000000000 l / 1 j I l @ O00000000000000r: i j , l I l O.752" min ' O.000' Thick a W Cap y Cap Chamel l O Fig. 4.2.2 REGION 1 STORAGE CELL GEOMETRY. 4-34
O 8.85" 1 0.s 12" " 90X I.D. l
! !r -BORAL SHEET ~~ ~~ ~"
0.0 5 (no in 1) Thick o, g . t 0.005" 0.025 g BIO /cm 2 j g, y g l ' I 00000000000000000 OO OO OO 00000000 I l l88ll08ll88l88 OO OO 00000000000000000 OO 00 j 00000000000000000 !
! i OOeOOeOOeOOeOOeOO 1.132"mini I
00000000000000000 vatercap ' j 000e000000000e000 1
! 00000900900900000 I= :
00000000000000000 ll j
,00000000000000000 I
j ~ j j i i I.232"min. 0.000' Thick U" # 0*P _1 Cap Channel , O Fig. 4.2.3 REGION 1 STORAGE CELL GEOMETRY. 4-35
O LATTICE SPACING _
! 9.03 .040" !! c 8.850" t 0.032" BOX I.D. rj I
i
'0 060" 2 0 M5" 00000000000000000 00000000000000000 08 00 u zx 8009000 9000 OO OO 00 OO 00 00 80 SO 0800000 0000 000000000000 00000000000000000 O 00*o0*o0*o0*o0800 00000000000000000 00000000000000000 i
00900900900900900 ! 00000000000000000 000e0000000009000 00000900900900000 00000000000000000 l _._i_.00000000000000000 _...___.t
, i
- SCRFLEX j _
l : 7 t/4" t 1/16' i l
' O.04t" t 0.007' Thick j c, 0.01302 0.000* 3 510/en8 f
'O Fig. 4.3 Region 2 storage cell geometry. 4-36
O V a i i eiiiii ie i 6 i iiii e i i iiiiii i iif iiiii
/
O
~
u l ~
~
L.1 / **
) - /
N f -
* / - N - e / - v = / - D ~ t c ! - w N
d W
/ / - O '
en r " Y i k N
" ~ /
l
/l - 3 / # . b m M K / 1 y ! / M / - = u - e I /
- t
% o \
s 8 / ", / Q
' l 3 - // m, -
E e l*l
/ / / / - A y / */ / - 5- u 0 - / e / - O e .# o d h * - A - / / - a M O
l
,I - ~
b U
/ / / /
g M M
~ r s / ! ~
D U
/ - @ c - ,/ ,
O J ]
- e4
/ - a - W -
e 3 a o.4 1 0
- T - u O
C
- O in " . sig 6
m
- j O
- U 'A
~ a ., , , ,,,,,,, . . ..,,,,, ,,,,,,,,,- , ,g o o o T ~5 e o w 0 m. O M. N. ed Q * = e4 ed g a @W
' e4 c Lo (DNI1Vd3eO 10H) 31INIdNI-M 4-37
i i i i i i i i i i i I i i i i 1.50 iiii i i i i i.40 i.30 w _ F
- w z -
t O " a 1.20 z - w _ l " x i.10 i.00 l I , , , , , , i i i i i ii i ' ' ' if i I 4 6 8 10 12 O 2 O waren one esTuse" asseme'Ies- 2"- Fig. 4.5 REACTIVITY EFFECT OF WATER SPACING BETWEEN l fuel assemblies. 4-38
5.0 THERMAL-HYORAULIC CONSIDERATIONS (G u) A primary objective in the design of the high density spent fuel storage racks is to ensure adequate cooling of the fuel assembly cladding. In the following, a brief synopsis of the design basis, the method of analysis, and computed results are given. Similar analysis has been used in previous licensing reports on high density spent fuel racks for Fermi 2 (Docket 50-341), Quad Cities 1 and 2 (Dockets 50-254 and 50-265), Rancho Seco (Docket 50-312), Grand Gulf Unit 1 (Docket 50-416), Oyster Creek (Docket 50-219), Virgil C. Summer (Docket 50-395), and Diablo Canyon 1 and 2 (Dockets 50-275 and 50-323). 5.1 DECAY HEAT CALCULATIONS FOR THE SPENT FUEL This report section covers requirement 111.1.5(2) of the NRC's "0T Position for Review and Acceptance of Spent Fuel Storage and Handling Applications" issued on April 14, 1978. This requirement states that calculations for the amount of thermal energy removed by the spent fuel cooling system shall be made in accordance with Branch Technical Position APCSB 9-2, "Residual Decay The Energy for Light Water Reactors for Long Term Cooling" (Ref. 1). i calculations contained herein have been made in accordance with this requirement. 5-1
-s 5.1.1 Basis C
The Byron Nuclear Power Station Units 1 and 2 reactors are both rated at 3411 megawatts thermal (MWt). Each core contains 193 fuel assemblies. Thus, the average operating power per fuel assembly, Po, is 17.6736 MW. The fuel discharge can be made in one of the following two modes: o Normal refueling discharge o Full core discharge Another mode exists where a normal discharge from one unit could be followed 17 days later by a full core discharge from the other unit. IO
' An equilibrium reload consists of 84 assemblies (with 18-month cycles). The four transitional reloads for each unit consist of 88 assemblies. The fuel transfer begins after 100 hours of decay time in the reactor (time after shutdown). It is assumed that the time period of discharge of this batch is 28 hours (three assemblies transferred to the pool per hour). The cooling system consists of a Seismic Category I spent fuel cooling circuit. The bulk 0
temperature analysis assumes a 105 F coolant inlet temperature to the spent fuel pool heat exchanger for these refueling cases. For the full core discharge, it is assumed that the total time period for the discharge of the full core is 64 hours (after 100 hours of shutdown time in the reactor). The discharge rate to the pool is assumed to be continuous and uniform. l 5-2
O For the 17 day back-to-back discharge mode, the following assumptions were made; 410 hours between the normal discharge and the full core discharge; 3288 hours between the normal discharge and the previous normal discharge; all the previous normal discharge were assumed double discharges (168 fuel assemblies) for conservatism. The fuel assemblies are removed from the reactor after a maximum postulated time of 4.5 years of cumulative operating time. Since the decay heat load is a monotonically increasing function of the cumulative reactor operating time, t o, it is conservatively assumed that every fuel assembly discharged has had the maximum postulatedo t of 4.5 years for the batch size of 84. The water inventory in the reactor cavity cooled by the residual heat removal (RHR) heat exchanger exchanges heat with the fuel pool water mass through the refueling canal. This source of heat removal is neglected in the analysis. Thus, the results obtained for the normal refueling discharge, full core discharge and 17 day discharge mode are conservative. The fuel pool cooling system (FC) consists of two independent trains, each consisting of one pump and heat exchanger. Either train is capable of j providing sufficient cooling for the pool. The following list identifies all relevant design data for the spent fuel pool heat exchangers: o Type Tube and shell g O 5-3 l
O o Quantity 2 o Performance data 6
- Heat transferred 15.833 x 10 8tu/hr Tube Side I
6
- Fluid flow 2.23 x 10 lb/hr - Pool water inlet temperature 1200F - Outlet temperature 112.90F - Fouling factor .0005 O Shell Side 6 - Fluid flow 2.72 x 10 lb/hr - Coolant inlet temperature 1050F - Outlet temperature 110.820F - Fouling factor .0005 The above data enables complete characterization of the thermal performance of a fuel pool heat exchanger.
O 5-4
5.1.2 Model Description [] Reference 1 is utilized to compute the heat dissipation requirements in the pool. The total decay heat consists of fission product and heavy element decay heat. Total decay heat, P, for a fuel assembly is given as a linear function of Po and as an exponential functior, of ot and 's P=Po f (to,13) (5.1-1) where: P = total decay heat per fuel assembly, linear function of P o
^ver^9e operati"9 aower per rue) essemb's O P o =
t = cumulative exposure time of the fuel assembly in the reactor o i = time elapsed since reactor shutdown s The appropriate uncertainty factor, K, was applied in accordance with NUREG-0800 (Ref. 1). Furthermore, the operating power, P o , is taken equal to the rated power, even though the reactor may be operating at less than its rated power during much of the exposure period for the batch of fuel assemblies. Finally, the computations and results reported here are based on the discharge taking place when the inventcry of foal in the pool will be at its maximum resulting in an upper bound on Jted decay heat rate. 5-5
(] Having determined the heat dissipatien rate, the next task is to evaluate the time-dependent temperature of the pool water. Table 5.1 identifies the loading cases examined. The pool bulk temperatura is determined using the first law of thermodynamics (conservation of energy). A number of simplifying assumptions are made which render the analysis conservative, principally: o The heat exchangers are assumed to have maximum fouling. Thus, the temperature effectiveness, S, for the heat exchanger utilized in the analysis is the lowest postulated value: S = .3875 for fuel pool cooler. S is calculated from heat exchanger technical data sheets. No heat loss is assumed to take place through the concrete floor, o No credit is taken for the improvement in the film coefficients of the heat exchanger as the operating temperature rises. Thus, the film coefficient used in the computations are lower bounds, o No credit is taken for heat loss by evaporation of the pool water. o No credit is taken for heat loss to pool walls and pool floor slab. The basic energy conservation relationship for the pool heat exchanger system yields:
" (5.1-2)
C t 01-02 where: C = Tnermal capacity of stored water in the pool t i
= Temperature of pool water at time, t t
O 5-6
rh Heat generation rate due to stored fuel assemblies in the pool; U Qi
=
is a known function of time, : from the preceding section. Qi
= Heat removed in the fuel pool cooler.
Q2 The pool has a total water inventory of approximately 63444,0 cubic feet when all racks are in place in the pool and every storage location is occupied. 5.1.3 Decay Heat Calculation Results The calculations were performed for the pool, disregarding the additional thermal capacity and cooling system available in the transfer canal, and the reactor cavity. IS For a specified coolant inlet temperature and flow rate, the quantity Q 2 shown to be a linear function of 1 in a recent paper by Singh (Ref. 3). As stated earlicr, Qi is an exponential function of .. Thus, Equation 5.1-2 can be integrated to determine t directly as a function of :. The results are The results show that the pool water plotted in Figures 5.1A through 5.20. j i never approaches the boiling point even with the most adverse heat load, under t. normal operating conditions. These figures also give Qi as a function of The first and third plots for each Four plots are generated for each case. case show temperature and power generation, respectively, for a period l l extending from t = 0 to r : 2n where 'n is the tctal time of fuel transfer. The second and fourth plots show the same quantities (i.e., 1 The temperature and power generation, respectively) over a longer period. O l 57 I l l
long-term plots are produced to show the temperature drop with time. ( Sumari:ed results are given in Table 5.2. Finally, computations are made to determine the time interval to boiling after all heat dissipation paths are lost. Computations are made for each case under the following two assumptions: o All cooling systems lost at the instant pool bulk temperature reaches the maximum value, o All cooling systems lost at the instant the heat dissipation power reaches its maximum value in the pool. Results are summarized in Table 5.3. Table 5.3 gives the bulk boiling vaporization rate for all cases at the instant the bofling commences. This rate will decrease with time due to reduced heat generation in the fuel. In m all cases, adequate time exists to take corrective action. gd 5.2 THERMAL-HYDRAULIC ANALYSES FOR SPENT FUEL COOLING This report section covers requirement III.1.5(3) of the NRC's "0T Position for Review and Acceptance of Spent Fuel Storage and Handling Applications," issued on April 14, 1978. Conservative methods have been used to calculate the maximum fuel cladding temperature as required therein. Also, it has been determined that nucleate boiling or voiding of coolant on the surface of the fuel rods occurs only at the locations where freshly discharged fuel assemblies are stored. O 5-8
5.2.1 Basis bm In order to deternines an upper bound on the maximum fuel cladding temperature, a series of conservative assumptions are made. The most important assumptions are listed below: o As state ( Tbove, the fuel pool will contain spent fuel with varying ti..a-after-shutdown (13). Since the heat emission falls off rapidi,/ with increasing i s, it is obviously conservative to assume that all fuel assemblies are fresh (t 3 = 100 hours) and they all have had 4.5 years of operating time in the reactor for cases 1 and 2. The heat emission rate of each fuel assembly is assumed to be equal (Ref. 2). o As shown in Figure 2.1 in Section 2, the modules occupy an irregular floor space in the pool. For the hydrothermal analysis, a circle circumscribing the actual rack floor space is drawn, it is further assumed that the cylinder with this circle as its base is packed with fuel assemblies at the nominal pitch of 9.03 inches (see Figure 5.3). (G The downcomer space around the rack module group varies, as shown s) o in Figure 2.1. The nominal downcomer gap available in the pool is assumed to be the minimum gap available around the idealized cylindrical rack; thus, the maximum resistance to downward flow is incorporated into the analysis. o No downcomer flow is assumed to exist between the rack modules. 5.2.2 Model Description In this manner, a conservative idealized model for the rack assemblage is obtained. The water flow is axisymmetric about the vertical axis of the circular rack assemblage, and thus, the flow is two-dimensional (axisymetric three-dimensional). Figure 5.4 shows a typical "flow chimney" rendering of 1 the thermal hydraulics model. The governing equation to characterize the flow l field in the pool can now be written. The resulting integral equation can be l solved for the lower plenum velocity field (in the radial direction) and axial l It velocity (in-cell velocity field), by using the method of collocation. l 5-9 l ___ _ --~-v -v ,--c
should be added that the hydrodynamic loss coefficients which enter into the (] v formulation of the integral equation are also taken from well-recognized sources (Ref. 4) and wherever discrepancies in reported values exist, the conservative values are consistently used. Reference 5 gives the details of mathematical analysis used in this solution process. After the axial velocity field is evaluated, it is a straight-forward matter to compute the fuel assembly cladding temperature. The knowledge of the overall flow field enables pinpointing the storage location with the minimum axial flow (i.e., maximum water outlet temperature). This is called the most "choked" location. In order to find an upper bound on the temperature in a typical cell, it is assumed that it is located at the most choked location. Knowing the global plenum velocity field, the revised axial flow through this choked cell can be calculated by solving the Bernoulli's equation for the flow circuit through this cell. Thus, an absolute upper bound on the water exit temperature and maximum fuel cladding temperature is obtained. It is believed that, in view of the aforementioned assumptions, the temperatures calculated in this manner overestimate the temperature rise that will actually occur in the pool. The maximum pool bulk temperature, t, is computed in Section 5.1.3 and reported in Table 5.2. The corresponding average power output from the hottest fuel assembly, q, is also reported in that table. The maximum radial Thus, it is peaking factor, Fxy, is 1.55 for the Byron Nuclear Power Station. conservative to assume that the maximum specific power of a fuel assembly, qA' is given by: 5-10
aaa e xx (5 2-1) .O where: F xy = 1.55 The maximum temperature rise of pool water in the most disadvantageously placed fuel assembly is given in Table 5.4 for all loading cases. Having determined the maximum local water temperature in the pool, it is now possible It is conservatively to determine the maximum fuel cladding temperature. assumed that the total peaking f actor, F g , is 2.32. Thus, a fuel rod can The produce 2.32 times the average heat emission rate over a small length. axial heat dissipation in a rod is known to reach a maximum in the central region, and taper off at its two extremities. For the sake of added conservatism it is assumed that the peak heat emission occurs at the top where the local water temperature also reaches its maximum. Furthermore, no credit is taken for axial conduction of heat along the rod. The highly conservative model thus constructed leads to simple algebraic equations which directly give the maximum local cladding temperature, tc* 5.2.3 Results Table 5.4 gives the maximum local cladding temperature,c t , at the instant the pool bulk temperature has attained its maximum value. It is quite possible, however, that the peak cladding temperature occurs at the instant of maximum value of qA, the instant when the fuel assembly is first placed in a stJrage location. Table 5.5 gives the maximum local cladding temperature at t = 0. Q 5-11
0 The local boiling temperature near the top of the fuel cladding is 240 F. However, the cladding temperature must be somewhat higher than the boiling temperature to initiate and sustain nucleate boiling. The above considerations indicate that a comfortable margin against the initiation of localized boiling exists in case 1. For full core discharge (case 2 and case
- 3) under the described assumptions, the maximum cladding temperature will give rise to localized nucleate boiling, but not to bulk pool boiling (5.4).
O O 5-12
REFERENCES TO SECTION 5
- 1. NUREG-0800 U.S. Nuclear Regulatory Commission, Standard Review Plan, Branch Technical Position ASB 9-2, Rev. 2, July 1981,
- 2. FSAR, Byron Nuclear Power Station.
- 3. Journal of Heat Transfer, fundamental Transactions of the ASME, August 1981, Vol.103, "Some Fundamental Relationships for Tubular Heat Exchanger Thermal Performance," K.P. Singh.
- 4. General Electric Corporation, R&D Data Books, "Heat Transfer and Fluid Flow," 1974 and updates.
- 5. 4th National Congress of the ASME, "A Method for Computing the Maximum Water Temperature in a Fuel Pool Containing Spent Nuclear Fuel," by K.P. Singh and A.I. Soler, Paper 83-NE-7, Portland, Oregon (June 1983).
1 - O 5 13 l
o o o i 4 l Table 5.1. I.IST OF CASES ANAI.YZED f No. of Total Time- - Fuel to Transfer Decay Time * } Assemblies Fuel into Before Transfer Discharged, the Pool Case No. Condition l h, hrs Begins, hrs N i
- j g Normal refueling discharge 84 28 100 4
- 2 Full core 193 64 100 discharge
. m 9 l 5 3 Back-to-back refueling 277 + 100 1 1 } t i 4 i t
- i Discharge is assumed to be into a pool containing fuel from 17 previous discharges of 168 assembiies.
4
- +
f 28 hours for first discharge (84 assemblies), 64 hours for second discharge d (full core). 17 days between the discharges and 4-1/2 months to the previous discharge. i t 4 4
l O O o Table 5.2 MAX IHilM 14X)l. lilli.K TEMI'ERATilRE, t , COI NCllll3Xi' 'IUI'Al. l'OWER, QI , AND ColNCIDi'XT SPECIFIC POWER, q, FOR Tile IKTITEST ASSEMill.Y Coincident Coincident Ilme to Maximum Ilme (After Coincident Iotal Iransfer Pool flulk Initiation Spectric Power Case 11 0 . of fuel Into Temp., t, of fuel Power, q, Q (10 6) llo . Assemblies Pool, hrs *f Iransfer), hrs titu/sec fltu / hou r flo t e s 28 138 37.0 55.15 .1985 Normal 1 84 , discliange w I 64 155 71.0 50.30 .1811 Fuii cone, 5 2 193 discliange 158 55.0 55.0 .I980 Abnormal 3 277 + discharge Discharge is assumed to be into a pool containing fuel from 17 previous discharges of 168 assemblies.
+ 28 hours for first discharge; 64 hours for the second discharge; 17 days between first and second discharge and 4-l/2 months to tiie previous discharge.
O Table 5.3 VAPORIZATION RATE FROM THE INSTANT ALL COOLING IS LOST : CONDITION 1 CONDITION 2 Loss of Cooling at Maximum loss of Cooling at ' Pool Bulk Temperature Maximum Power Discharge Case Time (Hrs) Vap. Rate Time (Hrs) Vap. Rate No. To (ib/hr) To (ib/hr) Boil Boil 1 9 35329.0 9 35788.0 2 4 54233.0 4 54827.0 3 3.8 57560.0 3.8 58395.0 0 5-16
Y Table 5.4 MAXIMUM LOCAL POOL WATER TEMPERATURE AND LOCAL FUEL CLA00!NG TEMPERATURE AT INSTANT OF MAXIMUM POOL BULX TEMPERATURE 3 Maximum local Maximum Local Case Water Fuel Cladding No. Temperature, OF Temperature, 0F Cdse Identified 1 194 239 84 assemblies 2 208 250 193 assemblies 3 227.4 275.1 277 assemblies O {. I l l . O I 5-17
10 Table 5.5 POOL AND MAXIMUM CLADDING TEMPERATURE AT THE INSTANT FUEL ASSEMBLY TRANSFER BEGINS Coincident Pool Cladding Temperature. F , Case No. Temperature, OF Bulk Local 1 236.6 122.6 186.0 2 236.6 122.8 186.2 3 236.6 122.8 186.2 O , h P i f i t I 1 I r l t O . 5-18 j
. . _ . _ , . . - . - - . ~ _ . . . _ . _ _ _ _ - - _ _ _ _ _ _ _ - . _ _ _ - - _ . . _ _ _ _ _ _ . . _ _ ~
m N E 3
, O I
v o o O
.8 O - 8 i
O o w W o E o 3 o H K W Q. 2 8 W
~
e r o nm N o .J U D E
.8 m < = x .; O O M O o o a. .o
- O 1 - .J 2 i
M W -
& D .a W
I 90 6 (
.8 w 3 1
(n s oI p 4 W r Z. -g W c Q. . M J 6 <
.o .
g 3 E 4 O 9 ", Z u *
,g l
O U, o l - I
.8n l .8~ .8 =
(d'Obbf3801Ib$dH31 l ' 5-19
d O
*CC a w m
g
)_( ~,,
o s
-s u k I
O n I O
. m z
C H nc T. w T E m w O H m N O
'w A 3 F >
C z
# O , O -
r o
~@ O )V o O T z =t C
0O- -- m M C C "I A Z m % M D s 2 y H
)n& z F O ~
m > O
- k.
O m s,. z
^
O e" O' O
-p O
O l? O O b V O O 0 421 0 021 0 611, 0 041 0 631 0 231 0 621 0.e41 j )F . GED( ERUT AREPMET I 9-20 l n ,- - a
*A '
._.. .-- - . .-----_---..-..~.;
O d . i _ LO' 1 E j i O t O. ,
= .,
t O
*I m.- t v. \
O W i m 1 i 1 o w*m N' - T E"~ ro I s g < a e, m
- v.
EN 10 ; x o G. O. .I O. r W l O. m_ O
. . . . I; 0.000 100.0 200.0 300.0 ^400.0 500.0 600.0 706:0. 800.0 ~ TIME (HRS) :
FIO. 5.1 C POWER DISCHARoED IN SPENT FUEL POOL -(0-700 HOURS) L o.onmat nEFUEtina.oSenan E>-
O O O i i 1 l l c l 1 to _ i O W l - O T N p I N o e-m
*7_
j ccru , tu x o a_ o Po-i 1 4 o S
!- ! .4 E.
W o & o 2 a
.o .W -W m y - G - . E M 4 o J E*
e 3 g o e
.o n = - a a vi o a- o w E
I &
, o- o- -g w .a oi E- W-M 3 a; 'O-- H E .a o 3 >- E ,
o g w a l _o W n
& j @ r 1
o 4 o N ! d ; O
.g
. o 6 o
.o
(') ? t.. [
-i o ;
o
.o N
> b t o , r o l ' .oe i i i e 2 t O , V
=
o . . O'Gdi D'Odi O'GYi O'OYi O'GET O ' 0 E'i O'Gdi 0'021*
- (d ,930) 38niveadH31 ,
o 1. 4 t
- 5-23 t b
-- . . - - - - . . . - - - _ _ . - _ - - . . _ . - . ~ . _ _ . - _ _ _ - . _ _ _ _ _ _ _ _ - _ _ _ _ _ _ - - . . _ _ ~ - . . .,.,,,,,,-I
"g4mAenwe.as A g 4 m AU.6-- .s, M .
o
.N A e .: < g N '
1
.D. '
o 14 > i%)~ ' O< 6
-: e -
P t o C-w- o
.o e '
c W E o 3 l
.o - e m_ < -E W
3 . - ,
- o. W
. .o > w cc @
E i
; M .J . -g 3 -0 g a , , .o -
I b Q m- - 0 W
>- 0 E
- o. c" g 0
O
"@ W W. {
Z
.O' .s H
3 g .J
- 3 i p th -
- I
.o' n
W a, r O E
- o. N.
.o 9 O
o. o
.n I o.
3N o -
.P .o L . e i O = !
i O. 1 0*Gdi 0*0dl O'Gil O'0/1 0'Gdi O'Odi 0'Edi 0
- 02f i (3'930)3801v83dH31 ,
5-24
L o
^ .S - us W
D o O
^
o I (Jx, R.- o _- o N e o I o w o o _o o .J o 0
-8 O A
J o W o 3
.o a m
H a
.2 W e W O . . A g e m <
T Z 2 0 tn e o C O O 5 W W
.S c u O t
z E E (~3 - < O (,,) H I O O M
.J e O J
o 3 C a i g w W o N o O l .o A l O cv. o e o *
.o O m -
Y-o o
.o N
o o
.o s -
o o
%)' o l 0*dG O'dG O'dP 0'dr O'dE 0'dE 0'd2 0'08 ,0i k (bH/019) d3MCd 1
5-25
~
l o
.m" a
M o O ('d3 o
- c " e-9" l
o O w a a o O O
-g o.
J W D o u.
.o
- H Z ^
. W W c.
(0 0 o E
.o Z N
g y) O
> O 3 4 W Q OS O n 'S m$ <
Z O
+ 0 0 40
\ o O a I J
.o 3 n g W b c
O o a.
.o l
- O N.
o d
.o 2 l
o
.o N
o
.oe e
b)
\>
S o O ' G'G 0-dg c g', g . 0'" - U ' 0# 0I * (bH/019f'd3M0d 5-26 L . .
ACTUAL OUTLINE OF POOL
^ 7"^t OuTelu.
O OF RACK ASSEMBLY
. EL S lES \_ %
RACK ASSEMBLY 4 ,
+++
O +++-- (
~ 'A rDEALIZED OUTLINE IDEALIZED OUTLINE * = ^88 man lO OF eOOL BOUNDARY FIGURE 5.5 IDEALIZATION OF RACK ASSEMBLY 5-27 _ _ _ _ __ __
O r - d> , ( h 7 A 6
< < r, / /
7 _[I C.7
. out 0
T ' in N' g E P1 a A e H E o 3 h <- O Heat Addition l 3 O l c' c V 3 {, ' O ' O __ _/ [L V _1
't t._/ [ T in "1
L7 f~~7 / / l
, ;.-,y
<O FIGURE 5.4 Thermal Chimney Flow Model 5-28
/7 6.0 STRUCTURAL ANALYSIS V
The purpose of this section is to demonstrate the structural adequacy of the spent fuel rack design under normal and accident loading conditions. The method of analysis presented herein uses a time-history integration method similar to that previously used in the Licensing Reports on High Density Fuel Racks for Fermi 2 (Docket No. 50-341), Quad Cities 1 and 2 (Docket Nos. 50-254 and 50-265), Rancho Seco (Docket No. 50-312), Grand Gulf Unit 1 (Docket No. 50-416), Oyster Creek (Docket No. 50-219), V.C. Summer (Docket No. 50-395), and Diablo Canyon 1 and 2 (Docket Nos. 50-275 and 50-323). The results show that the high density spent fuel racks are structurally adequate to resist the postulated stress combinations associated with level A, 8, C, and D conditions as defined in References 1 and 2. 6.1 ANALYSIS OUTLINE The spent fuel storage racks.are Seismic Category I equipment. Thus, they are required to remain functional during and after a Safe Shutdown Earthquake V (Ref. 3). As noted previously, these racks are neither anchored to the pool floor nor attached to the sidewalls. The individual rack modules are not interconnected. Furthermore, a particular rack may be completely loaded with fuel assemblies (which corresponds to greatest rack inertia), or it may be completely empty. The coefficient of friction, u, between the supports and l pool floor is another indeterminate factor. According to Rabinowicz (Ref. 4) the results of 199 tests performed on austenitic stainless steel plates submerged in water show a mean value of u to be 0.503 with a standard deviation of 0.125. The upper and lower bounds (based on twice the standard deviation) are thJs 0.753 and 0.253, respectively. Two separate analyses are performed for the rack assemblies with values of the coefficient of friction equal to 0.2 (lower limit) and 0.8 (upper limit), respectively. Analyses performed for the geometrically limiting rack modules focus on limiting values of the coefficient of friction, and the number of fuel assemblies stored. Cases studied are for the 12x14 rack, for the 7x14 rack, and for a special rack with 35 regular locations and 12 large storage cells. Typical simulations are: O 6-1
p 0 Fully loaded rack (all storage locations occupied), V u = 0.8; 0.2 (u = coef ficient of f riction) 0 Rack half full, u = 0.8, 0.2 The simulations were performed using a consolidated fuel weight of 3000f per cell. The case of nearly empty racks is found not to be critical. The method of analysis employed is the time-history method. The pool slab acceleration data were developed from the resporse spectra provided by Sargent and Lundy Engineers, Chicago, Illinois. The objective of the seismic analysis is to determine the structural response (stresses, deformation, rigid body motion, etc.) due to simultaneous application of the three statistically independent, orthogonal excitations. Thus, recourse to approximate statistical summatior, techniques such as the "Square-Root-of-the-Sum-of-the-Squares" method (Ref. 5) is avoided. For nonlinear analysis, the only practical method is simultaneous application, m Pool slab acceleration data are provided for two earthquakes: Operating Basis U Earthquake (0BE) and Safe Shutdown Earthquake (SSE). Figures 6.1 - 6.3 show the time-histories corresponding to the SSE condition. The seismic analysis is performed in three steps, namely:
- 1. Development of a nonlinear dynamic model consisting of inertial mass elements and gap and friction elements.
- 2. Generation of the equations of motion and inertial coupling and solution of the equations using the "component element time integration scheme" (References 6 and 7) to determine nodal forces and displacements.
- 3. Computation of the detailed stress field in the rack (at the critical location) and in the support legs using the nodal forces calculated in the previous step. These stresses are checked against the design limits given in Section 6.5.
A brief description of the dynamic model follows. O a 6-2
6.2 FUEL RACK - FUEL ASSEMBLY MODEL - (v3-Since the racks are not anchored to the pool slab or attached to the pool walls or to each other, they can execute a wide variety of rigid body motions. For example, the rack may slide on the pool floor (so-called "sliding condition"); one or more legs may momentarily lose contact with the liner ("tipping condition"); or the rack may experience a combination of sliding and tipping conditions. The structural model should permit simulation of these kinematic events with inherent built-in conservatisms. Since the Byron racks are equipped with girdle bars to dissipate energy due to inter-rack impact (if it occurs), it is also necessary to model the inter-rack impact phenomena in a conservative manner. Similarly, lift off of the support legs and subsequent impacts must be modelled using appropriate impact elements, and Coulomb friction between the rack and the pool liner must be simulated by appropriate piecewise linear springs. These special attributes of the rack dynamics require a strong emphasis on the modeling of the linear and nonlinear springs, dampers, and stop elements. The model outline in the remainder of this section, and the model description in the following section describe the detailed modeling technique to simulate these effects, with emphasis placed on the nonlinearity of the rack seismic response. 6.2.1 Outline of Model
- a. The fuel rack structure is a folded metal plate assemblage welded to a baseplate and supported on four legs. The rack structure itself is a very rigid structure. Dynamic analysis of typical multicell racks has thown that the motion of the structure is captured almost completely by the behavior of a six degrees-of-freedom structure; therefore, the movement of the rack cross-section at any height is described in terms of the six degrees-of-freedom of the rack base. The rattling fuel is modelled by five lumped masses located at H, .75H, .5H, .25H and at the rack base. Fuel mass connectivity is included to model the fuel assembly stiffness.
- b. The seismic motion of a fuel rack is characterized by random rattling of storage locations. Assuming that fuel assemblies in their individual all assemblies vibrate in phase obviously exaggerates the computed dynamic loading on the rack structure. This assumption, however, greatly reduces the required degrees-of-freedom needed to model the fuel assemblies which are represented by five lumped masses located at different levels of the rack. The centroid of each fuel assembly mass can be located, relative to the rack structure centroid at that level, so as to simulate a partially loaded rack.
O 6-3
- c. The local flexibility of the rack-support interface is modeled
(] U conservatively in the analysis,
- d. The rack base support may slide or lift off the pool floor,
- e. The pool flocr has a specified time-history of seismic accelerations along the three orthogonal directions.
- f. Fluid _ coupling between rack and assemblies, and between rack and adjacent racks, is simulated by introducing appropriate inertial coupling into the system kinetic energy. Inclusion of these effects uses the methods of References 4 and 6 for rack / assembly coupling and for rack / rack coupling (see Section 6.2.3 of this report).
- g. Potential impacts between rack and assemblies are accounted for by appropriate "compression only" gap elements between masses involved.
- h. Fluid damping between rack and assemblies, and between rack and adjacent rack, is conservatively neglected,
- i. The supports are modeled as "compression only" elements for the vertical direction and as "rigid links" for dynamic analysis. The bottom of a support leg is attached to a frictional spring as described in Section 6.3. The cross-section inertial properties of the support legs are computed and used in the final computations to determine support leg stresses.
O V j. The effect of sloshing can be shown to be negligible at the bottom of a pool and is hence neglected,
- k. Inter-rack impact, if it occurs, is simulated by a series of gap elements at the top and bottom of the rack in the two horizontal directions. The most conservative case of adjacent rack movement is assumed; each adjacent rack is assumed to move completely out of phase with the rack being analyzed. The effect of misalignment is also considered. The rack is assumed misaligned so that the rack has a gap that vaires linarly from 0 inch at one corner to 0.25 inch at the other corner along each side.
Two different configurations are consider:
- 1) All four adjacent racks are rotated in the same direction (configuration 1).
i
- 2) Three of four adjacent racks are rotated in the same direction and the fourth one is rotated in the opposite direction (configuration 2).
- 1. The form drag opposing the motion of the fuel assemblies in the storage locations is conservatively neglected in the results reported herein.
- m. The form drag opposing the motion of the fuel rack in the water is also conservatively neglected in the results reported herein.
O 64
- n. The rattling of the fuel assemblies inside the storage locations causes the "gap" between the fuel assemblies and the cell wall to change from a maximum of twice the nominal gap to a theoretical zero gap. However, the fluid coupling coefficients (Ref. 8) utilized are based on linear vibration theory (Ref. 9). Studies in the literature show that inclusion of the nonlinear effect (viz. vibration amplitude of the same order of magnitude as the gap) drastically lowers the equipment response (Ref.
10). Figure 6.4 shows a schematic of the model. Six degrees-of- freedom are used to track the motion of the rack structure. Figures 6.5 and 6.6, respectively, show the inter-rack impact springs ard fuel assembly / storage cell impact springs at a particular level. As shown in Figure 6.4, the model for simulating fuel assembly motion incorporates five lumped masses. The five rattling masses are located at the baseplate, at quarter height, at half height, at three quarter height, and at the top of the rack. Two degrees-of-freedom are used to track the motion of each rattling mass in the horizontal plane. The vertical motion of each rattling mass is assumed to be the same as the rack base. 6.2.2 Model Description The absolute degrees-of-freedom associated with each of the mass locations are identified in Figure 6.4 and Table 6.1. The rattling masses (nodes 1*, 2*, 3 , 4 , 5 ) are described by translational degrees-of-freedom q7-q16' Uj(t) is the pool floor slab displacement seismic time-history. Thus, there are sixteen degrees-of-freedom in the system. Not shown in Fig. 6.4 are the gap elements used to model the support legs and the impacts with adjacent racks. 6.2.3 Fluid Coupling An effect of some significance requiring careful modeling is the so-called "fluid coupling effect". If one body of mass (mt) vibrates adjacent to another body (mass m2), and both bodies are submerged in a frictionless fluid medium, then Newton's equations of motion for the two bodies have the form: (mi + M11) X1-M12 Xy = applied forces on mass mt 6-5
() '
-M 21 Xi + (m2 + M22) 2 = applied forces on mass m2 Xt , 5 2denote absolute accelerations of mass mg and m2, respectively.
M11, M12' M21, and M22 are fluid coupling coefficients which depend on the shape of the two bodies, their relative disposition, etc. Fritz (Ref. 9) gives data for M j for various body shapes and arrangements, it is to be noted that the above equation indicates that the effect of the fluid is to add a certain amount of mass to the body (Mit to body 1), and an external force which is proportional to the acceleration of the adjacent body (mass m2). Thus, the acceleration of one body affects the force field or another. This force is a strong function of the interbody gap, reaching large values for very small gaps. This inertial coupling is called fluid coupling. It has an important effect in rack dynamics. The lateral motion of a fuel assembly inside the storage location will encounter this effect. So will the motion of a rack adjacent to another rack. These effects are included in the equations of motion. For example, the fluid coupling is between nodes 2 and 2* in Figure 6.4. Furthermore, the rack equations contain coupling terms which model the effect of fluid in the gaps between adjacent racks. The coupling terns modeling the effects of fluid flowing between adjacent racks are computed 0 assuming that all adjacent racks are vibrating 180 out of phase from the rack being analyzed. Therefore, only one rack is considered surrounded by a hydrodynamic mass computed as if there were a plane of symmetry located in the middle of the gap region. Finaliy, fluid virtual mass is included in the vertical direction vibration equations of the rack; virtual inertia is also added to the governing equation corresponding to the rotationa? degree-of-freedom, g6(t). 6.2.4 Damping In reality, damping of the rack motion arises from material hysteresis (material damping), relative intercomponent motion in structures (structural damping), and fluid drag effects (fluid damping). In the analysis, a maximum of 4% structural damping is imposed on elements of the rack structure during SSE seismic simulations. This is in accordance with the FSAR and NRC guidelines (Ref. 11). Material and fluid damping are conservatively 6-6
O neglected. The dynamic model has the provision to incorporate fluid damping effects; however, no fluid damping has been used for this analysis. 6.2.5 Impact Any fuel aesembly node (e.g. 2*) may impact the corresponding structural mass node 2. To simulate this impact, four compression-only gap elements around each rattling fuel assembly node are provided (see Figure 6.6). As noted previously, fluid dampers may also be provided in parallel with the springs. The compressive loads developed in these springs provide the necessary data to evaluate the integrity of the cell wall structure and stored array during the seismic event. Figure 6.5 shows the location of the impact springs used to simulate any potential for inter-rack impacts. Section 6.4.2 gives more details on these additional impact springs. 6.3 ASSEM8LY OF THE DYNAMIC MODEL The cartesian coordinate system associated with the rack has the following O Q nomenclature: 0 x = Horizontal coordinate along the short direction of rack rectangular platform 0 y = Horizontal coordinate along the long direction of the rack rectangular platform 0 z = Vertically upward As described in the preceding section, the rack, along with the base, supports, and stored fuel assemblies, is modeled for the general three-dimensional (3-D) motion simulation by a sixteen degree-of-freedom model. To simulate the impact and sliding phenomena expected, 64 nonlinear gap elements and 16 nonlinear friction elements are used. Gap and friction elements, with their connectivity and purpose, are presented in Table 6.2. O 6-7
A descriptive model restricted to two dimensions (one horizontal motion plus vertical motion) of the simulated structure which includes gap and friction elements is shown in Figure 6.7. Note that only the too rattling mass is shown for clarity. The impacts between fuel assemblies and rack show up in the gap elements, having local stiffness Kg , in Figure 6.7. In Table 6.2, gap elements 5 through 8 are for the vibrating mass at the top of the rack. The support leg spring rates K6 are modeled by elements 1 through 4 in Table 6.2. Note that the local compliance of the concrete floor is included in K6 To simulate sliding potential, friction elements 2 plus 8 and 4 plus 6 (Table 6.2) are shown in Figure 6.7. The friction of the support / liner interface is modeled by a piecewise linear spring with a suitably large stiffness Kf up to the limiting lateral lead, uN, where N is the current compression load at the interface between support and liner. At every time step during the transient analysis, the current value of N (either zero for liftoff condition, or a compressive finite value) is computed. Finally, the support rotational friction springs KR reflect any rotational restraint that may be offered by the foundation. This spring rate is calculated using a modified Beusinesq O equation (Ref. 4) and is included to simulate the resistive moment of the support to counteract rotation of the rack leg in a vertical plane. This rotation spring is also nonlinear, with a zero spring constant value assigned af ter a certain limiting condition of slab moment loading is reached. The nonlinearity of these springs (friction elements 9, 11, 13, and 15 in Table 6.2) reflects the edging limitation imposed on the base of the rack support legs. In this analysis, this effect is neglected; any support leg bending, induced by liner / baseplate friction forces, is resisted by the leg acting as a beam cantilevered from the rack baseplate. The spring rate K6 modeling the effective compression stiffness of the structure in the vicinity of the support, is computed from the equition: 1 1 1 1 Y"T*T*T 5 1 2 3 6-8
where: K spring rate of the support leg treated as a i = tension-compression member = E SUPPORT *A SUPPORT /h (h = icngth of support leg) 2 K 2 = 1.05E c B/(1- v ) = local spring rate of pool slab (E = Young's c modulus of concrete, and B = length of bearing surface) K 3 = spring rate of folded plate cell structure above suppe-t leg (same form as K2 with E chosen to reflect the local stiffne,s of the honeycomb structure above the leg) For the 3-D simulation, all support elements (listed in Table 6.2) are included in the model. Coupling between the two horizontal seismic motions is provided both by the offset of the fuel assembly group centroid which causes the rotation of the entire rack and by the possibility of liftoff of one or more support legs. The potential exists for the rack to be supported on one or more support legs or to lif toff completely during any instant of a complex 3-D seismic event. All of these potential events may be simulated during a 3-0 motion. O 6.4 11st istcGRA11os or 1st tQuArious or Mo110s 6.4.1 Time-History Analysis Using Multi-Degree of Freedom Rack Model Having assembled the structural model, the dynamic equations of motion corresponding to each degree-of-freedom are written by using lagrange's Fonnulation. The system kinetic energy can be constructed including contributions from the solid structures and from the trapped and surrounding fluid. A single rack is modelled in detail. The system of equations can be l represented in matrix notation as: I (M] ('q'} = (Q) + (G) t where the vector (Q} is a function of nodal displacements and velocities, and (G) depends on the coupling inertia and the ground acceleration, j Premultiplying the above equations by [H]-I renders the resulting equation l uncoupled in mass. I We have: {'q') = (M]~I (Q) + (M]-I (G} l 6-9 L
(] v As noted earlier, in the numerical simulations run to verify structural integrity during a seismic event, all elements of the fuel assemblies are assumed to move in phase. This will provide maxir,um impact force level, and induce additional conservatism in the time-histcry analysis. This equation set is mass uncoupled, displacement coupled, and is ideally suited for numerical solution using a central difference scheme. The computer program "DYNARACK"* is utilized for this purpose. Stresses in various portions of the structure are computed from known element forces at each instant of time. Dynamic analysis of typical multicell racks has shown that the motion of the structure is captured almost completely by the behavior of a six-degree-of-freedom structure; therefore, in this analysis model, the movement of the rack cross-section at any height is described in terms of the rack base degrees-of-freedom (qi(t)....q6(t). The remaining degrees-of-freedom are associated with horizontal movements of the fuel assembly O masses. In this dynamic model, five rattling masses are used to represent fuel assembly movement in the horizontal plane. Therefore, the final dynamic model consists of six degrees-of-freedom for the rack plus ten additional mass degrees-of-freedom for the five rattling masses. The totality of fuel mass is included in the simulation and is distributed among the five rattling masses.
*This code has been previously utilized in licensing of similar racks for Fermi 2 (Docket No. 50-341), Quad Cities 1 and 2 (Docket Nos. 50-254 and 265),
Rancho Seco (Docket No. 50-312), Oyster Creek (Docket No. 50-219), V.C. Summer (Docket No. 50-395), and Diablo Canyon 1 and 2 (Docket Nos. 50-275 and O 50-323). 6-10
6.4.2 Evaluation of Potential for Inter-Rack Impact O Since the racks are closely spaced, the simulation includes impact springs to model the potential for inter-rack impact, especially for low values of the friction coefficient between the support and the pool liner. To account for this potential, yet still retain the simplicity of simulating only a single rack, gap elements were located on the rack at the top and at the baseplate level. Figure 6.5 shows the location of these gap elements. Loads in these elements, computed during the dynamic analysis, are used to assess rack integrity if inter-rack irrpact occur s. 6.5 STRUCTURAL ACCEPTANCE CRITERIA There are two sets of criteria to be satisfied by the rack modules:
- a. Kinematic Criterion This criterion seeks to ensure that the rack is a physically stable structure. Byron racks are designed to sustain certain inter-rack impact at designated locations in the rack modules. Therefore, physical stability of the rack is considered along with the
'oce142ee aer emeet decor et4om or O iocei4 zed 4ater-rec' 4=nects-the module is cermissible, so long as the subcriticality of the stored fuel array is not violated.
- b. Stress Limits The stress limits of the ASME Code, Section Ill, Subsection NF, 1983 Edition a?e used since this code provides the most appropriate and 4
consistent set of limits for various stress types and various loading conditions. The following loading combinaticns are applicable (Ref. 1). O l 6-11
. Loading Combination Stress Limit 0+L Level A service limits 0+L+T o 0+L+T e +E 0+L+T a +E Level 8 service limits 0+L+To+Pf 0+L+Ta + E' Level 0 service limits 0+L+F The functional capability d
of the fuel racks should be demonstrated, where: 0 = Dead weight inauced stressas (including fuel assembly weight) L = Live Load (0 for the structure, since there are no moving objects in the rack load path). f = Force caused by the accidental drop of the heaviest load d from the maximum possible height. (] P = Upward force on the racks caused by postulated stuck fuel f assembly E = Operating Basis Earthquake E' = Safe Shutdown Earthquake T = Otfferential temperature induced loads (normal or upset o condition) T, = Differential temperature induced loads (abnormal design conditions) The conditions Ta and Tocause local thermal stresses to be produced. The worst situation will be obtained when an isolated storage location has a fuel assembly which is generating heat at the maximum costulated rate. The surrounding storage locations are assumed to contain no fuel. The heated water makes unobstructed contact with the inside of the storage walls, thereby producing the maximum possible temperature difference between the adjacent cells. The secondary stresses thus produced are limited to the body of the rack; that is, the support legs do not experience the secondary (thermal) stresses. 6-12
/%
U 6.6 MATERIAL PROPERTIES The data on the physical properties of the rack and support materials, obtained from the ASME Boiler & Pressure Vessel Code, Section III, appendices, and supplier's catalog, are listed in Tables 6.3 and 6.4. Since the maximum 0 pool bulk temperature (except for the full core discharge case) is 150 , this is used as the reference design temperature for evaluation of material properties. 6.7 STRESS LIMITS FOR VARIOUS CONDITIONS The following stress limits are derived from the guidelines of the ASME Code, Section III, Subsection NF, in conjunction with the material properties data of the preceding section. 6.7.1 Normal and Upset Conditions (Level A or Level B)
' a.- Allowable stress in tension on a net section =Ft = 0.6 Sy or Ft = (0.6) (23,150) = 13,890 psi (rack material)
Ft = is equivalent to primary membrane stresses Ft = (.6) (27,500) = 16,500 psi (upper part support of feet)
= (.6) (62,400) = 37,440 psi (lower part of support feet)
I i
- b. On the gross section, allowable stress in shear is:
F* = .4 S (.4)#(23,150) = 9,260 psi (main rack body) Fy = (.4) (27,500) = 11,000 psi (upper part of support feet)
= (.4) (62,400) = 24,960 psi (lower part of
! support feet) ! O 6-13
n c. Allowable stress in compression, F a U 5 3 kl F, = S 1 (k1 )3 y ll - \ (kl rC )2 c-I c where: C = ( p, 2E )0.5 e 3y k1/r for the main rack body is based on the full height and cross section of the honeycomb region. Substituting numbers, we obtain, for both support leg and honeycomb region: f a= 13,890 psi (main rack body) F, = 16,500 psi (upper part of support feet)
= 37,440 psi (lower part of support feet) d._ Maximum allowable bending stress at the outermost fiber due to flexure about one plane of symmetry: = 13,890 psi (rack body)
Fb = 0.60 16,500S # psi (upper part of support feet F b=37,440
= psi (lowerpartofsupportfeet} - e. Combined flexure and compression:
I a mx bx , C,j fby , y { 0F x bx y by where: f, = Direct compressive stress in the section f = Maximum flexural stress about x-axis
= Maximum flexural stress about y-axis (by C mx = C my = 0.85 f
a Dx = 1 F. ,x
'a D
y =1F. gy where: E I'ex,ey
- 12, k1 2 bx,y)2 bx,y 6-14
(D v and the subscripts x,y reflect the particular bending plane of interest.
- f. Combined flexure and compression (or tension):
f f 0.65 y *F bx bx
+f by < 1.0 The above requirement should be met for both the direct tension or compression case.
6.7.2 Level D Service Limits F-1370 (Section III, Appendix F), states that the limits for the Level D condition are the minimum of 1.2 (Sy /F t ) or 0.7(S y/Ft ) times the corresponding limits for Level A condition. Since 1.2 yS is less than 0.7 Su for the rack material, and for the upper part of the suppcrt feet, the multiplying factor for the limits is 2.0 for the SSE condition for the upper section. The factor is 1.68 for the lower section under SSE conditions. Instead of tabulating the results of these six different Stresses as These dimensioned values, they are presented in a dimensionless form. so-called stress factors are defined as the ratio of the actual developed stress to its specified limiting value. With this definition, the limiting value of each stress factor is 1.0 for OBE and 2.0 or 1.68 for the SSE condition. 6.8 RESULTS 1 Figures 6.1; 6.2, and 6.3 show the pool slab motion in horizontal x, horizontal y, and vertical directions. This motion is for the SSE earthquake. Results are abstracted here for a 12x14 module (the largest module) and for a 7 x 14 module (largest aspect ratio) next to the cask pit. Results are also presented for the J1 rack which has a special configuration. O 6-15
= .-
A complete synopsis of the analysis of the 12x14, the 7x14 module, and the J1 module, subject to the SSE earthquake motions are presented in summary Tables 6.Sa through 6.51 which gives the bounding values of stress factors R$ (i = 1,2,3,4,5,6). The stress factors are defined as: R = Ratio of direct tensile or compressive stress on a net section 1 to its allowable value (note support feet only support compression) R = Ratio of gross shear on a net section to its allcwable value 2 R = Ratio of maximum bending stress due to bending about the x-axis 3 to its allowable value for the section
= Ratio of maximum bending stress due to bending about the y-axis R4 to its allowable value R = Combined flexure and compressive factor (as defined in 6.7.le 5
above)
= Combined flexure and tension (or compression) factor (as R6 defined in 6.7.lf above) ^> stated berore, the ei'owedie veioe or a4 (i =t.2.3.4.s.6) is t for the ost O condition, and 2 for the SSE (except for the lon'er section of the support where the factor is 1.68)
The dynamic analysis gives the maximax (maximum in time and in space) values of the stress factors at critical locations in the rack module. Values are also obtained for maximum rack displacements and for critical impact loads. Tables 6.Sa through 6.51 presents the results of the dynamic analyses. The fuel load in each case was assumed to be consolidated, with weight of 3000# per assembly. It is found that the results corresponding to SSE are most critical vis-a-vis tne corresponding allowable limits. The results given herein are for the SSE. The maximum stress factors (R g ) are below the limiting value for the SSE condition for all sections. It is noted that the critical load factors reported for the support feet are all for the upper segment of the foot and are to be compared with the limiting value of 2.0. f% G I 6-16 i a
Analyses (not included here) have been carried out to show that significant . margins of safety exist against local deformation of the fuel storage cell due to rattling impact of fuel assemblies and against local overstress of impact bars due to inter-rack impact. Benchmark analyses (not presented here) have also been carried out for the OBE condition. Comparison between OBE and SSE shown that the SSE condition controls the design of the fuel racks. 6.9 IMPACT ANALYSES 6.9.1 Impact Loading Between feel _ Assembly and Cell Wall The local stress in a cell wall is estimated from peak impact loads obtained ; from the dynamic simulations. Plastic analysis is used to obtain the limiting impact load that can be toleratec. Including a safety margin of 2.0, we find that the total limit load for tne number of cells (NC) is: QL = 9030 NC Q Impacts Between Adjacent Racks , 6.9.2 All of the dynamic analyses assu're, conservatively, that adjacent racks move completely out of phase. Thus, the highest potential for inter-rack impact is achieved. Based on the dynamic loads obtained in the gap elements simulating adjacent racks, we can study rack integrity in the vicinity of the impact point. The use of framing material around the top of the rack allows us to ' I withstand impact loads. It is shown that rack-to-rack impact loads can be accomodated with the load levels shown in Tables 6.51 and 6.5k without causing any permanent deformation in the active fuel region. Thus, impacts ! between racks can be accommodated without violating rack integrity. 6.10 WELO STRESSES j Critical weld locations under seismic loading are at the bottom of the rack at the baseplate connection and at the welds on the support legs. Results from the dynamic analysis using tne simulation codes are surveyed and the maximum loading is used to qualify the welds on these locations. 6-17
(3
# Baseplate to Rack Welds and Cell-to-Cely Welds 6.10.1 Section NF permits, for the SSE condition, an allowable weld stress t = .42 S u = 28,600 psi. The calculated weld stress, based on the highest load factor, is 7735 psi for the baseplate to rack welds.
The critical area that must be considered for cell-to-cell welds is the weld between gap channels and tubes. Where skip welding is used, this weld is continuous near the baseplate but is a skip weld (2" on 10.25" spacing) as we move up the tube. The critical shear stress in this area for the SSE condition is less than 14000 psi where account is taken of the skip welds or spot welds used in this region. Near the bottom of the rack where this shear stress dominates, there are no other stresses on the weld which need be considered. R Stresses in the channel-to-cell welds may also develop due to fuel assembly impact with the cell wall near the top of the rock. Inis will occur if fuel ("_T
,/ assemblies in adjacent tubes are moving out of phase with one another so tnat impact 10adt in two adjacent cells are in opposite directions which would tend to separate the channel from the tube at the weld. Our analysis shows that the maximum weld shear Stress in this area is less than 8000 psi.
6.10.2 Heating of an Isolated Cell Weld stresses due to heating of an isolated hot cell are also computed. The assumption used is tnat a single cell is heated, over its entire length, to a temperature above the value associated with all surrounding cells. No thermal gradient in the vertical direction is assumed 50 that the results are conservative. Using the temperatures associated with tnis unit, we show that the skip welds along the entire cell length do not exceed the allowable value for a thermal loading condition. () 6-19
O 6.11
SUMMARY
OF MECHANICAL ANALYSES V The mathematical model constructed to determine the impact velocity of falling objects is based on several conservative assumptions, sucn as:
- 1. The virtual mass (see Ref. 8-10 for further material on this subject) of the body is conservatively assumed to be equal to its displaced fluid mass. Evidence in the literature (Ref. 12),
indicates that the virtual mass can be many times higher. 2 The minimum frontal area is used for evaluating the drag coefficient.
- 3. The drag toefficients utilized in the analysis are the lower bound values reported in the literature (Ref. 13). In particular, at the beginning of the fall when the velocity of the body is small, the corresponding Reynolds number is low, resulting in a large drag coefficient.
- 4. The f ailing bodies are assumed to be rigid fcr the purposes of impact stress calculation on the rack. The solution of the imersed body motion problem is found analytically. The impact velocity thus computed is used to determlae the toaximum stress generated due to stress wave propagation.
4 O v With this model, the following anslyses are performed:
- a. DroppedfuelAccidentj A fuel assembly (weight = 1616 pounds) is dropped from 36 inches above the module and impacts the base. The final velocity of the dropped fuel assembly (just prior to impact)To is study calculated and, baseplate thus, the total eneroy at impact is known.
integrity,we assume that this energy is all directed toward punching of the baseplate in shears and thus transformed into work done by the supporting shear Stresses, it is determined that shearing deformation of tne baseplate is less than tha thickness of the basepiste so that we conclude that local piercing of the baseplate will not occur. Direct impact with the pool liner does not occur. l l The subcriticality of the adjacent fuel assemblies is not violated. I l
- b. Dropped fuel Accident 11 l
One fuei assembly drops fron 36 inches above the rack and hits the f top of the rack. Permanent deformation of the rack is found to be limited to the top region sur.h that the rack cross-sectional i geometry at the level of the top of the active fuel (and below) is not altered. The region of local permanent deformation does not extend below 6 inches from the rack top. An energy balance approach is used here to obtain the results. lO 6-19 i
i. t Jammed Fuel-handling Equipment {} c. A 4400-pound uplift force is applied at the top of the rack at the "weakest" storage location; the force is assumed to be applied onThe ; one wall of the storage cell boundary as an upward shear force.
- plastic deformation is found to be limited to the region well above the top of the active fuel, These analyses prove that the rack modules are engineered to provide maximum ,
safety against all postulated abnormal and accident conditions. 6.12 DEFINITION OF TERMS USED IN SECTION 6 S1, S2, 53, S4 Support designations pg Absolute degree-of-freedom number i ag Relative degree-of-freedom number i u Coefficient of friction Pool floor slab displacement time history in tne Ug i-th direction x y cocrdinatet horizontal direction (} z coordinr.te vertical direction , Kg impact spring between fu>l assemblies and cell Kg Linear component of friction spring l Kg Axial spring of support leg locations N Compression load in a support foot r Kp Rotational spring provided oy the pool slab Subscript i When used with il or X indicates direction (i = 1 x-direction, i = 2 y-direction, f = 3 z-direction) i I I r i 6-20
REFERENCES TO SECTION 6 () V i
- 1. USNRC Standard Review Plan NUREG-0B00 (1981).
- 2. ASME Boiler & Pressure Vessel Code, Section 111. Subsection W (1983).
- 3. USNRC Regulatory Guide 1.29 "Seismic Design Classificaticn," Rev. 3, 1978.
- 4. "Friction Coefficients of Wr.".er Lubricated Stainless Steels for a Spent Fuel Rack Facility," Prof. Ernest Rabinowicz, MIT, a report for Boston Edison Company, 1976.
- 5. USNRC Regulatory Guide 1.92, "Combining Medal Responses and Spatial Components in Seismic Response Analysis,' Rev. 1, February, 1976.
- 6. "The Component Element Method in Dynamics with Application to Earthquake and Vehicle Engineering," 5. Levy and J.P.O. Wilkinson, McGraw Hill, ,
1976.
- 7. "Dynamics of Structures," R.W. Clough and J. Penzien, McGre . Hill (1975).
- 8. "Mechanical Design of Heat Exchangers and Pressure Vessel Components,"
Chapter 16. K.P. Singh and A.I. Soler, Arcturus Publishers. Inc., 1984.
~s 9. R.J. Fritz, "The Effects of Liquids on the Dynamic Motions of Imme sed (d Solids," Journal of Engineering for Industry. Trans. of the ASME, February 1972, pp 167-172.
- 10. "Dynamic Coupling in a Closely Spaced Two-Body System Vibratino in Liquid Medium: The Case of fuel Racks," K.P. Singh and A.I. Soler, 3rd laternational Conference on Nuclear Power Safety, Keswick, England, May 1987.
- 11. USNRC Reguldtory Guioe 1.61, "Damping Values for Seismic Design e' Nuclear Power Plants," 1973.
- 12. "Flow Induced Vibration " R.D. Blevins, VonNostrand (1977).
- 13. "Fluid. Mechanics," M.C. Potter and J.F. Foss, Ronald Press, p 459 (1975).
I i O l 6-21 L -
O Table 6.1 DEGREES CF FREEDOM-Displacement Rotation e e e
~ Location U x
Uy U z x y 2 (Node) 45 1 pi p2 P3 94 96 2 Point 2 is assumed attached to rigid rack at the top most point. 2* p7 p8 3* pg p10 4* pii p12 pl3 p14 l 5* ! -- 0 1 s 915 pio
- j. where:
pg = gj (t) + Ui (t) i = 1,7,9,11,13,15 i
= qj (t) + U2 (t) i = 2,8,10,12,14,16 l
gj (t) + U3(t) i=3
=
Uj (t*; are the 3 known earthquake displacements, l l 1 O 6-22 l
ex U Table 6.2 NUMBERING SYSTEM FOR GAP ELEMENTS AND FRICTION ELEMENTS
- 1. Nonlinear Springs (Gap Elements)- (60 Total)
Node Location Description Number 1 Support 51 Z compression only element 2 Support S2 Z compression only element 3 Support S3 Z compression only element 4 Support $4 Z compression only element 5 2,2 X rsck/ fuel assembly impact element 6 2,2* X rack / fuel assembly impact element 7 2,2 Y rack / fuel assembly impact element 8 2,2* Y rack / fuel assembly impact element 9-24 Other rattling masses Bottom cross- Inter-rack impact elements p sj 25 section of rack (around edge) Incer-rack impact elements Inter-rack impact elements Inter-rack impact elements Inter-rack impact elements Inter-rack impact elements Inter-rack impact elements 44 Inter-rack impact elements 45 Top cross-section Inter-rack impact elements
, of risck Inter-rack impact elements (around edge) Inter-rack impact elements Inter -rack impact elements Inter-rack impact elements l
Inter-rack impact elements Inter-rack impact elements i 64 Inter-rack impact elements l l 1 (3
~J 6-23 i
i
. - = . . - . .. . _. . __
Table 6.2 (continued) NUMBERING SYSTEM FOR GAP ELEMENTS AND FRICTION ELEMENTS
- 11. Friction Elements (16 total)
Number Node Location Description 1 Support 51 X direction support friction 2 Support S1 Y direction friction 3 Support S2 X direction triction 4 Support S2 Y direction friction 5 Support S3 X direction friction 6 Support S3 Y direction friction 7 Support S4 X direction friction 8 Support 54 Y direction friction 9 S1 X Slab moment 10 51 Y Slab moment 11 52 X Slab moment 12 S2 Y Slab moment 13 S3 X Slab moment 14 S3 Y Slab momer*, (~ N 15 S4 X Slab moment 16 S4 Y Slab moment O V 6-24
Table 6.3 RACK MATERIAL DATA Young's Yield Ultimate Modulus Strength Strength Material _E (psi) Sy (psi) Su (PSI) 304L S.S. 27.9 x 10 6 23150 68100 Table Table Section III Table 1-3.2 Reference I-6.0 I-2.2 Table 6.4 SUPPORT MATERIAL DATA Young's Yield Ultimate , Modulus Strength Strength Material E (psi) Sy (psi) Su (PSI) 27.9 x 10 6 27,500 68,100 1 SA-351-CF3 (upper part of support feet) 27.9 x 10 6 62,400 90,000 2 SA-217-CA15 (lower part of support feet) l O 6-25 i
d Table 6.5a RACK
SUMMARY
M1'(7x14)
~ Fuel Assembly Rack / Rack Rack / Wall to' Cell Impact Impact Load Impact Load GB/8P GB/8P Run # Load (f) 4 M12- 4.035 x 10 4 8.343x10 0./0.
4 6.165x10 M18 4.354x10 4 2.062x10 5 0./0. 5 1.247x10 4 4 0,/0. M28a 3.30x10 -9.757x10 4 2.662x10 G8 = Girdle bar BP = Baseplate l ~ O 6-26
O-Table 6.5b RACK'
SUMMARY
M1 (7x14) STRESS FACTORS (Upper values for. rack base - Lower values for support feet) R Run i R3 R 2 R 3 R4 R 5 6 M12 .264 .054 .280 .347 .636 .715
.573 .142 .255 .251 .621 .806 M18 .240 .204 .455 .683 .933 1.057 .791 .559 1.001 .860 1.416 1.692 O M28a .104 .083 .246 .255 .394 .446 .332 .278 457 .495 .643 .771 l
O l 6-27 l
, - . - . - . _ ~ _ _ _ . . _ . , . _ , - _ _ _ . . . . _ . _ _ - _ _ _ . _ , _ . . . . - - . _ _ . . - _ - . _ _ . . . _ , , . _ . _ . , - , _ . . . .
.1 l
Table 6.5c , RACK:
SUMMARY
_M1 (7x14) Max. Max. Disp. Disp.. Max, Vert. Run # Remarks DX (in) OY (in) Disp. (in)
.5218 .7385 .075 M12 U = .2 Full Load of Fuel Assemblies 1.822 .322 .209 M18- V = .8 Full Load of Fuel Assemblies p = .8 ' ' 1% .4936 .2011 .0343 M28a -fuel load 1-positive x-half O
4 l l l l-1
~
O 6-28 i ln-_,.._,._.,_.-._._,_-__,._.__,..___._.~_,___
1 Table 6.5d RACK
SUMMARY
M1 (7x14) Max. Floor Max. Floor *
.Run # Remarks Load (#) Load (#)
4 feet Vertical / Shear M12 p = .2 Full Load 8.376x10 5 233593./44954, of fuel assemblies M18 p = .8 Full Load 7.'437x10 5 ,22164./184918. of fuel assemblies. M28a p = .8 50% 3.217x10 5 134546./83986. fuel load in positive x-half
.O
- vertical = vertical ioad Shear = shear load i
lO L 6-29 i
1 O Table 6.Se RACK
SUMMARY
C1 12x14; C3 RACK Max. Max. Disp. Disp. Max. Vert. Run i Remarks OX (in) DY (in) Disp. (in) C28 C1 rack, u = .8 .346 484 .076 Fuel centroid at x = 27.17" (1/2 full) C1 rack, u = .8 .346 425 .116 C18d 168 cells have fuel assemblies C1 rack, u = .2 .318 .363 .126 C12 168 cells have fuel assemblies C3 rack, u = .8 .1239 .1318 .144 C51 168 cells with fuel; adjacent racks per configuration #1 C3 rack, u = .8 .1375 .1321 .144 C53 168 cells with fuel; adjacent racks per configuration #2 O 6-30
4 Table 6.5f RACK
SUMMARY
C1 12x14; C3 RACK Max. Floor Max. Floor Run # Remarks Load (f) Load (#) 4 Feet Vertical / Shear 5 5 C1 rack, u = .8 6.221x10 2.297x10 C28 Fuel centroid at x = 27.17" 1.511x10 5 I' (1/2 full) 6 5 C1 rack, u = .8 1.013x10 3.831x10 C18d 168 cells have fuel assemblies 1.831x10 5 6 5 C1 rack, u = .2 1.131x10 3.585x10 l C12 168 cells have fuel assemblies .717x10 5 O 5 5 C3 rack, u = .8 9.957x10 2.724x10 C51 168 cells with fuel; adjacent 5 racks per .8319x10 l configuration #1 5 5 C3 rack, u = .8 9.767x10 2.724x10 C53 168 cells with fuel; adjacent 5 racks per
.9628x10 configuration #2
- Vertical = verticci load Shear = shear load O
6-31
l () Table 6.59 RACK SUKMARY J1 Max. Max. Disp. Disp. FMx. Vert. Run # Remarks OX (in) DY (in) Disp. (in) 35 cells .5722 .3601 .0395 J101 u = .8, filled with 3000# assemblies 35 cells .492 2.0788 .043 J201 u = .2, filled with 3000# assemblies 70 cells 1.214 .383 .185 Jill u = .8, filled with 3000# assemblies 70 cells .2893 .3348 .073 J112 u = .2, filled with 3000# assemblies
.8, 12 large .977 .628 .164 ' () J150 p =
cells with 10000# assemblies 12 large .5923 .3023 .015 J151 u = .2, cells with 10000# assemblies i l 1 l i l 6-32 i
-,,.,--__m . . , , _ . . . - - - , _ , . - - , - , _ . . _ _ , _ _ - _ . , . . . , - , . , - , - . _ . , _ , _ _ . _ , _ _ , , . , _ , , . . , . , _ , _ . , _ . _ _ _ . . . . , . - , _ , , . , , . . . . . . . . _ - _ _ , - . - _ . . , . .
,~ \/ Table 6.5h RACK
SUMMARY
J1 Max. Floor Max. Floor Run # -Remarks Load (#) Load (#) 4 feet Vertical / Shear J101 u= .8, 35 cells 2.057x10 5 118970. filled with 3000# 94444. assemblies 5 125000. J201 u= .2, 35 cells 2.392x10 filled with 3000# 24990. assemblies 5 226652. Jill u= .8, 70 cells 5.163x10 109011. filled with 3000# assemblies 5 195021. J112 u= .2, 70 cells 4.866x10 38439, filled with 3000# assemblies r~s 3.744x10 5 166930. (-) J150 u= .8, 12 large 86890. cells with 10000# assemblies 5 105420. J151 u= .2, 12 large 2.163x10 21084. cells with 10000# assemblies Vertical = vertical load Shear = shear load l l l l O 6-33 l c .- , ,.-,_ ,--.-- - .,..-,._. , - . - . . . ,., -, . - - . . - , , _ - - , - - - .- --
o - Table 6.51 RACK
SUMMARY
J1 Fuel Assembly Rack / Rack Rack / Wall to Cell Impact Impact Load Impact Load GB/BP GB/BP* Run i Load (#) J101 6.761x10 4 7.28x10 4 0./0. O. J201 5.361x10 4 6.827x10 4 0./0. 4 2.261x10 Jill 8.421x10 4 1.099x10 5 0./0. 4 2.749x10
.O 6.738x10 4 6.115x10 4
0./0. J112 O. J150 1.002x10 5 7.362x10 4 0./0. (at botton) 4 3.555x10 4 J151 9.583x10 4 2.863x10 0./0. (at bottom) 4 3.794x10
- BP = Baseplate GB.= Girdle bar O
6-34
Table 6.5j RACK
SUMMARY
J1 STRESS FACTORS (Upper values for rack base - Lower values for support feet) R R4 R R Run # R1 R 2 3 5 6
'J101 .092 .086 .242 .271 .448 .512 .293 .327 .604 .335 .741 .884 J201 .109 .028 .341 .212 .489 .559 .308 .082 .137 .151 .350 .454 Jill .233 .164 .375 .550 .684 .767 1.182 .558 .347 .524 .640 .969 J112 .216 .056 .320 .292 .530 .582 .480 .136 .250 .238 .539 .698 .155 .148 .324 .308 .427 .490 J150 .411 .320 .590 .549 .707 .855 l .089 .023 .222 .191 .309 .351 J151 l .260 .074 .136 .129 .295 .383 l
l O 6-35 i l
O Table 6.5k RACK
SUMMARY
C1, C3 RACK
~
Fuel Assembly Rack / Rack Rack / Wall
'to Cell Impact Impact Load Impact Load GB/BP GB/BP*
Run # Load'(#) C28 37570. 93410 0./0. O. C18d 70340. 101000 0 /0. O. 70600 67310 0./0. Q C12 46620 1.31x10 5 N.A. C51 69830. i-8654. l' t I 5 1.331x10 N.A. C53- 69870. 8727. GB = girdle bar BP = baseplate l l O I 6-36
Table 6.51 () RACK
SUMMARY
C1, C3 RACK STRESS FACTORS (Upper values.for rack base - Lower values for support feet) R 8 R R R 4 5 6 Run # R1 2 3
.325 .310 .433 .494 C28 .118 .076 .756 .807 1.171 1.403 .564 .474 .385 .264 .607 .679 C13d .193 .102 1.205 .967 1.402 1.675 .940 .673 483 .549 O C12 .215 .061 .250 .274 .431 .399 .991 1.286 .880 .241 .087 .266 .281 C51 .187 .052 .111 .528 .834 1.042 .669 .295 .366 .091 .266 .281 C53 .187 .053 .110 .571 .813 1.001 .669 .319 .473 i
I O 6-37 f
- o. -
0 0 5 1 O o. o e 4
't O
o. ks o. k s I o. o o t
's I
l o. 1 o o
\ 'o t }
o. o o
'o I
- o. )C I
h oE sS O 1 1 i 0 T 6
- o. (.
S o o - E fE W M .
/
I T G I T F o. S o A b _ s E l I I\ E o. S o S b s I o. K \ o C I b 4 A R L o. E o U ) b a F g N o. O o R b a Y B o. o O bt o o o. o o5o
.4. :)o g - o6 ? ' Y !a : o *9o 5 5 C 4 $' d 0
1 r o O I H
.) roO D % c "e
o _& D - _ e
.e c
Z - _- _ _c o g gO
~~ - T_ ,c-(g u; m --_m n --
y i .. (, i _ __=_ ._ Y _
^~ -- L ,. s -
G
.-_ -- t r <v -z
__ m _ g g o
^
s
.e _ k"'*% a ~ ~
t _
$0 t. .o ' . C m ^ ~ " ~ ~ ~ ._
4 .We g __,-=- .- _-
.- _n.
umme
-. __ e ammumme __
m . 2 h
~
S w m-__~&- _ ,
** / - .. . ==== - W _-
H N_ _
< W
_~ _ me
#:.?
e .
... 8 5 4 8
9 5 I 3 C O C-
- C C tr, -~
C 7~, 5; 1 , .~. cw e-- v; - i !
, v 1 I 5,C" 30C E-U 6-39
i t 4 W O , 8 f 8 l O. l a l
- e l SJ d d ku)x
< M C' Q - .- t E - ,I w
l W _S 4 l 5< Ch 1 9 J 5- - O ,
.s O
8 1 4 __.. cor'o oor'o oct'o co5$) g NOI$E373$h > n
- -----------------------.-..--.-..6_-4__0-___ _ __ _
0- Coupling Elements Z 2 j _7 A e e .
'~ a Qa a7 / p: -
9,,
* / $16
- li H
l ,
. = Gos 1
Ash Qu I
' ' I AY ='
o t 7 i
- $v / 7l , ~ ' / \
3/,S3 l ,l
'yS4 3 i2 f
1 i jfp =y--( o . AX i
~ $ 4
[ r k4
/
f 5 Support p S2 + mi f 1 m7 % X X3, YS - Location of fue! rod group mass centroid - relative tc q, centerline of fuel rack l O FIGURE s.4 Dynamic Model 6-41
TYP.TOPIMPACT
" ~ " '
O m
.l l } } ./
m W w
,W ^~t' r- i. H fff f! ! .
_.3$LCK STRUCTURE TYP. BOTTOM IMPACT _ELEhENT
) l l )a W ~
i >4
&w ^~f' & M ;f- ({ { {n, 1
l FIG. 6.5 G AP, ELEMENTS TO SIMUL ATE INTE_R-R ACK IMP ACTS 6-42
O impact Springs y
,i )
i! nn r ,_
' Mass
- 2*
- t;i r O f \ Fluid Dampers 4 ,
(not used in this analysis) Rigid Frame
=X FIGURE ~6.6 IMP ACT SPRINGS AND FLUID D AMPERS 6-43
Kx z n O Kw K, -f Y lwl 2 - [4 j 43 2 Gap Elements To Simulate Inter-Rack impacts Rack Centroid (4 for 2-D Motion ( Assumed At H/2) 20 for 3-D Motion) E H lQ [ Rigid Rack Baseplate Kw Kw 3 h) o Kg Kg h ld wj RbKg K K FIGURE e.7 Spring Mass Simulation For O Two-Dime n sio n al Motion 6-44
7.0 ENVIRONMENTAL EVALUATION n U 7.1
SUMMARY
Installation of high density spent fuel sto nge racks at the Byron Nuclear Power Station will increase the licensed storage capacity of the spent fuel pool from 1060 to a maximum of 2870 assemblies. Radiological consequences of expanding the capacity have been evaluated with the objective of determining if there is a significant additional onsite or offsite radiological impact relative to that previously reviewed and evaluated (Ref. 1). In addition, radiological impact to operating personnel has been evaluated to ensure that exposures remain as low as is reasonably achievable (ALARA). The decay heat loading and the radiological burden to the spent fuel pool water are determined almost entirely by refueling operations. The frequency of refueling operations and the conduct of refueling are independent of the increased capacity of the storage pool, except that the increased capacity should reduce fuel movement and allow continued normal operation. Since the fuel assemblies which will utilize the bulk of the storage capacity (and will g V ultimately fill all incremental capacity above that of the existing design) are aged, their contribution to either the peak decay-heat load or the increased radiological impact, in terms of increased doses, is negligible. A study perfermed by the NRC (Ref. 2) supports this conclusion. Consequently, the increase in the storage capacity of the spent fuel pool will neither significantly alter the operating characteristics of the current pool nor result in a measurable change in impact on the environment. 7.2 CHARACTERISTICS OF STORED FUEL \ l Because of radioactive decay, the heat generation rate and the intensity of gamma radiation from the spent fuel assemblies decreases substantially with l decay time. After a cooling time of about 4 years (Ref. 3), the decay heat I generation rate is less than 2% of the rate at 7 days, the nominal time at which depleted fuel assemblies are transferred to the spent fuel pool. The intensity of gamma radiation is very nearly proportional to the decay heat and O decreases w4th c->ing time 4n e s4m41er menner. 7-1
The bulk of the heat load is due to freshly discharged fuel; aged fuel contributes relatively little to the total heat load. Therefore, this expansion will not significantly increase the thermal dissipation to the environment. Since the intensity of gamma radiation follows the decline in decay heat generation rate, it is similarly concluded that there will be no significant increase in gamma radiation beyond the pool environment due to the expanded storage. It is important to note that the aged fuel in the expanded storage capacity will not contain significant amounts of radioactive iodine or short-lived gaseous fission products, since these would have decayed during the storage period. The Krypton-85 which might escape from defective fuel assemblies has been shown to do so quickly (Ref. 2) (i.e., within a short time af ter discharge from the core). Further, the residual Krypton-85 will be contained within the fuel pellet matrix and hence any leakage would occur at very low rates (Ref. 2). Cesium 134/137 (Ref. 2) is strongly bound within the fuel Any pellet matrix and its dissolution rate in water is extremely small. r~ Cesium dissolved in the pool water is easily controllable in the cleanup (Q system (demineralizer-ion exchanger resin bed) (Ref. 2). Thus, the planned release of gaseous storage expansion will not significantly increase the radionuclides. 7.3 RELATED INDUSTRY EXPERIENCE Experience with storing spent fuel underwater has been substantial (Refs. 2, 3, and 4). These references show that the pool water activity, normally low, experiences'small increase during refueling periods, which then decays rapidly with time. Typical concentrations (Ref. 5) of radionuclides in spent fuel pool water range from 10 4 t.Ci/mi to 10-3 aci/ml, with the higher value associated with refueling operations. References 2 and 5 also state that the increase in pool water activity during refueling can be attributed to: 0 Dislodging (sloughing off) of corrosion products on the fuel assembly during transfer and handling operations. O The possible short-term exposure of fuel pellets to pool water via a O cladding defect. 7-2
0 Mixing of the spert fuel pool water with the higher activity reactor coolant. Upon cessation of the refueling operations, the fuel pool water ard the reactor coolant system would be isolated fron each other, thereby terminating transports of corrosion products from the reactor ccolant system. Thus, deposition of crud is a function of refueling operations and is not impacted by the expanded storage. O Furthermore, it has been shown (Ref. 6) that release of fissicn products from failed fuel decreases rapidly after shutdown to essentially negligible levels. The dissolution of exposed fuel pellets (made of UO 2
) is very slow in water at fuel pool temperatures and the corrosion of the cladding (Zircaloy 4) at spent fuel pool water temperatures is virtually nil (Refs. 2 and 5).
Another me:.hanism available for the release cf the gaseous fission productsisdiffusionthroughtge00at the loo watercoefficient diffusion temperatures is (<l50 F), 7 p extremely small (Ref.7). Therefore, the small increase in activity of the spent fuel pool water is due to either crud transport, fission products release, or cross-flow from the reactor coolant system, and 15 only a function of refueling operations. The expansion of ruel pool storage capacity will not cause a significant increase in doses either onsite or offsite. The corrosion properties of irradiated Z!rcaloy cladding have been reviewed in References 2 and 4 and the c v clusion is drawn that the o corrosion of the cladding in spent fue'. pool water is negligible. V The minor incremental heating of poo's water, due to the expansion of storage capacity, is far too small to r.aterially affect the corrosion properties of Zircaloy ci ddding. 7.4 BYRON NUCLEAR POWER STATION EXPERIENCE At present, there is one spent (normal core offload) batch of spent fuel assemblies stored in the spent fuel pool. 7.5 SPENT FUEL POOL COOLING AND CLEANUP SYSTEM lt has been shown previously in Section 5 of this licensing report that the
- cooling system at Byron is adequate to handle the expected heat loads and l
maintain the pool temperature peaks within acceptable limits. it has been shown in Section 5 that the small increase in heat load due to the storage capacity expansion will neither significantly increase the thermal dissipation to the environment nor increase the propensity for corrosion of the cladding. l It has also been shown that the crud deposition in the spent fuel pool water O occurs dorins refuei4ns oute9es ene that the olemnee expension w411 not 7-3 l l l L
m
' increase long-term crud deposition. The fuel pool cleanup system (filter and h demineralizer) is designed to maintain fuel pool water clarity and is operated and maintained in accordance with tne Byron operating procelures. The cleanup system takes a surface skim from the fuel pool and cleans it through a process of filtration and demineralization to prevent crud buildup on the fuel pool walls at the wa ter-to-air inter face.
The spent fuel pool water is sampled and analyzed periodically to confirm proper operation o f the pool cleanup system. The spent fuel pool modification will not result in a significantly higher quantity of solid radwaste. 7.6 FUEL P0OL RADI ATION SH1 ELO!in 7.6.1 Source Terms The spent fuel gama source terms used for the fuel pool shielding evaluation were generated using the point reactor fission product inventory code RIBD. The following assumptions were used in the analysis: 0 Initial fuel enrichment = 4.2% 0 Net Reactor core power level = 3565 MWt 0 Average assembly discharge burnup = 45,750 MWD /MTU 0 Power level for average assembly = 17.67 MWt 0 Power level for peak power assembly = 29.16 MWt (peaking factor of 1.65) O Burnup for peak power assembly = 33,000 MWD /MTU (one 18-month cycle at maximum power) O Start of Pefueling Process = 100 hours af ter shutdown 0 Time after shutdown for dose rate calculotion = 9 days The peaking factor of 1.65 and burnup for peak power assembly of 33,000 MWD /MTU for one cycle were chosen to produce the highest possible gama source term attainable unoer operational conditions. The average assembly burnup and power level were chosen to represent a conservative gama source term for the 4n4 mum Permitted Period before O spent fuei. 18e 100-hour decer time is the 7- 4
f-refueling can begin per the Technical Specifications. gq V The photon energy production rates of an average spent fuel assembly is given in Table 7.1. 7.6.2 Radiation a Shielding For an equilibrium refueling cycle 84 fuel assemblies will be discharged into the pool starting no earlier than 100 hours after reactor shutdown at a rate of I fuel assembly per hour. For a full core discharge, 193 assemblies are discharged. To evaluate the adequacy of the shielding capability of the spent fuel pool walls, the radiation dose from all racks filled to capacity with assemblies decayed 9 days after shutdown. It was found that only the first row of racks closest to the wall contributes significantly to the dose rate on the other side of the wall due to self-shielding. Only the "first row" results were used. O The north and south perimeter walls of the spent fuel pool are 5 feet thick. The east and west perimeter walls of the pool are 6 feet thick. The wall separating the spent fuel pool from the transfer canal is 5 feet 6 inches thick. These walls, together with the water and fuel storage racks are incorporated in the dose rate calculations for the adjacent areas. The results of the calculations are provided in Table 7.2. Except for the area immediately adjacent to the Al rack, the maximum calculated dose rates through the pool walls are less than the currently designated radiation level limit of 20-mrem /hr for these areas. The area l adjacent to the Al rack is normally occupied only briefly during the local and integrated leak rate tests which are conducted every 5 years. Access to the area is controlled. The dose rates in Table 7.2 are considered an upper limit since they are c calculated for freshly discharged fuel. The dose rates will reduce by a factor of six, 60 days after the fuel is discharged into the pool. {i 75 l I
O The above discussion indicates that the shielding available arctnd the spent fut.1 pool is adequate fur installation of the high density storage racks. The radiation dose level at the side of the pool and on the spent fuel pit crane bridge due to the transfer of a peak power fuel assembly are 2.5 mrem /hr and 2.0 mrem /hr, respectively. These calculations assume a minimum of 10 feet of water cover over the active fuel. If the transfer of fuel assemblies into the spent fuel cask pit becomes necessary, the water cover will be less than 10 feet over the active fuel, however, the dose rate during this fuel
.novement will be much lower than for the transfer of the peak power assembly noted above since the radiation level of the fuel assembly will have had time to decay to a level which would more than compensate for the loss of water cover shielding.
Since the fuel transfer operation normally lasts less than 4 days (88 assemblies at 1 assembly per hour), the above radiation field does not create excessive operator exposure. U 7.7 RADIOLOGICAL CONSEQUENCES The design basis fuel handling accident (dropped assembly) in the Fuel Handling Buildin in Section 15.7 of the FSAR was reviewed for possible effects on radiological dose consequences. The review determined that the conclusions in the FSAR were still valid and that offsite radiological dose consequences were well within 10CFR100 limits. 7.8 RERACKING OPERATION installation of the fuel racks will include removal of the existing racks, making minor pool modifications, and cleaning and installing the new racks. The existing racks are bolted to the pool floor. The new racks will be cleaned prior to installation. The fuei handling This building overhead crane will be used to place the racks in the pool. effort wiil be performed efter the first refveiin9 for uait 1. The exist 4ns O 7-6
fuel racks will have been exposed to spent fuel and should only be nominally .m contaminated. Therefore, doses to individuals-involved in the reracking should be minimal. All pool modifications that can be completed prior to filling the pool with water have been done to minimize underwater work. Although divers may be needed for some tasks, all of the work associated with the installation will be sequenced to minimize potential radiation exposure of personnel due to the spent fuel located in the pool. AtARA considerations will be fully incorporated in the installation procedures for this condition. When the fuel handling building overhead crane is used over the pool, physical stops will be installed on the crane to preclude carrying racks over any stored fuel assemblies. Exact disposition of the existing racks has not been determined. They will be decontaminated and/cr ed and disposed of in accordance with the applicable Federal ano s. .e regulations. O 79 cosCoustoss Based upon the industry experience and evaluations discussed in previous sections, the following conclusions are made: Minor increases in radiological burden to the pool water, if any, can be adequately handled by the fuel pool cleanup system (filter and demineralizer), thereby maintaining the radionuclide concentration in the water at an acceptably low level. No appreciable increase in solid radioactive wastes (i.e., filte media and demineralizer resin) is anticipated. No increase in release of radioactive gases is expected since any long-lived inert radioactive gas potentially available for releast. (i.e. Kr-85) will have leaked from the fuel either in the reactor core during operation or during the first few months of residence in the pool. Further, Vol. 1, Reference 3 (pp. 4-16) has shown airborne in activity to be considerably lower than that allowable by Table 1 of 10CFR Part 20, Appendix B. Therefore, the planned expansion will not significantly increase the release of radioactive gases. O O 7-7
The existing spent fuel pool cooling system will keep the pool weter te pereture et en eccepteele ievel (see section s. O Thermal-Hydraulic Considerations). The existing radiation protection monitoring systems and program are adequate to detect and to warn of any unexpected abnormal increases in radiation level. This provides sufficient assurance that personnel exposures can be maintained as low as is reasonably achievable. Since the reracking will occur after the first refueling, procedural <:ontrols and necessary precautions will be taken to reduce radiation exposure to as low as is reasonably achievable, and hence, radiological impact will be minimized. Expanding the storage capacity of the spent fuel pool will not significantly increase the onsite or offsite radiological impact above that of the currently authorized storage capacity, nor is any significant increase in environmental radiological or nonradiological impact anticipated. O G 7-8
-y___ __ _ _
s
- \
REFERENCES TO SECTION 7 (~)$
- 1. FSAR, Byron Nuclear Power Plant.
- 2. NUREG-0575, "Handling and Storage of Spent Light Water Power Reactor Fuel, Vol.1, Executive Summary and fext, USNRC, August 1979.
- 3. NUREG-0800, USNRC Standard Review Plan, Branch Technical Position ASB9-2, Rev. 2. July 1981.
- 4. J.R. Weeks, "Corrosion of Materials in Spent Fuel Storaqe Pools," BNL-NUREG-2021, July 1977.
- 5. A.B. Johnson, Jr. "Behavior of Spent Nuclear Fuel in Water Pcol Storage," BNWL-2.36, September 1977.
- 6. J.M. Wright, "Expected Air and Water Activities in the fuel Storage Canal," WAPD-PWR-CP 1723 (with Addendum), undated.
- 7. ANS 5.4 Proposed Standard, "Method for Calculating the Fractional Release of Volatile Fission Products from 0xide Fuel," American Nuclear Society, issued for review, 1981.
- 8. "Licensing Repc.*t on High Density Spent Fuel Racks for Quad Cities, Units 1 and 2," Docket Nos. 50-254 and 50-265, Commonwealth Edison Compr.ny.
g June 1981. O 9. "Licensing Report for High Density Spent Fuel Storage Racks," Rancho Seco Nuclear Generating Station. Sacramento Municipal Utilities District, Docket No. 50-312, June 1982.
- 10. Finsl Safety Analysis Rtport, Limerick Generating Station Units 1 and 2 Section 9.1.
- 11. Safety Evaluation Report Related to the Operation of Limerick Generating Station Units 1 and 2, HUREG-0991, August 1983.
- 12. Source Term Data for Westinghouse Pressurized Water Reactors, WCAP-8253, July 1975.
7-9
- : +~
Table 7.1 PHOT 0N ENERGY PRODUCTION RATES OF AN AVERAGE SPENT FUEL ASSEMBLY (9 OAYS AFTER SHUTOOWN) Photon Energy Photon Energy Production Rate (MeV) ___ (MeV/sec)
< .05 4.16E+18 .05 .10 3.36E+17 .10 .30 2.24E+17 30 .55 1.40E+17 .55 .90 1.27E+17 .90-1.1- 4.50E+15 1.1-1.6 ;1.70E+16 1.6-2.0 .1.82E+14 2.0-2.4 3.13E+14 2.4-2.6 3.26E+14 > 2.5 9.42E+12 /~N U
TOTAL 5.01E+18 l , O 7-10 ? r
.4 J
Table 7.2 CALCULATED DOSE RATES IN AREAS ADJACENT TO THE SPENT FUEL POOL (9-DAYS AFTER SHUT 00VN) High Density Rack Location Dose Rate (mrem /hr) Floor el. 401 ft, 0 in., argas 54.0 adjacent to the north walls Floor el. 426 ft, O in., areas 2.3 adjacent to the edge of the pool. Floor el. 401 ft, 0 in...ppent fuel 4.0 pool _ heat exchanger area Fuel transfer canal ** 20.0 0
- A design water gap of 4-7/8 inches between the high density rack and the wall is used.
A design water gap of 4 inches between the high density rack and the wt11 is used. O 7-11
/s 8.0 IN-SERVICE SURVEllLAtlCE PROGRAM F0?. BORAFLEX NEUTR0!! ABS 0RBif4G MATERI At. 8.1 PROGRAM INTENT A sampling program to verify the integrity of the neutron absorber material employed in the high density fuel racks in the long-term environmen+. is described in this section. The program is conducted in a manner which allows access to the representative absorber material samples without disrupting the integrity of the entire fuel storage system. The program is tailored to evaluate the material in normal usa mode and to forecast future changes using the data base developed.
8.2 DESCRIPTION
OF SPECIMEtiS The absorber material used in the surveillance program, henceforth referred to '- as poison, is representative of the material used within the storage system. It is of the same composition, produced by the same method, and certified to the same criteria as the productinn lot poison. The sample coupon is of similar thickness as the poison u :d within the storage system and not less than 4 by 2 inches on a side. Figure 9.1 shows a typical coupon. Each poison specimen is encased in a stainless steel jacket of an identical alloy to that used in the storage system, formed so as to encase the poison material and fix it in a position and with tolerances similar to the design used for the storage system. The jacket has to be closed by tack welding in such a manner as to retain its form throughout the test perind and still allow rapid and easy opening without causing mechsnical damage to the poison specimen contained within. The jacket should permit wetting and venting of the specimen similar to the actual rack environment. 8.3 SPECIMEN EVALUATION After the removal of the jacketed poison specimen from the cell at a Q designated time, a careful evaluation of that specimen should be made to 8-1
determine its actual condition as well as its apparent durability for V continued function. Separation of the poison from the stainless steel specimen jacket must be performed carefully to avoid mechanical damage to the poison specimen. Immediately af ter the removal, the specimen and jacket section should be visually examined for any effects of environmental exposure. Specific attention should be directed to the examination of the Functional stainless steel jacket for any evidence of physical degradation. evaluation of the poison material can be accomplished by the following measurements: 0 A neutron radiograph of the poison specimen aids in the determination of the maintenance of uniformity of the boron distribution. 0 Neutron attenuation measurements will allodConsideration evaluation of the must continued nuclear effectiveness of the poison. be given in the analysis of the attenuation measurements for the level of accuracy of such measurements, as indicated by the degree of repeatability normally observed by the testing agency. A measurement of the hardness of the poison material will establish the continuance of physical and structural durability. The hardness O acceptability criterion requires that the specimen hardness will not O reduce the hardness listed in the qualifygg test d9cument for laboratory test specimen irradiated to 10 ra d s . The actual hardness measurement should be made after the specimen has been withdrawn from the pool and allowed to air dry for not less than 48 hours to allow for a meaningful correlation with the pre-irradiated sample. O Measuremen;. of the length, the width, and the average thickness and comparison with the preexposure data dill indicate dimensional stability within the variation range reported in the Boraflex laboratory test reports. A procedure will be prepared for execution of the test procedure and interpretation of the test data. l l lo l 8-2 l i b
1 l O m a. g .07 5 THK. FILLER - / FOR REGION 1. o 040"THK FILLER { FOR REGION 2.
- a. ( 304 L SST. ).
S / a: O
/
l TACK ' WELD _/
/ a 7 N g b ' /. g 3
- 7 ~:,
l sk 9 ?L-r
- 4 I- / l l ' i O ;%
I
/ k I so I
v ,/ /' lk; y NEUTRON ,
/ , jy / / '
ABSORBER s' s '
$ $ s' $%$ l' 3 g '- qf j 3 N .
S* 4 I, / . 06$THK. S.S. -30 4L . .O FIG. 8.1 - TEST COUPON S' . - - - -
9.0 COST / BENEFIT ASSES _SMENT - U A cost / benefit assessment has been prepared in accordance with the requirements of Reference 1, Sect.on V - Part 1. The assessment demonstrates that the installation of high density spent fuel storage racks is the most advantageous means of handling spent fuel . The material is presented merely for informational purposes. It is CEC 0's position that an environmental impact statement need not be prepared because the installation of high density fuel racks provides no significant impact on the environment. NRC precedent establishes that alternatives and economic costs need not be discussed when there is no significant environmental impact. However, for completeness, alternatives to reracking for additional spent fuel storage capacity are discussed in Section 9.3. 9.1 SPECIFIC NEEDS FOR SPENT FUEL STORAGE Disposal of Byron nuclear fuel is scheduled to be carried out by the Department of Energy in or af ter 1998 in accordance with Public Law 97-425 Nuclear Waste Policy Act of 1932. As Byron spent fuel- may not be accorded a high priority under the DOE Program, CECO is seeking to provide a spent fuel No storage capacity to support approximately 24 years of nominal operation. other contractual arrangements exist for the interim storage, or reprocessing of spent fuel from Byron Nuclear Power Station. Therefore, increased storage capacity in the spent fuel pool is the only viable option under consideration. Table 1.1, the Fuel Pischarge Schedule, indicates that with the high density spent fuel racks, loss of full core discharge capability (FCDC) will occur in 2009, 15 years beyond the current capability and 11 years beyond the scheduled repository fuel receipt date, per the DOE Mission Plan. 9.2 COST OF SPENT FUEL STORAGE The design and manufacture of the spent fuel storage racks will be undertaken by the organizations described in Section 1. It is expected that the total Project cost *in be bet ee" 52 5 "4 $2 9 mi" 'o"- O 9-1
O 9.3 At1ERNativES TO SeEnt FUEt 510a^oE CECO has considered the various alternatives to the pecposed onsite spent fuel storage. These elternatives are discussed below:
- a. Shipment of Fuel to a Reprocessing or Independent Spent Fuel Storage / Disposal Facility No commercial spent fuel reprocessing facilities are presently operating in the United States. CEC 0 has made contractual arrangements whereW spent nuclear fuel and/or high level nuclear waste will be accepted and disposed of by the U.S. Department of Energy. However, such services are not expected to be available before 1998. The existing Byron spent fuel storage capacity will not provide full core discharge capability beyond 1994. Spent fuel acceptance and disposal by the Department of Energy is not, therefore, an alternative to increased onsite pool storage capacity,
- b. Shipment of Fuel to Another Reactor Sit,e, Shipment of Byron fuel to another reactor site could provide short-term relief to the storage capacity problem. However, transshipment of spent fuel merely serves to transfer the problem to another site and does not result in any additional net long-term storage
- d. capacity. Accordingly, CECO does not consider the transshipment of spent fuel to be an appropriate alternative to high density spent fuel storage at the site.
- c. Not Operating the Plant af ter _t,h,e Current Spent Fuel Storage Capacity is Exhaustef As indicated in NUREG-0575, "Final Environmental Impact Statement on Handling and Storage of Spent tight Water Power Reactor Fuel," (Ref.
- 2) the replacement of nuclear power by coal generating capacity would cause excess mortality to rise from 0.59 - 1.70 to 15 - 120 per year for 0.8 GWY(e). Based on these facts, not operating the plant or shutting down the plant af ter exhaustion of spent fuel discharge capacity is not a viable alternative to high density storage in the spent fuel pool. The prospective 1986 expenditure of approximately $2.6 million for the high density racks is small
;ompared to the estimated value of replacement power equivalent to the plant's energy output: approximately $21 million per month in 1994 and $32 million per month in 2009.
A The subject of the comparative economics associated with various spent fuel options is the subject of Chapter 6 of NUREG-0575 (Ref. 2). Although the material presented is generic, it is of value in comparing the costs of the
- O 9-2
i various options._0f the options presented in that chapter, high density spent A fuel storage at the site is the most economic option at $ 3.50 per Kgu. The price of Independent Fuel Storage Facilities (IFSF), if available, would be
$54.35 per Kgu.
9.4 RESOURCE COMMITMENTS
-The expansion of the Byron spent fuel storage capacity will require the following primary resources :
0 Stainless steel - 284,815 lb/ unit 0 Boraflex neutron absorber - 22,000 lb/ unit o f which 10,925 lb is boron carbide (B C) 4 powder The requirement for stainless steel represents a small fraction of the total domestic production for 1986 (Ref. 3). Although the fraction of domestic production of B4 C required for the fabrication is somewhat higher than that for stainless steel, it is unlikely that the commitment of B4 C to this project will affect other alternatives. Experience has shown that the production of B C is highly variable and depends on need but could easily be expanded to 4 accomodate additional domestic needs. l O 9-3 l t
( ' O REFERENCES TO SECTION 9_
- 1. B.K. Grimes, "0T Position for Rev'iew and Acceptance of Spent Fuel Storage and Handling Applications," April 14, 1978.
2. NUREG-0575, "Final Environmental Impact Statement on Handling and Storage of Spent Light Water Power Reactor Fuel," Vols.1-3, USNRC, August 1979.
- 3. "Mineral Facts and Problems," Sureau of Mines Bulletin 671,19?O.
O O 9-4
10.0 QUALITY ASSURANCE PROGRAM V
10.1 INTRODUCTION
This chapter provides a general description of the quality assurance program that is implemented to assure that the quality objectives of the contract specification are met. 10.2 GENERAL The quality assurance program used on this project is based upon the system described in Oat's Nuclear Quility Assurance Manual . This system is designed to provide a flexible but highly controlled system for the design, manufacture, and testing of customized components in accordance with various The Oat Nuclear Quality codes, specifications, and regulatory requirements. Assurance Program has been accepted by ASME and has been approved by the CECO Quality Assurance Department and placed on CECO's Qualified Suppliers List. O The philosophy beh4ed Oet's Quel 4ty Assurence System is thet it sheii prov4de for all controls necessary to fulfill the contract requirements with sufficient simplicity to make it functional on a day-to-day basis. The system readily adapts to different designs and component configurations, making The following possible the construction of many varied forms of equipment. paragraphs provide an overview of the system and how it has been applied to Commonwealth Edison's speci fications. 10.3 SYSTEM HIGHLIGHTS The design control is organized to provide for careful review of all contract These requirements to extract each individual design and quality criterion. l criteria are translated into design and quality control documents customized to the contract requirement and completely reviewed and approved by responsible and qualified personnel . 10-1
The system for control of purchased material includes generating detailed p. V descriptions of each individual item of material along with specifications for, any special requirements such as impact testing, corrosion testing, monitoring or witnessing of chemical analysis, provision of over-check specimens, special treatments or conditioning of material, source inspection, and provision of performance documentation on any of the above. Material receipt inspection includes a complete check of all material and its documentation. Upon acceptance, each item of material is individually listed on a control sheet issued once a week to assure that only accepted material goes into fabrication. The fabrication control system provides that a shop traveller is prepared for The traveller is generated each subassembly and assembly in each contract. specifically to provide step-by-step instructions for fabrication, inspection, testing, cleaning, packaging, etc., which address all standard and special requirements of the contract specifications. Special attention is given to deployment of fabrication sequence and inspection steps to preclude the p possibility of missing poison sheets or incorrect sheets (incorrect 8 10 loading). All nondestructive examination procedures and test procedures are custom written to apply to CECO's requirements. The system provides for qualification and written certification of personnel performing quality-related activities including nondestructive examination and fat,rication inspection, welding, engineering, production supervision, and auditing. Other CECO requirements are fully covered in the Quality Aaurance Program, including document control, control of measuring and test equipnent, control of nonconforming material and parts, corrective action auditing, and other areas as specified by CECO. 10.4 SUMdARY Oat's quality assurance system provides the full measure of quality assurance - required by the contract. All special requirements of the specifications are 10-2 l
, , -, .- - - -- . ~ . - - - -- .-.-
covered, including source inspection of material and witnessing of material f testing by the engineer, furnishing of material certifications anel test reports within 5 days of shipment, and obtaining verification of qualification testing of poison materials. O i l l l l f O 10-3 l , _ __ --- -. . _ ..
. . i 5-
]- i I !9 , 1 i i 4 < i I j 1 1 l t l i : t i
,\?PENDIX A ;
I i BENCHMARK CALCtJLATIO!!S i i I G : s S
?
I e , F i I h i h i l r l t O A-1 I I l
- 1. INTRODUCTION AND SUf1 MARY The objective of this benchnarking study is to verify both the AMPX (NITAWL)-KENO (Refs. 1 and 2) methodology with the 27-
. group SCALE cross-section librarv (Refs. 3 and 4) and the CASMO-2E code (Refs. 5, 6, 7, and 8) for use in criticality calcula-tions of high density spent fuel storage racks. Both calcu-lational methods are based on transport theory and have been benchmarked against critical experiments that simulate typical spent fuel storage rack designs as realistically as possible.
Results of these benchmark calculations with -both methodologies are consistent with correspondino calculations reported in the* literature and with the recuirements of Reaulatory Guide 3.41, Rev. 1, May 1977. Results of these benchmark calculations show that the 27-group (SCALE) AMPX-KENO calculations consistently underpredict the critical eigenvalue by 0.0106 2 0.0048 ak (with a 95% proba-bility at a 95% confidence level) for critical experiments selected to be representative of realistic spent fuel storace rack configurations and poison worths. Similar calculations by Westinghouse suggest a bias of 0.012 i 0.0023, and the results of ORNL analyses of 54 relatively "clean" critical experiments show a bias of 0.0100
- 0.0013.
Similar calculations with CASMO-2E for clean critical experiments resulted in a bias of 0.0013 i 0.0019 (95%/95%). CASMO-2E and AMPX-KENO intercomparison calculations of infinite arrays of poisoned cell configurations show very good agreement and suggest that a bias of 0.0013 i 0.0018 is the reasonably expected bias and uncertainty for CASMO-2E calculations. Calculational Methods for Nuclear Criticality
- Validation of S a f e :.y . (See also ANSI N16.9-1975.)
A-2
The benchmark calculations reported here indicate that either the 27-group (SCALE) AMPX-KENO or CASMO-2E calculations are acceptable for criticality analysis of high density spent fuel storage racks. The preferred methodology, however, is to perform independent calculations with both code pack. ages and to utilize the higher, more conservative value for the reference design infinite multiplication factor.
- 2. AMPX (NITAWL)-KENO BENCHMARK CALCULATIONS Analysis of a series of Babcock & Wilcox (B&W) critical experiments (Ref. 9), which include some with absorber sheets typical of a poisoned spent fuel rack, is summarized in Table 1, as calculated with AMPX-KENO using the 27-group SCALE cross-section library and the Nordheim resonance integral treatment in NITAWL. The mean for these calculations is 0.9894 i 0.0019, conservatively assuming the larger standard deviation calculated from the keff values. With a one-sided tolerance factor V (K = 2.502), corresponding to 95% probability at a 95% confidence level (Ref. 10), the calculational bias is +0.0106 with an uncer-tainty of *0.0048.
Similar calculational deviations reported by Westinghouse (Ref. 11) are also shown in Table 1 and suggest a bias of 0.012
- 0.0023 (95%/95%). In addition, ORNL (Ref. 12) has analyzed some 54 critical experiments using the same methodology, obtaining a mean bias of 0.0100
- 0.0013 (95%/95%). These published results are in good agreement with the results obtained in the present analysis and lend further credence to the validity of the 27-group AMPX-KENO calculational model for use in criticality analy-sis of high density spent fuel storage racks. Variance analysis of the data in Table 1 suggests the possibility that an unknown factor may be causing a slight]y larger variance than might be -
However, such a expected from the Monte Carlo statistics alone. O A-3
Table 1
}
RESULTS OF 27-GROUP (SCALE) AMPX-KENG CALCULATIO!:S OF B&W CRITICAL EXPERIMENTS Westinghouse Experiment Calculated Calculated-meas. a k,ff Number keff I 0.9889 *0.0049 -0.008 II 1.0040 *0.0037 -0.012 III 0.9985 i0.0046 -0.008 IX(1) 0.9924 t0.0046 -0.016 0.9907 i0.0039 -0.008 L X l XI 0.9989 i0.0044 +0.002 l l l XII 0.9932 10.0046 -0.013 0.9890 i0.0054 -0.007 XIII XIV 0.9830 iO.0038 -0.013 0.9852 0.0044 -0.016 I) XV XVI 0.9875 t0.0042 -0.015 l 10.0041 -0.015 XVII 0.9811 XVIII 0.9784 10.0050 -0.015 XIX 0.9888 i0.0033 -0.016 XX 0.9922 iO.0048 -0.011 XXI 0.9783 iO.0039 -0.017 Mean 0.9894 0.0011(2) -0.0120 i 0.0010 Bias 0.0106 *0.0019(3) 0.0120
- 0.0010 Bias (95%/95%) 0.0106 *0.0048 0.0120 i 0.0023 Maximum Bias 0.0154 0.0143 (1) Experiments IV through VIII used B 4C pin absorbers and were (2)not considered representative of poisoned storage racks.
(3)Calculated Calculated fromfrom k eggindividual values andstandard deviations. used as reference. O A-4
factor, if one truly exists, is too small to be resolved on the V basis of critical-experiment data presentiv available. No trends aan, with absorbor sheet in kegg with intra-assembly water reactivity worth, or with s o '.u b l e poison concentration were identified.
- 3. CASMO-2E BENCHMARK CALCULATIONS 3.1 GENERAL The CASMO-2E code is a multigroup transport theory code utilizing transmission probabilities to- accomolish two-dimen-
' sional calculations of reactivity and depletion for BWR and PWtt fuel assemblies. As such, CASMO-2E is well-suited to the criti-spent fuel storace racks, since general caiity analysis of practice is to treat the racks as an infinite medium of storage cells, neglecting leakage effects. O CASMO-2E is closely analogous to the EPRI-CPM code (Ref. 13) and has been extensively benchmarked against hot and cold crit-ical experiments by Studsvik Eneraiteknik (Refs. 5, 6, 7, and 8). Reported analyses of 26 critical experiments indicate a mean k,ff of 1.000
- 0.0037 (lo). Yankee Atomic (Ref. 14) has also reported results of extensive benchmark calculations with CASMO-2E. Their analysis of 54 Strawbridge and Barry critical experi-ments (Ref. 15) using the reported buckling indicates a mean of 0.9987 * .0.0009 (lo), or a bias of 0.0013 t 0.0018 (with 95%
l confidence level). Calculations were probability at a 95% Strawbridge and Barry experiments repeated for seven of the Significantly large trends in keff with water cap and with ab- ' sorber sheet reactivity worth have been reported (Ref. 16) for GAM-THERMOS library. AMPX-KENO calculations with the 123-croup O l l A-5
selected at random, yielding a nean ke gg of 0.9987 i 0.0021 (le), thereby confirming that the cross-section library and analytical methodology being used for the present calculations are the same as those used in the Yankee analyses. Thus, the expected bias for CASMO-2E in the analysis of "clean" critical experiments is 0.0013 i 0.0018 (95%/95%). 3.2 BENCHMARK CALCULATIONS CASMO-2E benchmark calculations have also been made for the B&W series of critical exoeriments with absorber sheets, sinu-lating high density spent fuel storage racks. 51owever, CASMO-2E, as an assembly code, cannot directly represent an entire core configuration without introducino uncertainty due to reflector constants and the appropriateness of their spectral weighting. For this reason, the poisoned cell configurations of the central assembly, as calculated by CASMO-2E, were benchmarked against corresponding calculations with the 27-group (SCALE) AMPX-KENO V code package. Results of this comparison are shown in Table 2. Since the differences are well within the normal KENO statistical variation, these calculations confirm the validity of CASMO-2E calculations for the typical high density coisoned spent fuel rack configurations. The differences shown in Table 2 are also consistent with a bias of 0.0013 i 0.0018, determined in Section 3.1 as the expected bias and uncertainty of CASMO-2E calcula-tions.
- Yankee has attempted such calculations (Ref. 14) using CASMO-2E-generated constants in a two-dimensional, four-group PDQ model, obtaining a mean keff of 1.005 for 11 poisoned cases and 1.009 calculations for 5 unpoisoned cases. Thus, Yankee benchmark suggest that CASMO-2E tends to slightly overpredict reactivity.
O A-6
Table _2 RESULTS OF CASMO-2E.BENCHttARK (INTERCOMPARISON) CALCULATIONS k,III
-B&W Experiment No.Ill AMPX-KENO (2) CASMO-2E ok 1.1203
- 0.0032 1.1193 0.~0010 XIX 1.1149 i 0.0039 1.1129 0.0020 XVII 1.1059 i 0.0038 1.1052 0.0007 XV Interpolated (3) 1.1024
- 0.0042 1.1011 0.0013
.1.0983 i 0.0041 1.0979 0.0004 XIV 1.0992 i 0.0034 1.0979 0.0013 XIII i 0.0038 0.0011 Mean *0.0006 Uncertainty BWR fuel rack 0.9212
- 0.0027 0.9218 -0.006 (1) Infinite array of central assemblies of 9-assembly B&W criti-9).
(2) cal configuration (Ref.from AMPX-XENO corrected for bias of 0.0106 ak. (3)kInterpolated from Fig. 28 of Ref. 9 for soluble boron concen- , tration at critical condition. I i i I ( ! A-7 I L
REFERENCES TO APPENDIX A L]
- 1. Green, Lucious, Petrie, Ford, White, Wright, "PSR-63/AMPX-1 for Generating (code package), AMPX Modular Code System from ENDF/B,"
Coupled Multigroup Neutron-Gamma Libraries ORNL-TM-3706, Oak Ridge National Laboratory, March 1976.
- 2. L. M. Petrie and M. F. Cross, "KENO-IV, An Improved Monte Carlo Criticality Program," ORNL-4938, Oak Ridge National Laboratory, November 1975.
A Modular Code System for
- 3. R. M. Westfall et al., "SCALE:
Performing Standardized Computer Analyses for Licensing Evaluation," NUREG/CR-0200, 1979.
- 4. W. E. Ford, III et al., "A 218-Neutron Group Master Cross-section Library for Criticality Safety Studies," ORNL/TM-4, 1976.
- 5. A. Ahlin, M. Edenius, H. Haggblom, "CASMO - A Fuel Assembly Burnup Program," AE-RF-76-4158, Studsvik report (proprietary).
A. Ahlin and M. Edenius, "CASMO - A Fast Transport Theory 6. Depletion Code for LWR Analysis," ANS Transactions, Vol. 26,
- p. 604, 1977.
[}
- 7. M. Edenius et al., "CASMO Benchnark Report," Studsvik/RF-78/6293, Aktiebolaget Atomenergi, March 1978.
- 8. "CASMO-2E Nuclear Fuel Assembly Analysis, Application Users Manual," Rev. A, Control Data Corporation, 1982.
- 9. M. N. Baldwin et al., "Critical Experiments Supporting Close Proximity Water Storage of Power Reactor Fuel," BAW-1484-7, The Babcoch & Wilcox Company, July 1979.
- 10. M. G. Natrella, Experimental Statistics, National Bureau of Standards, dandbook 91, August 1963.
al., "Comparisons of Experiments and
- 11. B. F. Cooney et Calculations for LWR Storage Geometries," Westinghouse NES, ANS Transactions, Vol. 39, p. 531, November 1981.
R. M. Westfall and J. R. Knight, "Scale System Cross-section 12. Validation with Shipping-cask Critical Experiments," ANS' l Transactions, Vol. 33, p. 368, November 1979.
- 13. "The EPRI-CPM Data Library," ARMP Computer Code Manuals, Part II, Chapter 4, CCH3, Electric Power Research Institute, November 1975.
( A-8 i I
REFERENCES TO APPENDIX A (Continued)
- 14. E. E. Pilat, "Methods for the Analysis of Boiling Atomic Water ,
Reactors (Lattice Physics)," YAEC-1232, Yankee Electric Co., December 1980.
- 15. -L. E. Strawbridge and R. F. Barry, "Criticality Calculations for Uniform, Water-moderated Lattices," Nuclear Science and Engineerino, Vol. 23, o. 58, Septembe r 19 6 5.
- 16. S. E. Turner and M. K. Gurley, "Evaluation of AMPX-KENO Benchmark Calculations for Hiah Density Spent Fuel Storage Racks," Nuclear Science and Engineerina, 80(2): 230-237, February 1982.
1 0 h O , A-9 i e _..,....-.--,,.--__...,y,.,,
. - - . . , . . - , , . . _ - - - , _ . . . _ _ . - ,,.-,..._,_--,.,,,m., . . _ , _ , . , , _ , _ . . . - , _ . ,
.. . _ . - _ _ . - .- . .__ .. = - -_ _ _ . . _ . . . .
l 1 i APPENDIX B NRC QUESTIONS FROM DECEM3FR 22, 1986 MEETING { i 6 i i B-1 ;
i fs NRC Questions from December 2?, 1996 Meeting: U QUESTION #1 The information provided to BNL does'not address potential damage to fuel resulting from fuel to cell wall and baseplate impacts. Have these effects been evaluated? ANSWER #1 The fuel assembly has been analyzed by the fuel manufacturer and found to withstand lateral accelerations an order of magnitude greater than those developed due to rattling between the fuel assembly and the cell wall (less : than 1.0g). The testimony by a Westinghouse expert during the Turkey Point , ASLB hearings gives the details on the fuel assembly load bearing ability.
. The fuel assembly has a rugged inlet nozzle assembly which imparts considerable axial load bearing ability to it. Moreover, the forces developed due to bouncing of the assembly on the base plate are quite small compared to the loads associated with fuel drop accidents.
i 03EST10N #2 Has the pool floor evaluation included a check for potential local perforation . of the liner due to rack foot impact? How was this done? (]) ANSWER #2 Perforation of the stainless steel pool liner plate due to rack foot impact loading has been evaluated. A local bearing stress is calculated based on the maximum foot reaction. The net bearing area assumes a 9" x 9" shim plate distributing load over the leak chase channels. The critical bearing stress is 6275 which is less than the allowable concrete bearing stress of 6545 psi. The design compressive strength of the spent fuel pool basemat is 5500 ps! whereas the walls are 3500 psi. The latter value has previously been given in response to the second set of NRC questions. QUESTION #3 l Since the seismic impact loads may be sensitive to gap sizes, are there QA 3 procedures in place to assure that the proper inter-rack and wall to rack gaps are maintained during initial installing and subsequent loading /unloadinf of fue,? 1 ANSWER #3 The The rack to wall gaps will be maintained by temporary spacer blocks.
/~' nominal gap between racks is zero inches and will be checked to insure the actual gaps are below the 1/2 inches considered in the seismic analysis.
B-2
QA procedures are in the process of preparation and will be available to the NRC once they are approved. The procedures are designed to assure the assumed , i
-O rack to rack gaps are not exceeded and rack to wall gaps are maintained.
Since the load exerted by the fuel handling bridge crane is limited by the "load cell", the maximum load exerted on the rack during fuel loading and unloading is a small fraction of the weight of the module. Therefore, rack movement during fuel loading and unloading is not a credible phenomenon. 00EST10N #4 How was the conservatism of the single rack model demonstrated? The model appears to limit the amount of sliding and tilting of the rack between small gaps. This would not account for potential pileup of racks against the pool wall. Has this possibility been investigated? t ANSWER #4 The analysis of multi-rack motion has been carried out to determine the response of an array of four racks in the Byron spent fuel pool. The object of this analysis was to determine whether the peripheral rack module hits the wall. Another related purpose of this analysis was to obtain maximum values of rack kinematic displacements. The mathematical characterization of the modules in this analysis is not detailed enough to furnish refined data for The imposed seismic loading is the E-W SSE and the predicting stress levels.The four racks studied are labelled C2. Bl C4, and D4 in Fi: VERTICAL SSE. O 2.1. This array is the second row from the north wall. Realistic spring constant values for girdle bars are based on impact atAnthe ends of girdle bars with the full rack height providing flexibility. effective cross coupling gap was established based on our previous directions from the NRC for Diablo Canyon analyses. All racks were assumed to be full with consolidated fuel (weight - 3000#/ cell). This fuel weight is identical to the single rack design basis 3-D analyses. The table , Runs were made for .2 and .8 interf ace coefficients of friction. below provides results for comparison of key displacement data between the multi-rack and single rack analyses.The displacements increase over the single rack analysis but the units do not impact the mil. Friction 'tulti-Rack 3-D Single ^ Rack Coefficient Model 0.2 1.1027 inch .318 inch Maximum E-W Displacement .346 inch At top of rack 0.8 2.74 inch OUESTION #5 f ow has the conservatism of the rattling mass representation of the fuel The model does not appear to account within the rack cells been demonstrated? for the flexural rigidity of the fuel and the potential amplified response of O its resonant modes. Has this been investigated? 3-3
ANSWER #5 m The answer to this question is given in Section 6.2.1(a) of this report. (j QUESTION #6 Describe how i., - ek and rack to wall fluid coupling was handled by the model (note p.Il .' ,as missing from BNL copy of seismic report). ANSWER #6 The missing pige was provided at the meeting and the manner of handling fuel coupling was explained to be an extension of Fritz's classical paper. QUESTION #7 Uplif t of fuel from rack baseplate and resulting impact loads were not considered in the model . How was this effect evaluated? ANSWER #7 There is indeed no degree of freedom associated with uplif t of the fuel assembly mass located at the baseplate elevation. However, due to the low seismic accelerations involved, the resulting error is expected to be small . Furthermore, the fuel drop case considers drop of the fuel from 36 inches above the rack and is crictical compared to fuel uplif t due to seismic effects, p)
\
QUESTION #8 How were the rack to rack and rack to wall gap spring constants determined? How are the impact loads evaluateJ? ANSWER #8 The required rack to rack impact spring constant is at the girdleThe baractual and baseplate elevations where the inter-module impacts may occur. impact spring constant at the girdle bar location is estimated by treating The spring the problem of a laterailly loaded beam on an elastic foundation. constant vclue used in the dynamic analysis is, however, several times higher in order to maximize the computed impact forces. The impact loads in the spring elements are tracked at every time integration step. The maximum values are saved and printed out. A strength-of-materials type analysis is subsequently performed to est6blish that the locations subject to impacts will not fail as structural members. QUESTION #9 Have any studies been perforced to test the impact load sensitivity to gap size? ANSWER #9 Studies to test the impact load sensitivity to fuel as'.embly/ storage cell size were not performed. However, the fuel mas 2 in a rack is lumped together so it moves in phase resulting in an extremely anservative estimate of fuel impact effects. O i 3-4
QUESTION #10 ( ; How are the pool wall motions included in the analysis? Do the walls experience the same motion as the floor? ANSWER #10 Yes the walls experience the same motion as the floor since the walls are thick and rigid. QUESTION #11 1s the bouyant force of water considered in the calculation of frictional resistance of the rack feet. ANSWER #11 Yes, the bouyant force of water is included in the evaluation of the net vertical force on a support leg which is used it; the friction spring elements to compute the lateral friction force at each time step. QUESTION #12 Provide the How were the floor time histories used as input generated? corresponding response spectra. ANSWER #12 Time histories were generated based on the design basis spectra at the spent fuel pool floor level . The reponse sp'ctra generated frnm the time histories
.A match the design basis spectra closely (Ref section 6.1). Plots comparing
() the design basis spectra and spectra generated from synthetic time history are shown in Figures B-1, B-2, and ti-3. QUESTION #13 Clarify equation for P4 in Table 6.1 on page 6-28 of the licensing report. ANSWER #13 i
"4 are the "absolute" generalized coordinate,, qj are the associate P
attached to the pool floor and walls). QUESTION #14 to rack eccentricity l Explain the criteria used for the selection of fuel values (X33,YB) considered in the middle. ANSWER #14 The center of gravity of the assemblage of fuel assemblies in the horizontal I plane defines the coordinates Xe and Ym-(v3 3-5
LUESTION#15 Define the parameters given in the DYNAHIS computer run input and output. (] ANSWER #15 "Tnis was explained at the meeting. Some proprietary data was sent to BNL for their review. QUESTION #16 REck to rack and rack to wall gaps in the computer runs do not appear
~
consistent with gaps shown on the fuel pool layout drawing. E xplain. ANSWER #16 The rack to rar.k gaps for impact spring considerations should be looked up in the non-lincar spring tables. The marked output mentioned above explains the detail s , t l0 B-6
CALC N8. SARGENT 3 LUNDY pagJECT BYRON REY tieslutas PROJECT N8. 7500-03 8""P3"" 0 8d0 19 JUN 87 rnoc O SJH608 Fig. B-1 Response Spectra at Spent Fuel Pool Floor FREQUENCY IN CPS so.o 20.e 10.0 s.O t .0 10 05
' ' ' ' ' ' ' ' ' ' ' 4.0 20'- - - iii iisi i i iii 's e' si i i i i e 'i s i i iii iiis i ''''s i iii #i i iii iei 15 0 ~ -15 0 ~ . .~
10 0 10.C , e.0:
~ ~:.00 ~ ~
- .0c s.c: _
4.C
- -4 00 ~
f -DESIGN BASIS SPECTRA :
-s.00 3 0c
[ SPECTRA GENERATED FROM b a 2.C 3 l
/ [ SYNTHETIC TIME HIST 0hY -
2 00 E
/ : ~1 .50 ec15C~
E
/ n f E
1.00 E 1 - I h" r l 2 .GCE 5- 0.;; J%d' -
/ J / \ -
3 80 x_ 1 3 " g' - g isi i L __ e a ac . f f yf uvw y p; 3.50 u0 E0 ." g j / i[_ t s - D.40 4C g , 0 S4 ___ p/ g_3 30
~
3 3 20 0.C _ 3 16 0.1~ - . P4
'~ .
3 10 0 1C_ . 3 06 0.C
~
3 06 0.C0 ~ " 3 06 0.m 0.02 '.0S0 0.04 '0. 06 0.00 0 10 t, IE 0.20 0 30 0.40 0 80 0 80''.5' 1 0 1 20 PERIO IN SECON05 811SMIC RCCL SPECTRUM-55E b?ECTRA No. l NeoE ELE aTreN g I DIRECT!8N N-S AN04.E L8 CAT!8N SPENT FUEL P88L 0-7
CALC N8. SARGENT & LUNDY PR, JECT BYR8N REV DoDI E Dts - PR8JF.CT N8 7600-03 ORF'!NO 0 040 04 FEB 87 PROE O SJH604 Fig. B-2 Response Spectra at Spent Fuel Pool Floor FREQUENCY IN CP8 50 20 10 05 80.0 20 0 10.0 i i
,t, f ,f f , , ,1 ,i,, , f ,a i l ,1,1,7 i, ,1, ,i , i , ,t i,,5 :00 20.C _ ,i , , ,i , i, , , 1,1 1 1,i ,f ,I i, ,i i , ,1, , ,,,, i i : ~ ~
15 0 15.C _ 4 . - _ 10 0 10.C _ ::
- 00
- e.0C _ .
~ - :.00
- e.0C -
5 00 8.0C DESIGN BASIS SPECTRA 1 - f ~ 4 00
- 4 0C_ SPECTRA GENERATED FROM
- :
- [f SYNTHETIC TIME HISTORY
~ - s.00 3 0C: :
e . _ :.00
$;t.0C m
e 1.5C 2
/ / c$
5
- .50 0 E
_V__ s E
. 1 00 , t ..
E / / r,e L : 3 00 g
>- 0 0C s\,
E.0CE J0
/ / ,J,/ ._
E 3 00
/ s 3 50 - x-0 5C \ ~
0 4C _ rs
) ,- --/ _ 3 40 ' ~
3.30
^O 0 36 ./ D 3__ __ __
N 3.20 _ 0.tc _ :
~ ~
1 3 15 l 0.tC[ _
~ _
_ 3.10 l 0.tC _
~
3.00 1 0.C
~ ~ ~
3 06 l 0.C0 ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' '
' ' ' ' ' 3 05
! 0 05 0.15 0 20' ' ' ' .0 30 0 40 0 80 0 00 1 0 15 20 0 02' ' ' ' .03 0 0 04 '0.06 000010 PER180 IN $EC8N08 SPECTRR N8. SEISMIC RCCL 8PECTRUM-88E l ELEVRTION N00E ANGLE L8CRT!8N SPENT FUEL P88L 2 DIRECT!8N E-N B-3
l SARGENT & LUNDY EECI syren REv DestMttAs . PR8 JECT N6 7500-03 onnrlNo 0 040 04 FEB 87 rnoE g SJH604 l Fin. B-3 Response Spectra at Spent Fuel Pool Floor FREQUENCY IN CP8 50 t.0 10 15 so.0 20 0 10 0 iiiiiinni i i ' '
, iiiiiiiii i i e e > > > =0 0 to C _ i ,i i, i
i 'ii ' ' ' ' ii,. i i iii iii iiii siii iiii iiii y
. ~
t5 0 15 0 : , _ t0 0 10.C . -
-:.00 0 0C _ _ - -:.00 e.0c - :.00 5 0c -
_:.00 4 0C . .
, -DESIGN BASIS SPECTRA d,, -
3.0C ~ ~ SPECTRA GENERATED FROM e -
/ SYNTHETIC TIME HISTORY _~o ,
C t.0C . z -
/ 21.50 = :
e 1.5C ,
/ _
A - j -
- 1.00 g r\ /
W e
>- 0 .*9C- E # w / E 3..0 5 f()W .) -
M-M - - 3 00
- V 43 d 0 0C ,, 3 50 ~
S ((
- / . ~ 3 40 0 4C h CN
__ . / 3. 0 0 3C_ \ N -
%v'.l\sN '
1
- - 3 20 l 0.tc . ,
D t
~
3.15 0 15 \f m N
- _ 3.10 0 1C . " - 3 00 0.C:_ - 3 06 0.C0 - iI ii liII I I I I IIII I I I I I III iI I tftt QM g g g g ggai l l t i I tii 0.15 0.20 0 30 0 40 0 80 0 80 1 0 15 20 0.f.02 0.03 0 04 0 00 0.00 0.10 PERIOD IN SEC8HDS srectaa No.
O sEisaic ace' secciaun-sse N6DE ELEVRT!8N ANOLE L8 CAT 18N 2 PENT FUEL P86L '1 DIRECT!9N VERT B-9
.-}}