ML20116C953
ML20116C953 | |
Person / Time | |
---|---|
Site: | Prairie Island ![]() |
Issue date: | 03/31/1985 |
From: | NORTHERN STATES POWER CO., WESTINGHOUSE ELECTRIC COMPANY, DIV OF CBS CORP., WISCONSIN ELECTRIC POWER CO. |
To: | |
Shared Package | |
ML20116C944 | List: |
References | |
NUDOCS 8504290225 | |
Download: ML20116C953 (46) | |
Text
,, f. '
WESTINGH6USE PROPRIETARY CLASS 3 4
UPI LOCA PROGRAM PLAN USING SECY 83-472 METHODOLOGY Prepared by Northern State Power Wisconsin Electric Power Westinghouse Electric Corporation l
March 1985 I
{
i 8504290225 850419 FDR P
ADOCK 05000282 pon l
l
)
- i. t WESTINGHOUSE PROPRIETARY CLASS 3 UPI LOCA PROGRAM PLAN USING SECY 83-472 METHODOLOGY CONTENTS Section Title g 1 Sunnary 1 2 Introduction / Background 2 3 UPI LOCA Plan 4 3.1 Introduction 4 C. E 3.2 [ ] Code Description and Validation Program 5 3.3 Program Objective 14 3.4 Detailed Work Plan 15 4 Discussion of NRC Concerns 29 5 Program Plan Schedule, and Proposed Meeting 37 with NRC Staff 6 Conclusion 39 7 References 40 l
l C,E AppendixA-E ] Nomenclature 42 l
l
\
l .
g
. _ _ _ _ _ _ _ . _ - _ _ _ _ . _ . _ . . . _ _ _ _ . _ _ . _ _ _ _ _ _ _ _ _ - - _ m
\
WESTINGHOUSE PROPRIETARY CLASS 3 1.0
SUMMARY
This plan was prepared to address the NRC's concerns regarding the capability of the existing licensing evaluation model to accurately predict the perforiaance of the upper plenum injection (UPI) emergency core cooling system for a postulated Loss of Coolant Accident (LOCA). The objective of this plan '
will be to [
-. a.
_. This plan will cover cores with [
. ]'The thermal-hydraulic calculational tool which will be used in this program is the [ ]"i: ode. This program plan also includes a new method of assessing the LOCA margin for two-loop PWRs with UPI using the -
statistical methodology as described in SECY-83-472.
e 1
WESTINGHOUSE PROPRIETARY CLASS 3
2.0 INTRODUCTION
/ BACKGROUND After the core cooling hearings of 1972 to 1973,I } Westir.ghouse developed a series of computer codes which predicted the performance of the emergency core cooling systems for those PWRs with cold leg puriped safety injection and cold leg accumulators. These computer codes and modals were subsequently reviewed and approved for use on all Westinghouse PWRs with cold leg injection. It was j felt that the two-loop plants with upper plenum injection (UPI) would be covered by the approved evaluation models at this time.
The ECCS evaluation model for two-loop plants qquipped with UPI originally -
assumed that water injected into the upper plenum did not interact with steam and entrained water rising from the core and the UPI water was assumed to be delivered directly to the lower plenue through the cold leg injection -
locations. This approach was defended on the basis that a potential negative UPI impact of increased steam binding would be balanced by the positive UPI
. impact of improved core heat transfer.
In November 1977 the NRC issued a report (Reference 2) describing a model they devel> ped based on recent experimental data. Their report concluded that the then-current Westinghouse approach [ .
]Q8y calculating increased steam generation due to ton injection without giving credit in the hot channel for improved core heat transfer, the NRC model necessarily produced lower flooding rates and l hight'r peak cladding temperatures than the then-current Westinghouse model (all <sther factors being equal).
l This NRC model was then used to provide bases for continued operation of the two-loop plants, resulting in penalties in peak clad temperature to be applied to results of the old model (References 3 and 4).
e i
2
1 WESTINGHOUSE PROPRIETARY CLASS 3 The NRC' reviewed the Westinghouse revised model and found that it was acceptable as an interim basis for continued operation; however, they did not regard the model as a final accepted model and required that the twa-loop utilities continue to work to produce a final acceptable ECCS evaluation model or make a plant hardware change to eliminate UPI. .
In addition to the review of the Westinghouse UPI evaluation model, the NRC has also sponsored UPI ECCS experiments in the Japanese cylindrical core test facility (CCTF) to examine UPI phenomena, as well as better estimate PWR
, calculations using the TRAC code to examine UPI effects and to contrast the PWR response with UPI to a plant with cold leg injection. The tests and TRAC calculations did show thermal-hydraulic phenomena which were not accounted for explicitly in the Westinghouse evaluation model. However, there was no conclusion that the UPI emergency core cooling system was inadequate; rather, '
the NRC felt that the evaluation model which was used to estimate its perforiaance had an unacceptable level of uncertainty and that this evaluation 1 model should be improved. In the February 1985 letters to the two-loop utilities I ) the NRC identified ten areas where improved modeling was felt necessary.
Aprogramplanhasbeendeveloped,usingthe(' computer code, to
] c.f.
address each of these above areas. Thisplanwillalso[. .
. . . . . . . [willresultinanacceptablemodelwhichmeetstheAppendixK requirements.
. 3
a f WESTINGHOUSE PROPRIETARY CLASS 3 3.0 UPI LOCA MARGIN PLAN
3.1 INTRODUCTION
This program will utilize the Westinghouse [
]which is specifically designed to bee-
)
]presentedinSECY-83-472.
The SECY-83-472 document states that there are three areas of conservatism in the current licensing models: the required Appendix K conservatism, the conservatism imposed by the NRC staff to cover uncertainties, and conservatism imposed by the industry in many cases to reduce the cost for the analysis.
Basid on a review of the available experimental data and the best-estimate computer code calculations, the NRC staff has concluded that there is more than sufficient safety margin to assure adequate performance of the emergency core cooling system and that this excess margin can be reduced without an adverse effect on plant safety. Therefore, in the SECY-83-472 approach it is proposed that the licensee would utilize a realistic model of the PWR to calculate the plant response to a LOCA at the most realistic or most probable level (50% probability) and at a more conserv'ative 955 probability level. The calculation at the 95% probability would include uncertainties such as power i level, fuel initial temperature, nuclear parameters, and computer code uncertainties. The realistic PWR model t.nd the uncertainty analysis can be performed on a generic PWR model which covers a class of similar plants, that is, two , three , or four-loop PWRs such that all plants of that class are bounded by the generic plant model. The third calculation the licensee would have to perform would be to use a plant-specific realistic best-estimate model including the required Appendix K features, such as 1971 ANS decay heat + 20%
maximum stored energy, no rowet during blowdown, etc. This would be a new Appendix K calculation. The NRC staff will accept this new Appendix K calculation provided that the peak cladding temperature (PCT) calculated with the Appendix K calculation is greater than the PCT calculated at the 95 percent probability level but below the licensing limit of 2200'F.
4
. i 1 WESTIN6HdySE PROPRIETARY CLASS 3 The NRC staff interpretation of these results is that the required features of Appendix K have sufficient margin to cover all uncertainties inherent in a LOCA analysis combined at a 955 probability level. Such a series of calculations would provide an acceptable licensing basis by conforming to the Appendix K requirements and is estimated to result in 500*F peak clad temperature margin which can be used for operational flexibility, low leakage .
loading patterns to address pressurized thermal shock concerns, upratings, and more economical fuel designs. ;
3.2 [ bOEDESCRIPTIONANDVALIDATIONPROGRAM The[ ]c..ecomputer program has been developed to predict the thermal-hydraulic response of nuclear reactor primary coolant systems to small-and large-break loss-of-coolant accidents and other anticipated transients. It is derived from the merging of COBRA-TF for the reactor vessel and TRAC-PD2 for the loop components.
The COBRA-TF computer code provides a two-fluid, three-field representation of two-phase flow. Each field is treated in three dimensions and is compressible. Continuous vapor, continuous liquid and entrained liquid drop
' are the three fields. The conservation equations for each of the three fields and for heat transfer from and within the solid structures in contact with the fluid are solved using a semi-implicit, finite-difference numerical technique on an Eulerian mesh. The COBRA-TF vessel model features extremely flexible i noding for both the hydrodynamic mesh and the heat transfer solution. This flexibility provides the capability to model the wide variety of geometries i
encountered in vertical components of nuclear reactor primary systems. The flexible noding scheme allows representation of single hot assemblies with their unique hardware and power characteristics while permitting coarser noding for the remainer of the core. This reduces the run time of the calculations.
TRAC-PD2 is a systems code designed to model the behavior of the entire reactor primary system. It features special models for each component in the ,
i 5
= i WESTINGHOUSE PROPRIETARY CLASS 3 system. These include accumulators, pumps, valves, pipes, pressurizers, steam.
generators, and the reactor vessel. With the exception of the reactor vessel, the thermal-hydraulic response of these components to transients is treated with a five-equation drift flux representation of two-phase flow. The vessel component of TRAC-PO2 is somewhat restricted in the geometries that,can be .
modeled and cannot treat the entrainment of liquid drops from the continuous liquid phase directly, and thus, is replaced with the COBRA-TF vessel model.
The-TRAC-P02 vessel module has been removed and COBRA-TF has been implemented as the new vessel component in TRAC-P02. The resulting code is [ 3c,E The vessel component in[ ]C.E.has the extended capabilities provided by the three-field representation of two-phase flow and the flexible noding. The code has been assessed against a variety of two-phase flow data from experiments simulating important phenomena anticipated during postulated acpidents and transients.
Multiphase flows consisting of two or more fluids are separated by moving phase interfaces. The liquid and vapor phase properties are assumed to change discontinuous 1y across these interfaces. In general, the phases present can be any combination of liquid, solid, or gas. The flow pattern can take any one of a wide variety of forms, such as bubbly flow, droplet flow, gas-particle flow, and stratified flow. Exact conservation equations can be derived for each phase, and jump conditions relating variables on each side of the phase interface can be formulated, much as in single-phase shock wave theory. With appropriate initial and boundary conditions, these equations could theoretically be analytically solved for the exact motion of each phase and the phase interfaces. Except in a few simple cases, however, an exact analysis of multiphase flow is impossible because of its complex and essentially chaotic nature. The quantities of interest are the average behavior of each phase in the control volumes. Thus, most work in multiphase flow is done using average equations across the control volume.
C The average conservation equations used in the vessel module of [ ]a,E.re derived following the methods of Ishii . The average used is a simple Eulerian time average over a time interval, At assumed to be long enough to 6
__ _ . _ - _ _ . = . __ . _ _ _ _ - _ _ - _ _ _ _ - . . _ _ . - _ _ _ _ . --
WESTINGHOUSE PROPRIETARY CLASS 3 smooth out the random fluctuations present in a multiphase flow but short enough to preserve any gross unsteadiness in the flow. The resulting average equations can be cast in either the mixture form or the two-fluid form.
Because of its greater physical appeal and broader range of application, the two-fluidapproachisusedasthefoundationfor[ ]CE i
The two-fluid formulation uses a separate set of conservation equations and constitutive relations for each phase. The effects of one phase on an'other are accounted for by interfacial interaction terms appearing in the j equations. The conservation, equations have the same form for each phase; only the constitutive relations and physical properties differ.
c,E.
3 The three-field formulation used in the vessel module of [ 3isa straightforward extension of the two-fluid model. The fields included are i vapor, continuous liquid, and entrained liquid. Dividing the liquid phase into two fields is the most convenient and physically reasonable way of
. handling flows where the liquid can appear in both film and droplet form. In such flows the motion of the droplets can be quite different from the motion of the film, so a single set of average liquid phase equations cannot adequately describe the liquid flow or the interaction between liquid and vapor.
I t
Use of a third field, droplet flow, permits more accurate modeling of thermal-hydraulic phenomena such as entrainment, doentrainment, fallback, liquid pooling and flooding. These effects are significant when modeling the
- upper plenum injection system performance.
c Thephasicconservationequationsforthe[ ] v,E.
essel model describe the time-averaged behavior of a phase, which can be any phase in a multiphase
- flow. The phasic conservation equations are general within the assumptions 1
1 listed below.
- 1. Gravity is the only body force.
7
. i WESTINGHOUSE PROPRIETARY CLASS 3
- 2. There is no volumetric heat generation in the fluid.
- 3. Radiation heat transfer is limited to rod to drop and rod to steam.
- 4. The pressure is the same in all phases.
- 5. The viscous dissipation can be neglected in the enthalpy formulation for the energy equation.
Four mass conservation equations are required for the vapor, continuous liquid, entrained liquid, and noncondensible gas mixture. In vector form they -
are, respectively, as follows:
+
ep + V * (a p Uy ) = r"' + V . G[
a **
g e,pg +Y*(ae,U,)=-r,-S+V.G[
g
. g s,p g + V * (a,pg U,) = -r, + S= '
a =
g *g g + V * (e ggy p
e U ) = r g' + V
- G The individual terms in the equations become Rate of + Rate of mass = Rate of mass + Rate of mass + Rate of mass change gain by gain by gain by efflux due to of mass convection interfacial entrainment void drift transfer or chemical .
reaction l
8
. i
\
WESTINGHOUSE PROPRIETARY CLASS 3 Also, a gas mixture transport equation (for each species of noncondensible gas) is sc1ved explicitly at the end of each times step. The nomenclature used for these equations is given in Appendix A.
Two energy conservation equations in which the liquid and the entrained .
liquid are assumed to interact at a rate sufficient to nearly maintain equilibrium, are specified for the vapor-gas mixture and the combined liquid h(e,p ygh yg) + V * (s yp yghyg D y) = r"' hg +g gy + QQ - ? " (sy g fg) 3
E (*g + *e)p gh, + Y * (egggg e h U ) + V * (e,p g h,U,) = r"'hf+gh*O
- V * (a gT)
The individual terms in the equations become Time rate + Convection = Energy transport + Interfacial + Wall heat - Turbulent of change due to phase heat flux heat flux change transfer The use of a single energy equation for the combined continuous liquid and liquid droplet fields implies that both fields are at the same temperature.
In regions where both liquid droplets and liquid films are present, this assumption can be justified in view of the large rate of mass transfer between the two fields, tending to draw both to the same temperature.
Three momentum equations are solved in C08RA-TF, allowing the liquid and entrained liquid fields to flow with different velocities relative to the vapor phase. They are as follows:
g (*yygy)+V*(sp p U y ygUUy y) = s vP y + s yp gT tg av ty av
+ (r"'U)
+ V * (m y [g) g (a o U ) + 7 * (a p,U ggg g U)=
gg og VP + sgg eg t,, + tg -(r, U) - (S"'U) tv
+7* (e t )
n
WESTINGHOUSE PROPRIETARY CLASS 3 h(=,e,U,)+V -
(=,e,U,0,) =
., VP + e,p,g -I +i
.0)
- (r;e + (s, D) .
The individual terms in the equations become Rate of + late of = Pressure + Gravity + Wall change of momentum gradient force shear momentum change by .
convection
+ Interfacial + Interfacial + Interfacial + Momentum + Momentum drag between drag between momentum exchange exchange vapor and vapor and exchange due to due to continuous drops entrainment turbulence liquid .
C08RA-TF was developed for use with either rectangular Cartesian or subchannel coordinates. This allows a fully threewlimensional treatment in geometries amenable to description in a Cartesian coordinate system. For more complex or irregular geometries, the user may select the subchannel formulation (which neglects some of the momentum flux terms in the above equations) or a mixture of the two. The subchannel approach has been used by Wheeler et al., for bundle thermal-hydraulic analysis by the C08RA series of Codes .A comparison of the subchannel and Cartesian momentum equations is given in Section 2.3 of the Thurgood reportOI.
One of the important features of the C08RA-TF vessel model is that the governing equations form a complete set. No terms are omitted particularly in the momentum equations where wall shear, momentum exchange due to turbulence and all the interfacial drop terms are represented. The COBRA-TF vessel model w __ _ _ _ _ _ _ _ _ -
. 0 WESTINGH0'YSE PROPRIETARY CLASS 3 also has two energy equations such that thermodynamic non-equilibrium between the two-phases can be accounted for. .C ,
3 a.
c,E The currently planned program for [ .] code validation will [, ,
.. 3." System effects tests simulate most of the components of the PWR acting together during the transient and can be used to compare predicted event times and peak clad temperatures, as well as other local data. Single effects tests will provide data for specific scaled components of the PWR and can be used to C ' ] ct which for component models have a particularly important effect on[
3a c, E To examine [ ] blowdown performance, the L
] wS11 be predicted with [ ]."'khese tests are the best available system experiments to assess code performance during the blowdown
~
transient due to the relatively large scale of the facility and its nuclear core. The following quantities will be compared:
- CL The following are the most important phenomena affecting the prediction of peak clad temperature during blowdown:
cn.
Separateoffeettestswillbeanalyzedtovalidatethe[
~
fests will be used
[The{
11
WESTINGHOUSE PROPRIETARY CLASS .3 totestthepredictionofthe[ 3*
Predictions will be made of the[ ]Ind these will be compared to data.
I The[ .
t . )ca. will be used to assess the [ "
]cL The [. ]a.willbecomparedto[
'}ch-generated using a {
]&
Substantial validations exist for the TRAC PD-2 break model. Although uncertaintiesexistinthemodel,itisplannedin[
Multidimensionaleffectsmaybeimportantguringblowdown. The multidimensional capabilities of[ ] will be used to[
]9- The E
]CL. E.
To assess the predictive capability of [. ]C.Efor{
l
] The [ .
l hg There are integral systems effects tests which model a simulated PWR response for a reflood transient for a cold leg l
l injection transient; or a top-down flooding transient for[ ]a- The keyparameterstobecomparedinclude,[ ,
I i .
,] .1The [ _]R I are the largest scale system simulations that can be modeled and will provide usefulassessmentof[ ] cE.
12
WESTINGHOUSE PROPRIETARY CLASS 3 A
Inadditiontothe[ ] the [
]"bil be simulated with[ ft$assessthe These experiments are low pressure top
~
flooding experiments which can characterize the hot assembly thermal-hydraulic characteristics for a UPI situation. The key parameters to be examined in l these experiments [ ,
]Dhese experiments will help quantify the [ ] ]R since accumulator ECC condensation effects can influence the delivery and ECC bypass offacts; C 'illbeusedtomodelthe[
.]a.These comparisons should help reduce uncertainty in the [ i
]D c,G Itiscurrentlyplannedtoincludethe( _
v alidation [ .
3"Inaddition,the[ ]c"o'karisions with the [
Additionalverificationof{ ] with the [
. y l
l 13 l _ _ _ . _ - _ _ _ _ _ _ , _ . _ . _ _ - _ _ _ . _ . _ _ _ , _ , _ _ _ _ _ _ _ _ _ _ . _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ ___
~ ~
, O WESTINGHOUSE PROPRIETARY CLASS 3 ;
3.3 PROGRAM OBJECTIVE The objective of this program is to address the NRC concerns on the offactiveness of the upper plenum injection (UPI) emergency core cooling system to mitigate the consequence of a large break Loss of Coolant Accident .
(LOCA). This effort will [ ,
. ,, ]4Thisplanwillcover
]4.In addition the proposed plan will [ ,
~
3" 4
i 9
l Y
e 14
WESTINGHOUSE PROPRIETARY CLASS 3 3.4 DETAILED WORK PLAN Task 1. Set up a [ .
This task will [ . .
[The two-loop plant parameters. [,
J" a-bI 15
tfESTINGHOUSE PROPRIETARY CLASS 3 4
l ~
Task 2.[..
] Q., C.
The[
}"The {
4 In b .. .
]Q.There will also be J"
Therewillbeanumberof{
I l ]4 wo T options currently exist for the
{ i .
o,es
].4The[ *
]4Theendresultofthesecalculationswillbetodevelop a{ i
! ~
].4The work scope includes:
~ G 4
1 16
WESTINGHOUSE PROPRIETARY CLASS 3
_.. CL e
M M'
t 17
O O idESTINGHOUSE PROPRIETARY CLASS 3 Task 3. Run( . . . . . .... )
This task will use the two-loop [
Nny[ '
]"The results of these calculations will establish the C y
~
This task includes:
.w l
. em, 18
WESTINGHOUSE PROPRIETARY CMSS 3 l Task 4. Perform Plant Sensitivity Studies This task will determine the relative [ ,
The work scope includes:
a e
l e
f
~ ,
I 19 I
e WESTINGHOUSE PROPRIETARY CLASS .,3 Task 5. Perform Sensitivity Studies of Model Parameters In this task the effect on the calculated peak clad temperat*-e of the [
J The .
work scope includes:
~
Q_
e I
l ND 20 l - -.
WESTINGHOUSE PROPRIETARY CLASS 3 Task 6. Perforu {" :
]
This task will perform additional [
t'
[Thiseffortincludes:
9 6
4 L
i
(
i anse 1
21
i . .
WESTINGHOUSE PROPRIETARY CLASS .4 Task 7. Prepare Report on Best-Estimate Code This task will prepare a technical report for transmittal to the NRC on the [ ] c. EThe report will specifically
, describethe[ .
-]"The work scope includes:
q-. a n) suosit snis report to tne nut;.
i t
1 l
22
v tfESTINGH6USE PROPRIETARY CLASS 3 Task 8. PerformAnalysisof[ F This task will modify the [. .
.]O In addition, [ ,
]4This task will permit analysis of [
The specific work scope includes:
WP l
M Sumu '
23
WESTINGHOUSE PROPRIETARY CLASS 3 Task 9.[
. 3A Usingthetablesof[ , .
]a.Tothiscalculated[
\
I -
. s
]g The specific work' scope includes:
Q_
f s
l 24 1
WESTINGHOUSE PROPRIETARY CLASS 3 Task 10. [
3a RepeatTask9usingthe[.
]a. Add to the [
t
~
]
Thetotal[
]0.The specific Work scope includes:
- - a_
l i
I _
i i
l 25 l
O'
- tfESTINGHdDSE PROPRIETARY CLASS 3 Task 11. Perform a plant specific. { [
The [. .
]&The models~and assumptions which Will be added are [ .
]4,khiscalculationC. ~ ~
]Otherplant-specific [
, ']"Thespecificworkscopeincludes: -
GL k
1 l
l 26 l _. . - - _ _ _ . . _ - - - . . . _ - . - . - . - - - - - -- . - - - - - - - . - - - - - - - - - - - - - - - . - - - - . - - - - - - - - - - - - -
WESTINGHOUSE PROPRIETARY CLASS 3 Task 12. Perform a {, [
This task will calculate the [.. ,
. 3 The purpose of this [
show that the once burned fuel in a mixed core is less limiting than -
1 . . .
]aThehott ca.
]The specific work scope includes:
- O e
e M
eummun i
l l
l l
l i
. 27 l _. _ _ . _ _ . _ . . _ _
WESTINGHOUSE PROPRIETARY CLASS 3 Task 13. Prepare a Report Describing the [ 5
_] G'O A technical report will be prepared which describes the results of each of the [
~
t
].aThe report
{ . .. .
In addition, all (.
a,E
]The work scope includes:
CL d) Submit this technical report to the NRC.
28
- WESTINGHOUSE PROPRIETARY CLASS 3 1
4.0 OISCUSSION OF NRC CONCERNS The thermal hydraulic development effort presented in Section 3 will result in E
, 3 utilizing the ['
]$hamountof
. a g
In addition, it is felt that the [
] permit a more accurate calculation of -
the core behavior such that the [
]ENRCconcernsonthecurrentWestinghouseUPImodelidentifiedin the February 1985 letter to the two-loop utilities (5) are addressed below with a brief discussion of how the approach given in Section 3 will address the concern.
- 1. NRC Concern - Metal Heat Transfer to UPI Water Heat transfer from upper plenum metal to UPI water is calculated by a lumped thermal capacitance model. It is necessary to show that either (a) the i
coefficients used in the model are correct (b) sensitivity studies show that i
results are insensitive to the values chosen for the coefficients, or (c) the coefficients selected are conservative for all cases and conditions of interest.
,_ a, E, c.,
l f
l l
29
WESTINGHOUSE PROPRIETARY CLASS 3
%Ek s
~ ,
- 2. NRCConcern-UnnerPlenumInIectionFlowDistribution Westinghouse (M) assumes that safety inje: tion (SI) forces ECC water into the upper plenum and this water covers only part of the core. The percent core coverage is important because the M UPI EM uses different analyses for the covered and uncovered core regions.
Acceptance of the partial core coverage concept would require explanation of experimental data from the Semiscale and CCTF test facilities which indicate upper plenum pooling occurs in the test facilities and hence, may occur in UPI plants. It woyld also require determination and analysis of worst case coverage conditions and consideration of steam-water interactions.
a , E ,C.
t USEE en. n 5
I 30
WESTINGHOUSE PROPRIETARY CLASS.3 l
[ 3. NRC Concern - Decaw Heat and Stored Enerav
?
4 This section of the Westinghouse EM submittal says, in full, that i
"The decay heat calculation uses ANS plus 20 percent. The initial stored energy at the beginning of reflood is based on the stored energy at the
- l end of blowdown plus a calculation that adds the decay heat generated through the refill period to the core stored energy."
In evaluating the significance of this portion of the proposed model, the staff notes that the docketed EM submittals indicate that the refill period includes 14.45 seconds of pumped UPI before bottom of core recovery (80CREC).
- The proposed Westinghouse EM ignores this UPI water during refill, which should be conservative relative to vessel inventory. However, it is necessary -
to address the effect of this 14.45 seconds of pumped UPI on (a) steam generation and steam flow or steam binding, (b) changes in fuel and metal wall stored energy prior to SOCREC, and (c) reflood, this being done with (d) 1 and l ~
2 trains of SI operable. Alternately, (e) it is necessary to show that the proposed treatment is conservative and that injection of UPI water before 80CREC will not cause steam binding and delay SOCREC.
i
% E. C-e
- i 4
i h
i i
- 4. NRC Concern - Fuel Rod With Too Quench Front
! Operation of a safety injection train injects ECCS water into the upper plenum. The Westinghouse model assumes that this UPI water covers a fixed percentage of the core, and that it will flow downward into that part of the l core and cause a top-down quench.
4 31 l
j ..
WESTINGHdUSE PROPRIETARY CLASS 3 i
Westinghouse assumed that the time required for the top-down quench to reach the midplane of a fuel bundle was a specific linear function of the UPI mass 5
flow rate per fuel assembly. Adequate technical justification of the function
] describing the top-down quench would require (a) much more specific 3
identification of the experimental data used to develop the quench-time
] function, (b) showing that the data applies to the UPI situation, and (c)
- I showing that the flow per assembly is a valid correlation parameter.
j Westinghouse stated that if used a quench-time function which is a " bound to
- the [ experimental data) since more rapid quenching increases heat transfer to 4
l UPI water and results in more steam generation." This more rapid steam ~
i generation is presumably detrimental because it will impede reflood. However, l to justify use of this bounding quench function, it would be necessary to show
! that an overestimate of top-down quench speed (and/or of UPI flow per -
! assembly) will not cause an underestimate of peak cladding temperature.
i i
l a, E, c, i
i 5. NRC Concern - Core Heat Transfer Model The Westinghouse non-UPI EM uses very simple heat generation and heat transfer models to determine core exit fluid conditions. Westinghouse has made what appear to be relatively simple modifications to its non-UPI EM to account for
- heat transfer in the top-down quench region. No adequate justification is l given for using this simplified model under UPI conditions.
a, F, C.
~
J 32
WESTINGHOUSE PROPRIETARY CLASS 3 1- m et,E,C, l
J i
- 6. NRC Concern - Core Steam Generation It is necessary to clarify the location of the UPI water, and whether part of it is held up in a pool in the upper plenum, and it is necessary to demonstrate that the correlation for carryover from core to upper plenum is applicable to two-loop UPI plants.
The staff notes that water can be injected into the upper plenum at a rate which may be more than 10 times as large as the rate which water can be injected into the intact cold leg, and that these flows determine the peak cladding temperature after the accumulators are empty. The licensees must submit a reflood and refill model which considers these flows and relevant j .
thermal and hydraulic characteristics. The licensee must specifically show l how and why the water injected into the intact cold leg is or is not bypassed throughout the transient.
r-Q , E , C.
i
(
[
L
\
- 7. NRC Concern - Condensation
~
i The licensee must explain and justify (a) what happens if the falling subcooled UPI wter can condense all the rising steam in the covered region, I and (b) equations and assumptions coupling the covered and uncovered regions.
33
WESTINGHOUSE PROPRIETARY CLASS 3
- a,E.,C-
~
- 8. NRC Concern - Vertical Entrainment Steam may be generated by the UPI water falling through the core. This steam is assumed to entrain and carry part of the UPI water upward. Westinghouse has not shown that the experimental data on entrainment were applicable to the conditions in a UPI plant, or that there exists a reasonable technical basis for extrapolating the correlation outside the data base.
ct, E, c.
~
~7 i
i
" e 1
34
1 1
WESTINGHOUSE PROPRIETARY CLASS 3 j l
- g. NRC Concern - Morizontal Entrainment l 1
l Westinghouse assumed that 1.675 of the ECCS water which is injected into the upper plenum is entrained by horizontally moving steam which carries this 1.675 of UPI water into the hot legs. The Westinghouse entrainment model has not been proven valid because the model is based on (a) entrainment. test data which have not been demonstrated to be applicable to the conditions in the upper plenum of a UPI plant, and (b) air flow tests which have not been documented enough to be reviewed. Any new submittal will have to address the data from the CCTF experiments which imply that a frothy mixttre exists up to the het leg nozzle.
- G , E. , C. 1
- 10. NRC Concern - Total Steam Addition Due to UPI l No theoretical justification has been given for assuming that there is no j interaction of the steam-water mixture rising vertically from the botton l quench front and the steam-water mixture rising vertically from the falling ,
UPI liquid. Further, the quality of each steam-water mixture is calculated by l
a different method, and there has been no discussion of the effects of l difforent flow rates even though the two flow rates may be vastly difforent. ,
i 1
i 35 L
WESTINGHOUSE PROPRIETARY CLASS 3 t $
4 a
St>
e O
l 36'
WESTINGHOUSE PROPRIETARY CLASS 3 5.0 SCHEDULE AND PROPOSED MEETINGS WITH THE NRC STAFF Through the use of a more[ ,
) it is believed that the NRC technical concern on the[ ,
.e
]*Also,C...
]Ihe
]a,GA schedule for the program plan is given in Figure 1.
Meetings with the NRC staff are planned during this program and will provide theStaffwith{
].b These meetings will also serve to
. ]iey assumptions [
, . ]AThree meetings are proposed:
a, E.C-neue m
37
T C o
\'
9 4
g N
K
>=
. r w t N
W l i Q . 6. C e
38
e .
E STINGHOUSE PROPRIETARY CLASS 3
6.0 CONCLUSION
A detailed plan has been developed to address the NRC concerns on the performance of the upper plenum injection emergency core cooling system for a postulated large-break LOCA transient. This plan utilizes [ .
' .. J'.IhIandresultofthis effortwillbeto( , ,
0 l
i 1
I 4
39
WESTINGHOUSE PROPRIETARY CLASS 3
7.0 REFERENCES
- 1. Rulemaking Hearing: ' Acceptance Criteria for Emergency Core Cooling Systems for Light Water Cooled Nuclear Reactors" Docket R. M.- 50-1 (1973).
- 2. NRC Staff Report, Analysis Branch Division of Systems Safety " Safety Evaluation Report on ECCS Evaluation Model for Westinghouse Two-loop Plants " November 1977. ;
- 3. . . . . . .
et
- 4. I ~
] o.
- 5. Letter dated February 12.1985 J. R. Miller (NRC) to D. M. Musolf (NSP) titled, " Development of an Acceptable ECCS Evaluation Model Which Includes the Effact of Upper Plenum Injection Prairie Island Nuclear Generating Plant Units Nos. 1 and 2".
- 6. NRC Staff Report.
- Emergency Core Cooling System Analysis Methods,"
USNRC - SECY-43-472. November 1983.
c,E
- 7. Thurgood , M. J. , et al . , " [ ]- A Thermal Hydraulics Code for
- Transient Analysis of Nuclear Reactor Vessels and Primary Coolant l Systems," NUREG-CR-3046 Volumes 1-5, March 1983.
l S. ' TRAC-PC2: An Advanced Best Estimate Computer Program for Pressurized Water Reactor Loss of Coolant Accident Analysis' Los Alamos National
- Laboratory Report, LA-8709-MS, NUREG-CR-2054 (1981).
i
! 9. Ishii, M., 'Orif t Flux Model and Constitutive Equations for Relative j Motion 8etween Phases in Various Two-Phase Flow Regimes," ANL-77-47 I October 1977.
! 40 l
i
. .- . )
WESTINGH0USE pt0FRIETARY CLASS 3
- 10. Wheeler, C. L., e'.t al. " COBRA-IV-I: An Interim Version cf COBRA For Thermal-Hydraulic Analysis, Rod Sundle Nuclear Fuel Elements and Core,"
8WL-1962, March 1973.
- 11. Reeder, O. L., ' LOFT System and Test Description", NUREG/CR-024,7 .
July 1978.
- 12. McIntyre, 8. A. et. al. " Full Scale Controlled Transient Heat Transfer Tests-Data Report," EPRI-NP-1810, April 1981.
- 13. [ -A 14.[ . . .
- 15. Mirano, K., et al., Quick-Look Report on Large Scale Reflood Test-23",
JAERI - Memo 9767, October 1981.
- 16. Blaisdell, J. A., Waring J. P. and L. E. Nochreiter, "PWR FLECHT-SET Phase A Report, WCAP-8238, December 1973.
- 17. Nochreiter, L. E., et al., " Mixing of ECC Water With Steam: 1/3 Scale Test and Summary" EPRI 294-2, June 1975.
l 18, {
.t'~
l l
l l
41
7- -
QW e
APPENDIX A 46
,[
[ NOMENCLATURE O
42
APPENDIX A 1WO FLUID PHASIC CONSERVATION EQUATIONS APPENDIX NOMENCLATURE e internal energy ' '
f body force g accelgrationofgravity h enthalpy I phase. interface 1 identity tensor
- k thermal conductivity 1/Lj surface area concentration for the jth interface 1/Ls total surface area concentration M interfacial somentum exchange
~
g interfacial drag force I
M interfacial momente exchange due to phase change a mass Slow rate n unit exterior normal vector P pressu're 1 conduction heat flux 0 volumetric heat generation rate .
T A turbulent heat flux gj' interfacial heat transfer per unit volme i
5 phase interface stress tensor I turbulent (Reynolds) stress tensor T temper'a ture 1 stress vector t time at averaging time interval U fluid velocity 3 interface velocity V velocity of volme V A.1
., , vo . .-
\.
~ 5 t
o V arbitrary volme i position vector .
Greek Symbols SV boundary of V vj norwel displacement speed of interface j '
a void fraction 8 phase function
~
r not rate of interfacial mass transfer per mit volume '
A net rate of interfacial energy transfer per mit volme A second viscosity coefficient u viscosity p density I surface tension vector
- viscous stress tensor (stress deviator) 9 arbitrary fluid property subscripts 1 phase 1 ,
2 phase 2 Ig phase 1 interfacial limit 2g phase 2 interfacial limit '
I interfacial j Jth interface k phase k -
kg phase k interfacial limit Superscripts ,
1 interfacial surface average a material T turbulent '
t transpose A.2