ML19212A533
| ML19212A533 | |
| Person / Time | |
|---|---|
| Issue date: | 06/02/1978 |
| From: | Levine S NRC OFFICE OF NUCLEAR REGULATORY RESEARCH (RES) |
| To: | Case E Office of Nuclear Reactor Regulation |
| References | |
| RIL-027, RIL-27, NUDOCS 7912140487 | |
| Download: ML19212A533 (23) | |
Text
{{#Wiki_filter:NE% W b UNITED ST ATES 4 [ t NUCLEAR REGULATORY COMMISSION $.J WASHINGTON, D. C. 20555 1' c k, / [d-
- ~~
sun 21978 l4 J ; r _T l ) L' p 'i i c MEMORANDUM FOR: Edson G. Case, Acting Director o Office of Nuclear Reactor Regulation FROM: Saul Levine, Director Office of Nuclear Regulatory Research
SUBJECT:
RESEARCH INFORMATION LETTER - #27 - BEACON / MOD 2
1.0 INTRODUCTION
1.1 Identification of the Completed Research This Research Information Letter transmits the BEACON / MOD 2 computer code manual, describes its field of application and discusses the code's strengths and limitations. BEACON / MOD 2 is an advanced, best estimate computer code intended for eval-uation of short term thermohydraulic conditions within " dry" (full pressure) multicompartment containments, or within certain regions of the " pressure suppression" containment drywell. Developed by Idaho National Engineering Laboratory, this code was completed, documented and transmitted to the Argonne Code Center at the end of December 1977. The code 'Jsers Manual is given in Enclosure 2. 1.2 Relation of this Research to Regulatory Needs This particular research was initiated by RES, about three years ago, as a result of its survey of the existing containment codes and identification of the problem areas in need of important modeling improvements for Best Estimate analyses. Hence, this research was undertaken in anticipation of future NRR needs, rather than in response to a specific NRR Research Request. The areas identified by RES to require better resolution, especially for deterTnination of the margin of safety, involved the evaluation of hydraulic loads on the walls, barriers, and system components located either within the break compartment or in the adjacent compartments of a full pressure (dry) containment. Some regions of the pressure suppression (PS) contain-ment drywell also require a detailed examination of hydraulic loads and could be analyzed in the same fashion. 782190306 1570 064 M'/
- i. ' ' '
_o Edson G. Case, Those hydraulic loads are caused either by the direct impact of the fluid discharging from the postulated break, or by the creation of significant differences in fluid pressure in the compartments divided by the internal walls or barriers. Such loads reach their peaks with-in a relatively short time following the initiation of LOCA. 2.0 DISCUSSION 2.1 Modeling Requirements for Best Estimate Analysis of Hydraulic Loads The requirements listed below pertain to the locations in the dry and PS containments described above and to the early period of LOCA. (a) The ability to treat thermal and mechanical nonequilibrium is re-quired. Thermal nonequilibrium constrains the rate of phase change, hence density change, and thereby greatly influences the magni-tude of the local impact load of a flashing jet. It also affects the rate of pressurization within the compartment receiving the discharging fluid. Thermal nonequilibrium implies that the liquid ar.d the gaseous phases (vapor and air) do not have to be at the same temperature. Mechanical nonequilibrium implies that the fluid components need not be intimately mixed and could be traveling with different velocities. (b) Multidimensional analysis capability is essential for description of the evolution of a flashing jet and its impact load; for tracking of the propagation of pressure waves within the confined regions containing the ruptured pipe (such as a reactor cavity); and for determination of the manner in which the discharging fluid mixes with air. For example, pockets of the slowly circulating air could be formed, greatly affecting the composition of the fluid flowing through apertures between adjacent subcompartments. Fluid com-position influences the rate of flow through such apertures and, therefore, the buildup of pressure differentials between subcompart-ments. (c) The ability to describe the containment with a computing mesh consisting of multidimensional, one-dimensional and zero-dimensional (or lumped parameter) zones is required. This allows the opportunity to place the emphasis on the computational details where they are needed. 1570 065
~ , 3D Edson G. Case. None of the codes reviewed prior to initiation of this research effort had the above listed capabilities. 2.2 Analytical Techniques Employed in BEACON / MOD 2 The fluid flow solution routines are based on the first version of K-FIX code developed at the Los Alamos Scientific Laboratory. That version of K-FIX solves the conservation equations for water, water vapor, and air in two dimensions (axisymmetrical). In BEACON / MOD 2 the computational mesh was extended to include the two-dimensional cartesian mesh, a one-dimensional mesh for the description of interconnecting flow passages, and zero-dimensional (lumped parameter) zones for description of those regions of the containment which do not require the resolution of flow momentum. In addition, BEACON / MOD 2 includes the capability to consider " obstacle cells," within the multidimensional zones, to model internal flow barriers or system components experiencing hydraulic loads. Sources of mass, momentum, and energy can be prescribed for any cell. This feature allows, among others, consideration of the primary coolant system break fl ow. By sharing the momentum flux over different cell faces, the jet of fluid discharging from the break could be aimed in any desired direction to explore the consequences of jet impact loads. Multiphase, multicomponent description of thennal hydraulics allows for relative motion between liquid (droplets) and the gaseous phase (vapor and air) and for unequal temperatures. Consequently, both the mechanical and thermal nonequilibrium are allowed. The numerical technique of K-FIX is based on a Fully Implicit Exchange of mass and energy between the liquid and the vapor iihases. 2.3 Evaluation of Strengths and Limitations BEACON / MOD 2 offers considerable advantages over the existing containment codes for best estimate evaluation of hydraulic loads in multicompartment PWR type containments. It is particularly suitable for evaluation of the reactor cavity loads (for postulated breaks between the reactor vessel and the biological shield) in both PWR and BWR containments. Sample re-sults of such an analysis are given in Enclosure 1. The code has also shown a capability to describe the evolution of a two-phase (flashing) jet and the resulting pressure loads on the impacted barrier. This capability was studied by modeling two jet impact experi-ments - one featuring a very small,10 mm diameter jet; the other featuring 1570 066
Edson G. Case , a 100 m diameter jet. The tests were performed by KWU and Battelle Frankfurt, respectively. As shown in Enclosure 1, good agreements with test data were obtained by suitably adjusting the mass exchange (r ) and the momentum exchange (K coefficients. Although the same va18es of r and K gave good aghee)ments with test data, it is clear that a pureTy predkctive cap-ability was not demonstrated but, rather, the capability to handle the situation when certain basic physical models are empirically normalized. The interphase heat transfer in BEACON / MOD 2 is the product of an inter-phase heat transfer coefficient, an input constant selected by the user, and the differential temperature between the two phases. There are several interphase mass transfer options available in BEACON which are formulated for various flow regimes such as dispersed, bubbly, slug, etc. However, the MOD 2 version does not yet contain the logic to detennine the type of flow and to choose the proper formulation. The user may select the desired flow regime modeling as input; however, once selected, the flow regime is fixed for the remainder of the problem. The momentum exchange between phases is the product of a momentum exchange coefficient and the velocity difference between the two phases. The modeling options currently available for computing the interphase momentum exchange coefficient are: (1) Input constant, selected by the user (2) Dispersed flow momentum exchange assuming a constant bubble / droplet size (3) Dispersed flow momentum exchange assuming a constant bubble / droplet number density. The above brief description of the techniques employed in BEACON / MOD 2 to model the interphase exchange of mass, momentum, and energy illustrates some limitations in the present version of the code. Such exchange pro-cesses must be described with more physical models and withoJt the User intervention. The ongoing work is aimed at removing these limitations. Additional limitaticns involve (a) inability to consider de-entrainment of liquid droplets due to wall impingement; and (b) inability to consider 1570 067
-S y/ Edson G. Case wall heat transfer. Both of these limitations are being removed during the calendar year of 1978. Liquid droplets that impinge on walls or obstacles form falling liquid films which not only change the moisture content in the flowing stream - and thereby affect pressure drops be-tween compartments - but also affect the wall heat transfer. 3.0 RESULTS The product of this research is a computer code, hereby transmitted. Results of some of its applications are found in Enclosure 1. 4.0 RECOMMENDATIONS It is recognized that the code must be more extensively tested against experimental data. Nevertheless, the BEACON / MOD 2 code is recommended for calculations of the reactor cavity loads, for both PWR and BWR in-stallations, and for evaluation of jet impact loads. This recommendation stems from the fact that this code, due to its advanced thermal hydraulics modeling, does present very significant improvements over the presently avallable analysis tools. Further testing of this code will proceed as new (presently planned) test data sources yield additional information. This code version is not recomended for evaluation of inter-compartment pressure differentials because it does not handle wall heat transfer. However, BEACON / MOD 2A, scheduled for release to NRC the end of December 1978, will meet these requirements. ( n aul Levine,wM a.ar/ Director Office of Nuclear Regulatory Research Enclosures : - BEACON / MOD 2 Developmental Verification - C. R. Broadus, S. W. James, W. H. Lee, J. F. Lime and R. A. Pate, " BEACON / MOD 2, A CDC 7600 Computer Program for Analyzing the Flow of Mixed Air, Steam, and Water in a Containment System - Users Manual," December 1977, CDAP TR 002 See cc next page 1570 068
i_/ p ENCLOSURE 1 d BEACON / MOD 2 Developmental Verification The MOD 2 version of BEACON is applicable only to examination of the PWR and BWR reactor cavity loads and of the jet impingement loads. There is no test data source available as of now to verify the code's ability to correctly calculate the reactor cavity loads. However, LASL's K-FIX code on which BEACON / MOD 2 is based has shown the capability to properly track shock waves propagating through two-phase fluids, to cal-culate break discharge flow, and the PWR downcomer transients. There wera only two test data sources available for testing the code's ability to calculate jet impingement loads. Kraftwerk Union (KWU) has published test data on a small scale, two-phase jet impinging on a flat plate. This test was simulated with the BEACON code for a range of K, the drag coefficient between the liquid and D vapor phase and r and r, the scaling factors relating the evaporation or e c condensation rate and the amount of superheating of the liquid phase or subcooling of the vapor phase. Because there was little flashing in this test, information could not be obtained about r rr. But stagnation e c and indicated a high pressure was found to be a definite function of the KD degree of coupling between the phases. Another test on a somewhat larger (100-mm diameter) jet impinging on a flat plate was made by Battelle-Frankfurt in their test Program RS-50, Experiment C-12. As with the KWU test, simulation runs were made and the results were compared with the measured pressure distribution on the plate. From these simulations, a reasonable value of the r was obtained and the e was verified. previously obtained value of KD ( To match both tests, a high value of interphase heat transfer coeffi-cient was used. An analysis was made to determine the flow and pressure distribution in the reactor cavity following a recirculation line break in a typical BWR containment. The values derived from the jet experiments for r and g e 5 were used in the model. Peak pressures of slightly under 4.7x10 Pa were found in the cavity. 1570 069
.i j ..N, s E E e,==4 g a 3 ~N ,/ s I N si m n 5 to N n n,' \\.' m n s -W tii =i E i W = a es: 'y
- ma i f
+ Qm m th w cy:
- 3 U3 CD baJ CE:
~~ 1570 070 I
ii //* VL = ,d' d 9 h wE IE! W Ch. ( ) W m d g e x C,, = LaJ g g "a.
- s:
e o
== G S M S -o w o o M II ll =J e a g 3 II.. x
- ae 5,
E E s o 5 = t w a \\ W 2 CD at wa::
- D a
C/3 C/3 w iWj i i f o g o o o o o o La att ev=3 N M (ed 01)3HilSS3Hd 3 1570 071
3 n x. sr
- 's
\\x Y-_ :., 8 i <o 5 g E = m 31 / 2 e. a C.3 = M C:2
- 5 y
!E W
- 2C:) E g g x
=> w m a. sx -x cc w \\ W X 22 m2 x I LaJ a x ,,,a =e V3 'N \\ $gu> + xm3
- s::
N N, \\ \\,\\ s mei t _.,;[ i \\ \\ / N ix\\(g ///; g ,3,'\\ / a. \\ s.- j s,- a \\s\\. \\ @1 lA' \\ \\ \\ '/. l"- ~ ? y i ;, \\ j!! a. 1570 072
l~] ' M ~- l t-3,. y 7 N N = i l \\ E
- i. j N
tt! \\ )! N 1 = $l N 5 \\ t _f i N...<4 ,J = \\ 1,i \\.<<;,;,ll1 = i= =g \\ \\ ,,,,rl1
== :s i. .a J a om - \\ -m i w mw 5 \\ ,l 1113 is ex .*1 cn Hg t -4, .<,4 \\ Eis = = $ N x .lll11 C W, ~ \\ l.
- 11j, sd N
WE N m \\ 111,.t!1 = i o== y. N
- a. i
,l N i ! i i \\ \\ \\ '. E *in + N ssss e \\ 1 i s i s. ~ N N g i is N N
- is
\\ m \\ N' \\ l_ s l 1570 073 $.1 5 I!
a y g. 6 e ma use
- =
N N N N 5 O. g >== E 2 / W N g W Em If k\\ I E W Ww a Q. >< aug z E 4
- p., e9 E
a oc o .e-a .a 2 'S o y s = z act: g a s.a - if aW l -J I -J W M cm t 1570 074 Y
9 r-{ /; [~ BATTELLE-FRANKFURT EXPERIMENT C-12 1 l STEADY STATE CONDITIONS ...r; 15 I ~~ r, = 0.001 Cj K =10 4 3 a 12! 10 9 5 m. = , = 0.015 y => N 15 [5: 5 = 0.M o e ,M ~ DATA POINTS ki' f ? ..s. I l 0 O 100 200 300 RADIAL. DISTANCE FROM JET CENTIRUNE
i4 )i o l BATTELLEfRANKFilRT EXPERIMENT C42 2-PHASE JET 15 CALCULATED RESULTS 4 ,1 K =10 g Fe = 0.015
- 10(
a
- a
.O .g 8 K = 2.2 x10 g ...2 wg a 5 DATA POINTS l STEADY STATE CONDITIONS i 0 0 100 200 300 RADIAL DISTANCE FROM JET d 076
1 %k i EE o2o o >== w e m w 2 M Sh C g;; 9= 5 m~ m J i== N w / C/3 s -O b Cll2 cn 2ow n E x o o oE G d W o cn= ll w g w 2 Cw o
- = cc:
o g =c cz: CW \\ N m" w= O = D== 2 =s o o CO .ll2 l >= 4 ' w-Q 'N s. _ _. _ --=. L y La. i a O w, N x y 'p _. p 4 /,... CCh - f. //. 'k ,/ y, ,/ .. s,. - ./ / .s / w '** p _- ~.. A L.. 1570 077 ~ M
p 1 e REACTOR CAVIH CONRGURATION i ,3 1
- +
.~ _i f s -///////////// / ~ ~~ 1570 078
- .i
\\ -m. t e \\ G = = 3NVid AH13WWAS \\ \\\\\\\\\\\\\\\\\\\\\\\\\\\\ ^ N \\ zo E N 5 E E I \\ E E ll g N m E E . \\ E d w aa s.s ~E N 8 E E E g 5 N E C; o E N iii ts N k / ts \\/ oc u ~ \\ I l \\ i 3NVid AH13WWAS h e i v> 1570 077 E a
,M,',-' i' ~ 's % y" 4 [ - i I. ~ ~ REACTOR CAVITY BREAK NODE PRESSURE 9 ~ i i i O 8 O Fe=0.015, K =2.2x10 D n. g g MASS Fl.0W RATE = 5790.f % ~ 3 F w e gg: E . m OS 3 wa:: En. 0~ 0 0.1 0.2 0.3 0.4 a TIME (s) is ~
- i. ' <
O CD O
5
- ,,
- - - x -- -
w ..s f _........ '. -- Y... - ' ' ' ~ ' - - ~ ~ ' ~ ' " ' ' ' " ^ REACTOR CAVITY PROBLEM PRESSURE VS. VERTICAL DISTANCE FROM FLOOR 5 i k 4 f, p r O.300 s m=3 l g__ ~ ~ ::. %, 55:2 i m i e. 0.120 s i N 9 m i l w% 1 = 0.004 s l+-BREAK LOCATION l o ~ o -0 0 4 8 12 16 DISTANCE (m) 7 /
,.7. f.,,, ' 1 y. -- ~ REACTOR CAVITY PROBLEM PRESSURE VS. HORIZONTAL DISTANCE FROM BREAK 5 3 e 4 q i t ~ ac )' =. 3 ?. O i -0.300 s w l E2 ~ ~ ' _ _ _ _ _ _. - 0.120 s M m. xs_~ ~ _ _ _ 9 a:: a-0.004 s I w No =0 0 5 10 15 ~ DISTANCE (m)
s <Og s, N C#3 N N N N \\ I M l'" C ll a s., s, s s N \\ \\ l CL-i= N N N \\ \\ \\ -Ug-a Lad en lll> s N \\ \\ \\ k gg:. a o, N glL. I act
- =.
N 's N 's \\ \\ \\ l IL =c .s -s 's s s i L 1 1 1 5. s o Ck:: '\\ \\ ] l } } o-s. s. oas baJ s <.- ~ ~ ~, s I / / / J C2:: ~~ 1570 083 e
i--.-i-p <s -- ( e. O e A L l I s s s s i g 3;; as ~ j g s s N \\ \\ \\ l 'll E. _ .s s s s s s s i i i g s;.: p
- s. i 9
.,s s N N \\ \\ \\ I l 8 m = o 3 ea O N N \\ \\ \\ \\ i ) ~ I s, s, s N \\ \\ l 1 1 1 -a i O N, N \\ \\ k l l l T' O ?' h. s N \\ / / / J .- i 4-t N 4 / L \\ k .,g ' e. 1570 084 m-
\\
- . j'
\\ = I T I 1 I JE o +. s s j t, s L. en i \\ le ~ s ' \\ N ca s \\ l /. g .m .a s s s i m !m [,
- D ca d
' '. - <a W 4. g i w r w
- E
- i.. {
A
- =
/ [ et i-o / } / CC I / \\ l C \\ / / \\ o i / p- \\ $' O / / ,/ l cm /..h t / //.A s' -l g g ,//.g,:p' g g!; /.,- r j j s I . E d 6 Il ll ch, m
- L i
i e 1570 085 1 1 i
9 e . I: s ,e L s 4. h.h H g t.' - I k p 2 c. O !. h m.= F. m l QC t 9-== .J (#.5 I ie g r ,i M p' vs lCN y o im O 'o g ,/ s u. N o r m L O f Il 11 T> e w O I seg f 9 / glgg 'C 'N p f %-s O 1 / l '/\\ d% N x l,. i i s,y.s ~ x 9 6 r O 1570 086 P P 0, "
1 9 e 1}}