3 -

-

" GELEASED FOR ANNOUNCEMENT

. CIS
IN BUCLEAR SCIENCE ABSTRA ABSTRACT

 

 

OAK RIDGE NATIONAL LABORATORY
operated by

 

UNION CARBIDE CORPORATION
NUCLEAR DIVISION
for the
U.S. ATOMIC ENERGY COMMISSION

  

ORNL- TM-497
COPY NO. - {f

DATE -  August 16, 1966

 

ANALYSIS OF FILLING ACCIDENTS IN MSRE

J. R. Engel, P, N. Haubenreich, and S. J. Ball
CFsTI PRICTS

SV i e

———

HC, $¢;Lév; MN ,éfio

Whenever the MSRE is shut down, the fuel salt is drained
from the core. Then, during a normal startup, the graphite
and the fuel are preheated and the control rods are positioned
so that the reactor remains suberitical while it is being filled.
Certain abnormal circumstances could result in criticality and
a power excursion in the partially filled core. Various
postulated incidents were surveyed and the worst case was ana-
lyzed in detail. This case involved selective freezing in the
drain tanks to concentrate the uranium in the molten salt
fraction. TPhysical restrictions on the fill rate and safety
actions of control rods and gas control valves limited the
calculated power and temperature excursions so that any damage
to the reactor would be prevented.

NOTICE This document contains information of o preliminary nature
and was prepared primarily for internal use at the Oak Ridge National
Laboratory. It is subject to revision or correction and therefore does
not represent a final report.
 

e oy

e = e e

 

 

 

—~ Lecnumnce

This report was pupcrod os an account of Govommcm :ponsorod work., Neither the Unhed S!a!u,

_ nor-the Commission, nor any person acting en behalf of the Commission:

A, Makes any warranty or representation,’ uprus-d or implied, with respect to the cecum:y, h

complflonus, or usefulness of the ln!ormcfloa cmtctned in this report, or that the use of
any information, epporatus, mflrod or prccou disclosed in this report may not infringe
ptivately owned rights; or '

-B. Assumes any licbilities with nspoct h Ihc use of, or for domages nnulflng from the use of

any information, apparatus, method, or process disclosed in this report.
As’ vsed in the above, “person acting on behalf of the Commission® includes any cmployu or
contractor of the Commission, or employes of such contractor, to the extent that such employes .

or contracter of the Commission, or employee “of such contractor prepares, disuminuies, or

“provides acceas to, any information pursuant to llll omplflymm or contract whh the Commlulon,

or his amployment with such contractor.

 

 

 
Lo

PREFACE

. The analysi
written ea . s described
Publishingflihln 196k, Althizr; was done in 196
part of the Ms;ETGPOTt is being there has been 3 énd the report
analysis. Well safety Circuii Lesued as Writtean inordinate delwas
able to assess over two year ry was based on n because a si 'ay in
ss the credibilit s have now pass n the results of gnificant
y of some of thee:’ and now we aret:is
ssumptions invol etter
ved. We

Pla.n +
o do th.
ig ag
a, .
critical revi
lew Of
the s
afet
y syste
bt ’m.

RELEA )
SED :
FC .
MENT

IN NU
ICLEAR
SCIENCE ABST
TRACTS

 

 

LEG
Thnis ¥ O P red 28 8% aocount of Gwemmen gpon
1ates, e jpaiof any PeT on acting naif e
Mak rranty repres® on, eXP e8 or p\*.ed, ih resp® 1o th
racys &° ple a8, ©F ueefud a of the iniormation ine report {hat the use
of any ifoT ati appar L0, ethod oT 58 108€ thi report ¥ not nfringe
private‘nj owned vightBs
B. AsS y i jlitied pect use of mages resulting {rom h

ase of 39 info’ mations pparat Lnods pro 88 diecio® in this £,

g use o above ¢ acting 1f of W miseion" inciuded y
ployee or 5 the CO jon, OF mploye® of n comtT Ah at that
guch em loye nirac ommlest , ot e loyee guc a repared:
cnssem‘ma rovides c y informa parsud pie € b cont:a.ct

{th e jon, OF toyme g wit! ontractor
 
CONTENTS

Introduction.
Mechanics of Filling .

2.1 General Description . . .

2.2 System Volumes

2.3 Drain Tank Pressure During a Fill.
2.L  Amount of Gas in Tank During a Fill .
2.5 Gas Supply.

2.6 Course of a Normal Fill

2.7 Coast-Up of Fuel Level .

Nuelear Considerations . . . .
3.1 Effect of Fuel ILevel on Reactlivity .
3.2 Control Rods

3.3 Temperature Coefficients of Reactivity .

Survey of Possible Filling Accidents .

4.1 Filling with Control Rods Withdrawn .
h,2 Filling with Fuel at Low Temperature .
4,3 Filling with Concentrated Fuel .
Analysis of Maximum Filling Accident .

5.1 ©Specification of Accident.

5.2 Preliminary Digital Calculations .

5.3 Detailed Analog Simulation .

5.4 Discussion « o« ¢ ¢ o 0 0 0 . o

Conclusions .« o« v v o o o o o o o

O O &= & F o

. 11
Fig. No.

= w

10

11

12

13
1k

LIST OF FIGURES
Title

System Used in Filling Fuel ILoop . . .
Calculated Volume Calibration of Fuel Loop .
Calculated Volume Calibration of Fuel Drain Tank .

Drain Tank Pressure vs Liquid Ievel in Fuel Loop
at Various Salt Fill Rates . . . « . . . . . . . .

Amount of Gas Required in Drain Tank vs Salt Ievel
in Loop at Several Fill Rates

Salt ILevel in Fuel Loop During Normal Fill.

Coast-Up After Gas Addition is Stopped vs Fill Rate
at Several Initial Fuel Ievels . e s e e e .

Effect of Fuel Level on Reactivity .
Contrcl Rod Worth vs Position . . . . .

Liquid Composition Resulting from Selective Freezing
of Fuel Salt "A" in Drain Tank .

Bilock Diagram of Mathematical Model for Fill
Accident Simulation . . . . . . ¢ . 4 . . . ..

Cross Section of Typical Graphite Stringer with
Adjacent Fuel Stringers. « + o + o &+ o o o0 o o o « &

Net Reactivity Inserted During Maximum Filling Accident.

Power and Temperature During Maximum Filling Accident.

.10

.12
. 1h
. 15

. 21

25

26

. 29
. 30
ANATYSTIS OF FILLING ACCIDENTS IN MSRE

J. R. Engel, P. N. Haubenreich, and 5. J. Ball
1. INTRODUCTION

One of the features of the MSRE (and fluid-fuel reactors in general)
is that it can be positively shut down by draining the fuel out of the
core. The control rods provide a small shutdown margin to take the re-
actor subcritical whenever desired, but for any shutdown in which the
reactor is to be cooled dowh, the fuel must be drained. Draining and re-
filling the core is therefore an coperation which will probably be done
many times.

The normal procedure for a startup requires that the reactor and
the fuel be heated by electric heaters to near operating temperature
before the fuel 1s transferred from the drain tank to the core. The
rods normally are partially inserted so that the reactor 1s just sub-
critical at the fill temperature when the fuel fills the core. Criti-
cality is attained by withdrawing the rods after the fuel and cooclant
loops are filled and circulation has been started.

Tt is conceivable that criticality could be attained during a fill
before the core is completely full. This could result from one or more
of the following abnormal conditions: 1) the control rods are withdrawn
too far, 2) the temperature of the fuel and/or the core graphite is too
low, and 3) the uranium concentration of the fuel was increased (or the
poison concentration was decreased) while the fuel salt was in the drain
tank.,

If criticality is reached prematurely, the nuclear power will rise,
possibly causing damaging temperatures in the partially filled core. As
soon as the onset of such an undesirable situation is detected, the rocds
are inserted, the fill is stopped and the fuel returned to the drain tank.
Supplementing the effects of these actions will be the reactivity feed-
back from any changes in the fuel and graphite temperatures.

In order to evaluate the severity and consequences of various postu-
lated fiiling accidents it is necessary to have certain guantitative infor-

mation. This includes: 1) the filling rate (fuel level ve time),
10

2) the relation of kepp to fuel level for the particular abnormal situation —
being considered, 3) rod worth and the speed with which they can act in
the partially filled core, 4) how rapidly the fill can be stopped (level
ve time after action is taken to stop the fill}, and 5) temperature coef-
ficients of reactivity appropriate for this abnormal situation. This
information has been developed for a variety of cases and is presented
in Sections 2 and 3.
The relative severity of a number of postulated accidents is sur-
veyed in Section 4. In Section 5, the most severe of the postulated
accidents is analyzed in considerable detail. Conclusions are summarized

in Section 6.

2. MECHANICS OF FILLING

2.1 General Description

Figure 1 is a simplified flowsheet of the reactor fill, drain, and
vent systems showing only those features which are essential to a de- "
scription of the normal fill and drain procedures. All valves are shown
in the normal positions for filling the reactor from fuel drain tank
No. 1 (FD-1).

The reactor is filled by admitting helium, from a supply at 40 psig,
to the drain tanks to force the fuel salt up through the reactor drain
line into the primary loop. The fill rate is limited by a restriction
in the gas supply line and the maximum level in the loop is set by limit-
ing the pressure with PIC-517. Helium displaced from the loop by the
incoming fuel is vented from the pump bowl, at the high pcint in the loop,
through the auxiliary charcoal bed to the stack. (This vent route by-
passes the main charcoal beds to avoid the elution and release of xenon
and krypton which may be in those beds.)

When the fill is complete, salt is frozen in the drain line at
FV-103. The pump-bowl vent is switched to the main charcoal beds and,
after the drain-tank pressure has been vented through the auxiliary char-
coal bed, the pump-bowl and drain-tank gas spaces are connected through
HCV-54L. The system is now in readiness for operation, with the only

action required to drain being to thaw FV-103.
FHX

 

1l

 

ORNL-~DWG 63-7320

£
< s
FP PCv-522
533 XD o
P TO FILTER,
@ Y/ FAN, AND STACK

L\

‘~ \/ \/

2\ A

FFT
REACTOR FD-2 AUXILIARY
CHARCOAL
FFT BED
[} FD-2
Fvoos FFT r 5
FD-2 HCV- HCV-573
544 A~
“ .y HCV-8T72 -
[><] 7 Z % 50-psig
£ DISK
Fo-2
()
e
3
X [

PCV

517
40- psig
HELIUM
SUPPLY

Fig. 1. GSystem Used in Filling Fuel Loop.
12

2.2 System Volumes

The relation between the volume of salt transferred from the drain
tank and the liquid level in the fuel loop is shown in Fig. 2. The
datum level is the bottom of the reactor vessel {elevation 826.92 ft).
The midplane of the 65-in.-high graphite matrix in the core is at 3.75 ft
and the operating level in the pump bowl is &t 13.4 ft. The first 1.5 £t~
of salt transferred from the drain tank is required to fill the drain line.
The liquid volume-level relation in a fuel drain tank is shown in
Fig. 3. The volume of the tank above the salt level is the difference
between 80.5 £t° (the total volume of the tank) and the liguid volume.
The datum plane for Fig. 3 is 12.1 ft below that for Fig. 2 (81k4.82 tt).

2.3 Drain Tank Pressure During a Fill

The pressure difference between the drain tank and the pump bowl is
determined primarily by the fuel-salt density and the difference in the
levels in the loop and in the tank., Pressure drop in the fill line adds
to this pressure difference when fuel is flowing. The pressure in the
pump bowl is normally above atmospheric pressure because of pressure drop
in the gas vent from the pump bowl to the stack. (The pressure approaches
1 psig at zero flow because of a check valve in the vent line.) Figure 4
shows the relation between the actual pressure in the drain tank and the
level in the fuel loop during a fill. This figure is based on a total
salt inventory of 73.2 ft® and a salt demsity of 130 1b/ft>. The increased
pressures at salt fill rates of 1 and 2 ft3/min refiect the pressure drop

in the fill and vent lines.

2.4k Amount of Gas in Tank During a Fill

The amount of gas in the drain tank is a function of the temperature,
the pressure and the volume in the tank above the liquid. Figure 5 shows
the amount of gas in the tank (volume of gas at 32°F, 1L.7 psia) as a
function of loop liquid level for varicus salt fill rates. (The gas in the
tank was assumed to be at 1200°F. The pressure and the actual volume

were obtained from the information already discussed.)
13

ORNL DWG 66=7767

Tend ut T3Ad]

 

Salt Volume (ft3)

*

" Fuel Loop

ibration of

Calculated Volume Cal

. 2.

Fig
14

ORNL DWG 66-7768

 

(c33) susy uTsag wp sumrop 3Tsg

Salt Ievel in Drain Tank (ft)

Tank.

in

f Fuel Drai

1611 ©

alculated Volume Calibrati

C

3

ig.

E
15

URNL TWG €6-7767

 

 

 

 

 

 

 

 

i

 

~FILL RATE
{ctm)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

......

 

 

 

 

 

 

 

 

 

e e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1k

10

 

 

 

(81sd) UNSSTUI MNVI NIvEd

SALT IEVEL IN PRIMARY IOOP (ft)

Drain Tank Pressure vs Liquid Level in Fuel Loop at Variocus

Fig. 4.
Salt Fill Rates.
16

ORNL DWO 66="7710

 

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

  

 

 

 

 

 

 

 

 

 

. O )
' !
s . T, M
¢ . -4 |
: ; . .
) . : \ .
el e Ll Lol i
1 ! B e
i ! I :
i “ ..........
b i
. 4 - . -8 - . i
: i . . T ooy i
: ! | ! . LTI
Lo oo R
—— |<F1A§u b - e e ——
. i . TLTTTTIYTTI
1 + 4 - ~
i o ! . 1
. . i 1o
— T y T
! : “
3 “y ..
I ! .
AT i Ll
@ o ..m.
! i }
: ' HN
3l : L
o “ L : oot
Voo . P
: : Dol . Ll :
iyl el el il il il
it : : ! T TTTTTITTIATTTT
: : ; ! [
| : 1 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

   

   

 

 

 

 

0
0
O

(¢33 P3g) Husl uT=ld Tond Ul SBD

&

 

Toop (ft)

in Primary

Salt level

Amount of Gas Required in Drain Tank vs Salt Level in Loop

at Several Fill Rates.

Fig, 5.
17

2.5 Gas Supply

The gas addition rate, and therefore the salt fill rate, 1s Limited
by the gas addition system. This system was designed to avoid cver-
filling the fuel loop and to limit the severity of possible filling acci-
dents, the latter by restricting the fill rate. The first objective
wag achieved by limiting, with PIC-517, the pressure which can be put on
the drain tank to a value just sufficient to attain the desired level in
the pump bowl at tfie end of the f£ill. The fill-rate limitation was
attained by installing a capillary restrictor in the supply line and
limiting the primary ges supply pressure. (Details of this limitation
are described later in this report.)

During a normal fill, the controller PIC-517 will be set to give a
drain-tank pressure no greater than that required to just bring the salt
to the desired level in the pump bowl. Throughout most of the filling
operation, PCV-517 is wide open, because the tank pressure is well below
the setpoint pressure, and the controlling resistance is the capiilary.
Only when the drain-tank pressure approaches the setpoint does PCV-35LT
function to stop the gas flow.

2.6 Course of a Normal Fill

The f£ill rate, or level vs time, is determined by the combination of
the relations among level, fill rate, and amount of gas in the drain tank
and the characteristics of the gas supply system {gas addition rate vs
drain tank pressure). |

Figure 6 shows the predicted level in the fuel loop as & function
of time during a normal fill. For this prediction, done with the aid of
an analog computer, the gas supply was assumed to be at 40 psig. The
capillary restrictor was sized to give a salt fill rate of 0.5 £t°/min
when the level is at the core midplane with a gas supply pressure of
50 psig. The time required to complete the entire fill was 3'3/h hrs.
(The small delay at the beginning is the time reguired to £ill the drain

line.)
18

ORNL DWG 66-7771

 

12[

10

(33)

 

<« 0 =t

dooy Axsmrag ul ToAT] 3TeS

2ko

 

calt Level in Fuel Loop During Normal Fill,

Fig. 6.
19

2.7 Coagt=-Up of Fuel level

A filling operation can be interrupted at any time by any one of
three actions: venting the drain tank through the auxiliary charcoal
bed, equallzing draln-tank and loop pressures through line 521, and
shutting off the gas addition to the drain tanks. In an emergency alil
three of these would be done or attempted simultaneously. Either of the
first two actions would not only stop the fill but would allow the salt
in the loop to drain back to the tank. The third action, stopping zgas
addition, would be used if it were desired to hold up the fill at any
point.

Simply stopping gas addition does not immediately stop the salt flow
into the loop. This has important implications, for if the gas addition
is stopped while the level is rising through the active part of the core,
the level and the reactivity will continue to increase for some time
after the gas is shut off.

As can be seen from Fig. 5, the amount of gas in the tank at a
given loop level and an appreciable flow rate is capable of supporting
the salt at a considerably higher level when the fill rate has gone to
zero. The amount of level increase in the loop after gas addition is
stopped (coast-up) is a function of the initial level and flow rate.
Figure 7 shows this relationship for several initial salt levels in the
core. Analog calculations of the coast-up transient indicated that the

level approaches the final value with a time constant of about 1.4 min.

3. NUCLEAR CONSIDERATTIONS

The occurrence of a fill accldent presupposes the existence of an
abnormal condition which causes the reactor to be critical before the
core is completely filled with fuel salt. The basic nuclear character-
istics which must be considered are qualitatively similar for all of

the accidents examined.
20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ORNL DWG 66-7772

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fill Rate (cfm)

topped vs Fill Rate at

o
2

as Addition is

e
T

Coast-Up After

Several Initial Fuel Levels.

Fig. 7.
21
3.1 HEffect of Fuel ILevel on Reactivity

In all of the filling accidents, the system multiplication constant
is increased above unity by the continued flow of fuel into the critical
core. The shape of the reactivity curve as a function of fuel level,
along with the rate of fuel addition, determines the rate of reactivity
increase.

The effect of fuel level on multiplication was calculated for a
simplified model of the reactor which considered only the uniform-
channelled portion of the graphite-moderated region. Figure 8 shows
the effective multiplication constant as a function of the fraction (H/L)
of this region filled with fuel. The curve was normalized to a k of
0.997 for the full region; this 1s the target multiplication constant in
the MSRE for a fill under normal circumstances. The fractional levels,

O and 1 in this model correspond to actual fuel levels of 1.33 and

6.50 ft, respectively, in the MSRE, (See Fig. 2) Since the multipli-
cation constant in the MSRE rises above zero before the fuel level reaches
the channelled region and continues to rise until the entire vessel is
full, the curve in Fig. 8 indicates a greater differential change in re-
activity with level than actually exists. Thus, use of this curve shape
in the accident analyses leads to slight overestimates of their severity.
The flattening of the curve in Fig. 8 is due to the fact that the graphite
in the core is stationary. At low fuel levels, the graphite above the
fuel acts as a reflector which affects the multiplication substantially;

this added effect diminishes as the core fills.

3.2 Control Rods

The three control rods in the MSRE are clustered near the axis of
the reactor and enter the core from above. Since the core fills with
fuel from the bottom, the reactivity worth of a partially inserted rod
depends on the fuel level. Figure 9 shows the calculated fractional
worth of & rod as a function of the fraction of total insertion for a
full core and for a core with only 0.72 of the channelled region filled

with fuel. The total worth of the fully inserted rod was essentially the
22

ORNL DWG 66-7773

 

330y

H/L

Effect of Fuel Level on Reactivity.

8.

Fig.
 

 

70056

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ORNL~LR-DW(.

 

 

 

 

 

 

 

oo i

 

 

 

e bt o e

 

o e g -

 

 

 

 

 

 

tion.

 

 

 

 

 

031

e pe e p e el

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

T

23

 

 

 

 

T

 

 

 

 

 

 

 

 

oI

 

 

 

 

 

 

 

 

 

 

e g ]

 

—— -

 

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Control Rod Worth vs P

 

 

 

 

 

 

 

 

9.

 

bbwg v s rmy s

 

 

ig.

 

.

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 
24

game for both cases. The rod worths used in analyzing the filling acci-
dents were 2.9% 8k/k for a sirgle, fully inserted rod at all fuel levels
and 6.7% for all three rods.” The fractional worth of a partly inserted
rod at a given position was assumed to vary liunsarly with the fraction of
the core filled with fuel.

During a normal fill, the three control rods willi be withdrawn equal
amounts and positloned such that the reactor is Just subcritical
(k = 0.997) when completely full. (Protective interlocks will require
that the rodes be withdrawn a minimum distance to iznsure that negative re-
activity can be inserted in the event of premature criticality.) The
rods, poéitioned in this way, will control 4.3% Bk/k in the clean, full
reactor with the operating concentration of uranium in the fuel. They
will, however, control substantially less at the same position in a
partly filled reactor and can, therefore, insert correspondingly more
negative reactivity during a fill acecident.

Each control-rod drive is equipped with a magnetic cliutch so that the
rod can be dropped into the reactor if necessary. A O.l-sec release time
and an accelerating force of 0.5 g were assumed for the rods in the

accident analyses.

3.3 Temperature Coefficients of Reactlvity

The values of the fuel and graphite temperature coefficients of re-
activity depend oz the fuel composition., Table 1 lists the nominal compo-
sitions and the densities of three fuels bheing considered for use in the
MSRE. The temperature coefficients of reactivity shown are for the full
core.

Temperature coefficients of reactivity in the partially filled core
are quite different from the coefficients for a full core. This is due
in part to the difference in effective size and shape of the core. More

importantly, in the partially filled core an increase in fuel temperature,

 

*

Recent calculations indicate the worth of three fully inserted rods
ranges from 5.6 to 7.6%, depernding on the fuel composition., The usable
worth is from 5.2 to T.2%.
25,

Table 1

MSRE Fuel Salts Considered in Filling Accident Analyses

 

 

Fuel Type A B C
Salt Composition:  ILiF" 70 66.8 65
(mole %) BeFs 23.7 29 £9.2
7rF, 5 L 5
ThF, 1 0 G
UF, (approx.) 0.3 0.2 0.8
U Composition Us34 1 1 0.3
(atom %) Uz3s 93 93 35
y236 1 1 C.3
U=s8 5 5 el ki
Density at 1200°F (1b/ft>) 1445 134.5 ibe,7
Temp. Coeff. of Reactivity (°F~ 1)
Fuel -3.0 x 10" -5,0 x 10™°® =3.3 x 107°
Graphite -3.4 x 1075 -4,9 x 107 -3.7 x 107°

 

%99.9926% 147

with its attendant decrease in fuel density, raises the fuel level and
increases the effective height of the core. This is in contrast to a
full core in which a reduction in density expels fuel without changing
the core size. The coefficients are also affected by any difference in
fuel concentrations between the normal and abnormal situations.

The fuel temperature coefficient of reactivity for the partially
filled core was evaluated by comparing Kepe calculated for two cases,
each with Fuel B concentrated by selective freezing of 39% of the salt.*
Tn the first case, H/L was 0.60 and the fuel density was proper for 1200°F. "

 

*
See pp 20 - 21 for discussion of selective freezing.

*¥
This choice of salt type and H/L will be expiained later.
26

In the second, H/L was 0.6l and the fuel density was reduced to keep the
same total mass of fuel in the core., Neutron microscopic cross sections
were evaulated at 1200°F in both cases and the graphite density was un-
changed. These two cases simulated a rise in fuel temperature from 1200°F
to 1341°F in a time so short that the graphite temperature does not rise
appreciably. Use of constant microscopic cross sections implies that the
fuel temperature has no effect on the thermal neutroun energy distribution
when, in fact, it does. Results of these calculations gave a 8k/k of

0.061%, equivalent to a fuel temperature coefficient of -0.43 x 1075 °F 1,

L., SURVEY OF FILLING ACCITENTS

The relative severity of filling accidents can be described, quali-
tatively, in terms of the amount of excess reactivity available for ad-
dition to the core and the rate at which it can be added, particularly in
the vicinity of keep = 1. The amount of excess reactivity available
depends primarily on the circumstances postulated for the accident and
the composition of the fuel mixture. Three sets of circumstances which
can produce filling accidents have been considered; these are discussed
under separate headings below. The influence of fuel composition waes
examined for each type of accident.

The rate of reactivity addition involves, 1in addition to the factors
mentioned above, the rate of fuel addition. In order to restrict the
most severe accident to & tolerable level, it was necessary to 1imit the
salt addition rate under normal circumstances to 0.4 £t°/min. The normal
helium supply pressure to the drain tanks is 40 psig with an ultimate
1imit at 50 psig imposed by a rupture disc. The physical restrictions
which establish the normal fill rate limit the maximum rate tc 0.5 £1t°/min
with the salt level in the main portion of the core. All of the filling

accldents were examined on the basis of the 0.5 f£t°/min rate.

L.1 PFilling With Control Rods Withdrawn

The amount of excess reactivity that can be added in the MSRE by
filling the core with the contrcol rods withdrawn is limited to the amount

required in the fuel for normal, full-power operation. Although this
27

requirement varies somewhat with the fuel mixture it is not expected to
exceed 4% in any case and administrative control will be exercised to keep
the reactivity at or below this value. If the fuel were loaded with
sufficient uranium for 4% excess reactivity and all three rods were fully
withdrawn, the core would be critical at T4% of full. At this level,

the salt addition rate of 0.5 £t®/min corresponds to a reactivity ramp of
0.01% Bk/k per sec. Dropping the control rods after the power reaches the
normal scram level (150% of full power or 15 Mw) checks the excursion pro-
duced by such a ramp with no significant rise in fuel temperature. KEven
if only two control rods are dropped, sufficient negative reactivity is
inserted to prevent criticality from being attained again if the core is

completely filled.

4,2 Filling With Fuel at Low Temperature

In this accident it is postulated that the graphite has been pre-
heated to the normal startup temperature of 1200°F and fuel salt is
added at a significantly lower temperature. The amount of excess reac-
tivity available depends on the temperature coefficient of reactivity of
the fuel in the full reactor (see Table 1) and the degree of subcooling
of the salt. The heat capacity of the graphite in the core is 3.53 Mw-
sec/°F while that of the salt in the graphite-bearing regions is only
1.45 Mw-sec/°F. Therefore, if the fuel and graphite are allowed to come
to thermal equilibrium, the temperature rise of the salt is 2.4 times the
decrease in graphite temperature. Since the ratio of the graphite to the
fuel temperature coefficient is less than 2.4, heat transfer from the
graphite to the salt reduces the excess reactivity. |

The liquidus temperature of Fuel B, the salt with the largeét nega-
tive temperature coefficient of reactivity, is about 810°F. If salt at
this temperature were added to the reactor and heat transfer from the
graphite were neglected, the maximum amount of excess reactivity would be
1.9%. This is well below the 3.2% shutdown margin provided by the con-

trol rods.
28

4.3 Filling with Concentrated Fuel

The crystallization paths of all three salt mixtures being considered
for use ag MSRE fuel are such that large quantities of salt can be
solidified, under equilibrium conditions, before any uranium (or thorium)
appears in the solid phase. Selective freezing, therefore, provides one
means by which the uranium concentration in the liquid salt can be in-
creased significantly while the salt is in the drain tank. Since the re-
actor vessel is the first major component of the fuel loop that fills on
salt addition, approximately LO% of the salt mixture can be frozen in the
drain tank before it becomes impossible to completely fill the core.

The changes in liquid composition as selective freezing proceeds
depend upon the initial composition and the conditions of freezing.

Figure 10 shows the composition of the remaining melt for Fuel A as a
function of the fraction of salt frozen. The curves are based on the
assumption that only the equilibrium primery solid phase, 6 LiF+<BeFs<ZrF,,
appears.

The effect on premature criticality was evaluated for each of
the three salts with 39%, by weight, frozen in the drain tank as
6 LiF-BeFE’ZrF4.* Under these conditions the full reactor at 1200°F
had about 4% excess reactivity for Fuels A and C and 15% for Fuel B.

Fuels A and C contain significant amounts of thorium and 238U, respec-
tively, which remain in the melt with the 275U during selective freezing.
The poisoning effect of these species greatly reduces the severity of the
filling accident when they are present. The excess reactivities in this
accident exceed the shutdown margin of the control rods so it is necessary
to stop the filling process to prevent a second reactivity excursion after
the rods have been dropped. The accident involving Fuel B determines the
speed with which the fill must be stopped because the reactivity addition
rate for this case is 0.025% 8k/k/sec at k = 1 vs 0.01%/sec for Fuels

A and C.

 

*
The composition of the solid phase has little effect on the nuclear
calculations as long as it deoes not include fissile or fertile material.
N

ORNL DWG 66-7774

 

0.5

uotqowIy STONW

= o
o Q

U07398Ii STOW

Qd
@)
O

0.2

O.l

Welight Fraction of Salt Frozen

Liquid Composition Resulting from Selective Freezing of

"A'" in Drain Tank.

Fig. 10.
Fuel Salt
30
5. ANATYSIS CF MAYIMIM FILLING ACCIDRNT ~

Tt is clear from the preceding secticn that the most severe of the
fiiling accidents considered occurs when the mixture which remalns after
selective freezing of 39% of Fuel B in the drain tank is forced intc the
fuel loop. If it can be showz that this accident is tolerable, then all

of the other acciderts are also tolerable.

5.1 Specification of Accident

The accident which was analyzed in detail included a number of ab-
normal conditions in addition to filling the reactor with highly concen-
trated fuel. The conditions of the acecident and equipment performance,
both normal and abnormal, are described balow.

First, it was assumed that the gas addition system had failed tc the
extent that the gas supply pressure tc the drain tank was 50 psig, the
1imit imposed by the rupture disc, rather than the normal 40 psig. This
gave a salt addition rate of 0.5 ft°/min at the time the reactor first
became critical (at 55% full) and resulted in a reactivity ramp of
0.025% 8k/k per sec. The first corrective action was an automatic rod
drop initiated when the neutron flux reached 150% of design power (15 Mw).
A reiease time of 0.1 sec was assumed and the rods were allowed to fall
with an acceleration of 0.5 times the scceleration of gravity. However,
it was assumed that only two of the three rods actually dropped. Actileon
to stop the fill and initiate a drain was assumed to occur at the same
time asg the red drop. This action invelved automatic opening of the
equalizing valve, HCV-544, and the drain tazk vent valve, HZV-573, and
closing of the gag additiorn valve, HCV-572. It was assumed that orly one
of these valves, HCV-572, actually functioned and 1 sec* wag allowed for

the valve to close., This acticrn, coupled with the initial salt fill rate,

 

¥This time is not critical. Calculations using a 5-sec closing time
did not produce detectably different results.
31

gave a level coast-up of 0.2 ft which added a total of 0.52% reactivity
in excess of that compensated by the two inserted control rods. (If
either of the other two gas valves had functioned, the level would have

dropped and there would have been no second excursion.)

5.2 Preliminary Digital Calculations

The major portion of the transients associated with this acecident
was calculated with the aid of the ORNL analog computer. Howsver, startup
accidents typically begin at very low powers and the power varies over
several orders of magnitude. Since the useful range of the analog com-
puter covers only about two orders of magnitude for any variable, ex-
cessive range switching would be required to simulate the entire transient
on the analog facility. To circumvent this problem, the initial part of
the power transient, from source power to a power which began to affect
the fuel temperature, was calculated with the aid of a digital program,
MURGATROYD.1 MURGATROYD is a point model, nuclear kinetics program,
using six groups of delayed neutrons, with provisions for adding reac-
tivity in the form of steps and/or ramps. The digital program was
started at kgopp = 1 and a power of 1 watt and was used to calculate the
variation of power with time up to 10 kw. This portion of the transient
consumed 24 sec of reactor time and raised the fuel temperature 0.01°F.
The reactor period at 10 kw was about 0.7 sec. The results of this digi-

tal calculation were used as input to start the analog simulation.

5.3 Detailed Analog Simulation

Description of Model
In order to predict the excursions of power and temperature resulting

from the postulated fill accident, a mathematical model was constructed

 

10, W. Nestor, Jr., MURGATROYD — An IBM 7090 Program for the
Analysis of the Kinetics of the MSRE, USAEC Report ORNL-TM-203,
Oak Ridge National ILaboratory, April 6, 1962
32

which described the heat transfer, the nuclear kinetics, and the external
inputs to the system. Figure 11 shows a block diagram of the model used
as a8 basls for the simulation.

It may be noted that except for the heat transfer equations, all of
the computations are relatively straightforward for analog solution since
they involve ordinary differential equations. In spite of the preliminary
digital calculations, the power excursion exceeded the useful-range of the
analog computer and provisions for rescaling the power variable during
the solution were required.

As shown in Fig. 11, however, the fuel and graphite temperatures are
functions of both position and time, and thus are represented by partial
differential equations. To solve these equatlons on an analog computer,
one must use finite difference techniques. ©Since the accuracy of this
solution turns out to be a very important part of the simulation, it will
be discussed in detail.
oimulation of Local Fuel and Graphite Temperatures

Since the core consists of a large number of identical fuel charnels
and vertical graphite stringers, let us first consider the temperature in
a tTypical stringer cross-section, shown in Fig. 12. Due to the symmetry,
we can consider a basic heat transfer "element" as half of a fuel channel
cross-section and one-fourth of a stringer as shown shaded in Fig. 12.
Initially, considering the fuel and graphite as single regions having

mean temperatures T and T.,, heat would be transferred from the fuel to

¥ G’

graphite at a rate Q, where
QzK(TF—TG)

where K may readily be seen to be a function of the conductivities, the
geometry, and the conductance of a film at the Interface, It is important
to note however, the K also depends on the rate of change of the tempera-
tures. TIn general, the higher the frequency of the perturbation, the
greater the error in the computed heat transfer rate between the two

materials, where the approximate transfer rate is always lower than the
NUCLEAR
POWER

o,
ro g

 

ORNL DWG 66-7775

P(I‘, Z)-t)
|
i

ll9h%

i

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CRAPHITE
FUEL TG(I" Z, t)
Tp(r,z,t) | HEAT TRANSFER
o FROM FUEL
T T
Fo NUCTEAR G
l AVERAGE
TEMPERATURES
!
NUCIEAR y 1
KINETTCS EQUATIONS REACTTVITY
; ! COMPUTATION
{6 DELAYED-NEUTRON GROUPS) : _
A &
i
FUEL IEVEL
| COMPUTATICHN
HIGH POWER |——— EDESES T
IEVEL TRIP
(15 Mw) <
SHUT OFF ; Q X
——.w i m———
L FILL, VALVE FLN

 

 

 

 

 

i
¢
i

Fig. 11 Block Diagram of Mathematical Model for Fill
Accident Simulation
34

ORNL DWG 66~7776

 

 

 

 

 

¥ .
{ \ ‘
4 \ i
N e g
- ;
R T
f 4
; o | o TYPICAL
' ’ FUEL
ol ! . CHANNEL
!
SHADED PORTION g t
SHOWS TYPICAL v |
"RASTC'" HEAT . I
SFER —”aifi, . TYPICAT,
ETEMENT s | ,f”"""-“-~“\lm
X | 1 CRAPHTTE
| \_ ,} STRINGER
“ 2!1’

 

Fig. 12. Cross-Section of Typical Graphite Stringer
With Adjacent Fuel Channels

actual value.® Thus it is advantageous to use as fine an approximation
as possible in this case, since heat transfer from the hot fuel to the
graphite and the subsequent rise in graphite temperature is the major
mechanism for curbing the power excursion (because of the small tempera-
ture coefficient of the fuel in the partly full core).

In the simulation, the temperature distribution for & basic element
was approximated by five regions for the fuel and 15 regions for the
graphite, and using slab geometry with a corrected surface-to-volume
ratio. First-order central difference equations were used. The accuracy
of this simulation {assuming & negligible error in the geometry approxi-
mation) was found from a comparicon of the frequency response character-

istics of a distributed slab with those of a lumped-parameter approximation.

 

23, J. Ball, Approximate Models for Distributed-Parameter Heat
Transfer Systems, Preprints of Technical Papers, Fourth Joint Automatic
Conference, University of Minnesota, Mlnneapolls, Minnesota, June 19-21,
1963, AIChE, New York, 1963.

 
35

This comparison was made by means of a digital computer code which can be
used for general slab-geometry calculations of this type. The calculation
showed that for the approximation used, the simulated heat transfer rate
to the graphite would be within 10% of the exact solution values for
perturbation frequencies of up to 25 cycles per second. TFor the tempera-
ture changes involved in the incident being considered, the approximation
was more than adequate.

Two other simplifying assumptions were made:

1) Since, during a fill, the MSRE is a "stationary fuel' reactor,
the only means for axial heat transfer is by conduction, and this was
found to be negligible.

2) Radial conduction (between basic elements) was also assumed
ZET0.

Temperature Averaging

The entire reactor was divided into four major regions of fuel and
graphite. Average nuclear importances of temperature changes and frac-
tions of total power generation were assigned to the components of each
region. These assignments were based on the spatial variation of nuclear
importance and power density in the partly full core. Since the model
for calculating temperatures in a basic heat-transfer element was assumed
linear, the temperature changes in one region were proportional to those
in any other, and were directly related to the fraction of the power
generated in the region and inversely related to the volume of the region.
Thus, from the simulation of a single basic heat-transfer elemeut, a
direct computation was made of the mean fuel and graphite temperatures of
each region. The region temperatures were weighted with their respective
importances and summed to obtain the nuclear average fuel and graphite
temperatures for the reactor. These temperatures were used with theilr
respective coefficients of reactivity to compute the internal contribution
to the kinetic behavior of the reactor.

Other Aspects

The reactivity effects of the fuel and graphite temperatures were
combined with the other reactivity inputs to calculate the power translents.
The other reactivity inputs (see Fig. 11) included 1) the initial ramp

associated with the fill of the reactor, 2) a series of ramps to simulate
36

the dropping control rods, and 3) the decaying ramp associated with an
internal computation of the fuel-level coast-up after the gas-addition
valve was closed. The computer was set up to automatically insert the
appropriate reactivity term with the appropriate time delay during the
transient. The net reactivity was fed to a set of kinetic equations using
six groups of delayed neutrons for the actual power calculation.

Provisions were included for stopping the calculation at any time to
permit the necessary changes in the range of the power calculation.
Facilities were also available for recording any of the computed parameters.
Results of Simulation

The results of the fill-accident simulation are shown graphically in
Figs. 13 and 14. Figure 13 shows the externally imposed reactivity
transient exclusive of temperature compensation effects. The essential
features are the initial, almost~linear rise which produced the first
power excursion as fuel flowed into the core, the sharp decrease as the
rods were dropped, and the final slow rise as the fuel coasted up to its
equilibrium level. Figure 14 shows the power transient and some pertinent
temperatures. The fuel and graphite nuclear average temperatures are the
quantities which ultimately compensated for the excess reasctivity intro-
duced by the fuel coast-up. The maximum fuel temperature refers to the
temperature at the center of the hottest portion of the hottest fuel
channel. The initial power excursion reached 24 Mw before being checked
by the dropping control rods which were tripped at 15 Mw. This excursion
is not particularly important since it did not result in much of a fuel
temperature rise, After the initial excursion, the power dropped to about
10 kw and some of the heat that had been produced in the fuel was trans-
ferred to the graphite. The resultant increase in the graphite nuclear
gverage temperature helped to limit the severity of the second power ex-
‘cursion. Reactivity was added slowly enough by the fuel coast-up that the
riging graphite temperature was able to limit the second power excursion
to only 2.5 Mw. The maximum temperature attained, 1354°F, is well within

the range that can be tolerated.
ORNL DWG 66~7777

 

 

 

 

 

 

 

 

 

 

Accident.

37

I,l.j
4 | .
Ve e El D s
: i :
oy i :
1 : .
; . '
: 3
: ‘e
. +
. 1
SlLiIiId o --3 - il
. . +
‘ . '
. t
' i
: t
.- ey el
: ..
: M
{ i
. . . 1
T '
.
: '
.......... P .- e
' :
: .
. .
- '
1
i
.
e et

 

 

 

 

 

 

 

 

Eta bt s

 

 

TIME (sec)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(%) ZILIAIIOVEY LEN

ty During Maximum Filling

ivi

Net React

Fig. 13.
38

"
1

ORNL DW(

 

 

 

 

 

 

 

 

66-7778_ -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TIME (SEC)

Power and Temperatures During Maximum Filling Accident.

Fig. 1l4.
39

5.4 Discussion

The above analysis indicates that the reactor system is not likely
to be damaged by any of the f£illing accidents considered. The accident
that was studied in detail would probably be less severe than the calcu-
lations indicate because of the conservatism in estimating the effect of
fuel level on reactivity. The degree of fuel concentration by selectilve
freezing used in this study was chosen arbitrarily as that amount which
left just enough salt in liquid form to completely fill the reactor. It
is physically possibie, under ideal ccnditions, to achieve greater concen-
trations of uranium than that assumed. There is no valid basis for as-
sessing the amount of selective freezing that could occur in the drain
tank without being detected. There is also no assurance that freezing
in the drain tank would, in fact, leave all of the uranium in the remaining
melt.

In assessing the credibility of this accident, it should be recog-
nized that several equipment failures were postulated which compounded
the effects of filling the core with excessively concentrated fuel. These
are 1) the failure in the gas-addition system which allowed fuel to flow
into the reactor at an excessively high rate, 2) failure of one of the
three control rods to drop, and 3) failure of a particular two of three
valves to function to stop the fill. Flimination of any one of these
postulated failures would prevent the second reactivity excursion and re-
duce the effects of the accident to trivial proportions.

It has been pointed out that the severity of the filling accident
could be further compounded by postulating that the fuel loop vent valve,
HCV-533, is closed at the start of the fill and is manually opehed Just as
the initial criticality is achieved. Filling with HCV-533 closed would
allow the pressure in the primary loop to rise to about 5 psig, the normal
setpoint of the loop pressure controller, PCV-522. Analog calculations
of the fuel liquid level showed that, if HCV-533 were opened at this pres-
sure and all of the previously postulated failures occurred, the level in
the core would rise 1.5 ft, producing an excursion with temperatures that

would damage the reactor vessel. However, if either HCV-5L4L or HCV-573
40

opens on demand, the gas flow from the drain tank into the primary loop
or out through the vent line limits the pressure decrease in the loop and
keeps the accident within tolerable limits. It does not seem reascnable
to add two more improbable events, failure to have HCV-533 open at the
start and then opening the valve at the crucial moment, to the list of
conditions already postulated. If these two events are postulated, it is
probably reasonable to postulate that only one valve (rather than two)
fails to function on demand. This accident can be further mitigated by
reducing the setpoint of PCV-522 during the filling operation. A setting
of 1.5 to 2 psig would assure that the displaced gas from the fuel loop
goes to the auxiliary charcoal bed during a normal fill but would limit
the pressure in the fuel loop during a fill with HCV-533 closed. ©Such
action might increase the possibility of an activity release from the
normal charcoal beds but the consequences of this are minor, particularly

in view of the low probability of occurrence.

6. CONCIUSIONS

None of the filling accidents that were studied in detaill poses any
threat to the reactor system. On the other hand, it is possible to con-
ceive of a set of circumstances, however unlikely, that could produce
temperatures high enough to breach the primary contaimment. BSuch a breach
would probably occur as a melting of the control rod thimbles.* Even then,
the activity and/or salt release would be confined within the secondary
containment and the release would not approach the maximum credible acci-
dent for which the secondary contaimment is designed. Therefore, even the
worst conceivable filling accident does not represent a hazard to the

operating personnel or the environment.

 

*Since the reactor loop is less than half-full of salt during a
filling accident, it is not possible for a nuclear incident to generate
pressures of the magnitude required for a catastrophic rupture of the
reactor vessel,
-

N

\O 0= O\ o O

lo.

12.
13-17.
18.
19,
20.
21,
o2,
23,
2L,
25,
26.
27 .
28,
29.
30.

41

Internal Distribution

MSRP Director's Office 31. R. N. Lyon
Rm. 219, 9204-1 32, H. G. MacPherson
G. M. Adamson 33. R. E. MacPherson
R. G. Affel 34, H. C. McCurdy
S. J. Ball 35. H. F. Mcbuffie
S. E. Beall 3. A. J. Miller
E. S. Bettis 3T7. R. L. Moore
F. F. Blankenship 38. H. R. Payne
E. G. Bohlmann 39. A. M. Perry
W. H. Cook 40. H. B. Piper
J. L. Crowley 41. B. E. Prince
F. L. Culler Lo, H. C. Roller
S. J. Ditto 43, A. W. Savolainen
J. R. Engel 4h, D. Scott
E. P. Epler b5, M. J. Skinner
A. P. FPraas L6, TI. Spiewak
C. H. Gabbard h7. R. C. Steffy
W. R. Grimes 48. J. R. Tallackson
R. H. Guymon 49. D. B. Trauger
P. H. Harley 50. R. E. Thoma
P. N. Haubenreich 51. W. C. Ulrich
V. D. Holt 52, M. E. Whatley
T. L. Hudson 53. G. D. Whitman
P. R. Kasten 54-55, Central Research Library (CRL)
A. T. Krakoviak 56-57. Y-12 Document Reference Section
R. B. Lindauer 58-60. Iaboratory Records Department
M. I. Lundin 61. Iaboratory Records Department (RC)

External Distribution

D. F. Cope, Reactor Division AEC-ORO

C. B. Deering, Reactor Division, AEC-ORO

A. Glambusso, AEC, Washington

H. M. Roth, Division of Research and Development, AEC-ORO
W. L. Smalley, Reactor Division, AEC-ORC

Division of Technical Information Extension (DTIE)