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)