ML20207A196

From kanterella
Jump to navigation Jump to search
Forwards Draft Investigation of Steam Explosion Loadings W/ SIMMER-II, for Review.Meeting Suggested for 860115 to Discuss Results
ML20207A196
Person / Time
Issue date: 11/18/1985
From: Bassett O
NRC OFFICE OF NUCLEAR REGULATORY RESEARCH (RES)
To: Speis T
Office of Nuclear Reactor Regulation
Shared Package
ML20207A169 List:
References
FOIA-86-678 NUDOCS 8611070259
Download: ML20207A196 (52)


Text

'

i MEMORANDUM FOR: Themis Speis, Director Division of Safety Technology N0p 2 g g FROM: 0. E. Bassett, Director Division of Accident Evaluation 7

SUBJECT:

STEAM EXPLOSION ENERGETICS -

A recent letter to you from Dr. T. G. Theofanous on Steam Explosion Energetics suggested a " peer review process and feedback" and a presentation at a one day meeting some time after the middle of December. ,

The enclosed copy of the draft sumary chapter from LANL's report on "An Investigation of Steam Explosion Loadings wilh SIMMER-II" provides several ~

interesting examples for comparison with Dr. Theofanous's examples. One point of great interest will be the amount of venting calculated, given an explosion in the lower plenum. Currently these two reports calculate significantly different amounts of venting, but under different initial conditions b.ecause of the examples selected.

Giventheimportanceofthistopic,wesuggestasimultaneouslieerreview process for both reports and presentations by Dr. Theofanous and LANL~at a two day meeting with time allowed for peers to provide their' reviews and.

development of an overall consensus. Recall that some early results from the LANL work were presented at the Steam Explesion Review Group Keoting in Harpers Ferry. The final draft of the LANL report is expected to be ready for review by the first week in December. A simultaneous review process will sue contractor and staff time, and travel costs. .-

We suggest that three weeks be allowed for review of both reports beginning December 16, 1985. Considering the Christmas Holidays taken at many national laboratories, we suggest the two day rieeting be held about January 15, 1986.

For meeting arrangenent and coordination responsibilities, we suggest -

Mr. Cardis Allen and Mr. John Telford. -

Please let me know whether you agree with these suggestiors.

0..E. Bassa t, Director Division of Accident. Evaluation

Enclosure:

As stated -

cc: Harold Sullivan, LANL Charles Bell, LANL k M "

55 William Bohl, LANL distribution: circ; chron; g ; Curtis; lber; cal en; Ross; Minogue; Morrison; Bassett

... _h . ... __ _

NAME: TELF (D/im : CURTI

\  : RIS0N : BASSETT  :-  :  :

iiiir ii p ii T ii; w ii ; ii;R;;r T ii;J p ir T ; ~ ~ ; - - - i- - - -

8611070259 861030 PDR FOIA 2 #o' CURRAN 86-678 PDR

i CRAFT l

AN INVESTIGATION OF STEAM EXPLOSION LOADINGS WITH SI WER-II l

i by W. R. Bohl Contributors C. R. Bell T. A. Butler L. M. Hull J. G. Bennett P. J. Blewett

~~

D1 t -

,J< t t

g .-- - - - - ,

e 4 TABLE OT CONTENTS DRAFT ABSTRACT .......................................................................

LIST OT FIGURES ................................................................

LIST OF TABLES .................................................................

1.

SUMMARY

REPORT ......................................................... .

A.i Introduction ......................................................

S'. Modification of SIH4ER-Il .........................................

C. Calibration to SNL Steam Explosion Experiments ....................

D. Calibration to Los Alamos Shallow Pool Experiments ................

i

' E. . SIMMER-Il Reactor Case Calculative Results ........................

F. Containment Failure Probabilities .................................

G. Conclusions and Recommendations ...................................

II. SI2!ER-II MODIFICATIONS FOR THE MOLTEN-CORE /000LAhT IhTERACTION PROGRAM.

A. Introduction ......................................................

B. Summary of Modified Equations .....................................

1. AEOS Modifications ............................................
2. Heat Transfer Modifications ...................................
3. Vaporization-Condensation Modifications .......................

C. Discussion of AEOS Modifications ..................................

1

1. Review of Difficulties in the Manual Formulation .............
2. Exacination of Proposed Solutions to AEOS Problems ...........
3. Th e Va po r He a t Ca pa c i t y Mod i f i c a t i o n . . . . . . . . . . . . . . . . . . . . . . . . .

1 4. Modification of the Gas " Constant" at High Pressures .........

i 3

5. Definition of the Vapor Temperature for Single-Phase Cells and 30ther Modifications ...............................,..,.......

s

6. 3 Discussion of Results ........................................

' 7. Suggestions for Further Improvement ..........................

S. Conclusions ...........)......................................

( D. Modification of SIMMER-II Liquid-Liquid Heat Transfer for Water ...

l, E. Vaporization / Condensation Model Changes ...........................

c) F. Miscellaneous' Corrections for the Molten-Core / Coolant Interaction..

g . <

l Prograc ...........................................................

I III. LOVER HEAD MODEL FOR SIHMER-Il .........................................

A. Introduction ......................................................

B. Input .............................................................

C. Head Failure ......................................................

D. Lower Head Motion .................................................

l f \

Y i,

DRAFT E. Boundary Conditions and Edits .....................................

F. Correction Set ....................................................

G. Sample Problem ....................................................

IV. CORRELATION TO SNL STEAM EXPLOSION EXPERIMENTS .........................

V. ANALYSIS OF LOS ALAMOS EXPERIMENTAL DATA FOR SHALLOW POOLS .............

A. Introduction ...................................................... .

B. Experimental Correlation ..........................................

C. Hypothetical Cases with Nonuniform Interfaces .....................

D. Experimental Comparison with a Pressure Ratio of 50:1 .............

E. Experimental Comparison with a Nonuniform initial Interface .......

F. Conclusions .......................................................

VI. ANALYSIS OF STEAM EXPLOSIONS WITH SIMMER-II ............................

A. CASE 1 ............................................................

1. Water Surface Area Investigation .............................
2. Comparison to the ZIP Study Using the New Models .............
3. Case 1 Update ................................................

B. CASE 2 ............................................................

1. Explosion After 1 second of Mixing (Scoping Calculation) .....
2. Explosion After O.7 s of Mixing (Second Scoping Calculation) .
3. Ca s e 2 F i n a l i z e d . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . '. . . . . . . . . . .

C. CASE 3 ............................................................

D. CASE 4 ............................................................

E. CASE 5 ............................................................

Yll. PROBABILITY OF A STEAM EXPLOSION .......................................

V111. RESEARCH PRIORITIES ON STEAM EXPLOSIONS ................................

APPENDIX A: THE SIMMER-Il MANUAL ANALYTIC-EQUATION-OF-STATE TREATMENT .........

APPENDIX B: CFS NODE / VAPOR /VRBAEOS (EQUATION OF STATE CORRECTION SET). . . . . . . . .

APPENDIX C: REVISED SIMMER-Il INPUT DESCRIPTION FOR THE AEOS ..................

APPENDIX D: AEOS SIMULATION PROGRAM ...........................................

APPENDIX E: SUGGESTED VALUES FOR WATER AEOS INPUT .............................

APPENDIX F: STANDARD SIMMER-Il LIQUID-LIQUID HEAT-TRANSFER MODEL ..............

APPENDIX G: CORRECTION SET FOR MODIFIED LIQUID-LIQUID HEAT TRANSFER ...........

1 APPENDIX H: VAPORIZATION / CONDENSATION 00RRECT10N' SET .......................... I APPENDIX 1: SibNER-Il MANUAL TREATMENT OF THE VAPORIZATION /ODNDENSATION MODEL . l 1

APPENDIX J: MISCELLANEOUS 00RRECTIONS .........................................

1 APPENDIX K:

SUMMARY

DESCRIPTION OF THE HEAD FAILURE MODEL SUBROUTINE .......... '

APPENDIX L: DEFINITIONS OF VARIABLES IN THE PLUGW CORRECTION SET ..............

APPENDIX M: LOVER HEAD DYNAMIC ANALYSIS .......................................

,3 - - . .

l m l

Li

DNAFT APPEND 1X N: CONSIDERATIONS REGARDING RESULTS FROM THE SDOT FAILURE MODEL ......

APPENDIX 0: A LIMITED REVIEW OF SNL STEAM EXPLOSION EXPERIMEhTS ...............

APPENDIX P: SletER-Il INPUT FOR THE MD-19 EXPER]MEhT USING A UNIFORM MIXING ZONE APPENDIX 0: SletER-Il COMPARISON WITH THREE FIELD MIXING CALCULATIONS .........

APPENDIX R: SletER-Il INPUT FOR THE ANALYSIS OF THE MD-19 EXPERIMEhTS US]?Ki A NON-UNIFORM MIXING ZONE ........................................... .

APPENblX S: EFFECT OF CHANGING THE AXIAL AND RADIAL 00NSTRAINTS IN THE SIMULATION OF MD-19 ...............................................

APPENDIX T: EXPLOSION CALCULATION STARTING WITH A STANDARD SI4fER-II PREMIXED CONFIGURATION .....................................................

APPENDIX U: 514tER-11 PREMIXING WITH HIGH STEAM PRODUCTION RATES ..............

APPENDIX V: DETAILS ON THE UPPER BOUND SI41ER-Il STEAM EXPLOSION CALCULATION ..

i

W

DR/H

, AN INVESTIGATION OF STEAM EXPLOSION LOADINGS WITH S141ER-Il by W. R. Bohl ,

Contributors C. R. Bell T. A. Butler L. M. Hull J. G. Bennett P. J. Blewett ABSTRACT The S141ER-Il code was used to provide estimates of the maximum anticipated loads on the upper head of a pressurized water reactor (WR) following an in-vessel steam explosion.

The SleiER-II equation of state and heat-transfer models were upgraded for this purpose. A calibration to SNL steam (xplosion data and a compariton with Los Alamos shallow-pooli experiments was performed. A lower-head failure and motion model was developed. Analysis of parametric cases suggests that the upper bound on the conditional probability of alpha-mode failure, given core melt, should be 0.01 if the vessel upper head and bolts are near normal operating temperatures. Suggestions are made for activities that sould lead to increased confidence in the understanding and quantification of steam-explosion issues.

I. SU41ARY REPORT A. Introduction i

The purpose of this work was to provide a reasonable estimate of the maximum loads that might be expected at the upper head of a pressurized water reactor (PWR) following an in-vessel steam explosion. Using the determined range of loads, judgments were to be made on the potential for containment failure by missile production resulting free a stene explosion, commonly called alpha-mode failure. assuming core c:elt has occurred.

i

DEFT '

2 The calculations of loads were performed by the SI WER-II code. A previous

{

study of this type, herein called the ZIP study', was performed in 1950. To l provide some basic SlWER-Il capabilities lacking in the 1980 study, a brief program of code modification was undertaken as part of this work. These modifications are discussed in Sec. B. To calibrate the revised modeling to available steam explosion experimental data, Sandia National Laboratories (SNL) intermediate-scale thermite-water experiments were assessed, and one particular well-characterized experiment, MD-19, was calculated with SIMER-II. This discussion is in Sec. I.C. Evaluetion of the adequacy of the SI WER-II post-explosion expansion dynamics was investigated by comparing calculations with scaled shallow pool experiments. The experiments and the analysis performed are given in Sec. I . D. As the study proceeded, numerous scoping calculations on a full-scale PWR representation were performed; the results of these calculations are in Sec. 1. E. These reactor results were integrated into an overall scheme for judging the containment failure probability, presented in Sec. I. F. Finally, conclusions and recommendations from this study are given in Sec. I . G. To keep the volume of this summary chapter from beicg excessive, most of the above outlined sections will present s ucu::a ry ma t e r i a l onl y. Details are given in subsequent chapters of this report. i B. Modification of SI4fER-II The SI4fER-II code was modified to better represent a corium-water systec, and to implement a lower head failure and motion model suggested by the ZIP study. The details of these changes are given in Chaps. II and III and Appendixes A to N.

The four major areas of code modification for treating the corium-water system are an ic: proved equation of state (EOS), a more appropriate nonequilibrium heat transfer to low thermal conductivity liquid water, revised assumptions in the vaporization-condensation model regarding energies at which water is vaporized and condensed, and miscellaneous changes.

'M. G. Stevenson, " Report of the Zion / Indian Point Study, Volume II," Los Alamos National Laboratory report LA-5306-MS, NUREG/CR-1411 (1950).

t D f I 1 l J~ \1 ,

DR/fT The most significant EOS modification in terms of accuracy was to make the infinitely dilut e, wa t er-vapor -hea t capaci ty tempe ra t ure dependent , which better r'e pr e s ent s the increasing degrees of freedom of the water molecule and the hydrogen-oxygen system at high temperatures. For 3000 K infinitely dilute steam, the new heat capacity relationship gives 2.62 times the value for infinitely dilute 400 K steam.

The major modification to the EOS in terms of coding was to insert a relaxation constant so that as steam becomes more sup:rheated its internal energy becomes closer to the value for infinitely dilute water vapor. This is most important for stability because it limits the triple-valued nature of the vapor temperature as a function of density. It is also important for accuracy, removing a spurious pressure dependence for superheated steam, as well as allowing the inclusion of the -830 kJ/kg energy difference that appears between steam at the critical point and infinitely dilute steam' at the critical temperature.

A fitting procedure was applied to the new EOS relationships to compare them with both the Steam Tables' and the los Alamos SESAME tables.' Agreement is within 120% over the range of interest and generally is much better. Given the phenomenological uncertainties in steam explosions, these results cxceed the required accuracy. The main conceptual problem remaining (the one that relates only to the EOS) is a better variation of the gas parameter, R, near the critical point and at high steam density and temperature.

The SIWER-Il heat-transfer problem with the corium-water system relates to the low thermal conductivity of water. As a consequence of model development based on liquid sodium, the standard SIMER-II code transfers heat to bulk liquid. which then vaporizes. For water, both experiments and theoretical

. considerations support the view that steam-explosion liquid-liquid heat transfer occurs to the water surface. Some heat then is conducted into the wa t e r ; the remainder causes vaporization. Implementation of this revised concept allows 2

3. H. Keenan, F. G. Keyes, P. G. Hill, and J. G. Moore., St eam Table s (New York, John Wiley and Sons, 1978).

'Kathleen S. Holian, Ed., "T-4 Handbook of Ma t e r ia l Prope rt ies Data Bases ,

Vol. Ic: Equations of State " Los Alamos National Laboratory report LA-10160-MS (November 1984).

3 --

1 i

  • I l

DHAN

  • l s in:ula t ion of steam explosions in water-rich systems with nonequilibrium '

vaporization.

l The radiation heat-transfer regime occurring during the coarse premixing phase of a steam explosion was not modeled explicitly. Bounding cases that J

envelop a variety of heat transfer modes and magnitudes were evaluated instead. ,

Modeling a film boiling flow regime with radiation and film layer conduction in these chaotic environments is beyond the scope of the SIWER-II formalism at this time and would add only marginal in:provement to the overall treatment by itself. The lack of physics for fuel stream breakup and water interpenetration of the fuel are probably more important deficiencies in the SI WER-II code.

We found that inclusion of water subcooling into an effective heat of vaporization was necessary to obtain reasonable results from the revised heat transfer model. Past experience has suggested that the complementary process, putting vapor superheat into the effective condensation energy, also would be desirable. Although conceptually straightforward, these changes did result in somewhat complex coding modifications.

The miscellaneous changes involved programming the liquid water thermal conductivity as a function of temperature. . changing the maximum water droplet size algorithm, inserting an additional time-step control for the vapor energy work term, and correcting of two minor FORTRAN coding errors.

The necessity for treating structural effects can be introduced with the observation that only energetic s t e ar:: explosions are a concern in evaluating containment failure. These energetic explosions could cause a rupture in the lower vessel head area that could be expected to mitigate upwardly directed kinetic energy significantly. .

Finite element calculations indicated that vessel failure could be expected at the radius where the outermost vessel penetrations are . located. A circumferential split would be expected to proceed around the head at this radius and leave the inner portion of the head as a free body. A s ic:pl e ,

single-degree-of-freedom s pr i ng-c:a s s model was correlated to these finite-element results. This model accepts the average pressure. loading over the failure region and furnishes the t ic:e and downward head velocity when the head disengages froc: the vessel.

7 F*

l>

i Motion of the head segment into the lower vessel cavity was treated by a modification of the SIWER-Il plug e ject ion model. The failed portion of the lower head becomes a free body whose acceleration is computed from its mass and the integrated forces acting on it. A two-dimensional representation of the cavity was used in this study; the volume available to the venting materials

  • assumed the below core cavity was half filled with water.

C. Calibration to SNL Steam Explosion Experiments The selection of SIWER-II input parameters for use in reactor steam explosion analysis can best be accomplished with reference to the experimental data base. Two difficulties exist. First. any SIWER-Il model will be simplistic. The ability to propagate a mechanistic fragmentation (detonation) wave leaving behind a distribution of particle sizes does not exist in the code.

In fact, a generally agreed upon theoretical fragmentation mechanism does not exist for steam explosions. Second, the reported experiments have random aspects. Later experimental results have tended to contradict earlier experiments and especially contradict early simplistic theories. Also, the experimental information is sufficiently ambiguous that any unique simulation is l impossible. Because of these problems, one typical experiment was selected for computational simulation using a sicple SIWER-Il representation. Details for this calibration as well as a brief review of SNL steam explosion data are given in Chap. IV and Appendixes O to T.

The experiment selected was SNL test MD-19,* where 5.11 kg of iron-alumina thermite was dropped into 224 kg of ambient temperature water contained in a square Lucite box. Initially, a uniform interaction zone of a size reported by the experimentalists was used. It was reported to have a thermite / steam / water volume ratio of 0.04/0.4S/0.4S. A two-phase, equal volume frcction steam / water chimney was assumed to exist above the interaction zone. The bottom of the box was assumed rigid. The sides and top of the box were assumed to be at constant pressure. These lateral boundary conditions assume the Lucite has no strength.

Alternative boundary conditions are discussed in Appendix S.

'D. E. Mitchell, M. L. Corradini, and W. W. Tarbell, "I nt e rcied i a t e Scale Steam Explosion Phe nocie na : Experiments and Analysis,*' Sandia National Laboratories report SAND 81-0124 NUREG/CR-2145 (1981).

i l

' ~ ~

j Q

J\\r{u

~

DRVT Best results in the S141ER-Il simula t ron were achieved by having a water droplet size at a fractional multiple of the fuel droplet size. This was assumed to be reasonable because water surface tension is about an order of magnitude lower than that of iron-alumina thermite or corium. However, the fit was less than satisfactory, with a spurious high pressure tail resulting from ~

bulk heating of water in the interaction zone.

A second fit then was made using a previously calculated fuel / water distribution based on a three-velocity-field film boiling model.' Here the outline of the calculated mixing zone was correlated to the vague pre-explosion experimental configuration. The node with the peak thermite density had a the rmi t e /s t eam/wa t e r volume fraction of 0.27/0.37/0.36. The best fit achieved was with a fuel particle diameter of 300 um , a water droplet diameter of 75 Se, and a liquid-liquid beat transfer multiplier of 0. 2. (Note that the new liquid-liquid beat transfer model was used with the 0. 2 multiplier. This multiplier cnly acts on heat transfer to the water surface.) A comparison to the base pressure transducer is shown in Fig. 1.

The final exercise using SI41ER-Il was to perform some pre-explosion coarse-mixing calculations, with initially separated thermite and water.

Although the functional forms of some of the physics modeled in SlyfER-Il is less than desired, obtaining an appreciation of the qualitative errors involved was desirable. These results will be used in Sec. I.E in running bounding cases fo- the reactor meltdown problec:. The first calculation was to assume a uniform 15-me droplet size with the nominal (unmodified) S141ER-Il beat transfer model.

Steae production and the extent of thermite dispersion were underestimated.

This calculation was rerun with the heat transfer models of Sec. I. B. Excessive radial thermite dispersal resulted, and the interaction zone did not reach the bottom of the Lucite container in the required time. .The . configuration established with the first premixing calculation was exploded with the same models and parameters used in the fit shown in Fig. 1 to see the effect of different precixing on the characteristics of the explosion. The peak pressure increased to 20.8 MPa, there was a longer time at high pressure, but the width of the pulse at 8 MPa was about the sace. These results are consistent with

  • W. R. Bohl, 'A Calculational Advance in the Modeling of Fuel-Coolant Interactions," Proceedings of the L. M. F. B. R. Sa f e ty Topica l Meeting. Lyon France. American and European Nuclear Societies, pp. 111-557 to 111-566 (1982).

+

200 , , , . ,

DVFT 1 TEST -

DATA m

E

< 100 -

m f .

W .

O N -

.O .00200 . 00400 TIME IN SECONDS 20 . . . . . . . . ......

,3 S IMME R - ll  :

4 -

Q. .

2 ,.

l . r -

I w

C IO t -

D y) -

U) .

w -

& 5'---

c. -

l $ .

\

O l.O 2.0 .

TIME , Ms Fig. 1.

Comparison of water chamber base pressures in expe r ic:ent MD-19.

i

~

NAFT .

what would be expected with less local vapor volume for early expansion and therefore pressure reduction. Although the calibration exercise was unable to provide a strongly supported. unique set of calibration parameters, it did provide some confidence that the SI WER-11 treatment was qualitatively reasonable. Also, knowledge was obtained on the choices of parameters that would envelop the real behavior.

As a caveat, details from other FITS experiments, such as double explosions or significantly more efficient explosions through the use of rigid walls, cannot be modeled mechanistically with SIWER-11. Parametric studies coupled with additional model development are required to provide high confidence that consequences of all possible s t e ac: explosion phenomena are included in any study.

D. Calibration to Los Alamos Shallow Pool Experiments To investigate a basis for judgments on the ability of SIWER-Il to calculate post-explosion slug breakup, a series of shallow pool acceleration experiments were conducted and then analyzed with SIMER-II. Details of these experiments and their analysis are provided in Chap. V.

A schematic of the test apparatus is shown in Fig. 2. The depth ~of the pool and the height of the free space above the pool are scaled using the actual reactor vessel dimensions and the total amount of fuel available. As indicated in the schematic, a very thin (0. 006-ce) diaphragm supports the pool in a 102-mm-i.d. Plexiglas tube. The pool is 50 en d e e p , and there is 185 cc of free space above the pool. The bottom of the tube is separated from the nitrogen driver gas by a 0.1-ce Mylar diaphragm. The experiment is started by cutting the lower diaphragm; the pressure rises in the bottom tube and ruptures

'the thin diaphragm that supports the water. The increased pressure then induces the motion of the water that is of interest. Three sets of data are collected during th- test: the pressure history just under the top diaphragm (P1), the pressure history at the endplate (P4), and a high-speed movie. Although some random behavior was observed, the results of these tests are still interesting.

T,bree types of SIWER-II calculations were performed, s'ic:u l a t i ng bot h r e a l and hypothetical experieents in this geometry. These were (1) calibration cases relating to one experieental test, (2) hypothetical cases with postulated curvature in the thin Mylar diaphrage supporting the pool, and (3) calculative coc:psrisons with other shallow pool experieents using the calibrated paraceters.

5 DRAFT

,- m a V4 mm water Mylor V/g dophragm 4 m2 Pi .

Myiar \

dophragm um er JL y Ndophrogm nitrogen  : M er e m v//,/////////,

31////// win T

Fig. 2. ,

Shallow pool experimental apparatus.

For the calibratior. cases, the initial nitrogen driver gas pressure was

0. 56 MPa , and the evacuated space started at -5.6 kPa. Consequently, these tests are referred to as having a 100:1 pressure ratio. The high-speed movies indicated progressive slug breakup as a function of distance, with complete breakthrough occurring at a b'o u t the time of endplate impact. A reasonable simulation of the experimental pressure trace was obtainable in SIWER-II by using a small drop, diameter, 100 ue, and consequently maintaining close coupling -

between the liquid and vapor fields. Also, the value of o o, the vapor volume fraction where the two-phase to single-phase transition occurs, was adjusted to 0.025 to attempt to maintain two-phase pressures as long as possible upon impact. The best comparison with experimental data is shown in Fig. 3.

Two types of hypothetical situations were examined. First, the inner one-half of the lower water interface area was raised or lowered I ce to evaluate the efforts of surface shapes on pool breakup. Second, one of these nonuniform situations was scaled to up by a factor of 40 to evaluate the

. 200 -

PRESSURE RESPONSE FOR TEST SP7014 160-TEST i

120-DATA BD-5 -

4D-a y

CD

,.-. i- i - . .- . .

. .. i CD LO 2.0 30 4D SD 6.0 7.0 8.0 TlWE (ms) 20.0 -

PRESSURE RESPONSE FOR TEST SP7014 16.0 -

S IMM E R - il e

12D -

3 t s.:

i a: - .

f D

8.0-E 4D- '

.i kl'Q i

l OD , , , , , , , , ,

l 0.0 1D 2.0 3.0 4D 5.0 6.0 -7.0 8.0 l

TIME (ms)

Fig. 3.

Coeparison of P4 pressures in the shallow pool experiments.

I

.J scalability of the treatment. The nonuniform interface resulted in an instability growth that de finit ely made impact appear as more of a two-phase spray, but actual venting of the gas through the water was not achieved. A further calculation with half of the distorted lower interface actually blocked by a rigid disk was performed to assess an experimental concept for testing the effects of interface distortion. The results suggested that the qualitative difference of the pool behavior was not changed significantly by the disk.

Therefore, a distorted interface experiment was planned. Increasing the scale (and driving pressure) did produce the expected increase in peak impact pressure by a factor of 40, but the upscaled case resulted in a more coherent impact.

The two remaininE calculations made in comparison to expericents were (a) one with a 50:1 pressure ratio and (b) one with a depressed lower pool interface by a central disk. The impact pressure for the 50:1 case, which was obtained with the " calibrated" parameters, had peaks more than twice than measured. The character of the 50:1 calculation was much like that at 100:1.

Although slug breakup in the high-speed movie of the 50:1 experiment was similar to that for the 100:1, a significantly more diffuse impact is icplied by the '

pressure trace. With a disk to lower the center of the lower water interface by 3 cm, both the experimental and the calculation showed gas breakthrough before impact. However, now the calculative impact was significantly more diffuse than that measured.

The following conclusions were reached in this brief SI4!ER-Il exacination of shallow pool dynamics.

(a) The SI41ER-Il simulation of shallow pool breakup provides generally correct global behavior in terms of diffuse vs slug impact and overall impulse delivered by the fluid to the upper closure. It can not, and probably need not, match the very short duration pressure spikes measured in the experiments.

(b) With an initially uniform pool, the use of small drops in SI41ER-Il probably would exaggerate upper head loadings in the reactor case (be conservative).

L

(c) Even with exaggerated liquid-vapor coupling, the lack of a model for turbulent mixing leads SIWER-II to overly accelerate slug breakup in the case where small nodes are used to simulate a nonuniform interface in this experiment.

(d) Because the' node size at reactor scale is much larger than the turbulent mixing length, exclusion of a turbulent mixing model appears justified.

Also, the non-linear scaling of the SIWER-II liquid-vapor drag relationship was a major contributor in increasing impact coherence of '

the upscaled calculation.

(e) Because of unresolved scaling questions, .further analysis of these experiments was judged less cost effective than parametric calculations involving reactor meltdown sequences.

E. SIW ER-il Reactor Case Calculative Results The range of uncertainty in the expected conditions resulting in a steam explosion during a core meltdown is large'. If all possible initial conditions and modeling uncertainties were considered, a complete parametric stud'y of steam explosions and their effects would be an enormous effort. The effort involved in even one SIWER-Il reactor calculation is appreciable. With the limited scope of this study, the limited understanding of steam explosions that exists, and the limited extent to which this understanding currently has been programmed into SIWER-II, a massive parametric study was not appropriate. Rather, five scoping calculations were performed with SIWER-II to address the upper bound question. Details of these calculations, ,,a s well as some additional results that help to integrate these cases into the context of a meltdown accident sequence are given in Chap. VI and Appendixes U and V. The results are ,

summarized in this section.

The reactor structural representation for these cases is shown in Fig. 4.

The in-vessel configuration is the same as in the previous ZIP study. A covable

'H. Berman, D. V. Swenson, and A. J. Wickett, "An Uncertainty Study of PWk Steae Explosions, Sandia National Laboratories report SANDS 3-1435, NUREG/CR-3369 (1984).

m I

UJ j

b p

A EEEE!! +

======
=::

60 ______

i 55llll -

OUTLET PIPING l

______ _x 50 SOLID S T RUC T U R E I I II

, /

I I II m I I i

12. 5 4  :::::: =

INLET PlPING 40 DOWNCOMER b33333

=

l

_0 =eam =ae=m 0 '

_ss_e_e_s_s

~~

IIII CORE SUPPORT 20 i '

M!!! --+- FORGING '

i3 wow o / + 1,< /lJl I

'EKw < '.D iX4 l l <

5 //W a /%f'c'11 I

/ .

l

?) / c .  ? M A4 l l re n a MAN MOVE ABLE

~

W :W.V/A1l

l 5.20  :'M M'/"h'44 I I STRUCTURE

'AM M '/2E/4'1 I I

4/d%) A/N .1 1 l l M WM4 4J l l

'2 4 /

/

/ i/X'/: 4 'T/l l I v o w/ '

< '44 i =

l i > OUTLET TO KEYWAY I 5 10 15

+1.613  !

- 2.537  :

ris. 4.

PWR structural representation for 51WER-Il stearn explosion calculations R ""

j

DMFT lower structure has been added for the lower head failure model. The five SI W ER-Il cases can be summarized as follows.

In case 1. we used the same SI WER-Il corium and water geometry as in Cases I and 2 of the 1980 ZIP study. Thi s wa s a pJemixture comprising 20% of the total corium below a single-phase molten corium pool. The corium/ steam / water volume fractions were 0.50/0.25/0.25. The fuel and water were assumed to be fragmented to 300 um -diam globules within the premixture. In addition, 35 000 kg of steel particles were added between axial node 46 and node 60 to represent the mass, but not the strength, of the upper core structure for this case. The new SI MER-II lower head failure, equation of state, and heat transfer models were used.

The results of the calculation indicated a substantial reduction (compared with the ZIP study) in the likelihood of energetic missile production. The new models eliminated the spurious high-interaction-zone, single-phase pressures previously obtained. Lower head failure was calculated to occur at 3.5 es, and about 571 of the kinetic energy was directed downward as a result. Upwardly directed kinetic energy at head impact was -650 MJ. With the assistance of the downward venting of core material following lower head failure, the peak force on the head was reduced from 2.6 GN (giganeutons) to 1.0 GN. A plot of the integrated upper head loading is shown in Fig. 5. This head loading is near the threshold for failure if the head is near the normal operating temperature.

In case 2, we used more c:e chani s t ic but arguably conservative initial conditions for a large-scale coherent steam explosion. Starting with a l

cocpletely molten core at 3100 K, a 1.85-e-diam corium strean was allowed to l

pour into 18 000 kg of water in the lower plenum. Standard SIW ER-Il beat transfer was used with 20-ec-diam, prefragmented corium globules. Based on the results in Sec. 1.C. these assumptions. should under. estimate steam production and thereby permit extensive mixing. Mixing occurred around the edges of the corium streac, and was reasonably extensive because of large mesh size used (greater than 10 cc:).

An explosion of the c:ixed corium and water was triggered 0. 7 s after j initiation of the pour, when 40 000 kg of corium had entered the lower plenum and c6riue flow reversal was starting as a result of vapor production and plenum pressurization. This tic:e was chosen because the c:ixing configuration appeared to have the potential to produce caximum head loads. Examining details of t icing sensit ivity was beyond the scope of this study. but an explosion after

~P f

LJ b]n? \

b- 22.0 ;

11.0 ;i FORCE ON HEAD 10.0 ~

9.0i

? 8.0 '

c  :

3c 7.0 :~

v -

5 6.0 '

e O

c:

O 5.0:i A 4.0;:

3.0 -i

. 2.0i 1.0 -

0.0I .., . ,, ... .. ..... .... .... ......... ... ,

0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 TIME (s)

Fig. 5. Integrated upper head loading for case 1.

I s of mixing is discussed in Chap. VI. The SNL calibrated explosion parameters were used with the new SIWER-Il models. Steel particles were included to represent the upper internal structure. About 62% of the kinetic energy was directed downward. Upwardly direct'ed kinetic enerjiy at beat impact was -500 MJ.

The peak upper head force was 0.81 GN. as shown in Fig. 6.

In case 3. we used worst-case initial conditions as inferred from a recent SNL study'. Here 75% of the core (94 000 kg) was assumed to mix with 20 000 kg of water. To limit the height of the premixture to that of the original corium pool, a 19% initial steam volume fraction was assumed. The corium globule size was assumed to be finely frageented to 100 um in diameter. although the n e'w

' Ibid.

)

b 10.0 =:

.m i FORCE ON HEAD l 9.04: l 5

8.0i .

i 7.0 -!

p 6.0 [

i d

2:

5.0 -!

2 O I A 4.0 -2 i:

3.04

  • i:

2.0i :

5 1.0 -j 3

0.0 ...

.....c... mc ..c. ...c........c........c........c......c ......c,.m .,

0.70 0.71 0.72 0.73 0.74 0.75 0.76 0.77 0.78 0.79 0,80 0.81 TIME (s)

Fig. 6. Inte. grated upper head loading for case 2.

SIl44ER-Il models were used, no beat transfer f rom frozen corium t o , s t e am wa s calculated, limiting steam temperatures during the expansion to 922 K. No lower head failure was assumed, and the remaining core material and upper structure was assumed not to inhibit e x pay.s i on. The upward 1'y directed kinetic energy was 1830 MJ. However, as a consequence of the diffuse nature of the expansion, the l

peak upper head force was only 0.75 GN, as shown in Fig. 7. This result shows that a complex relationship exists between the explosion and its delivered loads to the system. Sic:ple correlations are not reliable.

In case 4, we took the initial conditions from case 2 but homogenized the 40 000 kg of corium and 10 000 kg of water in the lower plenum. A slow explosion" then was simulated where droplet sizes were increased an order of magnitude over the SNI. corre lat ed va lues. This reduces the heat transfer rate by 2 orders of magni tude. The idea was to simulate an incoherent, cultiple h

Ji k

f I

b- 9.0-l i FORCE ON HEAD 8.0 -'

1 7.0i 6.0i .

E 5.0 -'

r.o  :

u -

A 4.0 -'

3.0-2.0-'

1.0-0.0'................................................

0.00 0.02 0.01 0.03 0.04 , 0.05 TIME (s)

Fig. 7. Integrated upper head loading for case 3.

explosion environment that could be more representative of the reactor situation. In case 4, the lower plenum was not calculated to fail. Upwardly directed kinetic energy at beat impact was -760 MJ. The peak upper head force was 1.52 GN. as shown in Fig. 8. This case had "a significantly more coherent upper head impact than cases 1, 2, or 3 even though the explosion was relatively benign. Again, the complexity of relating explosion magnitudes and characteristics to loads on the head is evident. Here the most benign explosion produced the largest challenge to the head.

In case 5 we used the case 3 premixture but with corium and water thermally equilibrated before the expansion. The Low Alamos SESAME tables were used to give a starting pressure of 900 MPa. Refitted parameters for the SIW.iER-II EOS were derived to obtain these initial conditions. The remaining 25% of the core was placed on top of the premixed region. Both lower head failure and steel

-1S-D- 16.0 -

i FORCE ON HEAD 14.0 --

12.0 I m .

g 10.0 - i a . I B  :

5 a.0 -

co U -

Es.

6.0i .

4.0 -

2.0 ..

0.0 ', , ... ......,. .,. ..,,.. .. .......... .,

0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 ,

TIME (s)

Fig. S. Integrated upper head loading for case 4.

particles (for structure) were included in the calculation.' The maximum upwardly directed kinetic energy was 7250 MJ. The peak upper head force was 12.4 GN, which was obtained just before the calculation went unstable. The .

force plot is shown irI Fig. 9. This resulf' confirmed clearly that our calculated approach would produce the expected diasterous results for the upper limit assue:ptions and idealized physics.

F. Containment Failure Probabilities To obtain an estimate for the probability of containc:ent failure from an in-vessel steam explosion, given core melt, we must (a) judge how likely the initial conditons for a large scale steam explosion are, (b) e s t i c:a t e the in-vessel energetics and loadings produced by the spectrum of possible steam .

RA-

h U I V-II bw 16.0 -

FORCE ON HEAD 14.0 .

12.0 i p f

sn .

8 10.0 -

s .

D .

0 .

5 8.0 -

r.e U -

@ 6.0 '

h- .

4.0 .

2.0 ..

0.0 ' .. , .

0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 , 0.016 TIME (s)

Fig. 9. Integrated upper head loading for case 5.

explosions, and (c) evaluate the consequences of such loadings. A'brief study attempting to make such judgments is presented in Chap. VII. Here we summarize some conclusions from that study.

First, although th'e best-estimate probabili~ty is judged to be low (between 10 - 10), this is simply a guess. We currently lack the technology to construct a reliable probability density function for core melt consequences that would portray a best-estimate or most probable valu'e. Indeed. perforcsinE any "best-estimate" mechanistic core meltdown calculation, with models representing consensus phenomenology. was not possible because consensus phenomenology did not exist. The SIMER-II calculations performed were for addressing the upper bound question.

3

Second, any large-scale steam explosion or sequence of explosions leading to a containment challenge apparently must involve a sustained supercritical pressure. Although IDCOR aserts the maximum meaningful pressure during the expansion phase of an explosion is about half the thermodynamic critical pressure,' steam explosion experiments exist with expansion phase pressures greater than 34 MPa sustained for more than 1 s.' Thus, containment challenge from large-scale steam explosions cannot be ruled out on the basis of limited pressures.

Third, with the dissipation calculated by the SIMER-II model, to obtain sufficient energy a large-scale steam explosion apparently must involve efficient energy transfer from a considerable fraction of a molten core, 20% or more, in a few tens of milliseconds. This leads to such requirements as forming a large molten cor iuc: pool, having a coherent pour, not obtaining early triggering of an explosion, and limiting steam generation so as to permit the extensive premixing of corium and water. All of these are unlikely, but not outside the spectrum of reason. For example, triggering could be delayed by saturated water or elevated ambient pressures or the commonly quoted steady-state fluidization arguments, which limit mixing, could be defeated by the inertial effects involved in a large corium pour. Because of the edge-of-s pec t rue: character of these phenomena, a 0.1 probability was judged to be a proper upper limit to associate with obtaining initial conditions resulting in sufficiently energetic steam explosions. Further explanation as to the meaning of this probability value is in Chap. VII.

Fourth, obtaining the extreme conditions presented in case 5 appears to be I

impossible. Mixing 75% of the core requires water to rise through a molten corium pool in a file: boiling regime because of the necessity to accommodate the l vapor volume produced. In any. meaningful calculation attempting to attain this 1

l prec:ixed configurat ion, satisfaction of fluidization requirements resulting in 1

'Fauske and Associates, Inc. " Technical Report 14.1A Key Phenomenological l

Models For Assessing Explosive Steam Generation Rates," The Industry Degraded Core Rulec:aking System (IDCOR) (1983).

  • K. H. Wohletz, Los Alamos National Laboratory. Private Communication.

Some background on these experiments can be found in Michael F. Sheridan and Kenneth H. Wohletz, "Hyd r ovo l c a ni sc:: Basic Considerations and Review," Journal of Volcanology and Geothermal Research 17, 1-29 (1953).

l t

% me 1

a li i

corium dispersal would be assured. The assumptions of instantaneous temperature equilibration and pessimistic equation-of-state parameters are also unreasonable. Case 5 merely demonstrates that if arbitrary mixing is allowed, and the explosion efficiency is allowed to approach the thermodynamic limit, S144ER-II obtains the expected unacceptable result. -

Fifth, if the vessel head is to become a large missile, an integrated threshold head loading of approximately 1 GN is required if the head and bolts exist at the reactor operating temperature. Case 2 run with edge-of-spectrum initial conditions for a single large explosion, was below this level. Case 1 produced 1 GN, although its postulated premixing as well as its explosion / expansion sequence involved edge-of-spectrum considerations. Case 4 achieved more than 1 GN. However, true multiple explosions are not calculated.

The initial conditions are idealized (for example, with uniform material distributions radially), dissipation of energy in the upper core structure is neglected (with 900 M1 of fluid kinetic energy, this could be significant), and the time-dependence of the nonuniform upper head loading suggests formation of multiple missiles from the head spex, which could be stopped by the missile shields. In brief, case 4 is an edge-of-spectrum expansion calculation.

Because of the end-of-spectrum character of calculations exceeding 1 GN, if the head and bolts exist at operating temperature, a 0.1 probability is judged appropriate for obtaining vessel loadings of concern given initial conditions resulting in energetic steam explosions.

Sixth,- any attempt to establish a simple limit condition is fuzzy. An impulse limit is useful only if the impact involves single-phase liquid. Case 3 shows that high energies do not necessarily give unacceptable loads, if a highly two-phase spray is the impacting fluid. The 1 GN force limit is fuzzy because of the influence..of the s pa t i.a l loading distribution, the temperature of the vessel, and the response of the containment to any missile. Head loadings depend on the assumptions made in the explosion / expansion calculations. In this study, the presence of steel particles tended to concentrate loads, while lower head failure made upwardly directed impacts less severe. Thus, consistent analysis must be performed that links initial conditions, mixing characteristics, explosion characteristics, expansion behavior, head failure characteristics, and missile dynamics.

l

. G. Conclusions and Recommendations This study has resulted in significant progress over the ZIP analyzes. The SI MER-Il water EOS and heat-transfer mechanics are more credible.

The significance of lower-head f ailure can be evaluated consistently and shown to .

result in considerable mitigation for the constrained (by the corium pool) reactor meltdown case. The initial conditions used for a steam explosion can be related to a mixing calculation. Although SIMER-II models can only parametrize  !

the physics of steam explosions, the correlated models can now obtain a respectable pressure pulse for the experimentally observed water-rich system. l Further, parameters can be chosen to provide a calculational simulation of shallow pool behavior. The use of steel particles for upper structure reduces fluid breakup and dispersal which were nonconservative .in the ZIP analyses.

Because of the edge-of-spectrum character of both the calculative assumptior.s leading to a large scale stream explosion and those required to obtain significant head leading, the upper limit for the containment failure probability given core melt appears to be 10-8, if the upper vessel head and bolts exist near -550 F. For higher temperatures this probability must be increased. The importance of the high loadings in case 4, obtained , f rom slow heat transfer simulating incoherent explosions, cannot be ignored. A method to quantify a true best-estimate probability does not exist. Present guesses suggest values of 10~3 to 10" might be reasonable.

A detailed discussion of possible future research to obtain increased confidence on the steam explosion issue is presented in Chap. VIII. Beyond an extension of this study to better quantify the consequences of the caltulated upper head loadings, the proper priorities for additional research on the alpha-mode failure issue are judged to be (1) modeling of the meltdown process.

(2) large-scale steam explosion experiments with associated model development l for test analysis, test interpretation. and reactor application, and (3) smaller scale experiments to address specific mechanisms. The reasons for this judgment

(- are as follows, l

l 1. Understanding meltdown behavior would help place the alpha-mode of containment failure in its proper perspective among other safety issues.

mr.

l R *o;

l DWT 2. Both the state of any corium pool oefore contact with water in the lower plenum and the mode of contact depend on the meltdown sequence, which is uncertain.

3. The state of vessel internal structures as well as the temperature of the vessel head and head bolts depend on heat transfer during the core heatup (meltdown) phase of an accident.
4. Both SI41ER calculations" and the theories of Corradini, Fauske, and Theofanous lead to the conclusion that large-scale explosions will possess characteristics that differ from the smaller scale FITS tests. These differences are sufficiently qualitative to be' observable with current technology.
5. Ultimately, reduction of uncertainty in the vapor explosion field is believed only achievable by inclusion of increas.es in scientific knowledge, not in the performing of parametric studies.

! Finally, a judgment has to be made as to what level of residual uncertainty is tolerable. Otherwise, the steam explosion issues will tend to remain open-ended, i

i "This is shown when the calculations for test MD-19 are compared to Case 2 of Sec. V. Further discussion is in Chap. VIII on research priorities.

Norman A. Evans, "The Tf fect of Core Melt-Coolant Interactions on Severe Accident Risks in Light Water Reactors,' in Proc. International Meetint on Licht Vater Reactor Severe Accident Eva l u a t i or. , American Nuclear Society, 18.1-1 to 18.1-8 (1983).

l RAR

N' s  ? . l f.

t SANDIA NATIONAL LABORATORIES Albuguerque, NM 87185 Date: November 20, 1985 To: M. Berman, 6427 p

) _

From: M. F. Young, 6425 -

f

Subject:

Evaluation of Pressure Transducer Measurements in HEX and SHE Tests 1.0 Introduction This mem'o is an evaluation of the ur. essure traneducer measurements taken in the SHE (Straight Hish Explosive) and HEX (High Explosive) -tests performed at the FITS facility at SNLA. These tests were performed to get data on reliability and reproducibility of p'ressure gage readings obtained during the FITS tests in the Steam Explosion Program.

Several different gage types and mo6nting techniques were tested. An additional purpose was to check on the accuracy of calculations performed with the CSQII[1] two-dimensional hydrodynamics code to analyze the FITS data.

l The HEX series of tests were done using a wooden mandrel 11 inches long and 4 inches in diameter Wrapped'with 11 turns of primacord as the pressure source. This mandrel was l

placed in a 22 inch ID steel cylinder fille,d with approx-imately 2 feet of water in tests 1 through 9, and in a 22 l

l inch square lucite chamber, also containing 2 feet of water, in tests 10 and 11. The mandrel was positioned.either with the top of the mandrel level with the water surface, (high mount) or with the bottom approximately 4.5 inches from the bottom of the tank (low mount) . In all cases , the primacord was ignited above the water surface; for the low count cases, a single strand of primacord went down through?the water to the nandrel along the axis of the tank. Ch'arge weights rangem from 120 grains 1 to 1000 grains; all charges were PETN except for the 1000 grain charge, which was ilDX, a l somewhat hotter explosive than PETN. Table l.1 shows the actual test conditions in more detail.

e

1. 1 grain = 1/700 lbm or x 1

-t

'i M. Berman November 20, 1985 The SHE tests were begun with the intention of providing a more symmetric geometry that could be analyzed with CSQII,

.as the HEX tests were actually three-dimensional, with the primacord burning in a spiral around the mandrel. This burning pattern was shown to cause asymmetries in the pressure pattern around the tank [2] , as well as being a difficult source to model.with CSQII using the existing JCL equation of state for PETN. The SHE tests consisted of a single strand of PETN primacord along the axis of the steel tank, which was filled with 2 feet of water. The strand was again ignited above the water surface. The portion in the water was either 11 inches long (high mount) or 20 inches long (low mount) . These tests are also shown in Table 1.1.

2.0 Evaluation In general, although the SHE tests are geometrically symmetric, the consistency of gage response is not good in terms of pressure amplitude; the shapes of the pressure traces and arrival times are in good agreement, both between symmetrically placed gages and with CSQII results. Measured peak pressures also do not agree with CSQII predictions, i

being one-half or less of the calculated pressures. These results indicate that some problems exist with the gages, the principal candidates being tank vibrations coupling to the gages and gage resonance. Some characteristic times for the SHE and HEX test geometries are shown in Table 2.1.

The steel tank fundamental frequency, about 3kHz [3] , has a time constant of 330 s, which is considerably longer than the SHE.primacord burn time of 40 s; this,.however, does not rule out tank vibration as a potential problem, since the fundamental frequency is that of free vibration and we are mostly interested in the transient vibrational response of the tank induced by shocks for this problem. The gage resonance is reported to be 700kHz, and the rolloff fre-quency of the data acquisition system at the FITS site is 50kHz. This latter frequency limit is-due to the amplifier response. These frequency limits mean that any signal components with time constants shorter than about 20 s will be attenuated, and that gage resonance could be excited by signals with substantial components near 700kHz, or 1.4 s.

Unfortunately, this is close to the rise time of the main explosion pulse in the SHE experiments (10 s).

Also shown in Table 2.1 are the transit times for various paths in the steel tank, plus the digitizer holdup time. This latter time should be added in when calculating arrival times for the pressure transducer locations.

Figure 2.1 shows the pressure trace from the lower wall position in SHE01. The high frequency hash visible just

is it

, a+th .s. *

[t

\ '

li 3

}

l

s N t)

~ M.s Berman November 20, 1985 before the first large peak and leading into it is a result of the gage responding to longitudinal vibrations carried through the w411. This hash appears to be a mixture of fre-

'(  ; Jr quencies from 100-500kHz with amplitude *50 psi. Although this' response is not terrible by itself and can be easily A filtered out in this case, it does indicate that the gages will respond to wall vibrations. The strength and frequency a of wall vibrations is a function of the forcing pressure (the explosion) and the location at which the vibration is measured. Thus, if the gages respond in this case, they may also respond at'other frequencies and. amplitudes for dif-ferent explosion strengths. Figure 2.2 shows the overlaid pressure traces from the lower wall position in SHE05 and SHE06. SHE05 shows a pronounced 50kHz oscillation with amplitude $500 psi, probably gage resonance induced by the shock wave arriving through the water. Figure 2.3 is the

. same data af ter applying a lopass digital filter [4] to the l' original data (rolloff at 40kHz, 20dB down at 50Khz).

Although the 50kHz signal is removed, there still appears to be an oscillation at about 25kHz. This latter cannot be removed effectively because it is of approximately the same frequency as the pulse and reflections from the explosive source in the SHE tests.

. There is also some indirect evidence for gage resonance

s being a problem in the SHE tests in the form of many blown 0 gages, although CSQII (and calculations assuming an' inverse distance dependence for the pressure amplitude) indicate a maxim 6m wall pressure of less than 4ksi, and the. gages used are rated'at.either 5ksi, 10ksi, or 20ksi. The above implies

~

s resonance in the gages induced by the short rise time of the initial shock front in the water.

, , An example CSQII calculated pressure is shown in Figure l .' 2.4 for the low wall gage position. Comparisons between i s Figures 2.2, 2.3, and 2.4 indicate that the CSQII predicted r pressure is about 30% greater than that observed in the SHE tests on the first peak. The following reflection peaks are also higher.

The general conclusion that can be drawn from the SHE l test data is that, although the SHE tests are geometrically

[ symmetric, the source pulse is too fast for the instru-l mentation to respond, and in fact may cause destructive resonances in the gages. There may be some effect on the j gages due to wall vibrations, but this is not clear from the

! evidence at hand. Such an effect can possibly be measured by

! using blocked-off gages, although this procedure will remove l a possible coupling path from the wall to the gage through l

the water, as well as the direct source pressure. .,

1 S

~

i 10 l 'd 's l

,l i J . _ - - - _ - - .-

'.- M. Ber=an November 20, 1985 3.0 HEX Tests In the HEX tests, the total burn time is 580 s, giving a much wider overall pulse with a slower effective rise time of about 140 s. Since this wide pulse is actually made up of

~

11 overlapping narrow pulses from the 11 wraps of primacord, there are some irregularities present (see Ref. 2), and gages can still be blown, especially if they are located looking directly at the side of the mandrel. How uniform the pulse appears depends on the gage distance and orientation, both from the mandrel and from reflecting and free surfaces in the tank.

Figure 3.1 shows the data from HEXO7 at the low wall position after first being filtered with a 50kHz lopass digital filter to remove high-frequency noise. The pulse clearly shows the structure consisting of the individual pulses from the primacord loops. Gage positions farther away and closer to the mandrel axis see more of a smeared pulse; that is, the individual pulses from the wraps overlap more.

These smeared pulses tend to be more consistent from gage to gage and easier to interpret that those showing individual wrap pulses. Base-mounted gages generally see a very con-sistent pulse in the high-mount HEX tests. Figure 3.2 shows an overlay of three symmetric base gage positions from the HEXO7 test, showing good agreement. High wall positions in the high mount tests generally see two distinct pulses, one from the source and a reflection from the base of the tank.

Table 3.1 contains the test number, gage location, peak pressure of the main pulse, and the impulse of this first pulse from HEX tests 1 through 7. Table 3.2 is the corres-ponding data from HEX tests 10 and 11 in which lucite was used as the water chamber. Missing gages in the two tables had experimental problems, such as loose wires, destroyed during the test, or very low readings compared to the other gages,. indicating either blockage or an air bubble in the gage block. Some of the gages that are shown also exhibited problems (see notes in tables). The next table, Table 3.3, contains statistical data on the HEX tests for from two to four gages in symmetric positions. Gages not included either had problems or did not have a matching gage. The values shown are all for base gages, since at least one symmetric wall gage did not survive in each HEX test. The "*" numbers are the standard deviations, calculated using the estimated mean. The only two tests to be duplicated, 2 and 3, do not have the same pressures or impulses.

O _

^

. M. Berman November 20, 1985 l

The pressures in the steel tank HEX tests, 1-9, are seen ,

to have errors in the range *4 to *15%, depending on t? e test, with the corresponding impulses having errors in the range *2 to *11%. This level of accuracy is not apparent in the two HEX tests done in a lucite chamber, 10 and 11. An overlay of three symmetric base gage positions in HEX 10 is shown in Figure 3.3. I would say that the agreement between the gages is nonexistent. Apparently, the walls are moving sufficiently over the duration of the pulse to affect the pressure at these gage positions, although, since we have only two tests to look at, a high mount and a low mount, further tests will have to be done to verify this hypoth-esis. What can be seen from the lucite chamber tests is that the pressures at similar gage locations are an order of magnitude lower in the lucite (free boundary) tests, than in the steel tank (rigid boundary) tests (see Figure 3.4). This points to the importance of considering system geometry and confinement when predicting pressure response at a given location.

4.0 Conclusions and Recommendations As we have seen, the SHE tests have too fast a pulse for the instrumentation, and gages in the HEX tests can give different pulse shapes depending on gage location. A solu-tion to these particular problems is to use a true slow-burning propellant as the pressure source; the resulting pulse should then be symmetric, with a pulse duration of around Ims, and a rise time of around 250 s. This source will approximate the steam explosion source without the side effects of the primacord wraps in the HEX tests. The first test using this idea has been done in the steel tank and designated PR0pp. The charge used consisted of 15 inches of MS-B propellant contained in a 2 inch ID PVC pipe, ignited with a' pyrotechnic " rip" ignitor along the axis of the pipe.

Results from this test are shown in Figures 4.1-4.4. The pressure traces are seen to be uniformly consistent, regardless of position or mounting. The first sharp pulse occurring at O.5 ms is from the rip ignitor, not the main propellant charge. The main pulse traces, beginning at around 0.8 ms and lasting until 1.5 ms, show some structure differences due to echo and timing differences between gage location, but have envelopes that are very similar between all gage locations. The peak amplitudes, in particular, are

, virtually identical.

Another potential problem considered was the possibility of wall vibrations affecting gage response. For this reason, special isolation mounts were fabricated which hold the gages in frames weighted with lead bricks; the gages then poke through the tank wall and are sealed to the tank with

9

^ .* M. Ber=an November 20, 1985 Silastic compound. It was felt that these mounts would remove any possible wall-to-gage coupling, except locally through the water itself. Comparison of isolated gages with wall-mounted gages in Figures 4.1-4.4 shows no significant difference, although there is slightly more spread in the wall-mount gage readings, and they are slightly higher, than the isolated gages. More tests will help to resolve these differences.

The poor consistency in the lucite tests will neces-sitate more tests, using propellant charges and isolated gages in the chamber sides. Possibly, the location of gage mounting will have to be re-examined for these. tests; the British WUMT tests [5] , for instance , have a similar, but not identical, design, and seem to give consistent pressure data.

One further item that should be considered in-regard to steam explosions is the difference in behavior of explosive-type pressure sources as oppos=d to propellant-type sources resulting from the type of cor.i n.ement. The basic difference is one of scale;. explosives, at typical real-world scales, are said to detonate independently of confinement, although there is indeed an effect if the scal e is small enough; propellants are said to burn and definitely are affected by confinement at real-world scales, although the effect dis-appears at large enough scale. The point of the above comments is that steam explosions take place with a char-acteristic time typical of propellants, not explosives; thus at FITS spatial scales, using 2-50kg of thermite, the efficiency of burning is probably a strong function of confinement. We can'then expect efficiency to rise with increased confinement at this scale, as seems to be true in RC-2 {6] , but the effect should level out at some larger scale. We would, of course, like to know how large that scale may be, since we are ultimately interested in reactor size masses of corium and water.

5.0 References

1. CSQII Manual.
2. K. Schoenefeld, " Spiral Wrapping of PETN," Memo dated July 23, 1985.
3. W. Soedel, "A New Frequency Formula for Closed Circular Cylindrical Shells for a Large Variety of Boundary Conditions," Journsl of Sound and Vibration, Vol 70, pp. 309-317 (1980).
4. S. Stearns, Digital Signal Anal ysis , Hayden Book Company, Inc., (1975).

. M. Berman November 20, 1985

5. %UMT Data Reports, AEE Winfrith.
6. M. Krein, Estimation of the Conversion Ratio for the RC-2 Experiment," Memo dated May 17, 1985.

M. Berman November 20, 1985 I

Table 1.1 Independent Parameters for all Tests Prima-cord Total Charge Water Depth Charge Test (grains /ft) (grains) (inches) Location HEX 01 80 1000 24.0 High HEXO2 30 360 26.0 P.igh HEXO3- 30 360 26.5 High HEXO4 20 240 25.5 High HEXO5 20 280 24.0 Low HEXO6 10 120 25.0 High HEXO7 10 120 26.0 Low HEXO8 10 120 23.0 High -

HEXO9 10 120 24.0 Low HEX 10 20 240 25.0 High HEX 11 20 240 26.0 Low SHE01 100 92 24.5 High SHE02 100 167 24.5 Low SHE03 100 92 24.0 High SHE04 100 179 25.5 Low SHE05 100 92 24.0 High SHE06 100 92 24.0 High

~

4

.- M. Berman ~9- November 20, 1985 Table 2.1 Characteristic Times in HEX-SHE Tests Distance Time Description , (cm) 1psl 18" Primacord in air 2 45.7 70 11" Primacord in water (High) 27.9 43 20" Primacord in water (Low) 50.8 78 High mount SHE rise time -

10 Digitizer holdup time -

64 Tank center to wall 3 -

28 156 Wall to wall 56 311 High SHE to Base Center 34.3 191 Low SHE to Base Center 11.4 63 Base Center to Base Inner 15.3 85 Base Center to Low Wall 30.7 171 Base Center to High Wall 60.2 334 HEX burn time 365.8 563 HEX rise time -

141 High HEX to Base Center 46.9 261 High HEX to Base Outer 54.3 301 Gage Resonance -

2 (700kHz)

Amplifier Rolloff -

25 (40kHz)

Tank Fundamental Mode -

330 (3kHz)

2. Detonation velocity = 6500m/s
3. Water shock velocity = 1800m/s

M. Berman November 20, 1985 Table 3.1 Pressure and Impulse for Steel Chamber HEX Tests HEX Pressure Impulse Test WP/ Type / Loc (MPa) (MPa-s) 01 1/FCP/BC//D 29.7 (flattop) -------

3/K/BNE/0/T 29.7 (flattop) -------

4/K/BSE/0/T 31.4 (flattop) -------

5/E/BNW/0/T. 28.3 0.00644 6/PS/WNE/U 17.9/31.7 0.00726 7/E/BSE/0/D 30.0 0.00710 10/E/BSW/0 27.9 0.00605 11/K/WNE/L 37.2 (flattop) -------

02 1/FCP/BC//D ----- -------

2/E/BNW/0/D 4.4 (bubble) -------

3/K/BNW/0/T 16.6 0.00445 4/K/BSE/0/T 22.1 0.00440 5/E/BSW/0/T 20.0 0.00446 6/PS/WNE/U/T 6.2/9.0 0.00437 7/E/BSE/0/D 17.9 0.00427 11/K/WNE/L/D 55.24 -------

03 1/FCP/BC//D 11.0/28.3 (flat) -------

3/K/BNE/0/T 22.8 0.00565 4/K/BSE/0/T 23.4 0.00491 5/E/BSW/0/T 25.5 0.00607 6/PS/WNE/U/T ----- -------

7/E/BSE/0/D 22.1 0.00518 11/K/WNE/L/D 37.9 (flat) -------

04 1/PCB/BC//T 11.7 0.00382 2/E/BN/I/D 12.4 0.00293 3/K/BNW/0/T 13.8 0.00358 4/K/BSE/0/T 18.6 0.00368 5/E/WNE/L/D 10.3/24.1/13.8 0.00332

. 6/K/WNE/U/T ----- -------

7/E/BSE/0/D 16.6 0.00346 11/E/BS/I/T (bubble) -------

4. This gage, WP11, is a Kulite HKS-375, Serial No. 1090 10 5ksi range (34.5MPa) It flattopped in HEXO1, is here reading twice its range, flattopped again in HEXO3 and did not respond in HEXO4. Evidentally broken.

, M. Berman November 20, 1985 Table'3.1 (continued)

Pressure and Impulse for Steel Chamber HEX Tests HEX Pressure Impulse Test WP/ Type / Loc (MPa) (MPa-s) s 06 1/PCB/BC//T 8.3 0.00291 2/E/BN/I/D _13.1 -------

3/K/BNW/0/T 14.8 0.00329 4/K/BSE/0/T 14.8 0.00317 5/E/WNE/L/D 10.8/12.6/7.9 0.00268 6/E/WNE/U/T 4.5/12.4 0.00231 7/E/BSE/0/D 13.8 0.00267 11/E/BS/I/T (bubble) -------

07 2/E/BN/I/D 20.6 0.00375 3/K/BNW/0/T 19.7 0.00373 -

4/K/BSE/0/T 21.2 0.00386 7/E/BSE/0/D 19.4 0.00343 8/K/WSW/U 7.9 (low) 0.00063 9/E/WSW/L 9.0/12.1/11.0 0.00274 10/E/BN/I/T 19.3 0.00371 Nomenclature for WP/ Type / Location in Table 3.1 The numb'er is the WP channel number.

Under Type (first slash), PCB, FCP and PS are as shown, E means Endevco, and K means Kulite.

Under location (second slash), The first letter is either B (Base) or W (Wall). The letters after refer to the compass position (N = North, NW = Northwest, etc., and C = Center).

The letter after the next slash is:

0 = outer, I = inner, U = upper, L = lower.

The letter (if present) af ter the last slash refers to gage orientation for gages mounted in a gage block: T = tee, D=

direct.

i M. Berman November 20, 1985 Table 3.2 Pr.ssure and Impulse for HEX Lucite Tests HEX Pressure Impulse Test WP/ Type / Loc (MPa) (MPa-s) 10 1/PCB/BC//T 2.8 .000613 3/K/BNW/0/T 1.2 .000187 4/K/BSE/0/T 3.1 .000521 6/K/WW/U 1.6 .000305 7/E/BSE/0/D. 0.7 .000122 9/K/WE/L 2.2 .000130 11 2/E/BN/I 3.0 .000672 3/K/BNW/0/T 4.1 .000774 4/K/BSE/0/T 3.4 .000548 5/E/WW/L/T 0.8 -------

6/K/WW/U/D 19.3 -------

8/E/WE/U/D 1.3 -------

9/K/WE/L/D 4.3 -------

11/E/BS/I/T 11.0 .001544 Table 3.3 Pressure Statistics for HEX Tests HEX Gages Pressure Impulse Test Averaged (MPa) (MPa-s) 01 4,5,7,10 29.40*1.6 .00678*.00066 O2 3,4,5,7 19.1*2.4 .00439*.00009 03 3,4,5,7 23.4*1.5 .00545*.00052 04 3,4,7 16.3*2.4 .00357*.00011 06 3,4,7 14.5*0.6 .00304*.00033 07 3,4,7 20.1*1.0 .00367*.00022 2,10 19.4*1.7 .00373*.00003 10 3,4,7 1.7*1.2 .00028*.00021 11 3,4 3.8*0.5 .00066*.00016

-=v M. B e an November 20e 1985 d

k

?

co u

  • C

'*T T

~

d E o >>

o U2 -e o

LD

, y u 4

.J c m

o e si k o u

.N c e * $

Y o .

.)

s

- v2

, , m o~ S m

5 ( _

- j ,, o

'E ..,

tw

. ".b,

. y

'O

- a s

-. 4 r

&  ; o o

~ L 1

-W o m -

m D ' m 0

h u Jl L o.

c.

u a

T c

o* a

  • o (* o C

h' w o .-

s

,o o c 06 D

005$1 0050 0 dos 0obot g1sdj esnssesd -

e c

b assi en e in, it wee , nn, ,

e a vuste a=* ew

  • eher u.

i b

r SHC05 WPOS KUL!TE, LOWER SC CYLlHOER E

SHE06WPO4,D10 8Slin-HI WD02,rLUSH, LOWER NW HOLC IN CYL s originot. dato l U

q 8- ,

ui ..

9  ::

( LEGEND

!! o she0 w S '.

~

i .9...sht ..wp 0 .

9:  : )

b; ::

l R. :i.

li J

o.

n

.t. .

t.

- s '

8 .:

.  :' : }! '
\ :l :  :

!5 ge :n

l

- $ ^.

J! :6.* .
[-

r 4 4

. .I  :: -

7 t- a

U.!  :

1 Si

~
:  ::-  :* ..

1

- t.

"9  :  : i ; .

-6 g.

., 1.

LJ .

1  :

8 a .

i. 3 i i h 5 2 I N-  :-

9 ': e 7

9 8-

ll (','. '3 U.

3 o

m g

9 m.

. . ..m. u

. ./

r- .

s.o 7.o I

e.o 9.o Jhio.o 1.o 4.0 s.o a 8~~ ,.o i.o 2.0 ti.=n isI =10" ,

, o.a so i co

- L,n t

Figure 2.2 Overlay of innfiltered pressure traces from SIIEOS and SIIEOS, lower cylinder wall position.

. ' i e~m. xo<3 X , y emw r e.1 9g r

e w

o l

6s 0g 00 6 on ,--0 1 0

E Dw o.

l NS6.

EO9t '

. . . . v1 i S

Ge .'.

Eh h. s.

.'t a 0 d L Ls .

9 n T

C

. a o ?.

N -

5 I

0 0 E

.e E l

i L

R O . S H

C O

NM I P a.:.:.:.'.

.o 0

7 m

o r

L f C

Y R E

W s

E O l

.:....!. e S L o....'. .s o.

c a

R H r E

WU O

S de i. I: * . : : t s t L L F, r e ......:...:~  :

o.

5 i e r

C 2 t.  :..

...;...l i

e. u T 0 B i r
t. .

t s I

W s L - 0 e .

U I ..e;.*.. 4 rn K H-  : ; : i : : .....:..**. po S

O A l

. :. dt i

WiS @

o 93 0 ei rs 5 8 eo

- 0 D t p EH N E,

% l il S

i A0

2 fl 0 a f w W

6 r_ o r

0 E m0

3 ye ad H

S . l n ri m

m3 -

o. el vy Oc o 8C o.g_ - o. $ 9gu 9 g* oN o g. o 87 .

a.85 o 8* 3

- 2.ca eo3e* eta 2

. e r

u g

i a ta .t .

.e I .3 - I p :e - n'e F

.-l t !  ! .l l!' , ,4 ' , .; -) 4 , ,;: 23I;, i  ! ,1

W. Berman November 20, 1985 1.329E+01)

EULERIAN POINT 5

( 2.736E+01,

, , , , , , ,')

2M 160.0 o -

o 150.0 -

~

x -

140.0 -

130.0 -

120.0 -

}

110.0 -

h 100.0 -

i -
90.00 -

60.00 -

% 70.00 -

e { -

,. g 60.00 -

m 50.00 -

m w -

[ 40.00 -

30.00 -

20.00 -

I ,

10.00 - )

0.000 ,

750.

45C. 550. 650.

50.0 150. 250. 350.

~ '

X10 O TIME (SEC)

LONE pitN 1AS FALL LOWER CYUNDER WALL Figure 2.4 CSQ II calculation of SHE test pressure trace at lower cylinder wall position

, ;* ,,5 ' ' m r

8n' s' -

~y i , l 00 a u

1 d 0

_ - ~ 2 x i v

i d

- n i

g

  • n i

' w

- ' o h

,R . ' s 4E -

0 ,

6O ii 5 )

t

- N i 1 n

KI u o

TL -

m Y .- ) w 1 C - '

S o M D l

(

- E -

N 7 RL O O 1 O .- X 1 H -

C E i00. E I ,

I u 5 . i S m 8E i 1

( o r

N -

f .

D -

E ea p

. NR M cr EE I

aw

,W r, td 5O j

.- T r 0L eo rc P ,  : ua WT -

0 sm si

. 7C i

'5 er rp OE i p

_ XR  :

l o m

. EI -

j l r HD  :

af w

- s re ek wi p

I o, s

- 0

. . 1

- 0 .

- 3 0 0 0 e 0 0 0 0 r 0 0 0 u 0 0 0 g

- 4 3 2 1 i

F yDnD1CQ nQ_

( ((G _

;i:l' i il1 l  ; ))l ;iIi!jll1i;>I  : } 1J l

l November 20, 1985 M. Berman 8

O

. x O. O  !.E I I O

N *-*'

3 I 3 I ) o gg)M

..g.,

..........i...............,.......v'

- l  :

U e

He..st"3 . . .

a

..~.........:......J . . . . ...

,. . i.

: . e WMN .  :  :

....[.........i..........y.........'......'..

Ca-g, . gQ . . ,

v s, e Z y 8o,o w'.s

. . . . . . .. e U g .

... - ' o o r-. g z z'z l  :

< u

...... . .. .... . .g ... . . . . .Q e

.. %7%,o uo

- a Q' O .

t n o

y U- . 3 '

. . . . . . ... . .. . . . . . . . . . . . .. . . . . . . . . . . . . . . . O...

v Q- gXS-)g gg i  :

i  : .j . . . . . . . . . . . -

. . . . . . .. . . .... .. . .... ...... .. ..... . _...... ' c g (f)y _3 x w,sas . .

-

  • in o C

. ! . . . . . . . i . . . . . . . . .:, . . [. . . ....

2 U - M g.ts' . 4. . ..

....(.. .. ..

O Z -

O K .,, EEE x x.1 O e o

.. . . . . .. OO 4

@ O *O* i e

.,. O

  • U e O Q =e -

- . . . . . . . . ... . .... .... . .... ... . ... .. ..p.

i er 1 .. . .

O y

--) - -

O_ g k  :  : .:  : .. gr' #e t'" ..... ....

te C .s .........s..........i........

..e..........,.....

y o J .,3 - ~ ~ - -

y a_ ', . ..,, .. . . . . . . . .. . . .. . .. . .. ... .. ... .... g a

......p.-. .

- e r

a- U .8 c

..-...-i n .. . :. . ...... ... .. . a. . . . . . .. E-. l h i-- ( ) p U

l--a

[2j  :

.. .... ..s .. .. ..,

t c

g C

o. c L

g

-.. .. . . i m .c o

N U  : .

. . . t

..y O. t--A .....

...g........g...

.t........g..................i.

4 o

X D *

  • U O

. . . .... c y . . . . ' . . . . . . . . . . ~ . . ......................................1 o

j . . . .

E

...i...........l,. .. .. .. ... L I

t

. - .. q p I tC l , .

i a
..... .. p.... E

. g ...O. c

. . e I

I t

i e i J o I i i i i i '*

O O O O O O O O O o o o c; O o o o o l tn a sn a sn a tn tn. <

n n ~ ~ -. -

L

=

IISd) 380SS38d  ?

1 I

l s

l

x HEX 10 HIGH EXPLOSIVE CRLIBRATION EXPERIMENT BASE PLRTE PRESSURE GRGES WP03, WPO4, RNO WP07  ?

INote: Doto filtered and boteline odjusted.I Q i 5 500- i i 0 UTER BASE GnGES wo3, kut. , T.S. , em CORNER BASE _

450 .

'f. 'o Hf'04 kUL.t 1.5. ~

SE CORNLR DASE

~ -

pl I~ Y /,~...l

% e.,~0lNJ ~SE'CD4tjE'R,QAS(

400 - - .';,l l l= ;

350 - -

l 6 m  ; -

O- 300 - - .

  • fs'  ; e L 250 - -

.s , i a .- .

m m 200 - -

l ,

m s t  ;

n_ 150 - -

  • ,! ,; l g

100 - -

{ 2

'; *l a r.,

g 50 - -

) '

0-

. x m M i e

i 10.0 15.0 20.0 ._.

0.0 5.0 -4 e

TLme (s) x10 =

Figure 3.3 Comparison of three symmetric gage positions (outer base plate) for ilEX10 lucite chamber test.

r COMPARISON OF NW CORNER BRSE GAGES, OUTER TEE S10E ,

i 0

FOR EXPERIMENTS HEX 10 (LUCITE 1 RNO HEXO4 (STEEL) e inot.: ooto ritt- .4 ona bos. tin. oa f.t.d.: p

s 1

i 2500 i i

o WP03, KULITE HKS-375 1100-4-4'

_q E04 _ M4 LI_E. _H5MF5. J $9tt.4 i

2000 - -

th -

Q- 1500 - -

] '..

A

  • ! .  ?

L  ;

0 .'

$ 1000 - - ,

' e * .

L  :

o- ,

: z 500 - - n l

i.

.. . . _ - ____/ _ A _ ,_

.a i .

g O 20.0 i

5.0 10.0 15.0 '

0.0 x10 S TLme (s1 i

i in !!EXO7 (steel cylinder) and i

Figure 3.4 Comparison of out,er base position gages j llEX10 (lucite chamber). ,

I

M. Berman November 20, 1985 o --,

L 4 4 N x

~ - -m -

g

~

5 -

y i Qe - C m Ee -

~

S ad -

-~

e

-yyx -..

3 -

EE Ln kl aN m li; k

ni Jos .cs o- p 8

m I'

-o~ 0)  ?

a M23 ze _

E J $

J et H 3

y2 LJ ' ~

.I N ~on S

o Jo s ..

E a.

3 Mo 1* ;c O .

O DE '" E a .a O E

ce L- -

w' c

c. -

s C

a:

G.

~

e e

k

_ o I,

o' O o O o o O o g a W 9

- r.

g 'nd j 6Jossegg i

PR000HP08,KUL HK5-375 1508-9-2, 2' ~

FLUSH, LOWER SE HOLE CYLINOER tuote: Boto rtttered wtth IUILT.r0RI j.

500 - , ,

400 - -

't -

1 3. 0. 0. -

(

S v }

.h I -

ip 200 - -

'D L

tt 0

[ 1 100 - -

3 ls

.O a a a a a a

a e a r a t

15.0 20.[1 5 0.0 5.0 10.0 ILme {S1 x10" figure 4.2 PR000 propellant test, wall-mounted gage, lower SE wall.

t ,

  • T2 " Ly !3? abm*

00 .

1 i 0 2 x l

, l a

w E

N 0 N 0 ,

. r

' , 5 e 2L w

- Y 1

o l

8 CR1

- 0 e

0 E r. g 9 LT a g

0 Ol L 1

Hr E t n

H u 4E h

]

o 7Nt m 3Nt w '

- s o n

0 [

i SRd e ,

. t KEr 0 a e

l HWe t 1

l o

s O l. '

m i L L,it U f L ,

K, Ho S ot f i

T t s

e 2UD t 0L: t PFe t '

n a

W/o l 0TH I l

e 0X '

p 0E 7 0 o r

R ,

. p P '

5 0 0

0

^ R P

~ 3 4

~^ i e

r u

g 0 F 0

0 0 0 0 0 0 0 0 0 0 0 5 1

- 3 2 1

-. - @ 9.I L Wmr,LLi C

                              ,t            r

3 x ,yg , N. E %&* 8~ E*

                                                                                  '00   .

1 0 2 x l l a w W S 0 S r i5 e 1 1 w 1 L o l

       - Y     1                                                             '

9 CR e

       -       0                                                                                      g 8E. r                                                                   '

a g 0LT L 5 OI t n 1 HFE ' u H ) o n 5W h r 7St ' s n 3St w 0 ( i o

       -                                                                                 .           t SRd                                                                   ,

i0 a KEr e i i 1 e l o HWe t m i s Ol ' L o , L U L,f i T t s K, Ho So e t t 4UD t OL: ' n PFe a t l W/o J l e 0Tt ' p 0X I 0 o r 0E . p R .

                                                                            ! 5 P                                                                                               0 0

0 R p 4 4 e

                                                                       ~                               r u

g i

                                                                        -              0              F
                           -          _            -             -                     0 0   0            g           0              0                0 0   0            9           0              0 5   1            3           2              1
                               ). )

mf ,_ $' m'L

                             -    fO_

( C

!' >}}