INTRALABORATORY CORRESPONDENCE OAK RIDGE NATIONAL LABORATORY
February 1, 1973
MSR737
A CATALOG OF DYNAMICS ANALYSES FOR CIRCULATINGFUEL REACTORS
P.N. Haubenreich
Abstract
The effects of fuel circulation on reactor kinetics, mainly through the transport of delayedneutron precursors, have been the subject of a dozen or so significantly different investigations over the past 25 years. Each of the methods now available for calculating the dynamics of complex moltensalt reactor power plants involves approximations that limit its applicability. Good fundamental understanding exists, however, and mathematical formulations that can be developed into more realistic, practical computations have been established. This report briefly catalogs the various efforts on the dynamics of circulatingfuel reactors and lists references that give full details.
Introduction
Analysis of the dynamics of large moltensalt reactor power plants is an area in which sufficient progress has been made and enough reassuring results are available that the MSR program, faced with the necessity of contraction because of the budget, has suspended work. A great deal remains to be done, however, before another moltensalt reactor is built and operated. Penetrating safety analyses will demand more definitive calculations of the effects of credible reactivity events. The requirements on plant equipment, instrumentation, and controls to permit stable, convenient operation must be defined by dynamics analyses yet to be done. None of the existing analysis tools are completely satisfactory for these tasks.
A full evaluation of the situation with regard to dynamics analyses of moltensalt reactors will require a discriminating review of what has been accomplished up to now. This report is intended to be of assistance in such a review. It is essentially an uncritical catalog of the various calculations for the MSRE and later MSR designs, preceded by an account of earlier work on fluidfuel reactor dynamics which provided the basis from which these MSR calculations developed. Numerous references are made to sources of the fully detailed information that will be necessary for an evaluation.
When one considers the dynamics calculations for the MSRE and subsequent MSR designs, he finds that they tend to fall into two general categories. One category is comprised of what might be called control and stability analyses. The other includes the analyses of controlrod experiments and relatively large transients such as conceivably could occur under accident conditions. Although the purposes of these analyses are complementary, some of the requirements on the mathematical models and the computational tools are different. Furthermore, they were generally carried out by different people at ORNL. For these reasons the discussion of the MSR calculations proceeds along the lines of these two categories rather than chronologically.
Early Analyses of CirculatingFuel Reactor Kinetics
Fluid fuels are nothing newwithin little more than a year after the discovery of nuclear fission, Halban and Kowarski, in England, were conducting neutronmultiplying experiments in which the uranium was in the form of a U3O8D2O slurry. In 1943, as part of the Manhattan Project, preliminary designs were developed at Chicago for several types of reactors using D2O slurry fuel, intended as alternates for the Hanford reactors. (ref 1) In 1944 chemists and physicists at the Clinton Laboratories began to study a 10MWt aqueous homogeneous reactor that would both investigate thorium233U breeding and provide a high neutron flux for experiments. That same year (1944) the first actual fluidfuel reactor was built and operated at Los Alamos. This was the LOPO "water boiler", using an aqueous solution of enriched uranium. Almost from the outset it had been recognized that motion of the fuel could have a pronounced effect on reactor kinetics because of the displacement of the delayedneutron precursor atoms in the interval between their formation and the emission of the delayed neutron. For example, it was calculated that if the fuel in a bare tank were stirred vigorously the effectiveness of the delayed neutrons would be reduced by about half, and it was recognized that irregular stirring would cause fluctuations in the power level. (ref 2) Circulation of fuel through a heatremoval system outside of the core was expected to have other important effects on reactor kinetics.
The first actual reactor in which fuel was circulated outside of the core was the 1 MWt Homogeneous Reactor Experiment, which operated at Oak Ridge in 1952 and 1953. Close behind it was a moltensalt reactor, the 2 MWt Aircraft Reactor Experiment, which operated in 1954, also at Oak Ridge. Although contemporary and in close proximity, the HRE and ARE programs had quite different goals and proceeded along separate lines. Among many important differences between the reactors were the flow patterns in the core. The HRE core was a spherical tank, with a tangential inlet and central vortex outlet. In the ARE, on the other hand, the fuel flowed in tubes through a solid moderator, making 11 passes before going to the heat exchanger. In the HRE, effects associated with the rapid expansion of the aqueous solution fuel dominated the reactor's dynamic behavior during power excursions. In the ARE, with its moltenfluoride fuel, the damping effects due to fuel circulation and the emission of delayed neutrons assumed greater significance. As a result of these differences, the kinetics calculations for the HRE and for the ARE developed along separate and somewhat different lines.
In the HRE, the expulsion from the core of the expanding fuel resulted in a very large negative temperature coefficient of reactivity (around 1 x 103/°C). A variety of simple calculations indicated that large reactivity insertions could be tolerated without damage and in cases of interest from the standpoint of safety, the reactor would go well beyond prompt critical. The kinetics equations used in the safety analyses therefore used simplified representations of the delayed neutrons and concentrated the limited computational capabilities of those days on the description of the inertia and friction of the fuel solution in the relief pipe between the core and the surge volume. (ref 3) HRE dynamics experiments involved reactivity ramps at rates up to 0.7% dK/K per second and produced rather large power transients. V. K. Pare and Sidney Visner, in their analysis of these experiments, (ref 4) used a combination of equations developed by J. M. Stein and by K. Fuchs to describe the early and late stages of these transients. In the equations they used an effective value of beta (the delayed neutron fraction) equal to the measured steadystate value, which included the contribution of photoneutrons and the losses due to circulation.
The analyses of the dynamics of the ARE, which were concerned more with stability than with rapid transients, took into account the expected "slug" flow around the fuel loop. W. K. Ergen observed in 1953 that in an analytical model of such a circulatingfuel reactor, power oscillations tend to be damped by the circulation of the fuel apart from any effect of delayed neutrons. (ref 5) Ergen also recognized that "in a circulatingfuel reactor some delayed neutrons are emitted outside the reacting zone and this modifies the inhour equation in a manner which depends on ... whether fuel which is present in the reactor at a given moment returns to the reactor in a "slug" or gets mixed in the outside circuit with other fuel." Ergen derived the inhour formula for a circulatingfuel reactor with slug flow and flat flux in the core, (ref 7) but unfortunately, the details of his analyses were not widely reported because of the secrecy around the aircraft reactor work. In a brief, unclassified memorandum (ref 6) issued some years after the ARE operation, he reported that "the inhour equation, both with stationary fuel and with fuel circulating, was used to calibrate the regulating rod," but gave no details.
Meanwhile, at Brookhaven National Laboratory, studies had begun on the LiquidMetal Fuel Reactor, in which molten UBi solution circulated through a graphite core and an external heat exchanger. As part of the LMFR studies, J. A. Fleck, Jr. pursued a different, basic approach to the analysis of the kinetics of circulatingfuel reactors. Fleck developed a method for computing reactor periods by finding the eigenvalues of the system of differential equations describing the reactor kinetics, first with one (ref 8), then with several groups of delayed neutrons. (ref 9) In recognition of the utility of simpler relationships, he also derived approximate inhour relations for systems with slug flow, with either sinusoidal or flat power distribution in the core section of the loop. From his studies of these relations Fleck noted the important fact that "while the flux is increasing, there are more delayed neutrons in the system than would be predicted by the fraction [at steady power]. ... As the flux rises, the circulating fuel reactor simply loses less reactivity to the external loop than when it is just critical. The system thus behaves as though it contained an additional latent reactivity whenever the power rises." He concluded that the use of an effective delayed neutron fraction with a "conventional" inhour expression is permissible only when periods are large. In analyses of delayedneutron effects, Fleck neglected powertemperature feedback; later he did a study of the temperaturedependent kinetics of a circulatingfuel reactor with the "openloop" or constant inlet temperature boundary condition and neglecting delayed neutron effects. (ref 10) Results of Univac computations with Fleck's equations were encouraging, but there was a lingering thought that, because of the complexity of the kinetics at high power, the response of a complex, circulatingfuel reactor would not be accurately predicted and that further studies were required to assure stability against oscillations at high power."
The line of dynamics analysis that started with the HRE was developed further by P. R. Kasten for a second aqueous homogeneous reactor, the 5 MWt Homogeneous Reactor Test*, which operated at Oak Ridge from 1957 into 1961.
*This reactor was usually referred to in ORNL publications as the HRT and in other literature as the HRE2. An unusual feature of this reactor was that there were no safety rods or other fastacting reactivity controls: HRE experience and analyses of conceivable reactivity additions and peak pressures in the HRT indicated that such controls were not essential. (ref 12) The mathematical description of the HRT dynamics, developed by Kasten (ref 13) for safety studies, again stressed the hydraulics involved in the expulsion of the fuel solution from the core. Reactivity effects of temperature (fuel density) changes were approximated by simplified thermal equations. Kinetics of the chain reaction were represented by the equations for a stationaryfuel reactor, with the delayed neutron fractions reduced to account for steadystate losses due to circulation. Safety calculations on the Oracle determined reactivity steps or ramps corresponding to the maximum tolerable pressure surge, for comparison with rates of reactivity increase associated with physical events. (ref 14) Stability criteria were also established by Kasten, using different sets of linearized equations to examine highfrequency (nuclear power) and intermediatefrequency (loaddemand) effects. The stability of the coupled system of core, drain tanks, and condensate storage to fuel concentration effects was analyzed by Melvin Tobias. (ref 15)
The dynamics experiments on the HRT were minimal (because there were no control rods) but considerable effort was devoted to the analysis of some dynamic phenomena observed during operation. The small (up to 1%) fluctuations in reactor power were studied extensively by Jitsuya Hirota*. Hirota showed by noise analysis techniques (implemented only by recorder charts and desk calculator) that certain characteristics of the fluctuations could be related to system conditions, particularly to the status of the fuel solution. (ref 16) Larger excursions in reactor power that occurred on an average of a few times per hour when the reactor was operated above 1 MWt were analyzed by M. W. Rosenthal, S. Jaye, and M. Tobias. They developed an IBM704 program for computing the time variation of the reactivity associated with observed power excursions, based on conventional neutron kinetics equations and temperature feedback from a simplified thermal model. (ref 17) The indications were that the excursions were most likely caused by movement of separated uranium in the core. (A reactivity balance and physical evidence of local overheating were other indications of uranium separation by at least one of several possible processes.)
*Visiting scientist from Japan Atomic Energy Research Institute At about the time of the HRT startup, John MacPhee of AMF Atomics published a paper (ref 18) in which he considered some simple approximations for the effects of delayed neutrons. He derived transfer functions and computed the frequency response of a model in which temperature feedback was neglected, the core was perfectly mixed, and the delayedneutron precursors returned to the core after decaying during passage through the external loop in slug flow.
Reactivity effects produced by fluid motion in reactor cores were further elucidated by some original work done by Bertram Wolfe between 1959 and 1961. His initial interest was in "Borax"type excursions where the fuel is stationary but the moderator water is ejected at high velocity, "dragging" neutrons with it. Using perturbation techniques, he developed expressions for the reactivity effects (which turned out to be significant). Wolfe then applied the same mathematical treatment to circulatingfuel reactors to obtain expressions for the delayed neutron lifetime and the inhour relation for systems with slug flow and sinusoidal flux distribution in the core. (ref 19) His analysis did not include temperature effects on reactivity.
R. V. Meghreblian and D. K. Holmes, in their 1960 book on reactor analysis, devoted some 20 pages to an exposition of the kinetics of circulatingfuel reactors. (ref 20) Their presentation of the separate effects of delayed neutrons and temperature followed the general scheme of analysis used by Fleck (ref 9), but with a somewhat cruder model.
Control and Stability Analyses for MoltenSalt Reactors
MSRE
In the first investigation of MSRE dynamics, in 1960, an analog representation of the core and heat removal system set up by O.W. Burke (I & C Division) was used to determine temperatures following interruption of fuel circulation. (ref 21) Over the next year Burke analyzed various other transients to delineate the need for automatic control systems. (ref 22) Later in the design of the servo control system, S.J. Ball (I & C Division) used an analog simulation in which the core was represented by nine regions, each containing two fuel "lumps" and one graphite "lump," instead of the oneregion (threelump) representation used by Burke. Both Burke and Ball represented delayed neutrons by one group, with a total yield reduced to account for losses due to circulation. (ref 23)
When Ball and T. W. Kerlin (Reactor Division) made their more definitive analysis of the inherent stability of the MSRE in 1965, they used the same lumpedthermal model of the core, but more realistic equations for some other parts of the reactor system, including 6 groups of delayed neutrons with approximate representation of precursor losses and aroundtheloop feedback. (Their report (ref 24) includes a convenient comparison of the different models and parameter values used in the various MSRE studies.) Instead of an analog computer, Ball and Kerlin used a digital computer program, MATEXP, (ref 25) to solve the system of equations for the time response and another code to get the frequency response. The adequacy of their analysis was later verified when they measured the frequency response of the MSRE, using small input reactivity disturbances, and found that the results generally were in good agreement with their predictions. (ref 26)
Prior to operation of the MSRE with 233U, R.C. Steffy, Jr. and P.J. Wood, (Reactor Division) in 1968 predicted the frequency response of the reactor with the new fuel. Their model was the same as Kerlin's and Ball's except that a firstorder time lag had been added to represent the mixing in the reactor vessel plenums that had been inferred from comparison of observed and predicted frequency response. (ref 27) Results of extensive dynamics testing with 233U fuel later confirmed their predictions. (ref 28, 29) These theoretical and experimental analyses of MSRE dynamics were eventually summarized in two papers in Nuclear Technology. (ref 30,31)
In 1965, as part of the operator training program, Ball set up two reactor kinetics simulators at the MSRE, wired directly to the reactor controls and instrumentation. (ref 32) The reactor model used in these simulators was a somewhat simplified version of that being used in the stability studies. (The core was represented by one threelump region.)
The capability of the servo system for controlling the reactor with 233U fuel was proved in advance by Ball, Steffy and a team from the MIT practice school by an unusual experiment with the reactor during operation with 235U. (ref 33) The method consisted of modifying the flux signal to the servo system so that the response to rod motion resembled that which would result if the fuel were 233U instead of 235U. Prior to doing the test on the reactor, the method was verified by analog simulation, using a model similar to that used in the operator training simulator.
Throughout the years of MSRE operation, the neutron noise was extensively analyzed and techniques were developed by which it could be used for reactor diagnosis. D.N. Fry, R.C. Kryter, and J.C. Robinson (I & C Division) observed that changes in the amount of gas bubbles in the fuel salt could be inferred from noise measurements, (ref 34) using a mathematical model of system hydraulics and neutronics. (ref 35) (Their model, with its emphasis on hydraulics, differed somewhat from that of Ball et al.) Experiments with small induced pressure perturbations were used to test the model and to determine the best boundary conditions. (ref 36, 37) Comparison with the observed behavior indicated minor deficiencies in the hydraulic model, which Ulrich later refined. (ref 38) Noise studies were facilitated when C.D. Martin (I & C Division) programmed the online digital computer at the MSRE to collect the data and compute the neutron power spectral density upon demand. Following the observation that the neutron noise around 1 Hz was especially sensitive to the bubble fraction in the fuel, a continuous monitoring system was installed that displayed the magnitudes of the neutron noise and the pressure noise in the frequency range 0.5 to 2 Hz. (ref 39, 40)
MSBR
A preliminary study of the dynamics and control of a 1000MWe, singlefluid MSBR was made by W. H. Sides, Jr. (I & C Division) in 1969, using an analog computer simulation. (ref 41) For this study he used an abbreviated, lumpedparameter model of the thermal system in which the core was represented by one region (two fuel lumps and one of graphite). Two groups of delayed neutrons were included, with approximate representation of circulation effects. Later Sides went to an expanded lumpedparameter model in which the actual 2zone core was simulated by three core regions (each with two lumps of fuel and one of graphite) with salt flow in two parallel paths. (ref 42) Delayed neutrons again were represented by two groups, but the model allowed for the effects of variable fuel circulation rate on losses and feedback to the core inlet. Analog simulation with this model was used to study the integrated plant response to perturbations such as changes in load demand, salt flow, and reactivity. (ref 43)
After analysis of earlier results indicated the need for further refinement of the steam generator models, Burke set up a hybrid computer simulation of the referencedesign MSBR. (ref 44) This latest (1972) model featured a highlydetailed representation of the steam generator on ORNL's digitalanalog hybrid computer and a lumpedparameter analog model of the rest of the system like that used by Sides.
Analyses of Large Transients in MoltenSalt Reactors
MSRE
The analog system that was set up at the very outset of the MSRE design effort to examine system dynamics was incapable of representing extremely large power excursions, such as might result from some postulated accidents. Therefore C.W. Nestor, Jr. (Reactor Division) was assigned to develop a digital computer program for MSRE kinetics. The first result was a program, called MURGATROYD, (ref 45) which was a revised version of a program used previously for aqueous homogeneous reactor analysis, extended to include in a very approximate way the effects of separate fuel and moderator in the moltensalt reactor core. Using MURGATROYD, P.N. Haubenreich and J.R. Engel (Reactor Division) analyzed several conceivable reactivity accidents, including cold slugs, filling accidents and uncontrolled rod withdrawal. (ref 46) These results, plus an analysis of fuel pump stoppage effects that was done on Burke's MSRE analog, provided the basis for the discussion of nuclear safety in the preliminary hazards report (second addendum).
MURGATROYD gave unrealistic results in many cases because the reactivity effects of fuel temperature rise were simply taken to be proportional to a weighted mean of the core inlet and outlet temperatures. To eliminate this approximation, Nestor proceeded to develop another computer program, ZORCH, (ref 47) which included computation of timedependent axial distributions of fuel and graphite temperatures. The spatial distributions of temperatures were multiplied by a nuclear importance function and integrated over the core volume to obtain "nuclear average" temperatures that produced reactivity feedback. Neither in MURGATROYD nor in ZORCH was explicit account taken of the system outside of the core except for its effect on core pressure. (There was no provision for aroundtheloop feedback to the core inlet.)
Six groups of delayed neutrons were included in both MURGATROYD and ZORCH, with equations of the same form as in fixedfuel reactors but with fission yields reduced to allow for emission outside of the core. Both codes included a simple calculation of steadystate delayedneutron losses from input data on loop transit times, using an equation appropriate for a flat flux in the core. In preparation for accident analyses with ZORCH, Haubenreich computed "effective" values for yields of 6 groups of delayed neutron during steadypower operation, taking into account the effect on leakage probabilities of the initial energies and the spatial distribution of precursors in the MSRE core. (ref 48) The input data for ZORCH were then adjusted so that the kinetics calculations were done with these effective delayedneutron fractions. The analysis of various accidents with ZORCH was discussed in some detail in the MSRE nuclear analysis report (ref 23) and results were considered in the safety analysis report. (ref 49)
When the MSRE was started up in June 1965, preliminary analysis of the zeropower physics experiments revealed needs for improvement in the calculation of delayed neutronsthe decrease in reactivity due to loss of delayed neutrons was significantly less than had been calculated and the rodbump period experiments with fuel circulating were not adequately described by the conventional inhour equation with reduced delayedneutron fractions. To overcome these deficiencies B. E. Prince (Reactor Div.) developed an analytical model which described the effects of fuel circulation on delayed neutron source distribution while the power is changing on a stable asymptotic period and accounted for the delayed neutrons emitted in the fuel plenums above and below the core. (ref 50) When the inhourtype equation that Prince derived was used to calculate rod worths from the periods measured with the fuel circulating, values were obtained that agreed satisfactorily with those calculated with the fuel stationary. (ref 51) The reactivity loss due to circulation predicted by Prince's model was also in good agreement with the observations, reflecting the significant effect of the neutrons emitted in the upper plenum, which had been neglected in the earlier calculation.
In his development of the expressions needed for analyzing the MSRE experiments, Prince took a basic approach, starting with the general timedependent reactor equations written to include the transport of delayedneutron precursors. (ref 50) He carried through to a practical computation (by an IBM7090 program) only for the cases of a justcritical reactor or one on a stable asymptotic period. In the report describing this work, however, he pointed out the desirability of analyzing the transient delayedneutron modes and mentioned that he had begun some preliminary work on this subject. These extended studies, unreported until recently, (ref 52) show that, with inclusion of the abovementioned fluid transport terms in the general equations for circulatingfuel reactors, the application of the usual adjointweighting and integration techniques used to derive the pointkinetics equations from these general equations do not result in the usual set of timedependent ordinary differential equations associated with "pointkinetics" of stationary fuel reactors. However, a timedependent integrodifferential equation describing the kinetics of the neutron population can still be obtained. In reference 52, Prince describes how the kernels of this equation which are associated with weighted precursor and temperature distributions may be calculated, and discusses the potential advantages of using this type of formulation in future kinetics analyses of circulatingfuel reactors.
In 1967, when operation of the MSRE with 233U fuel was being considered, the relation developed to analyze the 235U experiments was used to calculate the reactivity increments necessary to produce a given stable period with the new fuel. It was found, as expected, that because of the smaller delayed neutron yields from 233U fission, the reactivityperiod relationship was more sensitive to circulation. (ref 53) An approximation that used effective yields of delayed neutrons in steadypower operation in the inhour equations for stationary fuel, although differing significantly from the exact relationship, suggested that ZORCH calculations with the same effective yields should be adequate for the necessary safety analyses. Therefore ZORCH was used to calculate the response of the 233Ufueled MSRE to two credible reactivity accidentsuncontrolled rod withdrawal and resuspension of separated uranium. (ref 54) Results of these calculations, which included effects of rod scram on level but not on period, were considered in the review of safety with 233U fuel. (ref 55)
Analyses of MSRE transients with 233U fuel were subsequently carried out independently by Burke and F.H.S. Clarke (I & C Division), using an analog simulation of the MSRE that had been set up primarily for use in determining the changes that might be necessary or desirable in the MSRE servo controller characteristics or the control rod trip points. (ref 56) Their delayedneutron representation included 6 groups, with approximate effects of fuel circulation on removing precursors from the core and feeding them back after decay in the external loop. In the precursor equations the core was represented as one wellmixed lump, but for calculation of thermal effects nine core regions (each with two lumps of fuel and one of graphite) were used. The servo controller and the rod trip circuits were realistically represented in this model.
MSBR
After the adoption of the singlefluid MSBR as the focus for the MSR program in 1968, O.L. Smith (Reactor Division) began to develop a mathematical model especially designed for analysis of transients such as would be of concern in safety analysis. Particular emphasis was placed on realistic treatment of the distribution of delayed neutron precursors along a fuel channel. He prepared a program for solution of his dynamics equations on a digital computer, and in 1969 published some preliminary results of studies of the sensitivity of MSBR dynamics to fuel and graphite temperature coefficients of reactivity. (ref 57) A report describing Smiths model was drafted and reviewed in 1970 but suggested revisions have not been made and the report has not been issued. The computer program is operable, however: in June 1972, Smith made the minor changes necessary to conform to the current ORNL computer systems and ran a trial case of a reactivity step in the reference MSBR.
References
1. J. A. Lane, "Homogeneous Reactors and Their Development," Chap 1 in Fluid Fuel Reactors, J. A. Lane, H. G. MacPherson and F. Maslan, Eds., AddisonWesley, Reading, Mass. (1958).
2. A. M. Weinberg and E. P. Wigner, The Physical Theory of Neutron Chain Reactors, p. 598, Univ. of Chicago Press (1958).
3. W. C. Sangren, Kinetic Calculations for Homogeneous Reactors, ORNL1205 (April 1952).
4. V. K. Pare and S. Visner, Experiments on the Kinetics of the HRE, ORNL2329 (1957).
5. W. K. Ergen, "Kinetics of the CirculatingFuel Nuclear Reactor," J. Appl. Phys. 25(6), 702711 (June 1954).
6. W. K. Ergen, The Physics of the FusedSalt Reactor Experiment, ORNL internal memorandum CF572130 (Feb. 1957).
7. W. K. Ergen, The Inhour Formula for a CirculatingFuel Nuclear Reactor with Slug Flow, ORNL internal memorandum CF5312108 (Dec. 1953).
8. J. A. Fleck, Jr., "Kinetics of Circulating Reactors at Low Power," Nucleonics, 12(10), 5255 (Oct. 1954).
9. J. A. Fleck, Jr., Theory of LowPower Kinetics of CirculatingFuel Reactors with Several Groups of Delayed Neutrons, BNL334 (April 1955).
10. J. A. Fleck, Jr., The TemperatureDependent Kinetics of Circulating Fuel Reactors, BNL357 (July 1955).
11. J. Chernick, "Reactor Physics for Liquid Metal Reactor Design," Chap. 19 in Fluid Fuel Reactors, J. A. Lane, H. G. MacPherson, and F. Maslan, Eds., AddisonWesley, Reading, Mass. (1958).
12. S. E. Beall and S. Visner, Homogeneous Reactor Test Summary Report for the Advisory Committee on Reactor Safeguards, pp. 100116, ORNL1834 (Jan. 1955).
13. P. R. Kasten, Dynamics of the Homogeneous Reactor Test, ORNL2072 (June 1956).
14. P. R. Kasten, Operational Safety of the Homogeneous Reactor Test, ORNL2088 (July 1956).
15. M. Tobias, HRP Quarterly Progr. Rep. Jan. 31, 1955, ORNL1853, pp. 2329.
16. J. Hirota, Statistical Analysis of Small Power Oscillations in the HRT, ORNL internal memorandum CF601107 (Jan. 1960).
17. M. W. Rosenthal, S. Jaye, and M. Tobias, Power Excursions in the HRT, ORNL2798 (Feb. 1960).
18. J. MacPhee, "The Kinetics of CirculatingFuel Reactors," Nucl. Sci. Engrg 4, 588597 (1958).
19. B. Wolfe, "Reactivity Effects Produced by Fluid Motion in a Reactor Core," Nucl. Sci. Engrg 13, 8090 (1962).
20. R. V. Meghreblian and D. K. Holmes, Reactor Analysis, pp. 590610, McGrawHill (1960).
21. MSR Program Quarterly Report, July 31, 1960, ORNL3014, pp. 4653.
22. MSR Program Semiannual Progress Report, Feb. 28, 1962, ORNL3282, pp. 1014.
23. P.N. Haubenreich, J. R. Engel, B. E. Prince, and H. C. Claiborne, MSRE Design and Operations Report, Part III  Nuclear Analysis, ORNLTM730 (Feb. 1964).
24. S. J. Ball and T. W. Kerlin, Stability Analysis of the MSRE, ORNLTM1070 (Dec. 1965).
25. S. J. Ball and R. K. Adams, MATEXP  A General Purpose Digital Computer Program for Solving Ordinary Differential Equations by the Matrix Exponential Method, ORNLTM1933 (Aug. 1967).
26. T. W. Kerlin and S. J. Ball, Experimental Dynamic Analysis of the MSRE, ORNLTM1647 (Oct. 1966).
27. R. C. Steffy, Jr. and P. J. Wood, Theoretical Dynamic Analysis of MSRE with 233U Fuel, ORNLTM2571 (July 1969).
28. R. C. Steffy, Jr., FrequencyResponse Testing of the MSRE, ORNLTM2823 (March, 1970).
29. R. C. Steffy, Jr., Experimental Dynamic Analysis of the MSRE with 233U Fuel, ORNLTM2997 (April, 1970).
30. T. W. Kerlin, S. J. Ball, and R. C. Steffy, "Theoretical Dynamics Analysis of the MSRE," Nuclear Technology 10, 118 (Feb. 1971).
31. T. W. Kerlin, S. J. Ball, R. C. Steffy, and M. R. Buckner, "Experiences with Dynamic Testing Methods at the MSRE," Nuclear Technology 10, 103 (Feb. 1971).
32. S. J. Ball, Simulators for Training MSRE Operators, ORNLTM1445 (April 1966).
33. MSR Program Semiannual Progress Report, Feb. 29, 1968, ORNL4254, p. 32.
34. D. N. Fry, R. C. Kryter, and J. C. Robinson, Measurement of Helium Void Fraction in the MSRE Fuel Salt Using NeutronNoise Analysis, ORNLTM2315 (Aug. 1968).
35. J. C. Robinson, "Analytical Determination of the Neutron FluxtoPressure Frequency Response: Application to the MSRE," Nucl. Sci. and Eng. , 42 (3) , 382 (Dec. 1970) .
36. J. C. Robinson and D. N. Fry, Determination of the Void Fraction in the MSRE Using Small Induced Pressure Perturbations, ORNLTM2318 (Feb. 1969).
37. J. C. Robinson and D. N. Fry, "Experimental Neutron FluxtoPressure Frequency Response for the MSRE: Determination of Void Fraction in Fuel Salt," Nucl. Sci. and Eng., 42(3), 397 (Dec. 1970).
38. W. C. Ulrich, An Extended Hydraulic Model of the MSRE Circulating FuelSystem, ORNLTM3007 (June 1970).
39. MSR Program Semiannual Progress Report, Aug. 31, 1969, ORNL4449, p. 36.
40. MSR Program Semiannual Progress Report, Feb. 28, 1970, ORNL4548, p. 7.
41. W.H. Sides, Jr., MSBR Control Studies, ORNLTM2489, (June 1969).
42. W.H. Sides, Jr., MSBR Control Studies: Analog Simulation Program, ORNLTM3102, (May 1971).
43. W.H. Sides, Jr., Control Studies of a 1000MW(e) MSBR, ORNLTM2927, (May 1970).
44. O.W. Burke, Hybrid Computer Simulation of the MSBR, ORNLTM3767,. (May 1972).
45. C. W. Nestor, Jr., MURGATROYD  An IBM 7090 Program for the Analysis of the Kinetics of the MSRE, ORNLTM0203 (April, 1962).
46. P. N. Haubenreich and J. R. Engel, Safety Calculations for MSRE, ORNLTM0251, (May 1962).
47. C. W. Nestor, Jr., ZORCH  An IBM 7090 Program for the Analysis of Simulated MSRE Power Transients with a Simplified SpaceDependent Kinetics Model, ORNLTM0345 (Sept. 1962).
48. P. N. Haubenreich, Prediction of Effective Yields of Delayed Neutrons in MSRE, ORNLTM0380, (Oct. 1962).
49. S. E. Beall, P. N. Haubenreich, R. B. Lindauer, and J. R. Tallackson, MSRE Design and Operations Report, Part V, Reactor Safety Analysis Report, ORNLTM0732 (Aug. 1964).
50. B.E. Prince, Period Measurements on the Molten Salt Reactor Experiment During Fuel Circulation: Theory and Experiment, ORNLTM1626 (Oct. 1966).
51. B.E. Prince, et,.al., ZeroPower Physics Experiments on the MoltenSalt Reactor Experiment, ORNL4233, (Feb. 1968).
52. B.E. Prince, Improved Representation of Some Aspects of CirculatingFuel Reactor Kinetics, ORNLTM in publication (Jan. 1973).
53. MSR Program Semiannual Progress Report, Period ending Aug. 31, 1967, ORNL4191, p. 58.
54. MSR Program Semiannual Progress Report, Period ending Feb. 29, 1968, ORNL4254, p. 37.
55. P.N. Haubenreich et al., MSRE Design and Operations Report, Part VA, Safety Analysis of Operation with 233U, ORNLTM2111, (Feb. 1968).
56. O.W. Burke and F. H. S. Clark, Analyses of Transients in the MSRE System with 233U Fuel, ORNL4397, (June 1969).
57. MSR Program Semiannual Progress Report, Period ending Aug. 31, 1969, ORNL4449, p. 66.
