ML19327C155
ML19327C155 | |
Person / Time | |
---|---|
Site: | University of Virginia |
Issue date: | 07/31/1989 |
From: | Freeman D VIRGINIA, UNIV. OF, CHARLOTTESVILLE, VA |
To: | |
Shared Package | |
ML19327C150 | List: |
References | |
NUDOCS 8911210027 | |
Download: ML19327C155 (107) | |
Text
V,
/
4 a
1 1
NEUTRONIC ANALYSIS FOR THE WAR REACTOR HEU TO LEU COM ER$10N PROJECT A Thesis Presented to the Faculty of the School Engineering and Applied Science University of Virginia l
l l In Partial Fulfillment of the Requirements for the Degree Master of Science by l David W. Freeman July, 1989
, ACKNOULEDGCIIENT jy This work was performed with support fron l Department of Energy Grants li
- DEFG0508EP.75300
'f l
in and 3 ft M
- DEFG0586ER75273 i}
3 0 5 39E5 E m m m er e x m m " ~ f" --- U l
l
I !
I i i
ABSTRACT ,
The University of Virginia Research Reactor is currently in the process of converting from high enriched uranium (HEU) fuel to
{
low enriched uranium (LEU) fuel. The conversion is in response to a mandate from the Nuclear Regulatory Commission. This thesis i addresses several neutronics concerns associated with the conversion project including: 1) control rod worths; 2) radial I peaking factors; 3) temperature feedback effects; 4) moderator void feedback effects; 5) reflector worth effects; 6) delayed neutron fractions; and 7) prompt neutron lifetimes. {
i Control rod reactivity worths for LEU fueled cores were found to be essentially the same as for the HEU fueled cores. Radial !
peaking factors for the LEU fuel were found to be slightly higher f than for HEU fuel due to the harder flux spectrum associated with the LEU fuel. Temperature and moderator void reactivity feedback ;
effects for both the LEU and HEU fueled cores were found to be very similar except for the doppler effect which was significantly I higher in the LEU fuel due to the higher uranium 238 loading. !
Reflector worths for several different water and graphite )
I reflector options were quantified for the LEU fueled cores, and f heavy water reflector tanks were analyzed. An effective delayed neutron fraction of 0.0074 was determined for both the HEU and LEU fueled cores. A prompt neutron lifetime of 67 micro seconds was determined for the LEU core which was about 15 percent lower than 1
y-I I 11 that determined for the HEU core.
The overall thesis results show that the LEU-fueled core will j behave very similarly to the HEU core currently being used.
I !
I l l
I ;
h
! i L
I 1 LI !
I ;
f I >
I l I
Lg
'I
.~. _ _ .
iii ACKNOWLEDGEMENTS
, i This report is dedicated to my wife, Donna, and my daughters. -
Sarah and Jodi, who constantly supported me throughout the l
project. !
Thanks go to Dr. Rydin for providing guidance and keeping the project on track. A special thanks goes to Vicki Thomas who typed this thesis and helped with the format, j Thanks go to Press Farrar who provided historical background and documentation for the reactor, f i
Finally, thanks are owed to Mary Fehr and Stuart Wasserman for -
laying the strong foundation of computer models which ! inherited.
I :
I I
o a
Li !
i
, . ~ . _ _ - - - . . , . . - . , , . _, .g. . .r...- . _ . _ , , _ . . . . . . , __, _ -,.. , ,, . . . .. __ __ ..- ,,,_
I l
I '
l l iv I l
l TABLE OF CONTENTS l
)
13 l hat ' i z; Abstract. . . . . . . . . . . . . . . . . . . . . . . . . . . 1 Acknowledgements. . . . . . . . . . . . . . . . . . . . . . . iii List of Figures ...... . . . . . . . . . . . . . . . . vii Li s t of Tables . . . . . . . . . . . . . . . . . . . . . . . .viii
1.0 INTRODUCTION
. . . . . . . . . . . . . . . . . . . . . . 1 l 4
- l.
2.0 BACKGROUND
INFORMATION . ............... 5 t
j l
2.1 Description of UVAR . .............. 5 !
2.2 Control Rods. .................. 8 l l
2.3 HEU Fuel. . . . . . . . . . . . . . . . . . . . . . 9 l E 2.4 LEU Fuel . . . . . . . . . . . . . . . . . . . . . 11 !
'E u 16 l
i 3.0 NETHODOLOGY. . . . . . . . . . . . . . . . . . . . . . .
3.1 Methodology for Developing Neutron Cross Sections . 17 ;
3.1.1 Core Region Cross Sections. . . . . . . . . 17 j
1 3.1.2 Control Rod Cross Sections . . . . . . . . . 19 .
l i
!. 3.1.3 Reflector Cross Sections . . . . . . . . . . 21
- l. .
i .
3.2 Diffusion Theory Modeling of the UVAR Core. . . . . 21 :
1 l' >
l 4.0 CONTROL ROD WORTH STUDY. . . . . . . . . . . . . . . . . 26 4.1 4 by-4 Cores and 4-by-5 Cores. . . . . . . . . . . 26 4.2 Core-Life Effects on Control Rod Worths . . . . . . 29 4.3 Reflector Effects on Control Rod Worths . . . . . . 30 4.4 Control Rod Position Effects. . . . . . . . . . . . 32 4.5 Control Rod Worth Study Summary . . . . . . . . . . 34 I
I ,
I k.I 5.0 RADI AL PEAKING FACYed . . . . . . . . . . . . . . . . . 36 5.1 Radial Peakir.g Factors for UVAR 4-by-4 and 38 4-by-5 Cores . . . . . . . . . . . . . . . . . . .
. 5.2 Control Rod Element Rotation Effects on Radial Peaking Factors . . . . . . . . . . . . . . . . . . 40 (g. 5.3 Control. Rod Insertion Effects on Radial
[g Peaking Factors . . . . . . . . . . . . . . . . . . 43 5.4 Core-Life Effects on Radial Peaking Factors . . . . 44 i
4 5.5 Core-Reflector Effects on Radial Peaking Factors . 47 5.6 Radial Peaking Fact .r Summary . . . . . . . . . . . 48 6.0 TEMPERATURE'AND VOID EFFECTS STUDY . . . . . . . . . . . 49 6.1 Tempercture Effects . . . . . . . . . . . . . . . . 49 6.1.1 Doppler Effects. . . . . . . . . . . . . . . 50 6.1.2 Moderator Temperature Effects. . . . . . . . 54 6.1.3 Power Defects. . . . . . . . . . . . . . . . 54 6.1.4 Gross bystem Heating Effects . . . . . . . . 55 6.2 MODERATOR VOID EFFECTS . . . . . . . . . . . . . . 55 6.2.1 Uniform Core Void Effects. . . . . . . . . . 55 6.2.2 Local Core Void Effects. . . . . . . . . . . 59 6.3 Temperature and Modcrator Effects Summary . . . . . 62 I- 7.0 REFLECTOR WORTH STUDY. . . . . . . . . . . . . . . . . . 63 7.1 Graphite and Water Reflector Options. . . . . . . . 63 7.2 Deuterium Tank Reflectors . . . . . . . . . . . . . 69 7.3 Black Boundary Reflectors . . . . . . . . . . . . . 72 7.4 Reflector Worth Study Summary. . . . . . . . . . . 72
l i
l i
! vi l
I
~
8.0 DELAYED NEUTRON AND PROMPT NEUTRON LIFETIME STUDY. . . 74 -
8.1 Methodology . . . . . . . . . . . . . . . . . . . 74 I 8.2 Model Verification. . . . . . . . . . . . . . . . 75 8.3 Prompt Neutron Lifetimes. . . . . . . . . . . . . 76 8.4 Effective Delayed Neutron Fractions . . . . . . . 79 8.5 Summary of Delayed Neutron and Prompt Neutron Lifetime Study . . . . . . . . . . . . . . . . . . 83 l
9.0
SUMMARY
. . . . . . . . . . . . . . . . . . . . . . . . 84 REFERENCES. . . . . . . . . . . . . . . . . . . . . . . . . 87 APPENDIX A MATERIAL CROSS SECTIONS FOR UVAR 2-D CORE MODEL ,
I i
r LI !
LI i I.
n I .
I l
3 I
vii ,
LIST OF FIGURES em Figure 2-1. UVAR Grid Pl ate . . . . . . . . . . . . . . . . 7
. Figure 2-2. HEU-18 Fuel and Control. Rod Elements. . . . . . 10 Figure 2-3. LEU-22 UVAR Fuel Element and Control Rod Element Designs . . . . . . . . . . . . . . . . . . . . 14 Figure 2-4. LEU-18 UVAR Fuel Element and Control Rod Element ;
Designs . ... . . . . . . . . . . . . . . . ... 15 l
Figure 3-1. UVAR 4-by-4 Core Model... . . . . . . . . . . . 24 >
Figure 3-2. UVAR 4-by-5 Core Model. . . . . . . . . . . . . 25 Figure 4-1. Alternate Control Rod Locations for UVAR 4-by-5 Core. . . . . . . . . . . . . . . . . . 33 Figure 6-1. Void Coefficient as a Function of Void Percent. 58 Figure 7-1. Reflector Options 1 through 6 . . . . . . . . . 65
-g' n' Figure 7-2. Reflector Options 7 through 13. . . . . . . . . 66 I
- I 3
I l B
I 1
9' ?
viii LIST OF TABLES Table 4.1. Control Rod Worths for the UVAR 4-by-4 and l 4-by-5 Cores. . . . . ... . . . . . . . . . .
27 Table 4.2. arison of Experimentally Determined and Combulated Cal Rod Worths for the 4-by-4 HEU 18 Core. . . . . . . . . . . . . . . . . . 28 Table 4.3. Core-Life Effects On Control Rods. . . . . . . 29 Tabl e ' 4. 4. Reflector Effects on Control Rod Worths. . . . 31 Table 4.5. Control Rod Worths for Alternate Rod Locations in a UVAR 4-by-5 Core . . ........... 32 Table 5.1. Radial Peaking Factors for UVAR 4-by-4 and 4-by-5 Core Configurations with Various I Table 5.2- .
Fuel Options . . . . . . . . . . . . . . . . .
Radial Peaking Factors for Control Rod 39 Orientation Scenarios. .. . . . . . . . . . . . 42 Table 5.3. Radial Peaking Factors for LEU-22 4-by-4 Graphite Reflected Core with Various Rod Insertions . . . . . . . . . . . . . . . . . . 44 Table 5.4. Core-Life Effects'on Radial Peaking Factors for the LEU-22 4-by-5 Graphite Reflected Core. 46 Table 5.5. Core Reflector Effects on Radial Peaking Factors . . . . . . . . . . . . . . . . . . . 47 p Table 6.1. Doppler Effects for the UVAR Reactor . . . . . 51 Table 6.2. Moderator Temperature Effects for the UVAR Reactor . . . . . . . . . . . . . . . . . . . 53 Table 6.3. Power Defects For 4-by-4 UVAR Cores. . . . . . 54 Table 6.4. Reactivity Effects Associated with Uniform Moderator Void Conditions With LEU-22, LEU-18, and HE'J-18 Fuel Loadings . . . . . . . . . . . 56 i
i k
^
%r
,i ix Table 6.5. Uniform Moderator Void Coefficients for LEU-22, LEU 18, and HEU-18 Fuels . . . . . . . . . . 57 Table 6.6. Local Void Reactivity Effects for Void Scenarios 1 Through 4. . . . . . . . . . . . 60 l
Table 6.7. Local Moderator Void Coefficients. . . . . . . 61 Table 7.1. Reactivity Effects Associated With Variations of a Full Graphite Reflected 4 by-5 Core. .. 67 Table 7.2. Reactivity Effects Associated with Variations of a Single Graphite Raf1ector Row Surrounding the 4-by-5 Core . . . . ... .......68 Table 7.3. Core-Life Effects on Relative Worth of Reflector Option 2 Using LEU-18 Fuel. . . . . 6f Table 7.4. Representative2 D 0 Tanks . . . . . . . . . . . 70 Table 7.5. Reactivity Effects Associated with Different-Sized D20 Tanks Bordering the West UVAR 1 Core Face. . . . . . . . . . . . . . . . . . . 71 Table 8.1. Energy Group Structure Used for the EXTERMINATOR Analysis. . . . . . . . . . . . . 75 Table 8.2. Delayed Neutron Fractions . . . . . . . . . . 75 Table 8.3. Comparison of k-effective Values Calculated by the EXTERMINATOR and 2DB codes. . . . . . . 76 Table 8.4. Prompt Neutron Lifetimes for the UVAR Core. . 70 Table 8.5. Total Effective Delayed Neutron Fractions for I LEU-22 Fueled Core . . . . . . . . . . . . . . 81 APPENDIX A Table A.1 Core Region Cross Sections for LEU-22 Fuel. . 91 Table A.2 Core Region Cross Sections for LEU-18 Fuel . 92 Table A.3 Core Region Cross Sections for HEU-18 Fuel. . 93 Table A.4 Control Rod Region Cross Sections. . . . . . . 94
l~
1.
l l
I x ,
Table ' A 5 ANL Developed Water and Graphite Reflector Cross Sections. . , . . . . . . . . . . . . . 95 1
Table A.6 UVA Developed Graphite Reflector Cross Sections. . . . . . . . . . . . . . . . . . . 96 3 <
l L5 .
1l. -
I i I ;
l.
I I .
-I I
I 5
I
.3
4
, 1 CHAPTER 1. INTRODUCTION This report presents severil studies designed to address l neutronics issues and concernv v:ociated with converting the University Of Virginia Research Reactor (UVAR) from highly-enriched uranium-(HEU) fuel to low-enriched uraniun (LEU) fuel.
The neutronics studies were performed using advanced neutronics.
- computer modeling liechniques. Studies presented in this report address the following topics: 1) control rod worths; 2) radial peaking factors; 3)' temperature and moderator feedback effects;
- 4) reflector worths; and 5) delayed-neutron and prompt neutron I lifetimes. Background information used to build the computer models and study methodologies are also addressed.
Several research reactor facilities in the United States use HEU fuel. The Nuclear Regulatory Commission (NRC) has directed all research facilities using HEU fuel to convert to LEU fuel.
HEU fuels are often enriched in excess of 90 percent uranium-235 while LEU fuels have an enrichment limit of 20 percent. The order
.to convert to LEU fuel was in response to security and safeguards issues raised over the HEU fuel.
Currently, the UVAR is in the process of being converted from HEU to LEU fuel. The conversion process is being performed by the UVAR staff and is non-trivial. The process includes: 1) a revision of the Safety Analysis Report; 2) Technical
)
1 ..
2 Specifications changes; and 3) other, as yet undefined, analyses i
deemed necessary by NRC.
To date, much work has been performed on the project by both UVAR staff and graduate students. This work is well summarized in the facility Progress Report (1) to the Department of Energy (DOE). The work has focused on the operating performance of the LEU fuel as compared to the currently used HEU fuel.
LEU fuel elements considered for use in the UVAR will use generic. LEU fuel plates described in Chapter 2. LEU fuel elements n using 18 plates (LEU-18), 20 plates (LEU-20), and 22 plates (LEU-
- 22) have been considered as potential replacements for the current 18 plate HEU fuel elements (HEU-18). The neutronics analyses have not-addressed LEU-20 fuel elements and have assumed that their characteristics were bounded by the LEU-18 and LEU-22 fuel element p
characteristics, in the fall of 1988, The University of Virginia (UVA) Reactor Safety Committee tentatively approved the use of either the LEU-18, LEU-20, or LEU-22 fuel elements. Work performed by Fehr (2) showed that the LEU-22 fuel element would be a superior replacement for the current fuel in terms of core lifetime. Based on this information, the UVAR staff decided to choose LEU-22 fuel as the preferred replacement element. For this reason, several of the studies presented in this analysis address only the LEU-22 I
fuel option.
I i
s I' 3 A key issue of concern was whether or not the UVAR control i
rods would provide an adequate reactivity insertion in I.EU cores for safe operation and whether their associated reactivity-effects would be significantly altered from those in HEU cores. This issue is addressed in Chapter 4 of this thesis and in Reference (3].
The focus of the HEU to LEU' conversion safety analysis is the thermal hydraulics analysis (4]. This effort is supported by information presented in the peaking factor study presented in i- Chapter 5.
Chapter 6 presents a study which addresses temperature and moderator feedback effects for the LEU fuel, and Chapter 7 addresses reflector worth behavior for the various fuel options.
NRC informally requested,. via a guidance handout at the 1988 I
g National Organization of Test, Research, and Training Reactors (TRTR) meeting, that delayed neutron fractions and prompt neutron lifetimes be calculated as part of the license renewal process.
The study presented in Chapter 8 addresses this issue.
It should be noted that no simple mechanism exists for i quantifying errors associated with the computer modeling techniques used throughout this thesis. Errors may be incurred in the modeling process, computer computations, and input data.- A default error value of plus or minus 5 percent has been chosen as I
representing what is believed to be a reasonabla estimate of the 1
e
L ,
,a
/ ,'
I 4 error involved with values presented in this thesis. Higher error estimates are quoted in this thesis when deemed appropriate. As discussed in Chapter 9, the computer models can be " fine-tuned" when experimental data from the new LEU core becomes available.
This' process should eventually provide information allowing a better estimate of computational errors to be made.
The current schedule for the conversion project includes the submission of a revised Safety Analysis Report (in which the author will be heavily involved) in August,1933. LEU fuel i element fabrication is projected for completion by the end of 1989. Assuming that NRC licensing concerns can be addressed ~in a timely fashion, .it should be possible to convert the UVAR to LEU fuel in.early 1990. -
p l'
s i
u L
5
2.0 BACKGROUND
INFORNATION This chapter presents background information which forms the basis for the UVAR neutronics models developed for analyses presented in this report. Presented below are descriptions of:
- 1) the UVAR reactor; 2) control rods; 3) HEU fuel; and 4) LEU fuel, u
2.1 Description of the UVAR Reactor r The UVAR reactor is a pool type research reactor constructed in 1959. Currently. licensed to operate at 2 megawatts, the UVAR i is used primarily as a research tool. Research activities, such as neutron: activation analysis and neutron radiography, use the reactor core as a powerful neutron source. The reactor assembly is comprised of fuel elements, control rod fuel elements, control rods,. and graphite reflector elements; all of which sit in the i reactor grid plate. The grid plate assembly is supported by an aluminum framework which is suspended from a movable bridge. Pool water (light water) serves as moderator, coolant, reflector, and radiation shielding.
1 The reactor sits under about 20 feet of water at the bottom L of a 75,000 gallon pool. The reactor is cooled by pool water
'downflow through the core elements in excess of 1000 gallons per minute. This water passes through a tube-in-shell heat exchanger and is returned to the pool. Heat is transferred to the atmosphore via a cooling tower mounted on the facility roof.
i
)
6 A full description of the UVAR is presented in the facility b Safety' Analysis Report (5). Certain key features of the UVAR, such as the grid plate and core components, are discussed below.
Grid Plate The UVAR reactor'is located on the reactor grid plate. The reactor is " built" by loading various core components into the grid plate to achieve a desired core configuration. Core configurations are widely varied and are limited only by shutdown margin and excess reactivity requirements.
m i The UVAR grid plate, shown in Figure 2-1, contains an 8-by-8 array of holes, approximately 21/2 inches in diameter, for positioning reactor components. For the purposes of this analysis, the grid plate has been subdivided into 64 equal grid cells measuring 3.189 inches by 3.031 inches (6). The grid cells j are depicted in Figure 2-1 with dotted lines.
Reactor Core Comoonents Reactor components include fuel elements, control rod fuel elements, graphite elements, and grid plate plugs. Reactor components are designed to fit into the grid plate holes and F extend no further than the grid cell boundaries.
Fuel elements, control rod elements, and graphite elements all have similar dimensions of about 3 inches in length and width, and 35 inches in height. The bottom of each element consists of a I
cylindrical tapered nozzle which fits snugly into the grid plate I
( g r
i I
,7 I
g-O!O!O!O!OlO!O!O O!OiO!O!O!O!OlO
! .O!O101010101010
..+..+ .+..+..+...+...+...
t oioioioio!O!O!o
_UlDidiUibidiUiU' OiO!OiOiO!O!O!O O!OiOiO!OiOiO!O UlUidiUidiUidlU Figure 2-1. UVAR Grid Plate
+'
8 holes . Graphite elements, used as reflectors, are generally placed on the outside faces of the fuel that makes up the core to
' enhance core reactivity.
f Grid plugs are used to prevent water flow through empty grid locations. UVAR Standard Operating Procedures (SOPS) require all unused crid l'ocations to be plugged in order to minimize non-element by-pass flow. A grid plug is a short. metal cylinder, approximately 3. inches in diameter, mounted on a tapered nozzle.
When inserted in the grid plate, the plug extends only a few-inches above the grid plate, which is below the active fuel region.
2.2 Control Rods The UVAR uses four control rods; 3 safety rods cr.d 1 regulating rod. Safety rod reactivity worths range from 3 to 5% ap and are used to scram the reactor. The regulating roo has a relatively low worth (0.3-1% ap) and is used primarily to compensate for small changes in reactivity. The control rods are oval and fit into the center hole of the coni.rol rod element.
Each type of rod is discussed below.
Safety Rods The UVAR safety rods are made of boron-statuless steel, clad in aluminum, and are about 1.5% natural boron by volume. The rods are about 25 inches long and have an oval cross section geometry with approximate dimensions of 2 1/4 by 7/8 inches. The rods are
T a
9
. vertically grooved to increase surface area. The safety rods are magnetically coupled to their drive mechanisms and drop by gravity into the core on a scram signal.
~ Reoulatina Rod The regulating rod, commonly referred to as the " reg rod", is made of stainless steel, clad in aluminum. It has the same dimensions as the safety rods but is not grooved. The rod is permanently attached to its drive mechanism and, therefore, does not drop on a scram signal. The reg rod has a relatively low worth and is generally located in the outer portions of the core, with the primary purpose of compensating for small reactivity changes associated with normal operations.
2.3 HEU Fuel The UVAR is currently fueled with highly enriched uranium (HEU). fuel elements, shown in Figure 2-2. Each fuel element contains 18 curved fuel plates and is referred to as "HEU-18" fuel throughout this. analysis. Each element initially contained 195 grams of U-235; however, the HEU fuel used in this analysis was partially burned, having an estimated 192.3 grams of U-235 per element. The element has approximate cross sectional dimensions of 3-by-3 inches, with an active fuel length of about 24 inches.
HEU control rod fuel elements are similar to the regular fuel elements except that the center 9 fuel plates have been removed to allow room for the control rod to move up and dowr.. The HEU control rod element is shown in Figure 2-2.
l
{I i l[h 10 I
r ! !r h -
i ,
9 J-
.. L_. ,
I E \ .
c// / m s < .
J -
N s x <
ll s a : . s s :
- , N s
- : : ;
. l..
l : : : s s l l l
- : : s s : : : '- -
. l. . s s -
s s l :
- l. ,
>///v, s Figure 2-2 HEU-13 Fuel and Control Rod Elements (5].
8-
g, ._.
D. 3 l'
l 11 Important parameters. associated with the HEU-18 fuel olements are. presented below:
I Parameter VAlyg
- 1. HEU Fuel Meat a'. width 0.02 inch
- b. length 2.375 inch
- c. height 23.5 inch
- 2. Clad Thickness 0.015 inch
- 3. ' Plate Thickness 0.05 inch
- 4. Water Gap 0.122 inch
- 5. U-235/ plate 10.7 gram
- 6. U-235/elenient 192.3 gram 2.4 LEU Fuel Generic LEU fuel plates have been designed by EG&G Idaho (in cooperation with DOE) and fabricated by Babcock and Wilcox in Lynchburg, Virginia. The plate is.made oversized so that it can be " trimmed" to fit different fuel element designs.
Detailed requirements for LEU fuel meat and fuel plates have J been published in TRTR specifications [7 - II). This information har, been reviewed for applicability to the UVAR neutronics work effort (12] and is incorporated in the LEU neutronics analyses ,
presented in this report.
3 L ,
12 I'
-The LEU fuel plate and fuel meat will have dimensions almost identical to the HEU fuel. The LEU fuel meat is 19,75 percent enriched uranium silicide powder (U3Si2) dispersed in aluminum powder which is " hot-rolled" to form the fuel core matrix, commonly referred to as the fuel meat. The fuel meat is bonded with, and encased in, an aluminum clad to form the fuel plate.
Specifications for the Generic LEU fuel plate are presented below[13):
PARAMETER yflyE
- 1. Type of Fuel U3 Si2
- 2. Fuel Meat Density 3.47 gram /cm3
- 3. U-235 Per Plate 12.5 gram
- 4. Plate Thickness 0.050 inch
- 5. Plate Width 2.775 inc.h
- 6. Clad Thickness 0.015 inch
- 7. Fuel Meat Thickness 0.020 inch
- 8. Fuel Meat Width 2.395 inch
- 9. Fuel Meat Length 23.25 inch Generic LEU fuel element designs de not exist because facility-specific elements are generally required. In response to i
the need for an element design to form a consistent basis for the
> thermal-hydraulic and neucronics analyses, preliminary UVAR element designs ware developed. These designs are presented in
5
- g .
I 13 references [13) and (14), and have been distributed to LEU Project Principals at both UVA and EG&G Idaho.
f The UVAR LEU elements use the generic LEU fuel plate described earlier. The elements will,use flat plate fuel and will >
have the same outer dimensions as previously-used UVAR 12-plate ,
elements. The LEU element designs are based on blueprints for previous HEU elerrients [15,16). The element designs currently address only two dimensions. The axial dimension in both the neutronics. and thermal-hydraulics analyses is- accounted for by an "
axial dimension function. .
Figures 2-3 and 2-4 present the UVAR LEU fuel element and
, control rod element designs for 22-plate elements (LEU-22) ard-18 plate elements (LEU-18). Applicable dimensions for the LEU fuel and control rod elements are presented in the figures.
l L
I s
I '
Lg l
LI og I
p I
--., 14 -
E
_ .s
+
l }
t.9968*
d 1 I i . 4 ;-
'E l
d _
. B.1575* -
8.87F8*
\
. l 2.9960*
l ~ .
E a
I I m l 'm R l l _I
,,,,,,. 2.s21e-e.arre*
Figure 2-3. 1.EU-22 UVAR Fuel Element and Control Rod Fuel Element Designs. (Figures by B. Hosticka, April 1988).
___ m______._.. ._._ _ _ . _ _ _ _ _
~ ^
- p. fj
- k. i (p' .i l'
15' .
g m
[
2.9950' I' l-m -
i j
i 1 i I
- l g l l' .
. L _
z I r.sais-m B.1975' j I .
\
I- a .ene'
[
L3 i
2.9960' ,
~
l 7 -
l- in ,
ca m
} 0 g ---
m' I a
~. 2.6210' r
- m D.1075*
B.8770' I
Figure 2-4. LEU-18 UVAR Fuel Element and Control Rod Element Designs (Figures by B. Hosticka, April, 1988).
L8
,1 ,
J l_
J 16 3.0 METHODOLOGY i Studies presented in this report were performed using several combinutions of various main-frame computer codes. The basic l methodologies supporting the studies are presented here. Study specific methodologies are presented, as appropriate, in the sections addressing the studies. Computer codes used in the analyses included: 1) LEOPARD; 2) THERM 0S; 3) GAMTEC; 4)
. EXTERMINATOR; and 5) 20B.
~The LEOPARD [17), THERMOS [18), and GAMTEC[19) codes are neutron interaction cross section generator codes and were used to develop energy-dependent cross sections for the analyses. Additionally, i,
the EXTERMINATOR [20] 2 D diffusion theory code was used to help develop control rod cross sections. .The methodologies for cros:.
section development are presented in Section 3.1, 2-D diffusion theory codes, 2DB and EXTERMINATOR, were used to develop detailed UVAR core models. The 2DB core model, as modified by the University of Michigan (referred to as 20B-VM),
was used for the bulk of the analyses. The EXTERMINATOR code was used to model the UVAR core for the purposes of calculating delayed neutron fractions and prompc neutron lifetimes. The 2-D diffusion theory codes are described in Section 3.2.
~ . _ -
p 17 v
3.1 Methodology For Developing Neutron Cross Sections l Cross sections used in this analysis can be classified as
- 1. Core-region cross sections
- 2. Control rod cross sections
\
- 3. Reflector cross sections A different methodology was, employed to develop each type of l
cross section.. These methodologies are discussed below.
3.1.1 Core-Region Cross Sections The LEOPARD computer code was used to develop multi-group cross sections for the UVAR core region. Variations in LEOPARD input parameters were employed to develop temperature and void dependent cross sections used in the temperature and moderator void study presented in Chapter 6. Additionally, the code calculates cross sections associated with fuel depletion at discrete burnup steps.
The LEOPARD code computes two and four group spectrum weighted macroscopic cross sections for an infinitely repeating unit fuel cell in slab geometry. The code allows both lattice and
., non-lattice portions of the unit fuel cell to be defined. The lattice region is defined as the fuel / clad / moderator repeating array. The non-lattice region contains any materials, such as structure and water spaces, which are not in the lattice. The
I 15 .
l 18 UVAR fuel cell model defines' the fuel meat, clad, and water gap array of the ' element to be the lattice region. The non-lattice
~
i region contains the element side plates, non-fuel bearing fuel plate edges, control rod guide plates, and control rod water i holes. .The fuel and control rod elements are slightly smaller than the grid cell dimensions shown in Figure 2-1, thus when P- loaded in the gridplate, a small amount of water exists between !
the edges of adjacent elements in all directions. This water is considered non-lattice, e
The LEOPARD models " distort" the actual fuel and control rod .
element dimensions somewhat while conserving materials and surface ;
area. The distortion consists of " compressing" the cell in the y [
direction and correspondingly " stretching" the cell in the x direction. This distortion was performed so that the lattice
. region of the cell extends completely across the grid cell in the
- l. x direction, thus eliminating the small amount of water that ltl L' . up .
actually lies between each element in the .x direction. This I distortion technique, developed by Fehr (2), allows for an uninterrupted fuel lattice region when several elements are loaded together in the 2-D core models.
L Burn-up cross sections account for fuel depletion, xenon and samarium buildup, and other fission products.
Two and four group cross sections calculated for the fuel lattice, non-lattice and control rod non-lattice regions are l~
l LI L - - - - - , .-
l l
B 19 presented in Appendix A.
3.1.2 Control Rod Cross Sections l Two group control rod cross sections were developed using a combination of various features provided in the THERMOS, GAMTEC, and EXTERMINATOR codes. The method of developing control rod cross sections is quite involved and time consuming to perform. "
Several. different modeling techniques were investigated by Wasserman (22]. The Wasserman study yielded methodologies for control rod modeling that produced results consistent with experimental data. The Wasserman study specifically addressed I 1
L control rods in the UVAR HEU-18 fueled core. An analogous -
l' methodology was developed to calculate control rod cross sections
~
for the LEU fuel. Control rod cross sections were used to calculate control rod worths (Section 4). A brief description of !
the methodology is presented here; however, the reader is referred I
to Wasserman (and Reference (3]) for a detailed presentation. j l- ,
Control rods are difficult to model using diffusion theory computer codes because the control rod is a strong absorber and li creates a steep local flux gradient. Diffusion theory cannot !
. adequately handle the steep gradient, thus transport theory is required. The methodology developed to model the control rod L-l' region consists of calculating cross sections for that region <
which, when used in a full 2-D diffusion core model, produce the correct results. These cross sections are called " effective" L
- I
l 5, ,
20 )
diffusion theory cross sections because they help the diffusion theory code " effectively" handle the local rod transport problem.
The effective cross sections were developed by first
- . calculating transport ~ thermal microscopic cross sections for the.
- control rod region. A 1-D THERMOS slab transport model of a control rod fuel element with control rod material in the central region was developed. The THERMOS model calculated region-smeared microscopic cross sections and region-dependent neutron absorption- i
. fractions. >
Fast' group microscopic cross sections were developed with a cylindrical GAMTEC model control rod fuel element. GAMTEC' I calculates fast group cross sections using a P-1 approximation.
.GAMTEC also calculates. resonance absorption fractions for the cell s model.
The thermal and fast microscopic cross sections were input into a 2-D EXTERMINATOR model of a rodded control rod element using x-y geometry. The area of the control rod region in the EXTERMINATOR model corresponds to the area used for the control rod water hole region of the 2DB models.
. The thermal absorption cross sections were manually adjusted in the EXTERMINATOR model until the thermal absorption fractions LW in the fuel and control rod regions matched those predicted by the '
THERMOS model. The fast absorption cross sections in the EXTERMINATOR model were also manually adjusted such that the L
.l I .
..m__... _ . _ . _ , _ , . . . _ . . . . . . , . _
2 ,
t 21 fraction of neutrons reaching thermal energies matched the resonance escape probability predicted by the GAMTEC P-1 approximation.
The methodology described above was applied to each fuel type 6
(i.e.,HEU-18, LEU-18, LEU-22). Resulting macroscopic cross sections are presented in Appendix A. !
3.1.3 Reflector Cross Sections
- Graphite and water are the predominant reflector materials for
, the UVAR reactor. Other reflector materials are addressed in ,
Section 7.
Two and four group water and graphite reflector cross sections ,
- were provided by Argonne National Laboratory (ANL) EPRI-CELL :
. computer models. These cross sections are presented in Appendix A.
Two group graphite reflector cross sections were developed by ;
Wasserman ht UVA using the GAMTEC code. The UVA cross sections
, were used in the 2DB core models and are presented in Appendix A.
- E iW The UVA graphite cross sections were found to be consistent with i the ANL cross sections.
3.2 Diffusion Theory Mdeling Of The UVAR Core
- The 2DB code was used to model the UVAR core for the purposes .
of calculating the k-effective eigenvalue solutions. Studies presented in this report focus on reactivity effects associated i
with various core changes as determined from changes in the 2DB calculated k-effective. The EXTERMINATOR code was used to develop I
I
. . = . . . . . ,
l 1
=
l I l
lEy 22 L
2-D UVAR core models for the purpose of calculating delayed neutron fractions and prompt neutron lifetimes. The EXTERMINATOR model is analogous to the 2DB model and is ad' dressed in Section 8.
Fehr (2) developed the basic 2DB core model for HEU-18 and for LEU fuels. The Fehr models were revised and updated somewhat to be consistent with the UVA LEU element designs described in Section'2. The reader is referred to Fehr for detailed modeling descriptions.
LThe 2DB core model consists of a 60-by-62 mesh in x-y geometry. The model extends about 10 to 15 cm beyond the UVAR !
.-i grid plate in both directions. The outer boundaries were defined 3.g-with a zero flux condition. The axial direction is accounted for
~
I. i i by spatially-dependent axial buckling terms developed from ANL 3-D 1 UVAR core models.
Cross sections described in Section 3.1 were used as input !
into the code. LE0r4RO-dev?le.Nd cro:s n dians were input I dire.:tly through the LINX computer code while other cross sections were input by hand.
Two basic core models were considered in this analysis; the 4-by-4 and 4-by-5 ccre model . The 4-by-4 core model, presented in Figure 3-1, consists of 12 fuel and 4 control rod elements in a 4-by-4 array. Figure 3-2 presents the 4-by-5 core model which consists of 16 fuel elements and 4 control rod elements in an 4-by-5 array. Individual fuel elements are identified by the number i
I .
, 7.. -- ___
Lo LI l
I 23 e
appearing in the lower left hand corner of the fuel element cells ,
shown in figures S-1 and 3 2. The cores are shown with an i options 1 full graphite reflector, which means that each grid plate !
location surrounding the core is loaded with a graphite element, j Core models shown in Figures 31 and 3 2 use " staggered" control rod locations. Staggered control rod locations are .
necessary because the emergency core cooling system sprays from l the east and west directions. If the control rod locations were not staggered, spray water to adjacent elements would be impeded, i The models are built with unit cells, each with dimensions corresponding to the grid cells described in Section 3.1. l The core model was ' built" by piccing appropriate material l
rross sections in desired model cells. Fuel region cross sections ;
developed with the LEOPARD code were input in cells to simulate fuel elements. LEOPARD control rod element cross sections were l input inte cells to simulate control rod elements. Control rod l cross sections and' reflector cross sections were utilized in the !
same manner.
I !
t
- s I-I
r; e .,
.I B "
I I s-
' hh
> d- <X>6D66666666666:
s Ig s <> <XXXXXX> <X> <> <> C
<X> s <><><XX> <XXX><X><XX><
rg XXMMMMXMMMMMMMMX> C
,3 <>C :::> <> :
- l x) ;; _
- :
c . . . .
- :q) c
<>C - _. . .
R . _
) (),;
X> ::~ ,
t . b 6
~ "
<> O I <>C X> :
8
' , n CX; r_. -
c -
p . : :.. : s~
,I xys-._-. : ;-
2
- v -
)c
<>C 8 .; > <>:
X> L.
3 l .
u__m . ._ a__ n - - -
- '> C
<>c ::::.:
X> <
"~
~
- 3 () :
I .:.. :: (> O
<>C ~ ~ ~ ~
~ ~ ~
... . . . i : (> :
X>2::~ - - -
~
- : :C)C l <>C
~ ~ ~
~
~2 x,
. .,c l '
g
- =
GRAPH TE y n' ' WA~~ER g.
Figure 3-1, UVAR 4-by-4 Core Model. (Figure by M. Fehr, January, 1989) 1 .
!l
L I l l l
'I -
25 I i l i
l P l I I I I I i
- >O QtO O' 'b3 Q M X AO O f 'N'%'M , 0 0 0 !
i OO:- -;
.-.- - _ - . _ ><; j
- )o<. ;7._ _
g._ _- - -
o , l
'I
\
co- -- -~
"-~
yo; ;
~:
C~~~' :
'> _ _ . : : 0 0 ,
) :
l
- )o ; i
\ I :)
OO
>O<..-._
~:::
~E'O D
i OO l
C OO q _ __,__ :
^
g:>D f I >O< a a 3 e a - ---
., ::: -: 00 '
1 co-- - --.- --
C .- _ > O :
- ~
- )o(
- _ _
R : .. _ :- _
g g :
,cogi : ' * "
>O I
_.. :. ~~.
- .__._ C R
~ '
0 0 OC :P> O '
I m :
. 1 ,a - - -
- >O<: '_-
R -
- oo i
co: -
Q --:_.. : _;.: :
) c, ;
I h. 'I
~
- >OC __ _ . .
- 3-0 0 '
OO : __-
i : : :- : - - -
]
- : : > O: l yoc.__._
p o ,
'O066066606660645i55 5 '0 0- !
I :>O000 O0000000000 0000000000 0000 000 000 0 0 0 000 0 ,
i g $333333333333333333333
{
l I 22:::: 35 33333 GRAPHITE 32; WATER lI l
t Figure 3-2. UVAR 4-by-5 Core Model. (Figure by M. Fehr, January, ,
1989)
I I
i i
26 l l
CHAPTER 4.0 CONTROL R00 WORTH STUDY ;
Control rod worths have been determined for 4 by 4 and 4 by 5 l l UVAR cores constructed with the LEU 22, LEU 18, and HEU 18 fuel options. The term cuntrol rod" encompasses both the safety rods f
[
(Rods 1, 2, and 3) and the regulating rod as described in Section !
- 2. The effects of core life and core reflector material on control rod worth have also been determined. Additionally, l control rod worth as a function of the positioning of the rods in i the core has been investigated, f 4.1 Rod Worths In 4-by-4 and 4 by-5 UVAR Cores l f
Control rod worths have been determined for both the 4-by-4 j I and 4 by 5 graphite reflected UVAR core models for each fuel
{
option. Rod worths were determined by replacing control rod water hole cross sections by the control rod cross sections :
presented in Appendix A. Control rod locations are shown in l i Figures 3-1 and 3-2 (See Chapter 3) for the 4 by 4 and 4 by-5 [
l cores, respectively. Control rod worths were initially developed f' l
in the Spring of 1988 and were documented in a UVA internal memorandum [23). The neutronics models have since been revised ;
l and the rod study.has consequently been updated. Results of the opdated study, not significantly different from the original l
study, are presented here. .
l >
l Table 4.1 presents calculated control rod worths for the UVAR l
, 4 by-4 and 4-by 5 cores for each fuel option. Rod worths 1
i s
I I
I 27 presented in Table 4.1 have an estimated computational error of the order of five percent.
I Table 4.1 Control Rod Worths for the UVAR 4 by-4 and 4 by 5 Cores.
I Rod RodWorths(hp)-
- 1. 4-by 4 UVAR Core a Reg. Rod 0.67 0.66 0.58 b Rod 1 3.75 3.69 3.77 c Rod 2 3.93 3.09 3.97 ;
d Rod 3 2.32 2.23 2,25 t t
!! 4 by 5 UVAR Core :
a Reg. Rod 0.39 0.39 0.33 l b Rod 1 2.46 2.38 2.47 i c Rod 2 3.21 3.18 3.21 i d Rod 3 2.79 2.75 2.79 l
Information provided in Table 4.1 shows that rod worths are !
l not significantly affected by the different fuel types.
The maximum control rod wortn is that of Safety Rod 2. This ;
l I is expected because Rod 2 is the most centrally located rod ar.d is '
l consequently exposed to higher flux and carries a higher adjoint i l
importanc than the other rods. This is also consistent with the l fact the peak flux value (determined in Section 5) is also found in the same region.
I Experimentally determined control rod worths for the 4 by 4 ;,
HEU-18 core are available [24) and are compared with the calculated worths in Table 4.2. Experimentally determined worths
i 1
LI ;
J
' I 28 were converted from dollars to % 6p using an effective delayed )
neutron fraction of 0.0074. (InformationpresentedinChapter8 {
shows that an effective delayed neutron fraction of 0.0074 is appropriate for the HEU core.) ;
Table 4.2 Comparison of Ex>erimentally Determined and Calculated ,
Rod Worths for tle 4-by 4 HEU-18 UVAR Core.
Experimentally >
I Rod 10 Reg. Rod Determined Worth (h o) 0.43 Calculated Worth ihr) 0.58 Difference 35%
i i
j 3.77
'I Rod 1 Rod 2 Rod 3 3.52 3.70 2.27 3.97 2.25 <1%
7%
7%
j The calculated control rod worths presented in Table 4.2 match well with the experimental values. It should be noted that .I an estimated 10 percent absolute error is considered arnropriate for the experimental values. The calculated worth of the l regulating rod is a bit high; however, because of the minor role j l
that the regulating rod plays in core safety (the regulating rod
- is non scramable), no further calculational refinement was considered necessary. (ANL has privately stated (25) that they also over-calculate the worth of the regulating rod. The reason l l
for this is not readily apparent, but may indicated a naed to f l
examine the cell boundary conditions or the cross section library for iron and other stainless steel constituents.) i i
1 I .
L 4
e k
29 4.2 Core Life Effects On Control Rod Worth Core life effects on control rod worth were determined for
. the single case of a 4 by 5 graphite reflected core using LEU 22 fuel. The results of this study are generally applicable to other control rods and fuel options.
The worth of control rod I was deteminet! at four different stages of core burnup: 1) 0 MW days; 2) 4 iG daysi 3) 194 MW-days; and 4) 394 MW days. The control rod vorths calculated at each stage of core burnup are presented in Table 4.3.
Table 4.3 Core Life Effects on Control Rod Worth Core Burnuo (MW-day) Worth ( b )
O 2.40 4 2.51 194 2.57 394 2.62 The results above show that control rod worth tends to increase slightly cver core life. The worth associated with 0 MW-day burnup is calculated for a " clean" xenon free core.
Experiments have shown that the UVAR reaches equilibrium xe,'on conditions after about 40 hours4.62963e-4 days <br />0.0111 hours <br />6.613757e-5 weeks <br />1.522e-5 months <br /> of continuous 2 Megawatt operat...i (i.e.,80MW-hrs). This corresponds roughly to the first bu nup step of 4 MW-days (96 MW hrs). The difference in rod worth between the O MW day and 4 MW day cases is most likely .a result of changes in spatial flux shape due to xenon rather than as a result
30 of fuel burnup. Control rod worth changes for the longer time steps (all at equilibrium xenon conditions) are due to fuel burnup.
Because control rod worths tend to increase with core burnup (and with xenon bulldup), the shutdown (safety) margin also increases in the direction of conservatism - a favorable trait.
Control rod worth curves form the basis for demonstrating compliance with NRC requirements and determining the worth of experiments. Because control rod worths change with core burnup, the rods should be periodically recalibrated to account for the burnup effect. Currently, UVAR procedures require rod recalibration every 1200 W hours (50 Mw days) for just this reason, As discussed earlier, control rod worths may change due to xenon buildup over operations extending over several days; however, no method e::tsts to account for this change. Because rod worths swing in the conservative direction (i.e., greater shutdown margins) this effect is not deleterious to safe operation.
4.3 Reflector Effects On Control Rod Worths Reflector effects on control rod worth were determined for three representative reflector options for the LEU 22 4-by-5 UVAR core. The reactivity worth of control rod I was determined for the UVAR using Reflector Options 1, 5, and 13 as described in Section 7.1.
mumm-mmmmmemu mm mmmm m i i : : i:-
)
31 1
Reflector Option I consists of a full grid plate of graphite surrounding the centrally located 4-by 5 core. Reflector Option 5 is similar to Option 1 except that the entire East core face is reflected by water. Reflector Optinn 13 consists of a f ull water reflector surrounding the 4 by 5 core.
The reactivity worths of control rod 1, calculated for the UVAR core using the various reflector options are presented in Table 4.4. Control rod worths were calculated relative to the unrodded core.
Table 4.4 Reflector Effects on Control Rod Worth a Reflector Control Rod 1 g Ootion Worth ( b )
- 1 2.47 5 2.75 13 2.84 The results presented above show that control rod worth tends ta increase as the core graphite reflector is replaced by water reflector. The rod worth for the all water reflector (Option 13) is about 15 percent higher than for the all-graphite reflector (Option 1). This implies that control rod worths should be reexamined (recalibrated) following significa.nt changes in core reflector. UVAR procedures address this effect by requiring control rod recalibration with removal or insertion of one or more graphite elements next to the core.
m .
32 4.4 centrol Rod Position Effects Control rod worths for control rods in the alternate locations shovn in Figure 4-1 have been calculated for LEU 22 I fuel.
- it has been suggested that rotating the UVAP. core by 40 degrees would lead to higher neutron fluxes in certain
/ experimental facilities. As discussed in Chapter 3. the UVAR core models used staggered rod positions to prevent interfering with q the emergene.y core cooling system (ECCS). If the core is rotated the rod pattern shown in Figure 4-1 is proposed to prevent interference with the ECCS.
Reactivity worths of the control rods in the alternate positions shown in Figure 4-1 for LEU 22 fuel are presented in Table 4.5 m
Table 4.5. Control Rod Wot U s For Alternate Rod Locations In UVAR 4-by-5 Core.
Rod Worth Rod InSer M M Ar)
Reg Rod 0.80 Rod 1 2.56 Rod 2 2.74 Rod 3 2.70
- Information provided in Table 4.5 shows that the safety rod
{ worths (Rods 1,2, and 3) show little variation, as expected, due to the symmetri al nature of their locations, i
l 33 l.
I s>'sp00<ysg': s sss sa ssgss; -
I -
4
-_. 200KyQQLgOO'X)OO . . . ._
30<
~
.>% - ._U'E I,5::::1.. :
.__...._.. ._.. . CO OO :_ ;>C<
I
- >Oc -
- : :;--- :: : OC i
oO;;;.; _-_..
.._;)o< -
- ) gg - .- .. . ...-
9 g l
5 po:
- >O<__.._._
- q:g .
~
-7 ) o ;
C C
~ ~
0 0 OO o R :. >Q I >O< OO 13 2 -- - - - -
co.. __._
3 g<
L>OI: __:
~ ~
. . . . : '. . OO g
OC:: :-- T1~q _.._
>O y >OC_ --
R E - f:_: : i. :0 0 CC__. . ._
l G ---r : -
. : > Of
>O;_ __.. _ _
- ': - -- - g D O CO -
- ::>O;
- >Of '
4_ _ _ _ _.. .
_: C 0 COI-'
~ ' ~ ' ' '
>O C- -- -
-Tp:
1 _i v r .: _ $ _ __ _ [..._[.' -
QQC s> O CODDOd66CD.
- >OO '
6A66 6606 s- s <
C0000 0< >O000 0 0 0 000 CQ '-O O 000' Os. C000000 00 0<
& c2222'3222223322MFJ
= xx ljijj GRAPHITE 's' WATER Figure 4-1. Alternate Control Rod locations for UVAR 4-by 5 Core.
i
)
1 34 Rod worths for Rods 1 and 3 in the alternate locations are about the same as calculated in the staggered rod positions (Table l 4.1) because their relative core locations are similar. Rod 2 in the alternate location model has about 1/2% 4 less worth than in the stiggered location model. This is because Rod 2 is more centrally located in the staggered model resulting in higher importance. By the same argument, the regulating rod worth in the alternate model is more than twice that calculated in the staggered model.
The sum of all control rod worths in the alternate model is about the same as for the staggered niodel.
4.5 Control Rod Worth Study Summary Based on the results presented above, the following summary l
statements can be made:
1
- 1. Control rod worths are not significantly affected by the choice of fuel option.
- 2. Control rod worths are dependent on the rod position in the core with centrally-located rods having the highest worths.
- 3. Good agreement was found between experimentally determined and calculated rod worths.
~ .
h-
\
35
- 4. Control rod worths change with xenon build up and core burn-up implying that control rods should be periodically l recalibrated,
- 5. ,ittrol rod worths change with core reflector changes implying the need to recalibrate rods following core reflector changes.
- 6. Control rod worths vary for tiie alternate rod-position model considered; however, the sum of the rod orths was almost identical to that calculated for the staggered-rod position model.
I I
I l
I I
36 5.0 RADIAL PEAKING FACTOR $
l Radial peaking factors have been determined for the UVAR core for the purpose of supporting the Thermal Hydraulics Analysis portion of the HEU to LEU Conversion effort. The sensitivity of radial peaking factors to variations in core configuration, control rod element orientation, control rod insertion, core t,urnup, and core reflector options have been evaluated. Axial peaking factors are not addressed in this analysis; however, experimental data pertai:iw to axial peaking factors for the HEU UVAR are available (26).
Peaking factors were determined using the 2DB computer code described in Section 3. The 2DB code providos a meth-centered flux map (and optional edge-centered fluxes) for each energy group. The 2DB code also provides average group flux values for specified edit regions. The " active" or " lattice" fuel region of the reactor core, where heat production occurs, was defined as a ,
single edit region. Peaking factors were calculated by taking the peak to-average thermal flux ratio found in the active fuel region of ihe reactor core.
Thermal flux peaking factors are considered representative of power peaking factors. Neglecting fast fission effects, pcwer production can be considered proportional to Effth, where Er is the thermal macroscopic fission cross section and th is the thermal nautron flux. If Lt is cr idered constant, power
h i
[.'
37 production is proportional to 4th. For all beginning-of-life 1
(BOL) unburned UVAR core models, Er was indeed constant over the a
fuel region. All analyses presented below use a BOL core model t except for the core burnup analysis presented in Section 5.5.
(The core burnup analysis accounts for a non uniform It f
- distribution.)
Thermal flux peaking factors were used in the original llEU UVAR thermal-hydraulic safety analysis and in the UVAR 2 megawatt upgrade safety analysis, and are considered appropriate for use in the LEU thermal-hydraulic analysis. The thermal-hydraulic analysis focuses on a single " hot channel" and relies on t.eaking factors to relcte the hot channel to the nominal channel which is normalized to core size and average power.
The calculated peak thermal flux in the fuel region was found to occur immediately adjacent to a control rod element water hole.
The model mesh boundary between the water hole and fuel region is located between the control rod guide plate and the first adjacent fuel plate (see Figures 2-3 and 2 4). Consequently, the fiux at
( the m:sh edge boundary between the fuel region and water hole is higher than the flux at the center of the first fuel region mesh.
Although the mesh edge flux is higher, the channel between the 61rst adjacent fuel plate and the control rod guide plate is not considered the " hot-channel" because it only receives heat from the single fuel plate. The next adjacent channel is the
- hot-i
- j. o 5
I I'
38 channel", rMeiving the maxirta heat flux t' rom the two surrounding fuel plates. Mesh centered flux values for the computer models are located in the " hot channel" and are therefore considered more h appropriate than edge centered values. ?iesh centered fluxes have been used throughout this analysis.
f The LEV peaking factor analysis focused on the four areas of study described below:
- 1) 4 by 4 and 4-by 5 graphite reflected core;
- 2) control rod element rotation effects;
- 3) control rod insertion effects;
- 4) core burnup effects; and
- 5) core reflector effects.
The results of these studies are presented below.
5.1 Radial Ptaking Factors for 4-by-4 and 4-by-3 UVAR Cores Radial ceaking factors have beer Jetermined for the UVAR reactor in 4-by-4 and 4-by-5 graphite reflected, unrodded BOL core configurations. Peaking factors for LEU 22, LEU 18, and HEU 18 fuel options have been evaluated for each core configurations Radial peaking factors determined for the above-defined core configurations and fuel options are presented in Table 5.1.
Information presented in Table 5.1 shows that peaking facturs are slightly higher in LEU fuel than in HEU fuel. This is primar-ily due to a harder flux spectrum downscattering in the control !
)
39 rod water holes. The LEU 18 peaking factors are about 4 to 6 percent hightr than the HEU 18 peaking factor.
Table 5.1. Radial Peaking Factors for UVAR 4 by 4 And 4 by-5 Core Configurations With Various Fuel Options.
Core Confiauratton/ Fuel Qqtiga Radial Peakina Fgip,r
- 1. 4 by-4 Core Configuration
,l a) LEU-22 1.72 b) LEU-18 1.67 c) HEU 18 1.61
- 2. 4 by 5 Core Cor figuration a) LEU-22 1,81 b) LEU-18 1.76 c) HEU 18 1.66 Peaking factors increase with core size and with the number of plates per fuel element. Peaking factors for the 4 b/-5 cores are about 3 to 5 percent higher than for the 4 by-4 cores. The peaking factor for LEU 22 fuel is also about 3 percent higher than for LEU 18 fuel.
Peak thermal flux values were found in locations immediately to the left of the Control Rod 2 water hole (see Figures 2.3 and 2.4). This location is ci"i tent t with the measured peak flux found by Sternberg [26).
l.imited c.xperimental data is available on U/AR HEU radial
i
- I 40 peaking factors. A radial peaking factor of 1.37 was experiment-ally determined for a partially burned UVAR 4 by 4 HEU-12 (12 flat plates per element) core via foil flux mapping (5). An estimated experimental error of 10 percent is considered appropriate for l
I this value. ;
i This experimental value was subsequently " corrected" (5) using a 1 D neutronics code to predict a value of 1.45 for the j HEU-18 core. The value of 1.45 is about 10 percent lower than the !
value of 1,61 estimated in this analysis, which appears to be i within the experimental and computational uncertainty. l The 4 by 4 core configurations are the limiting core configu- :
I rations (for each fuel type) for the Thermal Hydraulic Analysis t
- (4). This is primarily due to the higher avenge power density l
I and non element bypass flow rate associated rith the small cores.
l These factors outweigh the importance of the slight ~ increase in .
peaking factors associated with the larger cores. No core larger i l'E l than 4-by 5 was considered in this analysis.
5.2 Control Rod Element Rotation Effects on Radial Peaking Factors !
The control rod elements are asymmetrical in that there is one more fuel plate on one side of the control rod hole than on
<g Because the LEU tuel elements l the other side (see Section 2). ,
will use flat plate fuel, it should be possible to load the control rod elements (and fuel elements) into the grid plate in L3
,I
41 two different orientations, each 180 degrees rotated fron the other.
The curved plate HEU 18 fuel, currently used in the UVAR, cannot be rotated. This is due to the physical nature of the HEU-18 curved plate elements which interlock and consequently must all be loaded in the same orientation.
Radial peaking factors presented in Sectica 5.1 were calcula-ted for cores in which all control rod elements were loaded with the extra plate on the left side of the element. To characterize the effect of different control rod orientations on peaking factors, three different control rod orientation scenarios were i analyzed:
Scenario 1. All control rods loaded with the extra plate on the right side of the element.
Scenario 2. Rods 1 and 3 loaded with the extra plate on the right side of the element; Reg Rod and Rod 2 loaded with the extra plate on the left side of the element.
Scenario 3. All control rods loaded with the extra plate on the left side of the element.
Control rod locations are shown in Figures 3-1 and 3-2.
Table 5.2 pretents radial peaking factors determined for the different control rod ele,. Ant orientation scenarios. It should be noted that values presented for the HEU-18 fuel option in Scenario I
h t i
i i'
I 42 t
2 are fictitious because the curved piate fuel cannot be physically loaded in opposing orientations. HEU-18 peaking l factors for Scenario 2 are presented for comparison purposes only.
I Table 5.2 Radial Peaking Factors for Control Rod Orientation Scenarios j Radial Peakina Factors I Core Confiouration/ !
Fuel Loadina
- 1) 4 by-4 Core Scenario 1 Scenario 2 Scenario 3*
Configuration i l
a) LEU-22 1.72 1.66 1,72 !
b) LEU 18 1.67 1.62 1.67 c) HEU-18 1.61 1.57 1.61 l
- 2) 4-by-5 Core :
Configuration j I a) LEU-22 1.81 1.75 -
{
b) LEU-18 1.76 1.70 - I c) HEU-18 1.66 1.62 -
.
- Peaking factors for Scenario 3 were calculated for the 4-by-4 !
core configuration only, i Information presented in Table 5.2 shows that radial peaking factors for Scenario 2 are about 2.4 to 3.5 percent lower tnan those presented in Scenario 1. Scenario 3 values are virtually identical to Scenario 1 values. This is reasonabic because of the near symmetry of the core and the fact that all elements have the ;
same orientation in Scenarios 1 and 3.
I I- .
4 . ~ , -.-n . - . - - . , . . . ,. , --.-,-.,.n - , - . . . , . - . . - . ,--=
L L
8 i
I 43 The differences in values can be taken as an estimate of the overall precision associated with calculated peaking factors, i 5.3 control Rod Insertion Effects on Radial Peaking Factors Y The peaking factors presented in-Sections 5.1 ana 5.2 were l
calculated for unrodded cores. In order to gain insight into the i
effect of control rod insertion on peakir,9 factors, the LEU 22 4-by-4 graphite reflected BOL core was analyzed with various j control rod insertion scenarios. This analysis is not completely l
rigorous due to the limitations on diffusion theory when handling i interfaces at control rod boundaries and the approximations f involved in developing cross sections for the control rod regions (see Section 4). None the less, this analysis provides insight into the general behavior of peaking factors with control rod .
i inse-tion, Table 5.3 presents the radial peaking factors associated with the LEV 22 4-by-4 BOL graphite reflected core with various control ,
rod insertions.
Informnior provided in Table 5.3 shows that radial peaking l factors increase by as much as 3 percent with single rod insertions except for Rod 2. Because the peak flux location of the unrodded core is adjactnt to Control Rod 2, insertion of rods '
I other than Rod 2 push the flux up in the Rod 2 region causing an increased peak. Similarly, insertion of Rod 2 depresses the fluy in that region and ' pushes
- the peak location towards the Rod 1
'8 g
I :
i I i 44 :
l water hole. Consequently, the flux shape is flattened and the peaking factor decreases by about 2 percent.
i I '
Table 5.3. Radial Peaking Factors for LEU 22 ,
4 by 4 Graphite Reflected Core With i Various Rod Insertions.
Core Statjui Peakina Factor Peak location
- 1. Unrodded Core 1.72 Adjacent left side :
(No rods inserted) of Rod 2 l l
- 2. Rod 1 Inserted 1.77 Adjacent left side :
Only of Rod 2
- 3. Rod 2 Inserted 1.69 Adjacent Right i Only Side of Rod 1 i
- 4. Rod'3 Inserted 1.75 Adjacent Left side Only of Rod 2 ,
- 5. Re Rod Inserted 1.76 Adjacent left side On y of Rod 2
- I The Peaking Factors presented in Table 5.3 vary from 1.69 to 1.77 (about 5%) depending on which rod is inserted. Combinations l
of more than one rod inserted at a time were not evaluated because i the UVAR cannot be taken critical with more than ene rod fully l inserted due to regulatory restrictions on shutdown margin and I- e'A cess reactivity.
5.4 Cere-Life Effects on Radial Peakuig Factors The peaking factors presented in Sections 5.1, 5.2, and 5.3 wer? calculated for beginning of life (BOL), unburned cores, t
During reactor operation, the core is slowly " burned" resulting in I
e s
45 1
less U-235 and e buildup of fission products and poisons. The effect of core burnup on peaking factors was determined by analyzing a 4-by 5 LEU-22 graphite reflected core at various burnup times. Burnup times evaluated were: 1)noburnup;2)194 MW days burnup; and 3) 394 MW days burnup.
As discussed earlier, the validity of the thermal flux peaking factor is tied to the assumption that If remains constant over the fuel region. This assumption is valid for BOL cores; however, because cores burn unevenly, If is not constant through-out the fuel region in a burned core. In order to account for the non uniform It distribution in the core, the following approxim-ation scheme was used to calculate peaking factors for burned cores:
l
- 1. The peak value of (Ef4th)1ocal was determined by multiplying the local If determined for each fuel element by the associa-ted mesh thermal fluxes, in each case the location of the peak value of (Effth)1ocal coincided with the location of the peak thermal flux.
- 2. An average value I fwas determined by simply summing eier.ient specific Ef's over the 20 elements and then dividing by the 20 elements. The average thermal flux (Ith) in the active fuel region was provided as output by the 2DB code.
- 3. The peaking factor was then estimated as follows:
peaking factor = $8EEI Ef 5th <
)
L 46 Ideally, the denominator would be a flux-weighted average of It over each fuel region mesh point. Unfortunately, this value is not readily available from the information provided in the 2DB output *, The above described weighting scheme provides an adequate approximation for determining the general behavior of peaking factors with core burnup.
Table 5.4 presents radial peaking factors for the UVAR 4 by-5 LEU 22 graphite reflected core with 0,194, and 394 MW days burnup. Also presented in Table 5.4 is the min to max variation in Er over the core region for each burnup step.
Table 5.4. Core Life Effects on Radial Peaking Factors for the LEU 22 4-by-5 Graphite Reflected r. ore.
Burnuo (MW day) Peakina Factor 2;r variation 0 3.81 0%
194 1,77 3%
i 394 1.73 6%
- Certain edits were added to the 2DB code by the University of Michigan that appear to address power peaking factors; however, there are unresolved inconsistencies in these edits. Preference was oiven to the original coda edits for flux values which were used~throughout the peaking factor analysis.
I
)
1 47 Information presented in Table 5.4 shows that peaking factors decrease with core life. This is consistent with the fact that there is more local burnup in the peak flux region than in the rest of the core resulting in a flattened flux distribution.
5.5 Core-Reflector Effects on Radial Peaking Factors Reflector options for the UVAR core were analyzed for the LEU 22 and LEU-18 fuel options as discussed in Chapter 7. Peaking factors associated with the various reflector options have been determined. Table 5.5 presents peaking factors associated with three representative reflector options for LEU-22 and LEU-18 fuel.
Table 5.5. Core-Reflector Effects on Radial Peaking Factors
--- Ap(%Ak/k) ***
Reflector Ootion LEU 22 LEU 18
- 1. All Graphite Reflector 1.73 1.67
, 2. Single Row of Graphite 1.75 1.69 Surrounding Core
- 3. All-Water Reflector 1.76 1.73 Peaking factors were found to be relatively insensitive to the various reflector options and showed a maximum variability of about 3 percent.
l l
- l
\i 48
, 5.6 Radial Peaking Factor Sumary Peaking factors for various core configurations, control rod orientations, rod insertions, core burnup, and reflector options have been cbtained. Based on information presented above, the f
> I following conclusions have been made:
I i
- 1. Radial peaking factors for LEU fuel are slightly higher ;
than for HEU fuel. Peaking factors increase with core size and with the number of fuel plates per element. ;
- 2. Radial peaking factors may vary by as much as 3.5 percent l
- c. depending on control rod orientations.
- 3. Insertion of a single control rod can produce a change in l peaking factor ranging from a 2 percent decrease to a 3 l '
percent increase.
- 4. Peaking factors decrease somewhat with core life due to i
the higher burnup occurring in the peak flux region. I I' 5. Peaking factors are relatively insensitive to the core-1 reflector option, showing a maximum variability of about i 3 percent.
1 I
I I
I I
l I 49 j CHAPTER 6.0 TEMPERATURE AND NODERATOR EFFECTS STUDY Temperature and moderator void effects have been characterized for the UVAR reactor. Section 6.1 addresses temperature effects and Section 6.2 addresses moderator void effects. .
6.1 Temperature Effects ,
Temperature effects have been determined for the UVAR reactor ;
, for both a 4 by-4 and 4 by-5 graphite reflected core :
configuration. Temperature effects on core reactivity have been quantified separately in terms of: 1) doppler effect; 2) moderator temperature effect; 3) power defects; and 4) gross I system heating effect. Temperature effects for the 4 by-4 core l
configuration were determined in the Spring of 1988 and were f l
documented via a UVA internal memorandum [27). The temperature-effects analysis for the 4 by-5 core configuration was performed in early 1989 and incorporates the latest revised core models. '
Temperature effects presented in the following sections are determined from the change in core reactivity associated with that :
effect relative to the reactivity associated with a " reference" core condition. Reference core conditions were determined from a review of information used in the UVAR HEU-to-LEU thermal hydraulics wo P. where the fuel plate to moderator-temperature differential was attermined to be about 900 for 22 plate elements and about 1200 for 18 plate elements. The temperature gradient I
= - - - - -
=.a.:==. .: - - - - --- --
I I 50 across ti.e thin fuel plate was assumed to be negligible.
Temperatures used for the reference cores are provided below: ;
I Fuel Plate Moderator Element T g Temoerature (00) Temoerature (00) 22 plate 30 21 Ei- 18 plate 33 21 ;
Temperature effects for the 4-by 4 cores were determined for I LEU-22, LEU-18, and HEU-18 fuel options. Temperature effects for j the 4 by-5 cores were determined for LEU-22 and HEU-18 fuel options only. }
6.1.1 Doppler Effects [
Doppler effects for the UVAR core were determined by l adjusting the resonance temperature parameter in the LEOPARD code. ,
Table 6.1 presents reactivity worths and doppler coefficients associated with resonance temperatures of 75, 100, and 200 degrees I celsius. Differential reactivity worths were determined from the change in k-effective associated with each resonance temperature relative to the reference case k-effective. Doppler coefficients were determined by dividing reactivity worths by the associated change in resonance temperature.
L Doppler coefficients for each fuel option remain relatively stable over the 7500 to 2000C temperature range. Doppler I coefficients for LEU are about an order of magnitude greater than l
l' 0
)
'a.
I L.m 51 ,
l'"
Table 6.1. Doppler Effects for the UVAR Reactor ,
't;
'/ -- Resonance Reactivity Doppler Coeff ;
Temperature (0C) Worth Ao (Ao Der oc)
I. UVAR 4-by-4 Core Configuration (1988 Study)
LEU 22 FUEL
'a 75 -5.3x10-4 -1.2x10-5 4 -1.2x10-5
- l. g 100
-8.1x10 l: 200 -1.9x10- 3 1,1x10-5 ;
l' LEU-18 FUEL 75 -4.3x10-4 -1.0x10-5 4 5 ,
100 fI.
1 200 -6.7x10
-1.6x10- 3 -1.0x10
-9.5x10- 6 HEU-18 FUEL I .75-100
-3.8x10-5
-5.9x10*7 9,ox10-7
-8.8x10-7 200 -1.4x10-*- -8.3x10-7 II. UVAR 4-by-5 Core Configuration (1989 Study)
I L:
l LEU-22 FUEL L -
75 -5.1x10-4 -1.1x10-5 100 -8.4x10-4 -1.2x10-5 200 -1.8x10-3 -1.1x10-5 HEU-18 FUEL
. 75 -5.6x10-5 -1.3x10-6 100 5 -1.3x10-6
~
E00 -8.8x10
-2.1x10- 4 -1.2x10-6 eg i
. , . = ...
E' i l
I 52 l for HEU fuel.- This is consistent with the fact that LEU fuel has a much higher uraniui 238 loading and is thus more sensitive to resonance effects. Similarly, LEU-22 fuel has a slightly higher doppler coefficient than LEU 18 because of the higher U-238 U loeding per element. It should be noted that the values presented in Table 6-1 neglect doppler broadening of the uranium-235 fission
,I.
cross'section.
Recommended values ~ of doppler coefficients for the UVAR reactor are: e
- 1) LEU-22: - 1.1 x 10-5 % ap per OC
- 2) LEU-18: - 1.0 x 10-5 % op per DC
- 3) HEU-18: - 1,3 x 10-6 % ap per DC As estimated absolute error of 10 percent is believed to be g
assaciated with each of the above values.
6.1.2 Moderator Temperatura Effects j .
Moderator temperature effects were determined by adjusting
(
L. . the moderator temperature parameter of the LEOPApD code. Table 6.2 presents reactivity worths and moderator temperature coefficients for the UVAR core. Reactivity worths were calculated based on changes in %-effective relative to the reference case
.h (moderatcr temperature - 210C).
Moderator temperature coefficients are relatively stable over i the terrnrature range of interest and do not vary significantly l
betwet.. iuel types. Recommended moderator temperature
s
[-
p ,
53 coefficients for the UVAR are the values presented for the 5000
{ case for the 4-by-5 core because the UVAR normally operates with moderator temperatures below 5000. An estiniated 10%
absolute error is considered appropriate for the recommended coefficients.
I Table 6.2 Moderator Temperature Effects for the UVAR Reactor Moderator Reactivity Moderator Tempergture Temperature (DC) Worth (Ao) Coeff (op per DCi I. UVAR 4-by-4 Core Configuration (1988 Study)
LEU-22 Fuel I 50 75
-5.0x10-3
-1.1x10-2
-1.7x10-4
-2.0x10-4 LEU-18 Fuel 50 -5.3x10-3 -1.8x10 '4 75 -1.1x10-2 -2.0x10 '
HEU-18 Fuel I 50 75
-5.4x10-3
-1.1x10-2
-1.8x10-4
-2.0x10-4 II. UVAR 4-by-5 Core Configuration (1989 Study)
LEU-22 Fuel 50 -3.9x10-3 -1.4x10-4 75 -8,5x10-3 -1.6x10-4 95 -1.3x10-2 -1.7x10-4 HEU-18 Fuel 50 -4.2x10-3 -1.5x10-4 75 -C.8x10-3 -1,6x10-4 95 -1.3x10-2 -1.8x10-4 Na ' )
- : :._. . .= - -
j !
l 1.
L !
54 i
6.1.3 Power Defects ,
The power defect of reactivity is defined as the change in f reactivity associated with bringing the reactor from a " cold" zero power condition to normal (" hot") operating temperatures. The ;
UVAR coolant temperature rise across the core at 2 megawatts is
- about 7.200 (130F). The power defect for the UVAR was determined by first calculating the reactivity associated with a " cold" core in which the temperature of the core components and moderator were i 1
set at 17.400. This reactivity was compared with the reactivity of a " hot" core where the average moderator temperature was 210C.
The power defects calculated for the 4-by-4 core configuration for the various fuel options are presented in Table -
6.3.
Table 6.3. Power Defects For 4-by-4 UVAR Cores g Fuel Ootion __ag i
LEU-22 -6.8x10-4 LEU-18 -7.3x10-4 h HEU-18 -7.6x10-4 .
E' Absolute errors of the order of 5 to 10 percent are '
considered appropriate for the calculated power defect reactivities.
I-
- l :
a.. .= -- - -
I I '
55 6.1.4 Gross System Heating Effects Reactivity effects associated with gross system heating have been determined for a 2000 rise in all reactor components and moderator. This is consistent with actual operating conditions on hot humid days when the bulk temperature of the pool can swing from 680.: (2000) to 1040F (4000) during a single day nf operatior.. .
Gross system heating coefficients are presented below for the 4- ,
by-5 core configuration with LEU-22 and HEU-18 fuel.
Fuel Tvoe Temperature Coeff (65/0C)
LEU-22 -1.11 x 10-4 HEU-18 -1.13 x 10-4 The gross system heating coefficients for each fuel type are l
strallar and have an estimated absolute error of about 10 percent.
6.2 Modarator Void Effects 6.2.1 Uniform Core Void Effects ,
Ur.iform moderator void effects on core reactivity, and void coefficients, have been determined for UVAR 4-by-5 graphite reflected cores with LEU-22, LEU-18, and HEU-18 fuel loadings.
Uniform moderator voids ranging from 1 to 95 percent were l
analyzed.
Uniforr.: voids in the core-moderator region were simulated by adjusting core water densities in the LEOPARD fuel and control rod element models. Cross sections generated by the LEOPARD code were E
tg
!- )
'I I 56 I then used in the 2DB diffusion theory code eigenvalue calculation.
- Reactivity worths (Ap) associated with uniform core voids were determined by the change in k-effective relative to the no-void core condition.
Table 6.4 pretents reactivity worths associated with !
i different void fractions for LEU-22, LEU-18 and HEU-18 fuels. The ;
core moderator densities used to simulate the uniform core voids are also presented. .,
-E Table 6.4 Reactivity Effects Associated With Uniform ,
3 Moderator Void Conditions With LEU-22, LEU-18, and HEU-18 Fuel Loadings.
I % Void Moderator Densityfa/cm 3)
LEU-22 (n,)
LEU-18 (g,)
HEU-18 (g,) >
LI l 1 0.99 -0.193 -0.173 -0.152 1
2 0.98 -0.386 -0.348 -0.307 5 0.95 -1.16 -1.00 -0.887 10 0.90 -2.48 -2,18 -1.94 20 0.80 -5.77 -5.16 -4.66 L-. 50 0.50 -24.7 23.3 -24.3 90 0.10 -150. -161. -165.
I L 95 0.05 -211. -236. -253.
l I' P
h 57 Uniform moderator void coefficients were determined by ,
,' dividing the r? activity worths presented in Table 6.4 by their associated percent void. Void coefficients calculateo in this manner are presented in Table 6.5.
Table 6.5 Uniform Moderator Void Coeffic'~" '
LEU-18, and HEU-18 Fuels.
Void Coefficie,.
% Void LEU-22 }
t
-0.193
- i 1 .
p 2 -0.193 . >4 r 5 -0.231 -u, -
s.177 10 -0.248 -0.21b -0.194 20 0.289 -0.258 -0.233 l 50 -0.493 -0.466 -0.486
- l .
llml'.
90 -1.67 -1.79 -1.83
, 95 -2.22 -2.48 -2.66 i
p Figure 6-1 presents a plot of the void coefficient versus l percent void for each fuel option.
- Informau on provided above shows that in the 1 to 20 percent L-void region, the void coefficient for the LEU-18 fuel is higher than in the HEU-18 fuel . This is as expected because the harder fl u~ <
- r- . . - -iated with LEU fuel is more sensitive to o
l mo e 4 > % m;c Similarly, void coefficients for LEU-22 fuel ar. nigher cno. for LEU-18 fuel due to the smaller water gaps
- associated with LEU-22 (and consequently a harder spectrum).
L s
+-
c- . .s--m, , ._v ,.,_ . - , - e
- . ( p sa v .
n ..-
f.r
~
58
~
i ,
I' I
MODERATOR VOID COEFFICIENTS I o FOR 12U-82, LEU-19. AND LBU-82 FUBL
-o.s ! -
9
, j -o.4 -
l -0.4 -
. l -1.s -
-1.4 -
.I jx I ~
-1.s -
-1.e -
ta -n -
l- 8 -a.: - .
' O -L4 -
I h -a.e -
-Le 1 2 5 10 20 , 50 to DS PERCENT MODERATOR VDID D BEU-15 FUBL + LEU-18 FUEL o LBU-Et FURL t.
r Figure 6-1. Void Coefficient As A Function Of Void Percent fl P
- _N_._____-. . _ _ . _ _ _
I l g
I 59
~
A 20 percent total core void is significant and is considered a reasonable upper limit for void-effect studies. Void effects in this analysis were also considered for higher void f actior.s of 50, 90 and 95 percent. It should be noted that analyzing these '
higher void fractions " pushes" our current neutronics computer model.ing capabilities due to limited data libraries and the breakdown of diffusion theory with increasing void fractio.1s.
Therefore, results presented above for the higher void fractions are not considered rigorous but are presented to provide insight into the general nature of void effects in severely voided cores.
6.2.2 Local Core Void Effects Local void effects were analyzed for a UVAR 4-by-5 graphite reflected core with LEU-22, LEU-18, and HEU-18 fuels. Local void effects were determined by voiding 99 percent of the core moderator in local fuel elsaent regions. As mentioned earlier, severe void conditions are difficult to model with diffusion theory modeling techniques. On the othe hand, neutrons that ,
stream across a local void are not lost to the system, so that the predicted void effects may indeed be meaningful.
Four different local void scenarios were considered. Each scenario consists of voiding 99 percent of the core moderator in particular regions. The locations of the voided regions (see Figure 3-2) considered are as follows:
i Scenario 1 - West region of Element 11.
l L
1
4 -
{
[
):
I e
60 Scenario 2 - East regions of Elements 10 and 6, and West regions of Elements 11 and 7.
Scenario 3 - West region of Element 11 and entire Element 10 fuel region.
Scenario 4 - Fuel regions of Elements 10,11, 6, and 7.
Each scenario includes voiding the west region of Element
- 11. This is where the " hot-channel" is located (determined by the Peaking Factor study in Chapter 5) and is therefore the most likely region to vold.
Table 6.6 presents reactivity effects associated with the I-different void cases for LEU-22, dV-18, and HEU-18 fuels.
Reactivity worths are calculated relative to the unvoided core.
Table 6.6 Local Void Reactivity Effects for Void Scenarios 1 Through 4 LEU-22 LEU-18 HEU-18
.Cng (ho) (ho) (%)
1 -0.678 -0.730 -0.834 2 -2.88 -3.16 -3.59
\-
3 -3.49 -3.77 -4.23 4 8.67 -9.50 -10.7 I
' ei. '
-,_.m_ _ ___
'f' M
t 61 Unlike' the uniform moderator void effects discussed in Section ,
h,localvoidreactivityeffectsappeartobegreaterfortheHEU l fuel option which has a softer spectrum.
Local void coefficients have been developed by dividing the reactivity worths presented in Table 6.6 by the percent core water void. The' percent core water void represents the ratio of the ,
void volume to' the total moderator volume of the unvoided core.
Table 6.7 preser.ts local void coefficient-and percent core voids associated with each scenario for LEU-22, LEU-18, and HEU-18 fuel.
I Table 6.7 Local Moderator Void Coefficients (%6p/% void)
% of Core Scenario Voided LEU-22 LEU-18 HEU-18 l 1 1.3 -0.52 -0.56 -0.64 l 2 4.6 -0.63 -0.69 -0.78
~
3 5.9 -0.59 -0.64 -0.72 4 14 -0.62 -0.68 -0.76 Local moderator void coefficients appear to be relatively constant over the scenarios considered. Unlike the uniform void cases, higher void coefficients were consistently determined for h
the HEU soft-spectrum fuel option because these voids occurred in regions of high importance.
E 1 t LE !
u' _ . _ _ _ _ . _ . , _ .
o ... - . . . . .-
62 6.3 Temperature and Moderator Effects Summary Te'nperature and moderator effects were analyTed for the UVAR reactor. Based on the results preseated above, the following summary statements can be made.
- 1. Doppler effects in the LEU fuel were found to be about an order of magnitude greater than in the HEU fuel due to the higher U-238 loading.
- 2. Moderator temperature coefficients were about the same for the different fuel types and were larger than the doppler coefficients by an order of magnitude for the LEU fuel and by two orders of magnitude for the HEU fuel.
- 3. The power defects and gross system heating effects calculated for the different fue'. types were similar to each other and appeared to be dominated ay modentor temperature effects.
- 4. Uniform modcrator void coefficients wero.5 lightly higher in the LEU fuels thar in the HEU fual die to harder flux spectra as expected. This trend faltered as void fractions were increased past 50 percent.
- 5. Local moderator void coefficients appear to be relatively constant fo. the various scenarios and, unlike the uniform void coefficients, higher coefficients were found for the soft-flux spectrum cases (HEU-28) than for the 120 cases.
l 1
w _ _ . _ _ - _ - - - _ . . - - - - _ _ . . - _ _ - - - _ -
W
,)
ti3
, CHAPTER 7.0 REFLECTOR WORTH STUDY Reflector worths were determined for several different
~
s reflector materials and thicknesses placed on different faces of the UVAR reactor in a 4-by-5 core configuration. Ref'ector materials analyzed include graphite, water, deuterium tanks, and black boundaries (total absorbers). LEU-22, LEU-18, and HEU-18 fuel options were considered. The results of this study were first presented via a UVA internal memorandum (28). Due to model revisions incorporated after the completion of that study, reflector worths presented here have an estir..ated uncertainty of
= hpproximately 10 percent.
- 7.1 Graphite And Water Reflector Options 9 The UVAR reactor fuel and control rod elements are positioned on the UVAR grid plate, which has an 8-by 8 array of grid j locations. Grid locations,may ia occupied by a fuel element, control rod fuel element, graphite elemer.t. or grid plug.
For the purposes of this analysis, the 4-by-5 core configuration was assumed to be centrally located on the grid plate. Remaining grid locations, surrounding the core, were assumed to contain either graphite elements or grid plugs. Grid s plugs are short devices that do not extend into the active fuel region of the core, thus a grid location with a grid plug provides pool water reflector for the core in that location. The entire f 3
grirl plate assembly is, of course, surrounded by pool water.
3 4 d .
,s .- . ..
t 64 Several different reflector options were developed using d graphite and water reflectors. Reflector Options 1 through 6 are variations of a full graphite reflector as shown in Figure 7.1.
Option 1 consists of two full rows of graphite on the north, west, and east core faces, and a single row of graphite on the
' south core face. Reflector Options 2 through 6 were created by
[ replacing graphite with water reflector on one or inore of the core faces. . Option 7, shown in Figure 7.2, consists of a single row of graphite surrounding the core. Options 3 through 12 are variations of Option 7 as shown in Figure 7.2. Option 13 (not shown) is the all-water reflector option. Reactivity worths associated with Options 2 through 6, and Options 7 and 13, I calculated relative to Option 1, are presented in Table 7.1.
j The reactivity worth associated with going from a full
=
graphite reflector (Option 1) to a single row of graphite scrrounding the core (Option 7) is only about 1.5%, vnile the reactivity swing from Option 1 to the full water reflector (Option
- 13) is over 7 percent. This shows that the importance of the
--3 first row of graphite reflector far outweighs the outermost second row.
N Table 7.2 presents reactivity effects atsociated with Options
- 8 through 13. Reactivity worths are calculated relative to Option 7.
s
s-m 65 l
i el
' /'
l.
l jrap,htte, :
// -
Water ;
/
/
/
/
J
'7
/// / / /
n / ~
/.,a,e,/ ,//
- ' Option 1 Option a
^ =
Graphite *---. Pa t e r -e
/ ' '
vvaa 4.s cut / ,
e' /j ,,,, ,,, ,,,, ,,r 3 / / N/ [h/
=
/ /Y, ,i/
't, / d'
/c W# #
/ .
- I
,,1 : Water :
Option 3 Optton 4
]-
Y
[ Graphite !Graphita
/ / / / / /
4 99a4 4s3 cent j DvM est Coat I -
Water'
/
Water
& . 1 M
~
i /
// //
- 0,u n >
/ , . .. /////
0,u.a a figure 7-1. Reflector Options 1 through 6.
.g ,
~ ~
' ' i I
I o i f
,, Water ,
/ f Grapht'te/ '
/ '
g / .. ... . / 7 7 .. ... . '/
/ / / 7
/'
3 /j -
/ f y,
/ p U
//lllll / '
ff, kV*Ph\te',l , _
optio; 7 option 8
/ Graphite // c
/
en. ... con. /
,/
'// yu ,. , c ..
/
t t /
/
/
/
,/
s' t'
s/
,/ 1 u,
- . -/ *
- , [ '
//
Water
~_ Water _,
~
opttee s '
i onti.e to Wat9r ,
e Water- -
Grap'hlte / / ,
7/,Graphife'/
/ ...... .'
7 . ... .
/p l ,/ -
/ y/
s5 b
/
lff////~
optau :
a fffff/s i "(
option la i
L. Figure 7-2.
. Reflector Options 7 through 12.
L i8 .
\
+
i 67 k Table 7.1. Ret:tivity Effects Associated With Variations Of A Full Graphite Reflected 4-by-5 Core LEU 18 FUEL LEU-22 FUEL Option B.e.flector Description M' h" 1 Full graphite reflector: two. Reference Core rows of graphite on the east,
~
north, and west face; and one row of graphite on the south face.
2 Same as Option 1 with north -1.96 -1.99 face graphite replaced by water.
3 Same as Option 1 with south 1.17 -1.24
- f ace graphite replaced by water.
2 4 Same as Option I with north -3.35 -3.45 and south face graphite replaced by water.
5 Same as Ootion I witis east -2.27 -2.33 i
_ face graphite replaced by water. ;
-2.19 -2.28 .)
6 Same as 0) tion 1 with west face grap11te replaced by !
,-- water.
7 Single row of graphite -1.48 -1.55 surrounding entire core. ;
13 All-water reflector. -7.23 -7.42
- Reactivity worths calculated relative to Option 1.
M k i
.b
__),:
I
1 s a '
)i l
68 Table 7.2 Reactivity Effects Associated With Variations Of A Single Graphite Reflector Row Surrounding j .- the 4-by-5 Core.
LEU 18 FUEL leu-22 FUEL Drtiga Reflector Descrintion ( b ) . ," ( b )"
7 Single row of graphite Refercuce Core surrounding entire Core.
8 Same as 00 tion 7 with -1.36 -1.38 north face graphite replaceo by water.
9 Same e Option 7 with -1.12 -1.18 i south face graphite replaced by water.
l- 10 Same as Option 7 with -2.64 -2.71 north and south face graphite replaced by water. -
j 11 Same as Option 7 with -1.74 -1.78 .,
east face graphite l
replaced by water.
12 Same as Option 7 with -1,69 -1.74 west face graphite replaced by water.
13 All-water reflector. -5.75 -5.87 ;
- Reactivity worths calculated relative to Option 7.
l l-t q J--___-_---__
I I
69 Information in Tables 7.1 and 7.2 shows that reflector reactivity effects are slightly greater for LEU 22 fuel than for LEU-18 fuel due to the harder leakage spectrum of the latter.
The effect of burnup on the relative worth of a reflector variant was determined for Option 2 using LEU-18 fuel. Four time steps were considered ranging from 0 to 384 megawatt-days. The I results of this analysis show that the relative worth increases on the order of 10 percent over core life. Reflector worth increues with burnup because as the core burns, the flux shape is flattened resulting in higher leakage. These results ate considered to be representative of the behavior of reflector worths with burnup in general. The results of the analysis are shown in Table ',.3.
I Table 7.3 Core-Life Effects On Relative Worth of Reflector Option 2 Using LEU-18 Fuel.
Core Burnup Relative Worth Status Ao O Mw-days -1.96%
4 Mw-days -2.02%
194 Mw-days -2.09%
394 MW-days -2.16%
7,2 Deuterium Reflector Tanks Reactivity effects associated with various-sized deuterium !
tanks were anlyzed in response to current plans to install an I experimental D2 0 tank between the core and the neutron beamport j facilities. Preliminary designs have not yet been developed, so a 5
I.
c; ,
70 wide range of D20 tank sizes were analyzed. All D 20 tar.ks were assumed to be fabricated of one-quarter inch pure aluminum plate.
Two group microscopic cross sections were oeveloped for each tank using the GAMTEC code described in Section 3.
The tanks were assumed to be 59 cm high, consistent with the active fuel length, due to the 2-D limitations of our codes. It should be noted that the actual tank heights will probably be significantly shorter than the assumed 59 cm. Thus, results presented here.are conservative. Tank sizes were consistent with integer grid plate spacings. 2-D Tank dimensions are specified in Table 7.4.
Table 7.4 Representative D 20 Tanks 2-D Physical D2 0 Tank ID Dimensions Tank 1 8.1 cm x 7.7 cm Tank 2 8.1 cm x 15.4 cm Tank 3 16.2 cm x 15.4 cm Tank 4 24.3 cm x 15.4 cm Tank 5 24.3 cm x 23.1 cm Tank 6 8.1 cm x 30.8 cm Tank 7 16.2 cm x 38.5 cm Tank 8 16.2 cm x 80.8 cm i I:
L. _. _ .
l l;'7 ,
4 l;
71
- F.eflector Option 6 was used as the reference core to determine 2D 0 tank reactivity effects. Option 6 consists of a
. full graphite reflector on the north, east, and south core faces and water reflector on the west core face. Core reactivities associated with the D20 tanks were determined by placing the tanks on the west core face (thus replacing water reflector). Table 7.5 presents reactivity worths associated with each tank, Table 7.5 Reactivity Effects Associated With Different Sized D2 0 Tanks Bordering We:t UVAR Core Face.
I Reactivity IRDh Worth N o) 1 0.38 2 0.77 3 0.89 4 1.4 5 1.3 6 1.4 7 1.8 8 2.1 I
Information presented in Table 7.5 shows that a D 02 Tank may prcduce a sigr.ificant positive reactivity effect when replacing water reflector. To determine how this compares with graphite reflector material, graphite reflector was placed in the Tank 8 I
I
V 72 location. The reactivity associated with the graphite was 1.9%
ap, which is only slightly lower than the 2.1% ap associated with Tank 8. Therefore, while the D 2 0 tanks are significantly better reflectors than light water, they are only marginally better than graphite.
7.3 Black Bot ndary Reflectors A single case, using a 4-by-5 core model with LEU-22 fuel, was analyzed with a black (vacuum) boundary on the east face of the core. The associated reactivity worth was determined to be about -8% ap t' elative to replacing graphite, and about -5 to -6%
ap relative to replacing Water. It should be noted that diffusion theory codes cannot really handle black boundaries properly so that the worths presented are only approximations that are likely to be overestimates of the true situation.
7.4 Reflector Worth Study Summary Based on the information presented in the sections above the followir.g summary conclusions can be made.
- 1. Reflector effects are slightly more pronounced in the -
LEU-22 fuel than in the LEU-18 fuel due to the higher leakage associated with the harder flux spectrum.
- 2. The calculated k-ef?ectives associated with the BOL LEU-22 and LEU-18 full-water-reflected cores are 1.03 and 1.02 respectively. It may be reasonable to initially load the new core with a full water reflector and
?
73
. subsequently add graphite with burnup, as appropriate.
This should minimize core handling operations, core f configuration changes and rod calibrations. It is hoped that reflector worths determined in this study will provide input for such dec,isions.
- 3. Replacing water with a D2 0 tank can lead to a significaint <
positive reactivity insertion depending on the tank size.
Replacing graphite with a D2 0 tank results in a relatively small positive reactivity insertion.
- 4. Replacing an entire core face graphite reflector with a heavily absorbent medium or vacuum results in a negative reactivity swing on the order of 8 percent.
1 I i I
i
I: ,
74 :
CHAPTER 8.0 DELAYED NEUTRON AND PRONPT NELTTRON LIFETIME STilDY Effective delayed neutron fractions and prompt neutron lifetimes have been calculated for the UVAR 4-by-5 graphite and water-reflected core models using both LEU-22 and HEU-18 fuel .
options. Calculation of these parameters is in response to NRC guidance received at the TRTR-88 meeting in Oregon. A summary of
.I methodology and results are presented' below.
8.1 NETHODOLOGY
- l. Eftective delayed neutron fractions (peff) and prompt neutron lifetimes (f) were calculated using the EXTERMINATOR [20) 2-D diffusion theory code. A 60-oy-62 x-y nesh, identical to the mesh used in the 2DB models, was used.
g Four group cross sections for the core region were calculated from LEOPARD models for the fuel and control rod elements. Four group reflector cross sections for water and graphite were taken from ANL EPRI-CELL models. Composition-dependent axial buckling :
, g.
u terms were developed based on ANL 3-D modeling of the UVAR core.
l; The EXTERMINATOR code requires average reciprocal velocities for each energy group. The group structure for the LEOPARD code is not explicitly defined. Fortunately, the ANL EPRI-CELL models do specify a group structure consistent with the LEOPARD code.
This group structure is presented in Table 8.1. t I.. ,
Lg
ss@* 3 g#dkt %Y2 IMAGE EVALUATION b TEST TARGET (MT-3) ////gfTf 0fp!p#
\k//77;if()*>
4y,4<
+ %,;%;g l.0 'aM d
I E 9 RE l,l [!E Ils Ji i l.25 1.4 1.6 1 150mm >
- 6" *
- r % '
/
4?; y>>,3,pf y 77
[3+4 g//p AN awm . . .
. _.J
- @ Tk 3 kg* k' IMAGE EVAI.UATION 4)y 't/ TEST TARGET (MT-3) f,[g/f Y 36[>4 ,4 4 s W' 'k,,Y,#/
I.0 lf m LWi 5 Ef E l i.i p= Il!!M l.8 l l
1.25 1.4 1.6 1 4 150mm >
4 6" >
4&?% /A
+py>>p,,y, B'p,{
+g<g>
u J
d#C# l kt IMAGE EVALUATION ////p [Mif, 4
/// T[ %[# TEST TARGET (MT-3) p#4 3 k N'k l.0 l; m L24 p l} IL4 ;
j,l l m lll!&&
l.8 I I
l.25 1.4 1.6 4
4 150mm >
4 6" >
d*%
r%y >7///< .
~
o
\
t t
++& 4
+
'.k. . . .
_. .a
l (b
4+k*$1+#*f e gi h@*
IMAGE EVALUATION TEST TARGET (MT-3) 777/p j[f +4,#4 4)y,,,, /4,, s'
+ %
l.0 lg m saa u m g=n=
ba 5D bbs l1l llg l 91 14 g
4 150mm
- 4 6 >
++k;>%' i, /b
.,;p%,777,7
<g4%
t eye k . m .
- g. -
.x - __a......=,. - -
j l
t j I. ,
It- 75 l
I g Table 8.1. Energy Group Structure Used for the EXTERMINATOR 3 Analysis.
Energy Group E-maxieV) <v> (m/s) xn 5
l 1 1.0 E7 1.8 E9 0.739 2 8.2 ES 3.4 E8 0.261 8 3 5.5 E3 0.62 4.9 E6 2.6 E5 0.00018 0.0 4
I Delayed neutron fractions used in the analysis were based on ,
delayed neutron yields presented by Brady, et al (29] for six precursor groups. Information was available for nine and twelve precursor groups; however, EXTERMINATOR is limited to eight groups. The delayed neutron yields from Brady were normalized to a delayed neutron fraction of 0.0065. The delayed fractions (si) .
are presented in Table 8.2.
Table 8.2. Delayed Neutron Fractions
. . GrouD % Yield Bi 1 0.0351 0.0002282-2 0.1809 0.001176 3 0.1778 0.001156 5 4 5
0.3837 0.1565 0.002494 0.001017 '
6 0.0659 0.0004284 Total 0.0065 I
8.2 MODEL VERIFICATION The EXTERMINATOR UVAR model was verified by: 1) checking model inputs; and 2) comparing calculated k-effective values with I
e v , - . . , - . . -... ,-n... .__,.._,v.,..__,~,,.w..,,- -r,. _ _ , . - . - - _ . - - . - , - _ _ . - . -. - - . - - . _ _ _ . . . - , . , .
<l
- s. q
. 1 76 ,
! those calculated by our existing two group 2DB models.
Both two and four energy group EXTERMINATOR models were developed for the LEU 22 fuel core. A four group EXTERMINATOR l model for HEU-18 fuel was also developed. Table 8.3 presents a ;
I comparison of k-effective values calculated by the EXTERMINATOR and 2DB codes. As stated earlier, the model core is a 4-by-5 graphite reflected core.
!I Table 8.3. Comparison of k-effective Values Calculated By the EXTERMINATOR and 2DB Codes, k-eff k-eff Model EXTERMINATOR 2D1
. 1. 4 Group LEU-22 1.115 -
2 Group LEU-22 1.110 1.109 l
L 2, 4 Group HEU-18 1.118 -
i 2 Group HEU-18 -
1.115 L
( Information provided in Table 8.3 shows that the EXTERMINATOR and 2DB calculated values of k-effective agree closely. This provides confidence that the EXTERMINATOR models are properly set up.
L 8.3 PROMPT NEUTRON LIFETIMES I Prompt neutron lifetimes were calculated by the EXTERMINATOR code using first-order perturbation theory. The prompt neutron lifetime is the weighted integral time that the average neutron lives before being captured or leaked from the core. EXTERMINATOR l
l
. w I
i 77 computes the prompt-neutron lifetime by mumerically solving the following equation: '
,. Jcor,t&*(r,E)kd(r,E)dE]dV _
[
l fcore({d*(r,E')x(E')dE'{vr(r,E)d(r.,E))dV p f kefr ,
where d*(r,E) =' adjoint flux at location r and energy E, 1 = reciprocal neutron velocity at energy E, v(E) ,
d(r,E) = neutron flux at location r and energy E, xp(E') = fraction of prompt neutrons born at energy E',
u = prompt neutron production fu tor for uranium-235, and p Ir(r,E) = macroscopic fission cross section at location r and energy E.
g.
Prompt neutron lifetimes have been calculated for both LEU-22 and HEU-18 fuels. Both graphite and water reflected cores were evaluated. The prompt neutron lifetimes associated with the LEU-22 4-by-5 graphite reflected core with control rod insertion were also evaluated. The results of the prompt neutron lifetime '
analysis are presented in Table 8.4.
Table 8.4 shows that the prompt neutron lifetime for the LEU-
. 22 core is about 15% lower than that of the HEU-18 core. This is as expected because of the significantly higher U-238 loading of the LEU fuel. Calculated pronpt neutron lifetimes for the water I
h ,
l J
78 reflected cores are about 20 percent lower than for the graphite reflected cores. This is due to the higher leakage and absorption
. associated with the water reflector.
Table 8'.4. Prompt Neutron Lifetimes for the UVAR Core.
. . Core Descriotion 1 (usecs)
- 1. All-Graphite Reflector a) LEU 67.0 b) HEU-18 78.8
- 2. Water Reflected ;
a) LEU-22 53.0 I: b) HEU-18 64.4 <
- 3. Graphite Reflected LEU-22 Single Rod Inserted 67.0
- 4. Graphite Reflected LEU-22 All Rods Inserted 66.5 Rod effects wo. .k'.9rmined by inserting two group rod cross sections into the rod :.ste locations. The two group cross sections were input as EXTERMINATOR groups 3 and 4. Although using two group cross sections in a 4 group model is not 1
- . completely rigorous, this method does provide insight into the general behavior of prompt neutron lifetime with control rod l'
iI
'~
~ L. -. .- .: . . _ . , -
LI..a .
l l
l 79 insertion. The prompt neutron lifetime is not significantly affected by rod insertion.
8.4 EFFECTIVE DELAYED NELITRON FRACTIONS l
The EXTERMINATOR code calculates' effective' delayed neutron ;
fractions based on input delayed-group neutron energy spectra.
The effective delayed neutron fraction represents the importance :
, weighted production rate of delayed neutrons divided-by the total l importance weighted production rate of all neutrons. The l
EXTERMINATOR code computes the delayed neutron fraction using the following equation:
f i
fcoreI k Xdi(E')d*(r,E')dE'$sivrf(r,E)(r,E)dE]dV
~
p .
~
fcorekX(E')*(r,E')dE'{ur(r,E)(r,E)dE]dV I p f where xdi(E') = fraction of delayed neutrons in the ith precursor group born at energy E', and i
pj - delayed neutron fraction for the ith precursor group.
Other terms in the above equation are as defined in Section 8.3.
Several references were consulted to determine the delayed neutron energy spectra. Most references (e.g. [30], (31], and (32]) agree that delayed neutrons are born at energies between 200 and 650 kev. This range is clearly bounded by EXTERMINATOR group 3
. . - . . ~
73 80 2 (5.5 kev to 820 kev). Duderstadt [33] presents information that implies that.some delayed neutrons may be born at energies in EXTERMINATOR group 1. In order to check the sensitivity of the l ,
calculation:to reasonable assumptions of the delayed neutron spectrum, the following delayed neutron spectra were evaluated: -
E. j 0 , 0.1 , 0.2 1 0.9 0.8 x1 = 0 x2
= 0 x3
= 0 0 0 0 ;
i .
The above energy spectra are believed to bound the actual situation. ,
Four other delayed neutron energy spectra were also evaluated to provide insight'into the behavior of neutrens born at different
. energies. It should be noted that these spectra are not considered representative of the physical situation and were evaluated for purely academic purposes. These spectra are defined below:
'I.
- 1 0
0
, 0 0
, 0 0
0
, 0.739 0.261 0 0002 x4 = ys = 1 X6 = X7 =
I O O
1 0
The x3 spectrum is identical to the fission spectrum presented in Table 8.1.
Table 8.5 presents total effective delayed neutron fractions associated with each energy spectrum.
g
, = ..:.... u. . .: =
q l
1 81
- s Table 8.5. Total Effective Delayed Neutron Fractions for LEU-22 ,
Fueled Core Enerav Soectrurr Beff
- 1. Reasonable Spectrum Approximations ,
a) x 0.00736 9 -
b) x,1 - 0.00724 c) X3 0.00713
- 2. " Academic" Spectrum Approximations 3 a) x4 0.00620-3 b) X5 0.00807 0.00949 c) X6- ,
c) x7 0.00650 LL5 l The effective delayed neutron fraction appears to be relatively insensitive to reasonable spectrum approximations. A l shift of about 3% is associated with allowing twenty percent of the delayed neutrons to be born in group 1. Because the x1 spectrum was consistently supported in the literature, it was chosen as the preferred spectrum. Based on information provided l- in Table 8.5, a value of peff = 0.0074 is recommended for use with the LEU-22 fuel . An estimated absolute error of about 5% is l considered appropriate; however, it should be noted that the precision is of the order of 1%. ,
g-1 L
o r
_ c . _
LU c: ,
O 82 I L
- -)
j The academic spectrum assumptions show that peff increases as the energy of the delayed neutrons decreases, as expected. A
^
maximum value of peff = 0.00949 is obtained when all of the delayed neutrons are born in the thermal group. The x7 spectrum,
)
- which is identical to the fission spectrum (see Table 8.1),
produced peff = p - 0.0065, as expected. +
To gain insight into the leakage effect, peff was calculated for the LEU-22 fuel using an all-water reflector and a x1 ,
spectrum, peff associated with this case was 0.00752 which is.
- about 2% higher than that calculated for the graphite-reflected I case. This is reasonable considering that more prompt neutron leakage and absorption occurs with a water reflector. It should be noted that this value is well within the estimated 5% error associated with the above-recommended value of peff.
. The calculated value of peff with a single rod inserted was
~
0.00738 which is not significantly different from the recommended value. The calculated value of peff with all rods inserted is
+ 0.00741, which is also very close to the recommended value.
Finally, peff was calculated for the HEU-18 fuel core. The h
calculated value is peff = 0.00738, which is not significantly different from the recommended value for the LEU core.
I I
I: .
-r tg .g -,ew.+v-----. --
v -
)
N '
l 83 Jg 8.5 Summary of Delayed Neutron and Prompt Neutron Lifetime Study In summary, the following conclusions have been made based on t the information presented above:
- 1. The prompt neutron lifetime associated with the LEU-22 core is 67 psecs, which is about 15% lower than the value calculated for the HEU core. ,
- 2. Prompt neutron' lifetimes in the water reflected cores are about 20 percent lower than in the graphite reflected cores.
- 3. Prompt neutron lifetime is not significantly affected by -
" I control rod insertion.
- 4. A value of peff = 0.0074 1 5% is recommended for use with the LEU-22 core. Additionally, peff appears to be relatively insensitive to reasonable assumptions in the delayed neutron spectrum.
- 5. The value of peff does not change significantly with reflector options or insertion of control rods. j
- 6. The geff value associated with the HEU core is slightly
- higher, but not significantly different from the value
-' asscciated with the LEU core.
I I
I , . - - - - -
t
_-m,....+w. , -
t I
84'
9.0 CONCLUSION
S Several key neutronics issues have been addressed by the studies presentednin this thesis. Study-specific summaries have I been presented with each study.
1 RESULTS Control rod reactivity worths for the LEU fueled cores were found to be essentially the same as for HEU fueled cores. Control f
L rod worths were also found to change slightly with core burn-up and core reflector changes implying the need to re-evaluate 1
control rod worths periodically with burn-up and following l
.a g
reflector changes.
' - Peaking factors associated with LEU fuel were somewhat higher i than for HEU fuel. Higher peaking . factors and the relatively small water gap associated with the LEU-22 fuel have led Hosticka i
[4] to recommend adopting more conservative UVAR Limiting Safety System Settings. !
Temperature and moderator void feedback effects for both the LEU and HEU fueled cores were found to be very similar except for l the doppler effect which was significantly higher in the LEU fuel 1
due to the higher uranium-238 loading.
Reflector worths for several different water and graphite reflector options were quantified for the LEU fueled cores. The reactivity worths of deuterium tank reflectors were found to be
. - . ~ ,
t
- E 1 I
l 85 similar to graphite reflectors.
- An effective delayed neutron fraction of 0.0074 was determined ;
for both the LEU and HEU fueled cores. A prompt neutron lifetime of 67 micro-seconds was determined for the LEU fueled core which was about 15 percent lower than in the HEU fueled core.
CONCLUSION AND FUTURE WORK
,I '
The overall study results show that the LEU-fueled core will behave very similarly to the HEU core currently being used.
Future work on the project will include preparation and submission of a revised Safety Analysis Report and UVAR Technical Specifications to NRC, plus responding to NRC licensing concerns.
I-After the license renewal process is complete, an experimental data collection phase is envisioned for the initially loaded LEU core. Experimental data collected during this phase will provide information needed to " fine-tune" the UVAR neutronics models and will also address possible safety concerns such as
-I verifying bypass flow rates, fuel element flow distributions, and radial and axial flux maps. Additional experimental data to be collected on the clean core are recommended by Meem (34].
Additional investigation of radial peaking factors using h
transport theory models may be useful in verifying the values I presented in this thesis which were predicted by diffusion theory J
codes. Further investigations into the experimentally determined and calculated reactivity worth discrepancy for the regulating rod I
i, -
I W 86 ;
!- may be warranted and should include an evaluation of rod regton boundary conditions. New measurements of moderator temperature i coefficients and quantification of the corresponding effective core temperature rise are also warranted.
I-t I l il .
I I
I 3:
LI LI I-
I. :
I !
I l
87 !
l REFERENCES 1.. Rydin, R.A., D. Freeman, M. Fehr, B. Hosticka, R. Mulder, i Conversion of the University of Virainia Reactor to low l g
g' Enrichment Fuel, Second Year Progress Report, University of Virginia Department of Nuclear Engineering and Engineering
.l i
Physics, February, 1989. l
- 2. Fehr, M.K., Desian 0)timization of a low Enrichment University !
of Virainia Nuclear Reactor, Masters Thesis, University of I Virginia, School of Engineering and Applied Science, January i I. ; 1989.
- 3. ' Freeman, D., OA Packaae For Control Rod Worth. Shutdown Marain )
and Excess Reactivity Study, Draft 3 University of Virginia L
I~ Internal Report, August 9,1988.
l E3 4. Hosticka, b.. " Thermal and Hydraulic Analysis of LEU-18 and 22 Plate Fuel" Memo to the University of Virginia Reactor Safety l5 Committee, October, 1988.
- l h 5. UVAR Desian and ' Analysis Handbook, University of Virginia Department of Nuclear Engineering and Engineering Physics, Edited by J.P. Farrar, December 1983.
- 6. UVAR'Gridplate blueprint.
- 7. National Organization Of Test, Research, And Training 1
-~ I Reactors, Specification for Test Research Trainino Reactor LEU Silicide U3Si2 Fuel Plates,TRTR-6, Revision 4, May 20,1988.
1.
l_ 8. National Organization Of Test, Research, And Training Reactors, Specification for Low Enriched U Metal For Reactor l' Fuel Plates,TRTR-1;, March 12,1987.
- 9. National Organization Of Test, Research, And Training Reactors, Specification For Reactor Grade Uranium Silicide ,
I V3Si2 Powder,TRTR-14 Revision 2, July 1, 1987.
- 10. National Organization Of Test, Research, And Training
~
Reactors, Specification For Aluminum Powder For Fuel Plate l(
Core Matrix,TRTR-15, Revision 2, July 1, 1987.
^
- 11. EG&G Drawing 42264, Test Research Trainina Reactor LEU Fuel for UVAR 22-Plate Element, Revision 8, Dec. 5,1986.
I'
.I i
88 I 12. Freeman, D., " Review of Literature.Provided by L.K. Seymour's Letter of January 10, 1989", UVA Internal Memorandum, February 23, 1989. i
- 13. Freeman, D., " Preliminary Design of 18 and 22 Plate Fuel Elements, Memo to the University Of Virginia Reactor Safety
-Committee, April 20, 1988.
? 14. Freeman, D., "Preliminar'yDesign of 18 and 22 Plate LEU j Control Rod Fuel Elements," Memo to the University Of Virginia Resctor Safety Committee, May 6, 1988.
- 15. United Nuclear Corporation, "UVAR Fuel Element" Blue Print .)
- D-401580, July, 1969.
I
! 16. EG&G, "UVAR Control Rod Element" Blue Print #409339, October, 1 l
1978.
, 17. Barry, R.F., LEOPARD A Spectrum Deoendent Non-Soatial a Deoletion Code," WCAP-3269-26, Westinghouse Electric i Corporation, September, 1963. !
IP 18. Bennett, C. L., Purcell, W. l.. BRT: Battelle-Revised-Thermos, ,
Battelle Memorial Instiute, Reactor Physics Department, June, ,
,- 1970.
- 19. Carter, L.L., C. Richey, C. Hughey, GAMTECll: A Code for ;
Generatina Consistent Multiaroup Constants Utilized in Diffusion end Transoort Theory Calculations, Pacific Northwest !
Laboratory, March, 1965.
- 20. ORNL-4078, EXTERMINATOR-2: A Fortran IV Code For Solvina Multiaroup Neutron Diffusion Ecuations In Two Dimensions, Oak '
Ridge National Laboratories,1978.
on 2DB Code), Version #10, September, 1986.
- 22. Wasserman, S., Effective Diffusion Theory Cross Sections For UVAR Control Rods, Masters Thesis, University of Virginia, School of Engineering and Applied Science, June,1989. ;
- 23. Freeman, D., " Review Documentation of S. Wasserman's Work", !
University of Virginia Internal Memorandum, April 12, 1988.
- 24. Farrar, J.P., "Results of Testing For 18 Plate Texas A&M Fuel Elements In UVAR Reactor," Letter to Nuclear Regulatory Commission, June 17, 1975.
.I E
i I
i 89 ,
. 1
- 25. Private Communications with J. Matos, Argonne National +
Laboratories, Winter,1988. ,
- 26. Sternberg, H.I., Thermal Power Calculation and Correlation of I UVAR by Foil Irradiation and Heat Balance, Masters Thesis, University of Virginia, November 1981. -
- 27. Freeman, D., " Temperature and Moderator Effects for the UVAR Core." University of Virginia Internal Memorandum, May 31, 1988.
-I 28. Freeman, D., " Reactivity Effects Associated With Various UVAR Core Reflector Options", University of Virginia Internal L Memorandum, March 16, 1989.
- 29. Brady, M.C., R. Talmadge, W. Wilson, "Few-Group Analysis of Current Delayed Neutron Data", ANS Transactions, Volume 53,
, 1988.
1 1 30. Argonne National Laboratories, Rea: tor Physics Constants, ANL-5800, Second Edition, 1963.
I 31. LaMarsh, J. R., Introduction To Nuclear Reactor Theory, Addison-Wesley Publishing Company,.1965.
- 32. Keepin, G.R., Physics Of Nuclear Kinetics, Addison-Wesley-Publishing Company, 1965.
- 33. Duderstadt, J.J., L.J. Hamilton, Nuclear Reactor Analysis, John Wiley and Sons, Inc., New York, 1976.
- 34. Meem, J.L., "Possible Experiments On the New LEU Core", Memo to University of Virginia Reactor Stfety Committee, October 20, 1988.
1' l
I l
LB La I
z ;.: .
n' i
90 1
APPENDIX A. MATERIAL CROSS SECTIONS FOR UVAR 2-0 CORE MGDEL
'I This Appendix presents material cross sections for the UVAR 2-D Core model, Tables A-1 through A-3 present core region cross sections developed with the LEOPARD code. Table A-4 presents control roo cross sections for both the safety and regulating rods using the various fuel options. Table A-5 presents ANL calculated 1
water and graphite reflector cross sections and Table A-6 presents l
l UVA developed graphite reflector cross sections.
1.
1 I
I -
l' I
i L l~
F=w w
4 g: g gg g g ~ ' M M. : E E E E '. '
4-Table A-1 Core Region Cross Sections for LEU-22 Fuel (cm'I).
Description Group If Ia If Itr _Ic.g [(g-1 ) .g
- 1. Two Group Cross Sections
- 1. Fuel Lattice 1 2.32x10j 6.22x10 5.71x10-3 2.22x10'I 0.0 2 9.32x10' 1.26x10- 2.25x10~I 2.54x10[
1.26x10 1.13x10 2.57x10-2
- 2. Fuel Non-Lattice 1 0.0 #
3.83x10'2 0.0 1.89x10-I 1.75xt0-I 0.0 2 0.0 1.31x10- 0.0 6.15x10'I 6.02x10-I 1.38x10 -2
- 3. Control Rod Element 0.0 # t 1
0.0 I Non-Lattice 2 0.0 4.08x10'2 1.67x10- 2.70x10f 2.32x10f 0.0 0.0 1.61x10 1.59x10 3.76x10-2 II. Four Group Cross Sections
- 1. Fuel lattice 1 6.91x10'# 3 3 I
.i
' 2 1.11x10 # 1.91x10 # 1.68x10~I 8.46x10-2I 0.0 4.05x10 3 8.83x10'2 9.89x10 8.21x10-2 3
6.47x10'2 1.84x10- 1.57x10-2 3.04x10'I 4.09x10'0 2.09x10'I 3.05x10 4 9.32x10- 1.26x10"I 2.25x10^ 0 9.44x10-2 1.26x10 1.13x10 8.57x10-2
' 2. Fuel Non-Lattice 1 0.0 3.47x10~4 0.0 1.42x10-3 2 0.0 # 1.00x10-I 0.0 1,25x10'# 0.0 2.28x10-I 1.80x10-l 4.15x10 3 0.0 0.0 2.26x10 d
4 0.0 7.35x10'2 1.31x10' 1.79x10~I 4.77x10-2 0.0 6.15x10~I 6.02x10~1 4.58x10-2 e
a
- 3. Control Rod 1 0.0 ,
3.77x10~# 0.0 1.76x10-I Non-Lattice 2 0.0 3.88x10-5# 7.69x10^-2 0.0 0.0 3.16x10^ 1.93x10 9.89x10-2 3 0.0 0.0 4.77x10^ 3.53x10~I 4 0.0 8.67x10'2 1.67x10- 0 1.22x10~I I 0.0 1.61x10 1.59x100 1.23x10 3 -
s
,, .v . - y- ,-- - - , , os -_ _au_ ___ _ _ _ _ _ - _ . _ ' m__m- 1 d
M M M Wm MM M M M W W M M M Imj 1m TW i
r I
i Table A-2. Core Region Cross Sections For LEU-18 Fuel (ca~I).
Description Group bf Ia Y f Itr g,o (o-1) g
- 1. Two Group Cross Sections
- 1. Fuel Lattice 1 1.93x10-3 5.41x10-3 4.75x10-3 2.61x10'I 2.26x10-I 0.0 2.92x10 -2 0 0 2 7.90x10-2 1.09x10 1.91x10~I 1.40x10 1.29x10 l
- 2. Fuel leon-Lattice 1 0.0 0.0 1.93x10'I 1.77x10~I 0.0 2 0.0 3.92x10'2 1.37x10- 0.0 6.81x10'I 6.67x10-I 1.51x10-2
- 3. Control Rod Element 1 0.0 4.15x10'8 0.0 2.29x10~I 0.0
, Mon-Lattice 2 0.0 1.70x10-2 0.0 2.67x10f' 1.65x10 1.63x10 0
3.81x10-2 l 11. Four Group Cross Sections
-# 3
- 1. Fuel Lattice 1 5.68x10-# 9.80x10 1.57x10 #
1.71x10'I 8.23x10-2 0.0 i 2 3.31x10-4 3
6.87x10-4 8.09xt0 3.09x10-I 2.05x10-I 8.76x10-2 -I ~
3 1.59x10-2 1.30x10-2 4.3Bx10~I 3.25x10-1 1.04x10 5.36x10 2 -
0 0 4 7.90x10 1.09x10-I 1.91x10'I 1.40x10 1.29x10 9.65x10-2 8
- 2. Fuel Non-Lattice 1 0.0 0.0 1.44x10-II 9.99x10-2I 0.0 3.57x10 4 2 0.0 1.22x10- 0.0 2.30x10'I 1.78x10 I 4.34x10-2 2
3 0.0 7.52x10-4
-2 0.0 2.39x10- 1.89x10~
4 0.0 1.37x10 0.0 6.81x10-I 6.67x10-I 5.13x10 4.97x10-2
- 3. Control Rod 1 0.0 3.23x10 0.0 1.74x10-I 7.59x10-2 0.0 '
Non-Lattice 2 0.0 3.89x10-5 0.0 3.12x10-I 1.89x10'I 9.76x10-2 3 0.0 8.75x10~# 0.0 4.79x10~l 3.54x10-I 1.23x10'I 1.65x100 0
4 0.0 1.70x10-2 0.0 1.63x10 1.24x10-I N
,,.,%---c. ,- , - -
..- - *w ,y ,,...__e.-.__m__ --
-a _
W W W W' W W W W W WM W: 6 M-M M M ,W -
-- I Table A-3. Core Region Cross Sections for HEU-18 Fuel (cm-I)
DescriptM Group If _Ia dIf Itr __Ig.g .89-1).c I. Two Group Cross Sections
- 1. Fuel lattice 1 1.60x10-3 3 3.90x10-3 2.58x10
-I 3.01x10 2 0 2.42x10-I 0 0.0 2 6.96x10-2 9.64x10- 1.68x10-1 1.40x10 1.30x10 3.03x10-2
- 2. Fuel Non-Lattice 0.0 -#
1 3.89x10 0.0 1.98x10-I 1.81x10-I 0.0 2 0.0 1.34x10-2 0.0 7.21x10-' 7.08x10-I 1.60x10-2
- 3. Control Rod Element 0.0 -8 d 1 4.13x10 0.0 2.58x10 2.22x10-I 0.0 Non-Lattice 2 0.0 1.63x10-2 0.0 0 0 1.50x10 1.48x10 3.47w10-2 II. Four Group Cross Sections
- 1. Fuel tattice # #
5.87x10 " 8.23x10 -2 1
2.21x10'# 6.02x10 1.69x10-1 0.0 2 4.24x10-# 6.96x10-4 3.05x10-I 2.85x10 2.02x10-I 3 4.65x10-3 8.6Bx10i 1.12x10-2 4.29x10'I 3.21x10^
8.63x10-2 4 6.96x10-2 -1 0 0 1.03x10-1 9.64x10-' 1.68x10 1.40x10 1.30x10 - 9.89x10-2
- 2. Fuel Non-i.attice 1 0.0 #
0.0 3.55x10 # 1.46x10-3 1.00x10-I 0.0 2 0.0 1.18x10- 0.0 2.34x10^ 4.5;x10-2 3 0.0 7.40x10-# 0.0 2.52x10-I 1.80x10^I 1.99x10- 2 4 0.0 1.34x10-2 0.0 5.37x10 2 7.21x10^ 7.08x10-I 5.21x10
- 3. Control Rod 0.0 #
1 3.76x10 5 0.0 1.69x10-I 7.% x10-2 0.0 -2 Non-Lattice 2 0.0 0.0 5.26x10 # 3.01x10-I 1.90x10-I 8.86x10 3 0.0 8.55x10- 0.0 4.44x10-I 3.32x10-1 1.11x10-I 4 C.0 1.63x10 -2 0.0 1.50x10 0
1.48x10 0
1.11x10
-I 4
%_ _._,.--.-__<% m w - - - - - _ --
r-"
, - , 4 -d' Ye- e- **
w.- _--v.. ma _ __.-._._____2s--=.2m -
---e --d
6~M M U ]
4 Table A.4 Control Rod Region Cross Sections Group a tr q-q (q-1)-g Rod /Fual
- 1. Safety Rods -l I
0.0182 0.270 0.573 0.0
- 1. LEU-22 1 1.53 0.0229 2 0.0171 2.00 0.0186 0.271 0.577 0.0
- 2. LEU-18 1 1.57 0.0233 2 0.166 2.05
. 0.0'. 78 0.270 0.573 0.0
- 3. HEU-18 1 1.% 0.0229 2 0.171 2.06 II. Regulating Ro<1 0.00268 0.284 0.552 0.0
- 1. LEU-22 1 0.0163 2 0.0384 1.23 1.37 0.00220 0.284 0.552 0.0
- 2. LEU-18 1 1.39 0.0164-2 0.0425 1.25 .
0.00234 0.281 0.546 0.0
- 3. HEU-18 1 1.37 0.0160 2 0.0368 1.24 ,
i I g-g represents in-group scattering 2 (g-1)-g represents one group downscatter ,
. p
.e s.
g_.
~ W fM 'M - W . W W W W W W W W W W W W- W WE Table A.S. M.t Developed Water and Graphite Reflector Cross Sections (cm-I).
Dascription Group 1
a I tr e-c (9-1)-c I. Water Reflector i 1. Two Group 1 4.58x10-# 2.65x10~I 0 8.39x10-I 4 3.0
- 2 1.97x10-2 2.30x10 2.63x10 4.91x10-2
- 2. Four Group # I 1
3.60x10'6 1.59x10-I 0.0 1
2 3.12x10-I 2.70x10 8.83x1 I I 3 9.68x10 1.00x10-3 1.07x10~I 5 5.69x10*I 0 1.38x1 1.50x10' 4 1.97x10~2 2.30x10 2.63x1 1.40x10
II. Graphite Reflector
- 1. Two Group 1 1.06x10-# I I
' 2.69x10~I 3.64x10~I 0.0 2 3.28x10-3 4.93x10~ 5.21x10- 7.23x10-3
, 2. Four Group 1 8.90x10~ 1.44x10~I 1.85x10-I 0.0 2 2.73x10'I 3.06x10-I 3.75x10-I 3.10x10-2 3 1.75x10'# 3.51x10~I 4.24x'O~I 1.90=10-2 1.5Bx10 -2 4 3.28x10~ 4.93x10-' 5.21x10~I 1
i 1
?
3 u- r .,n.. = - . . * , + - e~. - - - +<w- - - - - - -- u s- - - - - - - - -J -
M M .' O : M M M MMM M ' W ;l
~
M -M M -M M -- 8 M: 8 1:'
i 4
f Table A.6. WA Deveicped Graphite Reflector' Cross Sections (cm*I)
Descriotion Group a Itr be-o -I(o-11-o i Graphite Reflector 1 1.0'x10 2.63x10'I 4.26x10~I '0.0 2 3.ux10-3 6.25x10'I 7.26x10-I .5.92x10-3 I
6 P
l ,
i l
i .
l l
t i
f n
.-.r. ~ # ..e_ , . , . - - . , ,y, , - - , .,
, , - . . .- r- _ _ _ _ _ _ _ ___m _. _m___ _ _ _ _ _ _ _ _ _ . _ - = - - -