ML19210C813

From kanterella
Jump to navigation Jump to search
Acceptance Criteria for Mark I Containment Long-Term Program
ML19210C813
Person / Time
Site: Hope Creek  PSEG icon.png
Issue date: 11/06/1979
From:
NRC COMMISSION (OCM)
To:
Shared Package
ML19210C809 List:
References
NUDOCS 7911200143
Download: ML19210C813 (46)


Text

{{#Wiki_filter:. NRC ACCEPTAtlCE CRITERIA FOR THE MARK I CONTAlfiMEtlT L0f4G TERf1 PROGRAft 4 1365 005 7911200 / e

TABLE OF CONTENTS 1. INTRODUCTION 2. SUPPRESSIM 000L HYDRODYNAMIC LOADS 2.1 ContairA :: 1'ressure and Temperature 2.2 Vent System Pressurization and Thrust Loads 2.3 Net Torus Vertical Pressure Loads 2.4 Torus Pool Swell Shell Pressures 2.5 Compressible Flow Effects in Scaled Pool Swell Tests 2.5 Vent System Impact and Drag Loads 2.7 Pool Swell Impact and Drag on Other Internal Structures 2.8 Froth Impingement and Fallback Loads 2.9 Pool Fallback Loads 2.10 Vent Header Deflector Loads 2.11 Condensation Oscillation Loads 2.11.1 Torus Shell Loads 2.11.2 Downcomer Loads 2.11.3 Vent System Pressures 2.12 Chugging Loads 2.12.1 Torus Shell Loads 2.12.2 Downconer Loads 2.12.3 Vent System Pressures 2.13 Safety-Relief Valve Discharge loads 2.13.1 SRV Discharge Device 2.13.2 Discharge Line Clearing Transient 2.13.3 Air-clearing Quencher Discharge Shell Pressures 2.13.4 SRV Discharge Line Reflood Transient 2.13.5 SRV Air and Water Clearing Thrust Loads 3.13.6 SRV Discharge Line Temperature Transient 2.13.7 SRV Discharge Event Cases 2.13.8 Suppression Pool Temperature Limits 2.12 Submerged Structure Drag Loads 2.14.1 LOCA Water Jet Loads 2.14.2 LOCA Bubble Drag Loads 2.14.3 Quencher Water Jet Loads 1 % 5 006 i

2.14.4 Quencher Bubble Drag Loads 2.14.5 LOCA Condensation Oscillation Drag Loads 2.14.6 LOCA Chugging Drag Leads 2.15 Secondary Loads 2.16 Differential Pressure Control Requirements 3. STR'JCTURAL ANALYSES AND ACCEPTANCE CRITERIA 4. REFERENCES 13bb ii

_j. 1. INTRODUCTION The purpose of the Mark I Containment Long Tem Program is to perform a complete reassessment of the suppression chamber (torus) design to include suppression pool hydrodynamic loads which were neglected in the original design, and to restore the original intended design sat'ety margins of the structure. This reassessment will be accomplished by a Plant-tjnique Analysis (PUA) for each BWR plant with a Mark I containment, using load specifications and structural acceptance criteria that are appropriate for the life of the plant. The following acceptance criteria have been developed from the staff's review of the Long Tem Program Load Definition Report (LDR), the Plant Unique Analysis Applications Guide (PUAAG), and the supporting analytical and experimental programs conducted by the Mark I Owners Group. These criteria specifically address the dynamic loading conditions. Unless otherwise specified, all other loading conditions and structural analysis techniques (e.g., dead loads and seismic loads) will be in accordance with the plant's approved Final Safety Analysis Report (FSAR). Similarly, references to original design criterb or original loading conditions shall be defined as those criteria or loading conditions which were found acceptable by the staff during the operating 1; cense review of the FSAR. For ease of reference, LDR refers to " Mark I Containment Program Load Definition Report," NED0-21888, PUAAG refers to " Mark I Containment Program Structural Acceptance Criteria Plant Unique Analysis Applications Guide," NED0-24583, and other supporting topical recorts are referred to by their report numbers. A complete set of the references used in these criteria, listed in numerical order, is presented in Section 4. 2. 513PRESSION POOL HYDRODYNAMIC LOADS 2.1 CONTAINMENT PRESSURE AND TEMPERATURE RESPONSE The pressure and temperature transients for the drywell and wetwell shall be detemined by the use of the analytical models and assumptions set forth in Section 4.1 of the LDR. These techniques have, in the past, been found to provide conservative estimates of the containment response to a LOCA, by comoarison to the staff's CONTEMPT-LT computer code. Plant-specific results for each break size shall be presented in the PUA, along with the input conditions, in sufficient detail to allow the staff to perfom confimatory analyses to assure proper application of these models. 1365 008

. The timing and duration of specific loads are based primarily on the plant-specific containment response analysis for the pool swell-relt.ted loads, while the condensation periods are non-mechanistically maximized. However, tht duration of the generic SBA condensation loads are assumed to be limited by manual operation of the Automatic Depressurization System (ADS) at 10 minutes into the accident. Therefore, as part of the PUA, each licensee shall identify procedures (including the primary system parameters monitored) by which the operator will identify the SBA, to assure manual operation of the ADS within the specifiad time period. Longer time periods may be assumed for the SBA in any specific PVA, provided(l) the chugging load duration is correspondingly increased, (2) the procedures to assure manual operation within the spec 9ied time period &re identified, and (3) the potential for thermal stratification and asymetry effects are addressed in the PUA. 2.2 VENT SYSTEM PRESSURIZATION AND THRUST LOADS The vent system pressurization and thrust loads shall be defined in accordance with the procedures set forth in Section 4.2 of the LDR, with the following exception. In order to assure the proper transition between vent clearing and bubble breakthrough for those plants that propose operation with a differential pressure control, the vent clearing time shall be derived from a containment analysis assuming no drywell/wetwell differential pressure and this time shall be applied to the vent system transients calculated from a containment response with the prc)osed drywell/wetwell differential pressure. In addition, for clarification, in the equation for F2V in Section 4.2.lc of the LDR, P3 shall be replaced by P2. 2.3 NET TORUS VERTICAL PRESSURE LOADS The downward and upward net vertical pressure loads on the torus shall be derived from the series of plant-specific QSTF (Quarter Scale Test Facility) tests, in accordance with Section 4.3.1 of the LDR. However, based on our review of the pool swell tests conducted by the Mark I Owners Group and confinnatory tests performed for the NRC by the Lawrence Livermore Laboratory, we will require that the following mar-gins be applied to each loading phase: UP = UP + 0.215 (UPmean) mean 2 DOWN = DOWN + 2 x 10-5 (DOWN mean mean where "mean" refers to the average of the QSTF plant-soecific test results (1bf). These margins shall be applied to the QSTF "mean" load function prior to scaling the load function to full-scale equivalent conditions. The margin for the downward loading function 1365 009

. shall be derived in terms of a fraction of the load at the time of the peak downward load, and that fraction shall be applied to the entire downward loading phase. The margins spe 'fied above may be reduced or omitted where minimum conservatisms ( .e., smallest parameter deviation from the nominal plant condition over the range of tested conditions) in the QSTF tested conditions for a specific olant can be demonstrated by the application of the QSTF sensitivity test series (NEDE-23545-P). The sensitivity tests may not be used to adjust the mean torus vertical pressure loads. If the plant configuration is changed to the extent that a single QSTF test series no longer represents a conservative configuration of the plant, then a new series of QSTF tests shall be performed. Application of the sensitivity tests and interpolation between plant-specific test series will be reviewed on a case-by-case basis. For those plants that use drywell/wetwell differential pressure control as a load mitigation feature, an additional structural analysis shall be performed assuming a loss of the differential pressure control to demonstrate the capability of the containment to withstand this extreme condition, as specified in Sections 5.3, 5.4, and 5.6 of the PUAAG. For this analysis, a single plant-specific QSTF test run may be used to define the loading function; however, the downward and upward loading phases shall be increased by the margins specified above for the base analysis. 2.4 TORUS POOL SWELL SHELL PRESSURES The spatial distribution of the torus shell pressures during pool swell shall be defined from the plant-specific QSTF test results and the azimuthal and longitudinal distribution factors defined in Section 4.3.2 of the LDR. However, the QSTF results shall be adjt.sted to incorporate the rargins specified for the nat torus vertical pressure loading function as follows: 1. During the downward loading phase, the average pool pressure shall be increased by the equivalent differential pressure, as a function of time, corresponding to the margin for the downward load. 2. During the upward loading phase, the torus airspace pressure shall be increased by the equivalent differential pressure, as a function of time, corresponding to the margin for the uoward load. 3. The pressure distributions shall be maintained such that the integral of the torus shell pressures will equal the net vertical pressure function with the margins included. 1363 010

4 2.5 COMPRESSIBLE FLOW EFFECTS IN SCALED POOL SWELL TESTS The QSTF plant-unique test series are based primarily an a " split-orifice" vent flow scaling relationship. Prel.minary calculations performed by EPRI and GE indicate that compressibility effects, which could not be accurately scaled in the testing program, could result in a higher loading condition at full-scale conditions than that derived from " scaled-up" test data. The original intent of these analyses was to provide justification for the scaled flow distribution in the EPRI 1/12-scale, three-dimensional puol swell test program. The loading functions predominantly affected by this finding are the torus downward and upward vertical pressure loeds, the torus pool swell pressure distribution, the vent header pool swell impact timing, and the vent header deflector impact timing. Based on our review of the preliminary analyses perfomed by EPRI and GE, which were presented to the staff in a meeting on July 24, 1979, we concluded that there is suf11cient margin in the loading functions to justify proceeding with ienlementation of the Mark I LTP, while this assessment continues. We will require, however, that the Mark I Owners Group complete the assessment of compressible flow effects and justify the adequacy of these load specifications prior to the issuance of the staff's Safety Evaluation Report, which is currently scheduled for December 1979. In the event that the adequacy of the load specifications cannot be demonstrated, these loading conditions will have to be reassessed. The vent header and vent header deflector impact timing do not, however, appear to be sufficiently conservative. Based on our review of the material presented thusfar concerning the EPRI 1/12-scale three-dimensional pool swell tests and the compressible flow effects analyses, we conclude that the downcomer orificing used for the " split-orifice" tests do not provide a prototypical pool swell response. Therefore, the vent header and vent header deflector impact timing shall be derived from the " main vent orifice" tests (using the same longitudinal load distribution methodology in the LDR), until a flow distribution analysis, acceptable to the staff, can justify somc less severe loading conditions. 2.6 VENT SYSTEM IMPACT AND DRAG LOADS 2.6.1 Vant Header Impact and Drag Loads The load definition procedures set forth in Section 4.3.3 of the LDR are acceptable, subject to the following clarifications: \\ bY

. 1. The experimental data of local vent header pressure in each of the Mark I plants shall be obtained from the QSTF plant-unique tests. 2. The specification, for each Mark I plant, of the pressure inside the vent header relative to that in the torus airspace at the time of water impact on the vent header shall,be determined from the QSTF plant-unique tests. 3. The plant-unique header impact timing (i.e., longitudinal time delay) shall be based on the EPRI " main vent orifice" tests as described in Section 2.5. 2.6.2 Downcomer Imoact and Draq Loads The load definition procedurc3 set forth in Section 4.3.3 of the LDR are acceptable, subject to the following clarifications. A pressure of 8 psid is to be applied cnifomly over the bottom 50* of the angled portion of the downcomer, starting from the time at which the rising pool reaches the lower end of the angled section and ending at the time of maximum pool swell height. The pressure is to be applied perpendicular to the local cowncomer surface. The structural analysis for the downcommer impact shall either be dynamic, accounting for the approximate virtual mass of water near the submerged parts of the downcomer, or a dynamic load factor of two shall be applied. 2.6_3 Main Vent Impact and Drag Loads The impact and drag loads on the main vent shall be evaluated in the following manner: 1. Subdivide the submerged portion of the main vent pipe into six equally wide segments (see Figure 2.6-1). If this sub-division results in AL < 0.3D fewer segments may be used such that al s 0.30. 2. Determine the velocity and acceleration histories at Points I tarough 7 in Figure 2.6-1 from the QSTF data and aopropriate corrections for longitudinal variations along the torus (at Point 7, only the initial impact velocity is required). 3. Using the velocity components nomal to the vent pipe, calculate the impact and " steady" drag pressure using the method in Section 2.7.1 (Cylindrical Structures). At Point 7, only impact force is to be considered. 1363 012

D 4 p 4 3 y Highest pool

levation

~ V ces a % pool elevation at V a cos a Itia in vent a V = pool impact velocity at station a = pool acceleration at station Figure 2.6-1 Schematic Diagram Illustrating the Methodology for Main Vent Impact and Drag 1363 013

. Using the acceleration components normal to the vent pipe, calculate the. acceleration drag pressure using the equation Pa" gFstatic buoyancy Where Pa is the acceleration pressure averaged over the pro-3 jected area (psi), o is the density of water (lbm/ft ), C is is the acceleration dua to gravity (ft/sec2).(ft/se:2) the diameter in feet, the cross flow acceleration and gc 4. Sum the pressures due to impact, viscous drag and acceleration drag and multiply by D to obtain force per unit length at Stations 1 through 7. 5. To obtain a smooth loading history for the main vent as a whole, the linear intercolation method suggester "or the vent header deflector in Section 3.5 of NED0-24612 may be used. 2.7 POOL S_WELL IMPACT AND DRAG ON OTHER_ INTERNAL STRUCTURES The impact and drag loads for internal structures above the suppres-sion pool (except the vent header, downcomers, and vant header deflectors), as specified in Section 4.3.4 of the LDR, shall be modified such that the structures are classified as either cylindri-cal (e.g. pipes), exposed flat surfaces (e.g., "I" beams), or grant-ings. The following load specifications for each of the three structural classifications shall be used to replace the methodology in the LDR. Any structures that cannot reasonably be classified as one of these geometries will be reviewed on a case-by-case basis. The longitudinal velocity distribution shall be based on the " main vent" EPRI pool swell tests, as discussed in Section 2.5. 2.7.1 Cylindrical Structures For cylindrical structures, the pressure transient which occurs upon water impact and subsequent drag is depicted in Figure 2.7-1. The parameters in Figure 2.7-1 shall be defined as folicws: 1. The maximum pressure of impact P will be detemined by max 2 P = 7.0 x 1 o V max 7 344 9e wht.re Pmax is the maximum pressure averaged over the projected 3 areu ind), o is the density of water (lbm/ft ), V is the impact velocity)(f'

c) and g is the acceleration due to gravity c

(ft/secz. 1363 014

8-Pmax I U t E I P p D

  1. y-

/ / / t i TIME Figure 2.7-1 Pulse Shape for Water Impact on Cylindrical Targets 5 red < x 10 1.5 / __. _.. __ ; i. ;_

i. j. i ! j j v1

.{ ,j..i l-

. 1

.1-~.7 f_. ! Mi.. [ l. C 5' - Re 1 5 x 10 D 0 'N N, r 8.. ' 1.0 ~ ' 7... j a.. ..l _. a. i ! ! + I l l l l l l l i l l I l 1 5 ! I l._4._ _ - f i i i 0 2. u. 6.

s. 10 20.

uo. 60. eo. FROUIE NUMEER (N = V /Dg) F 1365 015 Figure 2.7-2 Drag Coefficient for Cylinders Following Impact

. 2. The hydrodynamic Mss per unit area for impact loading shall be obtained frorr.ne correlation (cylindrical target) depicted by Figure 6-8 in NEDE-13426-P. A margin of 35% will be added to this value to account for data scatter. 3. The impulse of impact per unit area shall be detennihed by: .b V h I ^A~ Wg) p c where I is the impulse per unit area ( hydrodynamic mass per unit area (.lbm/ftgsi-secl, M /A is the p H } and V is the impact velocity (ft/sec). 4. The pulse duration will be determined from the following equation: T = 2I /P p 5. The pressure due to drag following impact shall be determined by: C o D cy-p = D I 144 g c where P is the average drag pressure acting on the projected D area of the target (psi), Cp, is the drag coefficient a3) defined by Figure 2.7-2, and o is tMe density of water (lbm/ft 2.7.2 Flat - Surface Structures For flat-surface structures, the pressure transient which occurs upon water impact and subsequent drag is depicted in Figure 2.7-3. The parameters in Figure '.7-3 shall be defined as follows: 1. The pulse duration (t) is specified as a function of the impact velocity: t = 0.0016W for V 1 7 ft/sec r = 0.011 W for V > 7 ft/sec V where W is the width of the flat structure (feet) and V is the impact velocity (ft/sec). 1365 016 2. The pressure due to drag following impact shall be determined by: PD* O 2 144 9j \\ e where PD is the average dra area of the structure (psi)g pressure acting on the frental ,CD is the drag coefficient (C0 = 2, flat strips normal to flow, independent of Reynolds number), and o is the density of water (lbm/ft ), 3 3. The hydrodynamic mass per unit area for impact loading shall be obtained from the correlation (flat targets) in Figure 6-8 in NEDE-13426-P. A margin of 35% shall be added to this value to account for data scatter. 4. The impulse of impz::t per unit area shall be detennined by: b V I = p A 144 g. where Ip is the impulse per unit area (psi-sec), Mg/A is the hydrodynamic mass per unit area (lbm/ft2), and V is the impact velocity (ft/sec). 5. The maximum pressure (Pmax) shall be calculated from the impulse per unit area and the drag pressure as follows: P max D 2.7.3 Gratings The static drag load on gratings in the pool swell zone of the wetwell shall be calculated for gratings with open areas greater than or equal to 60% by forming the product of the pressure dif-ferential (figure 2.7-4) and the total grating area (not only the area of the metal bars). The pressure differential curve in Figure 2.7-4 is based on a velocity of 40 ft/sec. If the maximum pool velocity in the area where gratings are located differs from 40 ft/sec, the force on the grating will be calculated as follows: F = AP x A Y grating max 40 To account for the dynamic nature of the initial loading, the load shall be increased by a multiplier given by: 'I36'J 017 F /D = 1 + [ l + (0.0064 Wf)2 )<~2; for SE Wf < 2000 in/sec.

11 - 1 y mm a 3 P D 2 .Ji m max J j .cW z Q. >== a< c.r we<z I w$p P / D / / / / \\ l TIME r Figure 2.7-3 Pulse Shape for Water Impact on Flat Tarcets -_..__y__ i ... 1 _ _ _t.__ x_ 2 -N.-_-._ _3 ...; __......4..... ._.i.... -..r. ._c 21 i ...... _ _ + _.-~~ s _~.. v = 40 ft/sec --. .g .p.. _. 4._ u _ __ a._ _._.i. ._.7__...... .9.. ...1. 4 ..____.2.. t 5 .... __,5-x ..__r._. x. Z s + _-. ' _.. W. = x. + m m x-- .__..r_______.___ ...i... i_.__ , s a 1o g._..._m. _......... _ ___.4._.___ 7 r...._._..... _... 4...... _._. m p_._ _ _ _ _._. _.. I _ ._a m. .p. .-.-- -.3.__ _.._._ -. _. 4 g 5 .........._.-...[ 4__.. ...__._t.. ..__4... ... _.. 4 ---_y_4--...

4..

..-..._..._J_.. o ~ -' " '4 -~ ' ~~t ... [. ~ ~ t'~ ~~~ I'~~'~~~ ~ f~~~~ ' ' ~ ~ ~ - ~ ~ ~f 0.5 0.6 0.7 0.8 o.g 1.o \\56,3 >jQ OPEN AREA FRACTION Figure 2.7-4 Pressure Drop Due to Flow Across Gratings

- 12 where: F static equivalent load SE = width of grating bars, in, W = f = natural frequency of lowest mode, Hz O = static drag load If Wf > E000 ft/sec (not expected for gratings) the force on the bars of the gratings will be calculated by the method outlined above for flat-s*x iaced structures. 2.7.4 Load Apolication These load specifications correspond to impact on " rigid" structures. When performing the structural dynamic analysis, the " rigid body" impact loads shall be applied; however, the mass of the impacted structure shall be adjusted by adding the hydrodynamic mass of impact, except for the gratings. The value of the hydrodynamic mass shall be obtained from the appropriate correlation in Figure 6-8 in NEDE-13426-P. When the impact loading is primarily impulsive and calculations have already been performed in accordance with the LDR methodology, simple adjustments may be made to the LDR analyses. Under these conditions, a parabalic pulse shape, as proposed in LDR, is acceptable provided corrections are made to account for the 35% margin in the impluse and with additional corrections for the drag force imediately following impact. Fcr structures with a natural frequency less than 30 Hz, loading can be treated impulsively (i.e., independent of pulse shape) when the conditions fall into the region above the straight line shown in Figure 2.7-5. The following corrections must be applied to the previously calculated stresses: 1. The celculated stresses will first be multiplied by a factor of 1.35 to account for the data scatter in the impulse data. 2. The calculated stresses will then be multiplied by ai additional factor to account for the presence of drag following the impact. This factor is determined as follows: 1363 019

. TYPICAL IMPACT STRUCI'UEES: -a% a r% f V[ \\ / \\ / \\ i i i i I l ) k /[ i / \\ / N / s / N/ wmdL [ 4-- W --> 4-W -+ .L w :au~--,l-1.-- - i.: _ - -. y. n 7-- =- r. : ::+ - - f,= =,:nj:.a: 1 -= +- - ~ ~ - - - - __ = .s-..- - +- - - 1


._--.--..~_._._~.b-_.'___1_..,

I .. - _'; l -i -+-- ..__:)_.---+ ..u ---4 / .=.=.-----.----------f_ - L _- - _. L: :...: t _ : _ 1 -

L :. ;. :j :.

l----- [ + - 7:

g... : 4. _... - t - - - = r.:._ z= =u :_;--- : - - - -. - -

-+--- -.---.__s.-.__-. }..._:)=_:-------*--_--*------- m

t... i.-
.. -- +

+ --+- - ---. - : j. : ry [d - "~ 5. IMPULSIVE LOAD ~ -~ g 10 L__: '..- \\ I. t o ...___=.f..=;&'_--_._ _.. _. _ _. _.... _. REGION -- 3R-5. >=_.._=. ;. h r __ s. -'- = } -. M_--g_..._ _ __. _ ,j V .~. ..I,,.. ..4_. .4___ .1.' ' :b: *: _":c--- = l _=p_A - -- -- C

:

= _r__.._u-._.__._ a ~ ~ ~#_. :- T= "~ ~ ~ *' " - '===- r- -== = t - __ t.a rn...:. -- *~'~] =' "~-- "E F T / 5 2} h s E.._~ $ -- =u =- -- =5- .._.r .-E -.. - u_ ---4 g _ _ :.=_ U _r-iEjif " - =:n- :=:l:-

-. L - - ---..

g _ _ _7--- y 5 - - - ~L ug<x.--: = t - = r_.=r -- = -r. = r - - : -.-- - - m - - -. t# f - : l. = g : :_-- =h- :.;._-- ;: :: __ g - -.;_-.__;---.._-._ ~ 5 N 5 0-o 't h D**D D'T 7 (feet) d as .o o e Figure 2.7-5 Impulsive Impact Loading Region for Structures with a Natural Frequency Less Than 30 Hz 1365 020

_ 14. Calculate the drag pressure,(Pimpactvelocitybytnemethodsoh9)corr a. dr 5ections 2.7.1 or 2.7.2. b. Form the ratio: Pdrag/#sax where Pmax is the amplitude of the parabolic pulse used in the original stress analysis multiplied by 1.35. c. Determine the dynamic load factor (DLF) from Figure 2.7-6, corresponding to the two cases: (1) parabolic pulse without drag and (2) parabolic pulse followed by drag, d. Multiply the calculated stress by the factor DLFwith drag /DLFw/o drag 2.8 FROTH IMPINGEMENT AND FALLBACK LOADS Froth is generated b vent header and (2) y (1) impact of the rising pool surface on the bubble breakthrough, as described in Section 4.3.5 of the LDR. The following load specification was derived from the high-speed film records of various pool swell tests and an analysis of pool acceleration following vent header impact. The impingement loads for Region I and Region II and the froth fallback loads, as described in Section 4.3.5, shall be defined as follows: Y Af P = f 144 gc where: P = froth impingement pressure (psi) f 3 of = froth density (1b /ft ) m V = froth imoingement velocity (ft/sec) g = gra"itational constant (ft/sec2) e Region I: The froth velocity shall be based on a source velocity equal to 2.5 times the maximum pool surface velocity prior to vent header impact, which is corrected for subsequent deceleration due to gravity starting at the 45 tangent on the bottom of the vent header, as shown in Figure 4.3.5-1 of the LDR. The froth density shall be assumed to be 20% water density for structures or sections of structures with a maximum cross-sectional dimension of less than 1365 021 D""D j'T }$ oJu,3..kUL 'o o P M max dra [ s / E TIME ....i V L.... J..i . _ _.. j.I_.. L _ I. .. [. [ 2. 0 -.- p i i i i P o.b - ~ [ _f I i I '.._.o V V / X ) l ~~ ~~~ ~~ hh/ fy . _L _.. /$/XV/i i 3 !.._. i / i/ __ a VI I i / /I /l I l I S l l l .i // l i 2

  • o _'

. [_ l i /f/ i I i 1 g I ,/i I i i S / / 1 i i 3 l [ ![ [b 4x l l { j / /7 i " P" 3 _.e / 7. _.. 7. T = Natural Period I N t t O. 0.1 0.2 03 0.4 05 t/ TN Figure 2.7-6 Effect of Drag Following a Parabolic linpact Pulse 1 % 3 022

, or equal to one foot, and a proportionately lower density for structures greater than one foot; i.e., o = (0.2/x) o, where x is the dimension in feet. w The load shall be applied in the direction most critical to the structure within the 90* sector bounded by the horizontal opposite the vent header to the vertical upward as shown in Figure 2.8.1. The load shall be assumed to be a rectangular pulse with a duration of 80 milliseconds. Region II: The froth velocity shall be based on a source. velocity equal to the maximum pool surface velocity, directly beneath the structure under consideration, which is corrected for subsequent deceleration from the elevation of the maximum velocity. The froth density shall be assumed to be 100% water density for structures or sections of structures with a maximum cross-sectional dimension less than or equal to one foot, 25% water density for structures greater than one foot, and 10% water density for structures located within the projected region directly above the vent header. The load shall be applied in the direction most critical to the structure within the + 45* sector of the upward vertical. The load shall be assumed to be a rectangular pulse with a duration of 100 milliseconds. Fallback: The froth fallback velocity shall be based on the freefall velocity from the upper surface of the torus shell directly above the subject structure. The froth density shall be assumed to be 25% water density, with the exception of the projected region directly above the vent header which is 10% water density. The load shall be assumed to directly follow the froth impingement load, with a duration of one second. 2.9 POOL FALLBACK LOADS The proposed load definition procedures set forth in Section 4.3.6 of the LDR for suppression pool fallback loads on internal structures following pool swell are acceptable. 2.10 Vent Header Deflector Loads The load definition procedures set forth in Section 4.3.9 of the LDR are applicable only to the four deflector types shown in Figure 4.3.9-2 of the LDR, and are generally acceptable, subject to the following constraints and/or modifications: 1365 023

..d J

  • o VENT HEADER q.a Tomus Figure 2.8-1 Direction of Load Appli;ation for Froth Region I

.',63 024

. 2.10.1 QSTF Oeflector Load

  • An individual plant may choose to use deflector load data taken directly from the QSTF plant-unique tests. This technique is subject to the following requirements:
1. If the QSTF deflector load measurement does not have a sufficiently fast response time to resolve the initial impact pressure spike for the deflector types 1 - 3, inclusive, the loading transient shall be adjusted to include the empirical vertical force history of the spike shown in Figure 2.10-1.

This impulse need not be applied for the type 4 deflector.

2. The QSTF plant-unique loads shall be adjusted to account for the effects of (a) impact time delays and (b) pool swell velocity and acceleration differences which result from uneven spacing of downcomer pairs. The longitudinal load variation shall be evaluated at the instant when the undistrubed pool surface would have reached the local elevatio a of the center (half-height elevation) of the deflector. The three-dimensional' load variation shall be based on the EPRI " main vent orifice" tests, as discussed in Section 2.5.
3. In applying the load to the deflector, the inertia due to the added mass of water below the deflector shall be accounted for.

The added mass per unit lengtn of deflector may be estimated by: % " 19 c Vw where: 4 = hydrodynamic mass per unit length (lb/ft) I = total impulse per unit length associated with the impact transient (lb-sec/ft) V = impact velocity (ft/sec) = deflector width (ft), as shown in Figures 2.104 w through 2.10-5 9c = gravitational constant 1363 025

. T] mw m W ju d gg F = Vertical upward force on aeflector per unit length 8* ... a ___

  • ua..._ __:.. t.:. _: t __. _r.. g d = Diameter of cylinder in

.7 .... y. 7 deflector types 1 - 3 ._.f p p .1

p... -l.

.. _ p~ ^ V = Impact velocity ' - ~ ~ - ---[ :_ -- r : h -..- p~ ~- ~ ~ ' o = Water density .k(: .[----~ ~ ~ [ ; + f~~:i:~~[ -- ' 6 t = Time from begining of impact

=.

=a. % _... I. _._.... F

=.. = -

-n.__. 4 r-l pVd U :l_k l !ir ~ E _ _._' E Nc3=, 2 ( -- *

  • I
{ -- * ~ ~ ".! CC 2

I. ..!._.i.. J - * ~ ~ -

  • - - I--- -

k ~! t < _.:.. :fi - Uh:::1- .{.._... .,__ g _ . g '.. ( _. i. . } }. ..[...._.p. I i,....:. 1..___.. .._._____L.'-:_' Z:." 1_:r:. _ L :. x - 2 _d.2- :--:-- u

hn___"]

_ct.-[_-

E.: :... :==:. _:22: n+

. r. __._. m.Nn _ _ : t u.___: =d.....__:-t.


r-2.____

-- - * - - f r___ t _. : In.. t:: 1- :t: _:-._ __.-t ri _f=ricrf:. A. \\u mr.-- :k_.. -.. ___i _.._.. r--.._ r: $ _ __ _ __3

: r:.

rx - h:_t:c

mr_

c

l

--j if.p fi F %;;.T- }E.lh_Fj._i((-{-.21-{ - --: - pu -u....._ r t. c r:- p: ;; r : :r:. %; =:t_.z :__.a= :2. t---- i ...,.__t_ __\\: i _ ___.__.. &....:- : f.- ] . N --- -l

:. [.

q. 4 j p :.- j o. 0.1 0.2 Vt / d Figure 2.10-1 Impact Force Transient for Addition to the Empirical Data for Deflector Types 1 - 3. 1365 026

... 2.10.2 Analytic Deflector Loads The deflector load definition which is based on empirical expressions for impact and drag forces together with plant-specific definition of the pool swell velocity and acceleration transients, as described in Sections 4.3.9.1, 4.3.9.3, and a.3.9.4 of the LOR, is acceptable, with the following modifications: 1. The impact transient and " steady drag" contributions to the l'oad shall be computed from the correlations shown on Figures 2.10-2 through 2.10-5, for deflector types 1 - 4. respectively. For times past the periods shown, the last value shall be extended for the duration of the transient. 2. The three-dimensional load variation and timing shall be based on the' EPRI " main vent ortfice" tests, as discussed in Section 2.5. 3. The gravitational component of the acceleration drag shall be included in F, as defined in NEDO-24612. A 4. In computing the deflector response to the load, the added mass of the water shall be accounted for, as described in Section 2.10.1.3 above. 2.11 CONDENSATION OSCILLATION LOADS The following criteria have been developed in consideration of the fact that the " condensation oscillation" loads (i.e., high vent flow rate with low air content) have been derived from a single FSTF test run (M8). The condensation oscillation regime is a hamonic phenomena and, therefore, statistical variance or load magnitude uncertainty cannot be established from one test run. Although we conclude that the M8 tested conditions are conservative and prototypical for the Mark I design, a reasonable measure of the uncertainty in the loading function is necessary to assure the margins of safety in the containment structure. However, based on our assessment of the phenomenological studies conducted by the industry and the NRC Office of Nuclear Regulatory Research, we believe that the following load specifications are probably conservative and fom a sufficient basis to proceed with implementation of the Mark I Long Tem Program. We will require that the Mark I Owners Group confim the condensa-tion oscillation loads (i.e., torus shell loads, downcomer lateral loads, vent system pressure, and submerged drag source) by perfoming a sufficient number of additional large break, liquid blowdown tests in FSTF to establish the uncertainty in the load magnitudes. 1365 027

. o d-e mm 9 m g 9-o J u S . lk lni _, 6o n i V 8. 6. ._..__'I F .a -. ..._ g 2 aVd 2 4.- (C )= from Figure 2.10-2b D 2. n .\\ j it t L., I l i l ) O. --L--- 0.136 05 1.0 Vt/d Figure 2.10-2 Impact and Steady Drag Force Correlation for Type 1 Deflector \\b0 h 2. 1 [' w.. .i. __f....- . g i g .[7 ~] - p- ~ j._. i oVd 5 i < 5 x 10 F--. - 1 { t, ./ u 6 l ' +- - f - -l .} p __.. } }.. 4. g j g (C ). Vl '.l ll I !!l D l ..i ..l.., l,, g, ~ ~ ' ' - T.' - j --! oVd > 5 x 105 Fj j t u l j l , s --j l l : ! i p j

i b

I ~ !' l ll l lf; } - r- - -- '[- l 1 1 1 g. l.4 i, lj I-l ~ _.. i

1. ;

r--- '-- l l l f l l t l 4 1 1 l T' r-l j .+. ...y...-- l---{- I i l lI l' ~ b 0. 1' 5. 10, 50. 100. 2 V /gd Figure 2.10-2b (C )= for Type 1 Deflector D 1365 029

. -w-(F 3> i d=0 707w () 1 ^ D"*D D'9' b. V S 6' h. p 2 pVw M* 2 I /1 j .\\ .i / / l / I

0. 7..

e .. _. - - f -- i, \\ / t ~ y j \\ / \\ / / .0. _ L _ \\ .207 0. 0.136 05 0.854 1. Vt/d Figure 2.10-3 Impact and Steady Drag Force Correlation for Type 2 Deflector Obb b w ---. I d=0.414w L n V 2.90 F 2 p v, _g, 2 i. 0.7 /, .'.. i - s _... i / l \\ f r t s _ _ /_... O. t g .207 0. 0.136 1.0 1.60 2.0 Vt/d Figure 7.10-4 Impact and Steady Drag Force Correlation for Type 3 Deflector 3363 031

25 - w 90 db V mm '3' D _w & ,m1 A 4 J~: 5 ..!!! !. ' f'. f. : f.; f :. f f[< ' f i. : h. ; d.. -. 2 '- U -45. n~ '~ ~ ~ '.'.;' ~ ..s.- ... l : 1;.,;.. 1: ': ; - ; l...;.... g 3 g_... ;.- ....?.. . : ;;*::p: ;;.-+ 2 .r.r ;!. E;E :L:E;i..A ip...iij..gi i aun...eg--- ._..... f..* : ~* ~ ~ f*: - f.I:k". :.! h;. [{! -]- - . C '*f--- '5E[ g -., g -... -.7...f:...f... j r j. q: y. __.g,...._...__ ..., _ _-..; } :. _: t--- c' -- p -- ;;.- j...: -~.4 ..p. .+ y.. n____ ---~+-.g_ - t - -- -

p..;_

2 ~6.-.-.%-. r._., b ',f -[ IM ~II.E.lN"E M Eli"I- "~~~~~~8 2 1{ .... - =... = =.,--, - -+ _ I 2. {' t L.... -, -.7... f,. - ..[- -.} _.a g .. {.. .-} - ' - -+- ig:p t.:j:- :{:--g ; --[--j. l2.g.... s. y- ----

I:Y! :li I ijF-.j!2jh:t}$ili.iii ll;,.y= i.~

. f' ' [ f *~, " ,k :'... -.t--.--. ;.- t. . + -. -- l :- - I[ fj 4 !! :I.E+ y- [. T F :q(=),='J..-... t .j.

j..

I-g, .._4..._ ~ ;. - )' : i. -. d-.:[- lh. j.._[.._ 1 _..:_.;ji--J..j. i....[. h. ,1 7 T._I_E g O. pp--.p-p -. p-.: -.. I .} I-O. 1. 15 2. 2Vt/w Figure 2.10-5 Impact and Steady Drag Force Correlation for Type 4 Deflector }Nb 2.11.1 Condensation Oscillction Torus Shell Loads The load definition and assessment procedures set forth in Section 4.4.1 of the LDR for the condensation oscillation loads on the torus shell are acceptable, subject to the following confirmation: 1. Provided that the " rigid wall" load derivation technique described in NEDE-24645-P is demonstrated to be conservative, in response to Question 7 in our request for additional information (D. Eisenhut, NRC, to L. Sobon, GE, dated July 30,1979). 2. Provided that sufficient justification can be provided to exclude a condensation oscillation asymetric loading condition, in response to Question 2 in our request for additional infomation. 3. Provided that the uncertainty in the load magnitude is demonstrated to be less that the demonstrated conservati!ms in the load specification, by the testing program described above. For clarification, the load specification set forth in Section 4.4.1 of the LDR shall be used in conjunction with a coupled fluid-structure analytical model. The condensation oscillation loading for the IBA is a continuous sinusodial function with a peak amplitude and fre' Jency range of that specified for the " pre-chug" load. We will require that the conservatism in the IBA condensation oscillation loads be demonstrated as part of the response to Question 7 in our request for additional infomation for the specific flow regimes of interest. 2.11.2 CONDENSATION OSCILLATION 00WNC0!1ER LOADS 2.11.2.1 Untied Downcomer Loads The condensation oscillation downcomer lateral loads for " untied" downcomers shall be defined as described in Section 4.4.3 of the LDR, based on the methodology in NEDE-24537-P. However, in computing the dynamic load factors DLF = P) g max 1363 033

. where: P = maximum static equivalent lateral load for max plant-unique downcs.ner P) = maximum static equivalent lateral load in FSTF DLF = plant unique downcomer dynamic load factor DLF) = FSTF downcomer dynamic load factor the plant-unique loading condition shall be derived as follows. The plant-unique DLF shall be calculated using a damping value consistent with the requirements of Regulatory Guide 1.61, " Damping Values for Seismic Design of Nuclear Power Plants," and a natural frequency detennined from the structural analysts of the downcomer - vent header system. The plant-unique driving frequency shall be specified as that frequency in the range 4 - 8hz which produces the maximum structural response. The natural frequency and damping values for the FSTF DLF shall be conservatively established from a " pluck" test of an untied downcomer in FSTF, with a nominal water level of 3 feet 4 inches and an amplitude in the range of the response level. The driving frequency for the FSTF DLF shall be assumed to be 5.5hz. 2.11.2.2 Tied Downcomer Loads The condensation oscillation downcomer loads for " tied" downcomers, as described in Section 4.4.3 of the LDR, are unacceptable. We will require that a load specification be derived from the maximum dynamic load components on each downcomer in a tied pair. The load definition and structural analysis technique shall be confinned by comparisons of the predicted structural responses to the measured strains in the FSTF vent header and tie-bar. The FSTF natural frequency and damping values shall be conservatively established by performing a " pluck" test for a tied downcomer pair in FSTF, with a nominal water level of 3 feet 4 inches and an amplitude in the range of the response level. 2.11.3 Condensation Oscillation Vent System Pressure lead The load definition procedures set forth in Section 4.4.4 of the LDR for the oscillatory pressures in the vent system durir.g the condensation oscillation period, are acceptable subject to confirmation by the additional testing as described above. 's 2.12 CHUGGING LOADS 2.12.1 Chugging Torus Shell Loads The load definition and assessment procedure set forth in Section 4.5.1 of the LDR for the chugging concensation loads on the torus shell are acceptable, provided the " rigid wall" load derivation technique is demonstrated to be conservative in response to Question 7 in our request for additional information (D. Eisenhut, NRC, to L. Sobon, GE, dated July 30,1979). This load specification shall be used in conjunction with a coupled fluid-structure analytical model. 2.12.2 Chugging Downcomer Loads 2.12.2.1 Untied Downcomer Loads The chugging lateral loads on untied downcomers shall be defined as described in Section 4.5.3 of the LDR, which is based on the methodology in NE0E-24537-P, with the following exceptions: 1. The load specification for comparisen to the ASME code primary stress limits shall be based on the maximum measured resultant static equivalent load in FSTF. 2. The fatigue usage analysis for each downcomer shall be based on a statistical loading with a 95% probability of non-exceedance. 3. The multiple-downcomer loading to assest statistical directional dependence shall be based on a non-exceedance probability of 10-4 per LOCA. 2.12.? 2 Tied Downcomer Loads For tied downcomers, the strains in the downcomer itself shall be evaluated exactly as in the case of the untied downcomers, using tied downcomer data. The strain in the tie bar shall be evaluated by assuming that one of the two tied downcomers is subjected to a dynamic load of triangular shape, with an amplitude of: RSEL (RSEL = Resultant Static F = max , ft Equivalent Load) d where RSEL is the maximum measured RSEL for an untied downcomer during chugging, f is the lowest natural frequen.y of vibration of an untied downcomer for the specific plant, and the duration of the load, td, shall be assumed to the 3 mil'iseconds. The load direction shall be taken as that (in the horizontal plane) which result in the worst loading condition for the tie bar and its attachments to the downcomers, s 4, 2.12.3 Chugging Vent System Pressure Loads The load definition procedure set forth in Section 4.5.4 of the LOR for the oscillatory pressures on the vent system during the chugging period are acceptable. 2.12 SAFETY-RELIEF VALVE DISCHARGE LOADS 2.13.1 Safety-Relief Valve Discharge Device The acceptance criteria set forth below for the quencher discharge loads and submerged structure drag load source strengths are applicable only to the "T" quencher configuration described in Section 1.1 of NEDE-24542-P. For plants using other types of quencher discharge devices, the SRV discharge load definition, submerged structure drag load source strength, and pool temperature limits will be evaluated en a plant-specific basis. 2.13.2 SRV Discharge Line Clearing Transient The load definition and assessment procedure, described in Section 5.2.1 of the LDR, for the pressure and thrust loads on the SRV discharge line and quencher, which is based on the methodology presented in NEDE-21864-P and NEDE-23749-1-P, is acceptable. 2.13.3 SRV Air-Clearing Quencher Discharge Shell Pressure Loads 2.13.3.1 Methodology for Bubble Pressure Prediction The load definition procedures describad.n Section 5.2.2 of t'le LDR and the methodology in NEDE-21878-P for predicting the quencher bubble pressure are acceptable, with the following exceptions: 1. The load definition procedures set forth in Section 5.2.2 of the LOR are acceptable for SRV discharge line water-leg lengths less than or equal to 13.5 feet. In the event that the water-leg length for a particular discharge ifne exceeds 13.5 feet, the load prediction for a 13.5 foot water-leg length shall be Wed. 2. The proposed methodology for predicting bubble pressures due to SRV subsequent actuations is not acceptable. The pressure amplitude predicted for the SRV first actuation shall be used in conjunction with the bubble frequency range for subsequent actuation, as specified below, for structure, equicment, and piping assessment in response to events containing SRV sub-quent actuations. 1365 036

. 2.13.3.2 Methodology for Torus Shell Pressure Prediction Based on the predicted air bubble pressure-time histories, as discussed above, the torus shell pressures at various locations in the suppression pool shall be calculated by the load definition procedures described in Section 5.2.2.3 of the LDR in conjunction with the appropriate pressure attenuation model. For quenchers located on the torus center-line, the pressure attenuation model described in Section 2.4 of NEDE-21878-P in conjunction with the bounding factor presented in Section 3.2 of NEDE-21878-P shall be used. The load adjustment and attenuation factors proposed for the "off-center" T-quencher configuration presented in a meeting with the staff on May 30,1979, are acceptable. We will require, however, that this load specification and its bases be documented in a supplement to the LDR. 2.13.3.3 Multiple - Discharge Loads The torus shell loads due to multiple SRV actuations shall be calculated as follows: 1. The peak values of bubble pressure due to a single valve actuation shall be combined by linear superposition (ABSS method) with the aporopriate pressure attenuation model, as discussed above. All bubbles shall be assumed to oscillate in-phase with the frequency ranges specified below for both first and subsequent actuations. 2. In the event that the combined peak torus shell pressure exceeds 1.65 times the local predicted peak bubble pressure due to a single valve actuation, the resultant torus shell peak pressure for the design assessment may be taken at the lower value. 2.13.3.4 Frequency of Pressure Wave Form The pressure wave form predicted by the methodology described in Section 5.2.2 of the LDR within the following uncertainty ranges (stretched or compressed time scale) that will produce the maximum structural, equipment, or piping system response shall be used fo the design assessment: 1363 037

. 1. First Actuation - the frequency range shall be 0.75 times the minimia - ' *1cted frequency to 1.25 times the maxinium predicted .ce' 2. Suu<cquent Actuation - the frequency range shall be 0.60 ti ss the mimimum predicted frequency to 1.40 times the max" predicted frecuency. 2.13.4 SRV Discharce Line Reflood Transient The transient analysis technique to compute the plant-specific reflood heights in the SRV discharge line following valve closure, as described in Section 5.2.3 of the LDR and based on the methodology in NEDE-23898-P and NEDE-21864-P, is acceptable. 2.13.5 SRV Air and Water Clearing Thrust Loads The load definition and assessment procedure for the quencher and quencher support thrust loads, described in Section 5.2.6 of the LOR, is acceptable. 2.13.6 SRV Discharge Line Temoerature Transient The cransient analysis technique to compute the maximum temperature loads on the discharge line and quencher device, as described in Section 5.2.7 of the LDR, is acceptable. 2.13.7 SRV Discharge Event Cases The kind and number of SRV discharge events shall be based on the plant-specific system configuration and a conservative assessment of plant operational history. The following load cases shall be considered for the design assessment: 1. A first actuation, single valve discharge shall be considered for all event combinations involving SRV events. Single valve subsequent actuations shall be considered for the SRV, SBA, and IBA event combinations, as detennined from a plant-specific primary system analysis. 2. Asymetric SRV discharge, both first and subsequent actuations, shall be censidered for SRV, SBA, and IBA event corbinations. The degree of asymetric discharge for each event combination shall be determined from a plant-specific primary system analysis designed to maximize the asymetric conuition. \\.)0 3. ADS valves discharging on first acutations shall be considered for the SBA and IBA event combinations, followed by subsequent acutations determined from a plant-specific primary system analysis. 4. The maximum number of valves that will actuate for the SRV event combinations shall be determined from a plant-specific primary system analysis for the design basis transients, which assumes that all valves actuate at their set-coint pressures. All first actuations shall be assumed to occur in phase, followed by subsequent actuations determined from the primary system analysis. All of the event combinations above include the earthquake events (OBE and SSE) in combination with the SRV discharge events. 2.13.8 Suppression Pool Temoerature Limits As part of the PUA, each licensee is required to either demonstrate that previously submitted pool temperature analyses are sufficient or provide plant-specific pool temperature response analyses to assure that SAV discharge transients will not exceed the following pool temperature limits. 1. Local Temperature Limit The suppression pool local temperature shall not exceed 200 F throughout all plant transients involving SRV operations, for any quencher device that has (a) the hole diameter equal to and (b) greater than or equal hole spacing than that of the generic Mark I T-Quencher. 2. Local and Bulk Pool Temperature The local to bulk pool temperature difference shall consider the plant-specific quencher discharge geometry and RHR suction and discharge geometry. The analysis of the plant-specific local to bulk pool temperature difference shall be supported by test data from either the existing Monticello pool temperature data or in-plant tests. The " local" temperature is defined as the temperature in the vicinity of the quencher device during discharge. For practical purposes, the water temperatures observed in the sector containing the discharge device at shell locations on the reactor side of the torus downstream of the quencher centerline at the same elevation as the quencher device may be considered " local" temperatures. The ' bulk" temperature, on the other hand, is the tertperature calculated assuming a uniform distribution of the mass and energy discharged from the SRV. i365 039 3. Suppression Pool Temperature Monitor System The suppression pool temperature monitoring system is required to ensure that the superession pool is within the allowable limits set forth in the plant Technical Specification. The system shall meet the following design requirements: a. Each licensee shall demonstrate that there is a sufficient number and distribution of pool tempera-ture sensors to provide a reasonable measure of the bulk temperature. b. Sensors shall be installed sufficiently below the minimum water level, as specified in the plant Technical Specifications, to assure that the sensor properly mnitors pool temperature. c. Pool temperature shall be recorded in the control room. A sufficient number of temperature monitors shall be provided to permit the operator to establish the bulk pool temperature. Operating procedures and alarm set points shall consider the relative accuracy of the measurement system. d. Instrument set points for alarm shall be established, such that the plant will operate within the suppression pool temperature limits discussed above. e. All sensors shall be designed to seismic Category I, Quality Group B, and energized from onsite emergency power supplies. 2.14 SUBMERGED STRUCTURE DRAG LOADS 2.14.1 LOCA Water det Loads The load definition and assessment procedure described in Section 4.3.7 of the LDR, which is based on the " Moody Jet Model" (NEDE-21472-P), is ace ptable subject to the follow-ing constraints and/or modif(ations: 1. The plant-specific jet discharge velocity, V (t), and D acceleration, aD(t) = dV /dt (t), from the QSTF plant-D specific test series shall be used as the driving sources for the jet model. 2. Forces due to the pool acceleration and velocity induced by the advancing jet front shall be computed for struc-tures that are within four downcomer diameters below the downcomer exit elevation, even if the structure } ig 4l} 1s not intercepted by the jet. The ow field shall be computed by modelling the moving jet front as a hemispherical cap ce..tered one downcomer diameter (D) behind the " Moody" jet front positions, containing the same amount of water as the " Moody" jet, and moving with the velocity of the " Moody" jet front. The formulas for the hemisphere radius (R ) and the trajectory s of the hemisphere center (x ) are: e R (t)

  • D (f +(3x (t)/0) 3) 1/3 for yf(t) > D f

3 7 R (t) = 9*f(t)l/3 for x (t) < D 3 ( ) f 2D x (t) = x (t) -D for x (t) > D c f f x (t) = 0 for x (t) < 0 c f where x (t) is the position of the " Moody" jet front as a f function of time, as computed in NEDE-21472-P. Using formulas I and 2 in NEDE-21472-P and assuming an average constant acceleration of the particles contained within one downcomer diameter behind the " Moody" jet front, the cross-sectional area in this region can be approximated by: T (1 + 1/(1-x/x )b) = nD2 A (x,t) f where x (t) is the " Moody" jet front position as computed f in NEDE-21472-P. The volume contained in this portion of the jet can be obtained by integrating A(x,t) from (xf -D) to x for xf greater than D, and from x=0 to x = xf for xf f less than D. When the jet is modelled by a more realistic hemispherical c;tp, while conserving the total volume of the fluid, the cap radius and position is given b. the equations above. The equivalent uniform velocity and acceleration at the location of the structure (x,y) shall be obtained from the time dependent potential 4)(x,y,t) induced by the jet front: 9 1365 041

  • j (*'#'

- [T) r [d"s\\- 1 f"s (x - x*) U "s c T ( (dt/ 2\\j F where'r = {(x-x )2 y )b and y is the transverse distance 2 ofthestructurhfromthejetaxis,and(x-x)isthe. distance fromthestructuretotheeffectivejetfronkcenteralongthe jet axis. The potential is the superposition of the expansion and motion of the sphere as given in any standard hydrodynamics text (e.g. Milne Thompcson, Theoretical Hydrodynamics, Fourth Edition, pp. 455-566). The local uniform floi; velocity is U, (x,y,t) = 4 3 as in NE00-21471,while the acceleration is a(x,y,t) = BU. at This calculation need only be performed for r> R and x > x. s g If either of these conditions are not satisfied, the methodology in the LDR will bound the load and is, therefore, acceptable. 2.14.2 LOCA Bubble Drag Loads The load definition and assessment procedures described in Section 4.3.8 of the LOR, which are based on the methodology in NED0-21471 and experimental confirmation in NEDE-23817-P, are acceptable subject to the following constraints and/or modifications: 1. Flow Field a. 05TF plant-specific test results (NEDE-21944-P) will be used, b. Model E in NEDE-21983-P will be used for the method of images simulation of the torus cross-section. After contact between bubbles of adjacent downcomers, c. the pool swell flow field above the downcomer exit elevation will be derived from the OSTF plant-specific tests. } } h 'j O R

2. Drag Load Assessment Drag forces can be computed for circular cylinders as given in a. NEDO 21471, but a conservative drag coefficient of CD = 1.2 must be assumed, independent of the Reynolds number, b. Drag forces on structures with sharp corners (e.g., rectangles and "I" beams) must be computed by considering forces on an p equivalent cylinder of diameter Deq = 2 2 Lmax, where Lmax is the maximum transverse dimension. Lmax is defined as the dia-meter of a circumscribed cylinder about the cross-gection of 2 the structure. For example, Lmax equals (a2 + b 72 for a rec-tangular cross-section of sides a and b. c. Long slender structures must be considered in segments of length (L), which oc not exceed the diameter (D Jr Deg). Al ternatively longer segments may be used as long as the equivalent unifom flow velocity and acceleration are evaluated cor.sarvatively for every point on any such segment. d. Interference effects due to the proximity of walls shall be consider $d for each structural regment that has its center less than 1.5 diameters from a boundary. Interference effects between neighboring structures shall be considered whenever the centers of the segments are less than 3D, where D = 1/2 (D1 + D ), the average diameter of the two structures. 2 For structures near walls, the multiplier (1 + Aw shall be used to increase the acceleration drag and the mu)tiplier (1 + D ) shall be used to increase the standard drag. Aw w and Dw that bound theory and experiments are given below as functions xw (xw = r/D - 1/2, where r is the distance from the segment center to the boundary). 0.05 < x,< l.0 A = 0.05/x g g D,= 0.12/x, x,< 0.05 A = 1.0 D,= 2.4 For structures with neighbors that are less than 36 away and within 30' of being parallel, the multipliers (1 + A ) and t (1 + D ) shall be applied to the acceleration and standard I drag components respectively. Bounding expressions for AI and DI are given below as functions of x1 (xI = rl2 /U - 1 where rl2 is the distance between segment centers). }3bb I = E [\\ uj + O \\ 0,05,< XI<2 A 0.2 D2 2/ DI = 0.2/XI where D 's the diameter of the structure under consideration and D2 's tM diameter of the neighbor. If more than one neighbor must be considered, the AI and DI values may be sunrned over the neighbor structures. For x! less than 0.05, the two neighbor structures shall be considered as an effective single structure. The effects of wall proximity and neighbor structures may be superimposed in order to compute overall multipliers as follows: (1 + Aw+ AIk) k (1 + Dw + IDIk) k In situations where interference effects must be considered, but the correction techniques outlined above are not applicable, a detailed interference effe-ts analysis shall be performed. 2.14.3 QUENCHER WATER JET LOADS The load definition procedure described in Section 5.2.4 of the LDR, which is based on the methodology in NEDE-25090-P, is acceptable, subject to the appropriate documentation of the confinnatory tests discussed in NEDE-25090-P. 2.14.4 0UENCHER BUBBLE DRAG LOADS The load definition and assessment procedures described in Section 5.2.5 of the LDR, in NEDE-21878, and in NED0-21471-2, are acceptable gubject to the following constraints and/or modifications: 1. Flow Field The determination of the charging, formation, and rise of the a. oscillating bubbles is subject to the same conservative factors that are used for the quencher torus shell pressure loads, as described in NEDE-21878-P. b. Drag loads on the quencher anns and the SRV discharge line shall be computed on the basis of asymetric bubble dynamics. Either a full 1800 phase shift shall be considered for full strength bubbles on opposite sides of these structures, or a more detailed assessment of the asymmetry of the bubble source strengths and phasing must be obtained from the experimental information in NEDE-21878-P. Model E in NEDE-21983-P shall be used for the nethod of c. images representation of the torus cross-section. 1363 044

. 38. 2. Drag Load Assessmen a. Drag forces for circular cylinders shall be computed on the basis of acceleration drag alone, under the condition that UmT/D < 2.74, where Um is the maximum velocity, T 15 the period of bubble oscillation, and D is the cylinder diameter. For UmT/D > 2.74, the standard drag shall be included with the drag coefficient CD = 3.6 in order to bound the relevant experimental data. b. The constraints specified for the LOCA bubble drag load assessment also apply to the quencher bubble drag loads, with the exception of the drag coefficient. 2.14.5 LOCA CONDENSATION OSCILLATION DRAG LOADS The load definition and assessment procedures described in Section 4.4.2 of the LDR and the methodology described in NEDO 25070 are acceptable subject to the following constraints and/or modifications: 1. Flow Field a. An average source strength shall be established by considering equal source strengths at all eight downcomers in equation B-4 in NEDO 25070. A maximum source strength shall be defined as twice the average source strength. For each structure,the loads shall be u mputed on the basis of both the average source at all downcomers and for the maximum rource applied at the nearest downcomer. b. The fluid-structure interaction effects shall be included for any structural segment for which the local fluid acceleration is less than twice the torus boundary acceleration. This may be accomplished by adding the boundary acceleration to the local fluid acceleration. 2. Drag Load Assessment a. The constraints and mcdifications specified for the quencher bubble drag loads apply. b. These loads may be applied quasi-statically to structures, only if the highest significant Fourier components occur at frequencies less than half the lowest structural frequency. \\h 2.14.6 LOCA CHUGGING DRAG LOADS The load definition and assessment procedures described in Section 4.5.2 of the LDR and the methodology in NEDO-25070 for the pre-chug drag loads are acceptable subject to the constraints in Section 2.14.5 for the condensation oscillation drag loads. The application for the post-chug drag loads is subject to the following constraints and/or modifications. Flow Field a. The maximum source strength history shall be obtained by using the maximum measured pressure (not necessarily at the bottom center) in a Type 1 chug in equation B-4 of NEDO 25070, with f(r) based on the single nearest downcomer. For each structure, the phasing between the two nearest downcomers that maximizes the local acceleration shall be established. The local acceleration shall then be computed on the basis of the two nearest downcomers chugging at maximum source strengths at the above established phase relation. b. The fluid-structure interaction effects shall be included for any structural segment for which the local fluid acceleration is less than twice the torus boundary acceleration. This may be accomplished by adding the boundary acceleration to the local fluid acceleration. 2. Drag Load Assessment a. The constraints and modifications specified for the quencher bubble drag loads apply. b. Unless the lowest structural natural frequency times the duration of the " spike" in the source strength is greater than 3, the loads shall be applied dynamically. Either sufficient Fourier components will be included to bound the " spikes" or the load shall be applied in the time domain using the source time history. The term " spike" refers to the short-duration high overpressure peak, such as that exhibited in Figure 6.2.1-20 of NEDE-24539-P. 2.15 SECONDARY LOADS The following loading conditions may be neglected for the PUA: a. seismic slosh pressure loads b. post-swell wave loads c. asymetric pool swell pressure loads d. sonic and compression wave loads } 3[3 O e. downcomer air clearing loads

i 40 - l 2.16 OIFFERENTIAL PRESSURE CONTROL REQUIREMENTS Those licensees that use differential pressure control (aP) as a pool swell load mitigation feature for the LTP, shal; demonstrate confomance with the following design criteria as part of the PUA: 1. There shall be no unacceptable change in the radiological consequences of an accident as a result of the inclusion of the AP systm. 2. Steam bypass of the suppression pool via the aP system shall be eliminated by appropriate system design, or sucn bypass shall be demonstrated to be acceptable by calculation. 3. Design and installation of the AP system shall be commensurate with other operational systems in the plant. 4. When the AP system involves tha addition of containment isolation valves, the additional valves shall be included in the plant's Technical Specifications and the valve design and arrangement shall conform to the requirements of General Design Criterion 56 in Appendix A to 10 CFR 50 and the regulatory positions in Standard Review Plan Section 6.2.4. Subsequent to the PUA, a license amendment shall be submitted to incorporate the following Technical Specification requirements for the AP system: a. Differential pressure between the drywell and suppression chamber shall be maintained equal to or greater than "X" (where X is the plant-specific differential pressure and values less than one psid will not be credited for load mitigation), exceot as specified in b and c below. b. The differential pressure shall be established within 24 hours after placing the plant in the RUN mode, during plant startup. The differential pressure may be reduced below "X" psid 24 hours prior to a scheduled plant shutdown. c. The differential pressure may be ; educed to less than "X" psid for a maximum of four hours during required operability testing of (specify here those safety-related systems for which operability tests either release significant amounts of energy to the suppression pool or cannot te performed with the aP established), d. In the event that the specification in a above cannot be met, and the differential pressure cannot be restored within six hours, an orderly shutdown shall be initiated and the reactor shall be in a cold shutdcwn condition within the subsequent 24 hours. 1365 047 e. A minimum of two narrow range instrument channel shall be provided to monitor the differential pressure. Error in the aP measurement shall be no greater than + 0.1 psid or the allowable aP shall be increased to offset the error in the measurement. The instrument channels shall be calibrated once every six months. In the event that the measurement is reduced to one indication, operation is permissible for the following seven days. If all indication of the differential pressure is lost, and cannot be restored in six hours, an orderly shutdown shall be initiated and the reactor shall be in a cold shutdown condition within the subsequent 24 hours, 3. STRUCTURAL ANALYSIS AND ACCEPTANCE CRITERIA The staff finds the general analysis techniques and proposed structural acceptance criteria set forth in the " Mark I Containment Program Structural Acceptance Criteria Plant Unique Analysis Applications Guide," (PUAAG), NEDO 24583, Revision 1, dated July 1979, acceptable. The proposed criteria will provide a sufficient basis for demonstrating the margins of safety required or steel structures and piping in the ASME Boiler and Pressure N sel Code and for concrete structures in the American Concrete Institute Code. Revision 1 to the PUAAG was presented to the staff in a meeting on June 29, 1979. We will require that this revision be formally submitted to complete the documentation required for this program. 1365 048

- 42 4. REFERENCES (Listed by Report umber) NEDE-13426P " Mark III Confinnatory Test Program One-Third Scale Pool Swell Impact Tests, Test Series 5805," General Electric Proprietary Report, August 1975. NEDO '671 " Analytical Model for Estimating Drag Forces on Rigid Submerged Structures due to LOCA and Safety Relief Valve Ramshead Air Discharges," General Electric Topical Report, September 1977. NEDE-21472 " Analytical Model for Liquid Jet Properties for Predicting Forces on Rigid Submerged Structures," General Electric Proprietary Report, September 1977. NEDE-21864-P " Mark I Containment Program Final Report Monticello T-Quencher Test, Task Number 5.1.2," General Electric Proprietary Report, July 1978. NEDE-21878-P " Mark I Containment Program Analytical Model for Computing Air Bubble and Boundary Pressures Resulting from an SRV Discharge Through a T-Quencher Device, Task Number 7.1.1.2," General Electric Proprietary Report, January 1979. NEDO 21888 " Mark I Containment Program Load Definition Report," Revision 0, General Electric Topical Report, December 1978. NEDE-21944-P " Mark I Containment Prnram Quarter feale Plant Unique Tests, Task Num M r 5.5.3, Series 2," Volumes 1 - 4, General Electric Proprietary Report, April 1979. NEDE-21983-P " Mark I Containment Program Submerged Structures Model Main Vent Air Discharges Evaluation Report, Task Number 5.14.1," General Electric Proprietary Report, March 1979. NEDE-23545-P " Mark I Containment Program 1/4 Scale Pressure Suppression Pool Swell Test Program: LDR Load Tests - Generic Sensitivity, Task Number 5.5.3, Series 1," General Elec,tric Prdprietary Report, December 1978. NEDE-23749-1-P " Mark I Containment Program Comparison of Analytical Model for Computing S/RVDL Transient Pressures ano Forces to Monticello Data," General Electric Proprietary Report Addendum, February 1979. NEDE-23817-P " Mark I Containment Program 1/4 Scale Test Report Loads on Submerged Structures due to LOCA, Air Bubbles and Water Jets, Task Number 5.14," General Electric Proprietary Report, September 1978. 1363 049

. NEDE-24537-P " Mark I Containment Long Term Program - Development of Downcomer Lateral Loads from Full Scale Test Facility Data - Task Number 7.3.2," General Electric Proprietary Report, May 1979. NEDE-24539-P " Mark I Containment Program Full Scale Test Program Final Report, Task Number 5.11," General Electric ' Proprietary Report, April 1979. NED0 24583 " Mark I Containment Program Structural Acceptance Criteria Plant Unique Analysis Applications Guide, Task Number 3.1.3," General Electric Topical Report, December 1978. NEDO 24612 " Mark I Containment Program Vent Header Deflector Load Definition, Task Number 7.3.3," General Electric Topical Report, April 1979. NEDE-24645-P " Mark I Containment Program Analysis of Full Scale Test Facility for Condensation Oscillation Loading," General Electric Proprietary Report, May 1979. NEDO 25070 "enalytical Model for Estimating Drag Forces on Rigid Sutmerged Structures Caused by Condensation Oscillation and Chugging Mark I Containments," General Electric Topical Report, April 1979. NEDE-25090-P " Analytical Model for T-Quencher Water Jet Loads on Submerged Structures, Task Number 5.14.2," General Electric Proprietary Report, May 1979. 1363 050

.}}