ML20107C019
| ML20107C019 | |
| Person / Time | |
|---|---|
| Site: | 05000447 |
| Issue date: | 02/23/1982 |
| From: | Pfefferlen H GENERAL ELECTRIC CO. |
| To: | Boehnert P Advisory Committee on Reactor Safeguards |
| Shared Package | |
| ML20107C011 | List: |
| References | |
| FOIA-84-237 HCP-10-82, MFN-016-82, MFN-16-82, NUDOCS 8502210130 | |
| Download: ML20107C019 (10) | |
Text
G E N E R A L h E l.E C T R I C uuctema powEn SYSTEMS DIVISION Md68[$0B)9b-iN
" ^"
""'^
M h-2 February 23, 1902 U. S. Nuclear Regulatory Commission Advisory Committee on Reactor Safeguards Washington, DC 20555
' Attention:
Paul Baehnert Gentlemen:
SUBJECT:
RESPONSE TO REQUEST FOR INFORMATION At the January 22, 1982 meeting of the Fluid Dynamic Subcommittee of the ACRS, Dr. Theof anus asked for information related to full scale data on Mark III LOCA loads.
This information is addressed in a draft of Attach-ment D to Response 38.32 for GESSAR.
I have attached a copy for your information.
The attachment contains material which is the type General Electric maintains in confidence and withholds from public disclosure and is therefore identified as proprietary in the attached affidavit.
As such, we request that it be withheld from public disclosure in accordance with the provisions of 10CFR2.790.
Also enclosed for your information is a copy of a 1972 ASME paper " Liquid Surface Motion Induced By Acceleration and External Pressure" by F. J. Moody and W. C. Reynolds.
This work by Moody served as the basis for General Electric's current pool swell models.
Very truly yours, j0ll.,*1 (
H.C.Phefferlen, Manage 7 PWR Licensing Programs Nuclear Safety and Licensing Operation HCP: hjr/CO2024 Attachment cc:
L. S. Gifford (GE-Bethesda)
J. A. Kudrick (NRC), w/o Attachment M. B. Fields (NRC), w/o Attachment H. Faulkner (NRC), w/o Attachment R. Villa (GE), w/o Attachment 8502210130 840831 PDR FOIA c
SHOLLY84-237 PDR s
'8t GESSAR II A 000 238 NUCLEAR ISLAND Rev. 4 ATTACHMENT D TO RESPONSE 3B.32 - GENERAL ELECTRIC'S 2D MARKER-AND-CELL POOL SWELL ANALYTICAL MODEL
1.0 INTRODUCTION
General Electric's pool swell analytical model evaluates the suppression pool surface profile and velocity in a Mark III con-
- tainment during the pool swell transient following a postulated LOCA.
The code computes the two-dimensional potential flow field, representing the suppression pool and governed by the Laplace equation, with the Marker-And-Cell (MAC) techniques.
Favorable comparisons have been obtained between.model predictions and GE Mark III Pressure Suppression Test Facility (PSTF) test data.
2.0 MODEL DESCRIPTION the tes l swet \\ amatyi seat medet simw\\ates the swo-dimensicaal
(\\pW field af the Sylametry piame (plane AA in 7tgura D-1) in.t Mark Ill.eentainment unit cell,.which is defir.ei as the regio.
bounded by the inner and outer containment walls radially, and by two adjacent symnetry planes circumferential1y.
Thus, a unit cell.is.a 9-deg sec':or of the containment.
Listed below are the individual components of the pool swell model.
Some of these are not part of the pool swell code, and the details can be found elsewhere.
i 3BO.3.2.32D-1
,A._
t
I.
GESSAR II
'?A7000 238 NUCLEAR ISLAND Aev. 4 l
SIMULATION A%% :
PLANE p
SYMMETRY x
PLANES 1
CONTAINMENT OUTER WALL CELL lW W e,'
I f
e VENT CONTAINMENT
- j.
l':!
- l;3 INNER WALL 4:::
^1 I
/
-^E e
L >,
y WEIR WEIA ANNULU:1 WALL Ah 22A7000.REV4 Figure D-1.
Unit Cell Plan View
-3B0.3.2.32D-2
GESSAR II 22A7000 238 NUCLEAR ISLAND Rev. 4 r+
RESPONSE 3B.32 - GENERAL ELECTRIC'S 2D MARKER-AND-CELL POOL SWELL ANALYTICAL MODEL (Continued) 2.1 COMPONENTS OF POOL SWELL MODEL The pool swell analytical model is composed of the following components:
e dryaell response model e
vent clearing model e
vent clearing jet model e
vent flow model e
suppression pool hydrodynamic model e
bubble model The drywell pressure response and vent clearing time, velocity and acceleration are obtained from the containment analytical model (Reference D-1), and are not part of the pool swell code.
The drywell response following vent clearing is used as input by pool swell model as the driving force in the vent flow and bubble charging model.
The vent clearing jet model employs the transient, accelerating water jet model described in Referen.ce D-2 to obtain the jet penetration, ip.
The assumption is made that t.he initial shape of the bubble is cylindrical'with a diameter d and length Ap, y
as shown in Figure D-2.
The nose is rounded slightly to avoid possible computational difficulties in the pool hydrodynamic model.
The following sections describe the analytical models in the code.
3BO.3.2.32D-3 y
GESSAR II 22A7000 238 NUCLEAR ISLAND p,y, 4
RESPONSE 3B.32 - GENERAL ELECTRIC'S 2D MARKER-AND-CELL POOL SWELL ANALYTICAL MODEL (Continued)
/H/f// /
)
O*
'"""f=
=
Figtire D-2.
Initial Loca Bubble Profile 2.2 VENT FLOW MODEL The assumptions in the vent flow model are 1.
Air flow only - this removes any condensable steam from the bubble and results in conservative bubble' pressure response.
2.
Fanno flow - That is, the flow is one-dimensional, compressible, adiabatic, frictional and follows the perfect gas law.
Vent flow is needed to provide bubble charging and evaluate-bub-ble pressure.
2.3 SUPPRESSION POOL HYDRODYNAMIC MODEL It is assumed that the flow is two-dimensior.al, incompressible, inviscid and irrotational.
Thus, the governing equation is the Laplace e;uation for the velocity potent.ial, I:
27h=0 3BO.3.2.32D-4
}
f GESSAR II 22A70go 238 NUCLEAR ISLAND g,y, 4
RESPONSE *3B.32 - GENEPAL ELECTRIC'S 2D MARKER-AND-CELL POOL SWELL ANALYTICAL MODEL (Continued)
The boundary conditions are At rigid walls:
f=0 h norma], to walls At free surfaces:
r (n)2" n + 1 (4)2 Si+e o
+
+
r P
=
constant Using finite differencas (Figure D-3), the Laplace equation is transformed-into a set of simultaneous algabraic equations.
Given the initial values of velocity potential on the free sur-faces, these simultaneous algebraic equations are then solved Eby overrelaxation for the velocity potentials in the interior region, which in turn permits velocity components and potential rates to be calculated at the free surfaces.
Markers used to describe the position and shape of each free surface are then moved with their given velocities, providing the boundary con-ditions for flow field computation at-the new time step.
The process is repeated until finished.
2.4 BUBBLE MODEL The' assumption of a perfect noncondensable gas in the bubble allows _the determination of bubble pressure through the relationship.
MB-B P.
B VB 3BO.3.2.32D-5
[
.n..
- - ~ - -
~
r.'.
GESSAR II 22A7000 238 NUCLEAR ISLAND Rev. 4 POOL SURF ACE
'C o
a
,, N I,,
o MARK.
ERS 4s
.)
,s es es V
%f
,f Af
%W h,
6
( 'r l3 d)
E3 l3 l) l)
I) f s
,) [ CELLS BUBBLE SURFACE f
'e r
),,
, s,
.)
s (H
l3 l3
.' )
', )
l}
T T
o e
e c
- )
c o
Z-ny
+
,w s.
er J
w w
w w
w i
' Y 0
O - REGULAR MESH POINTS O -IRREGULAR ME5H POINTS i
2247000, REV 4 Figure D-3.
Marker-and-Cell Representation of Flow Field 3BO.3.2.32D-6
I
GESSAR II 22A7000 238 NUCLEAR ISLAND Rev. 4 RESPONSE 3B.32 - GENERAL ELECTRIC'S 2D MP.AKER-AND-CELL POOL 4
SWELL ANALYTICAL MODFL (Continued) where
't P
=
bubble pressure B
t bubble air mass M
=
m dt
=
B vent tyc drywell temperature (assumed to remain constant)
T
=
B R =
gas constant for air V
=
bubble volume B
To define V the bubble growth is divided into two phases.
In B,
the early_ phase the bubble is relatively small compared to the cell size, thus the bubble growth is predominantly three-dimensional and the cell walls-(i.e., symmetry surfaces dividing adjacent bubbles) do not have a strong influence on the bubble expansion.
The bubble volume during this phase is given by the expression (Refer to Figure D-4).
V B
- B
- A B in which the two-dimensional area, A is given by Green's B,
theorem i
C an.f. K ~ iS an CO.*f
- r 10.'. '. 007.503".I
.00 If dffOr*IEi 10EI g
3BO.3.2.32D-7
,.,.mn.,
,-n-,
n--=
~~
. o'
- GESSAR II 22A7000 2 38 NUCLEAR ISLAND Rev. 4 4
RESPONSE 3B.32 - GENERAL ELECTRIC'S 2D MARKER-AND-CELL POOL SWELL ANALYT1,AL MODEL (Continued) z h
i m
4L an BUBBLE AREA Ag i
i U
?
Y x
Figure D-4.
Two-Dimensional Bubble Profile For the special case where the two-dimensional bubble profile is an ellipse with semi-axes a and b, n ab A
=
B and V reduces to B
V
=
B F.B n ab 3BO.3.2.32D-8
m- -
- f.**',
GESSAR II 22A7000
.n.-
238 NUCLEAR ISLAND Rev. 4 RESPONSE 3B.32 - GENERAL ELECTRIC'S 2D MARKER-At:D-CELL POOL i
SWELL ANALYTICAL MODEL (Continued) which corresponds to the volume of an ellipse rotated about the y axis if K is set equal to 1.0.
B The bubble growth enters the second phase when the lateral expan-sion is constrained by the side walls (that point being when the lateral dimension equals G) and is forced to expand upward.
The bubble volume is then calculated as l
t 4
U l
A B
7 B7 B For.the special case where the two-dimensional bubble profile is an ellipse, the bubble volume reduces to YB
- B "
- A B
qal k which in the volume of an ellipsoid with semi-axesp a, b, a.nd w/2 if K is set equal to 1.0.
B For the bubble constant KB, a value of'1.25 is recommended based on data comparison and physical argument.
3.0 MODEL-DATA COMPARISCNS The pool swell analytical model predictions he.ve been compared to PSTF test data.
A total of 11 test cases were used in the model-data comparisons.
Only air blowdown-tests were used in the com-parisons.
The parameters compaired are the-pool surface elevation history, the surface velocity as a function o:! elevation, the pool surface and bubble profiles and the wate: ligament thickness.
1 3BO.3.2.32D-9
UNION OF CONCERNED SCIENTISTS W6 Conrintimt Aseriue N.W.. S.1101. Washington, DC 20036. (202) 296 5600 30 March 1984 FREEDOM OF INFORMATION Mr.
J. M. Felton, Director ACT REQUEST Division of Rules and Records
)
Office of Administration h sg g 3 fg U.S.
Nuclear Regulatory Commission Washington, D.C.
20555
/
RE:
Freedom of Information Act Request for Documents Concerning NRC FIN A-3320, " Hark III LOCA Pool Dynamics"
[Sholly FOIA 84-19]
Dear Mr. Felton:
Pursuant to the Freedom of Information ' Act, please make available at the Commission's Washington, D.C.,
Public Document Room copies of documents in the following categories:
1.
All documents submitted to the NRC under a research program-at Brookhaven National Laboratories, NRC FIN A-3320, " Hark III LOCA Pool Dynamics".
This request
- p specifically includes the Form 189 application, all periodic progrets
- reports, letter
- reports, draft reports, final reports, correspondence, memoranda, and other documents pertaining to the work done in this program.
2.
All correspondence among NRC, utilities (and their consultants and contractors),
and Brookhaven National Laboratory in which any draft reports prepared under NRC FIN A-3320 were transmitted to any one outside the NRC and any comments upon such draft reports.
Should you or your sta f f have any questions regarding this request, please contact me at (202) 296-5600.
Your cooperation in responding to this request is appreciated.
Sincerely, hn 0
N-Steven C.
Sholly Technical Research Asso jete i
- t
_1 _w y p i [ ( pf. 6 O Main Office: 26 Church Street. Cambridge, Maseethusetts 02238. (617) ?17 5552