ML20003F925
| ML20003F925 | |
| Person / Time | |
|---|---|
| Issue date: | 12/16/1980 |
| From: | Minogue R NRC OFFICE OF NUCLEAR REGULATORY RESEARCH (RES) |
| To: | Harold Denton Office of Nuclear Reactor Regulation |
| References | |
| RIL-107, NUDOCS 8104230895 | |
| Download: ML20003F925 (22) | |
Text
s. 2.
e
,....,'o, UNITED STATES
+
['
- ~j NUCLEAR REGULATORY COMMISSION y.,
2 a wAsm NGTON, D. C. 20555
.j.
DEC 161980 MEMORANDUM FOR: Harold R. Denton, Director Office of Nuclear Reactor Regulation FROM:
Robert B. Minogue, Director Office of Nuclear Regulatory Research 107
SUBJECT:
RESEARCH INFORMATION LETTER RELAP4/M007 COMPUTER PROGRAM
- 1. 0 INTR 000CTION The subject of this Research Infor: ration Letter is the RELAPt/M007 computer 7
program that was released to the National Energy Software Center in April of ~~-
1980 (Ref. 1).
RELAP4/M007 is a best-estimate computer program designed to calculate thermal-hydraulic phenomena in a nuclear reactor. Although this code was originally designed to address PWR large-break loss-of-coolant accidents (LOCAs), capability exists for analyses of certain small-break LCCAs and non-LOCA transients / accidents in PWRs, together with a limited ca:: ability to address accidents and transients in BWR reactor systems. The relen a of RELAP4/M007 computer program represents the completion of the RELAF4 series.
The RELAP4 series, and the RELAP4/MCD7 in particular, was developed in response to a number of requests, including:
NRR: Provide capability for a statistical LOCA study and for sensitivity studies to assess data pertaining to LOCA rule changes (Ref. 3).
ACRS: Iaprove blowdown and reflood understanding and provide quantitative p praisals of ECCS (Refs. 4 and 5).
Two Research Information Letters have been issued (Refs. 6 and 7), describing the RELAP4/M006 code (Ref. 2) and results of its assessment. The major conclu-sions cited in these letters are:
1.
Good results are generally obtained for blowdown regime.
2.
Accuracy in the reflood regime is not as good as desired in all cases.
3.
Refill is not adequately calculated.
4.
A continuous integral LOCA calculation cannot be performed.
The letters also identified certain weaknesses and deficiencies in the RELAP4/
M006 code, indicating the need to improve accuracy of calculations for the reflood and refill periods of LOCA.
I
(
In RELAP4/M007, many of thesa deficiencies have been corrected and further improvrments in modeling have been made. This document contains a summary of the mo4e important models and features incorporated in the RELAP4/M007 code, the results of the code checkout (developmental assessment) calculations, and recommendations for applying the program to various reactor accidents or transients.
It should be noted that predictive capabilities of this code have not been independently assessed.
"i
_ _ ~
/*
t f
Harold R. Denton 2
- 2. 0 DISCUSSION 2.1 MODEL IMPROVEMENTS 2.1.1 Improved Vertical Slip Model The two major areas of improvement include:
(1) a set of flow regime dependent relative velocity (slip) correlations based on the modified Bennett mass flux-versus-void fraction map and (2) an improved junction void fraction cal-culation based on the volumetric flux weighted void fraction rather than the previously used, volume-weighted void fraction. The new vertical slip model should improve calculation of the gravity-dominated vertical two phase flow and, consequently, prediction of the refill hydraulics in large-break LOCAs
~ ~
and the phase separation in small-break LOCAs. Guidance for the selection of junctions that should be modeled, using the new vertical slip model, will be provided in the user's manual (Ref. 1),
2.1.2 Thermal Nonecuilibrium Model An explicitly formulated nonequilibrium model was incorporated in RELAP4/M007 to improve prediction of phenomena associated with ECC injection by allowing subcooled liquid and saturated vapor to coexist within user-specified control volumes (nodes).
This model was introduced to facilitate calculation of the refill stage. The energy of subcooled water is tracked in the calculation, and the resulting interphasic heat transfer rate is calculated using the flow regime dependent constitutive package. A multiplier on the interphasic heat transfer rate is available to allow for scoping and sensitivity studies.
Guidance for the selection of the multiplier will be provided in the user's manual (Ref. 1). The nonequilibrium condition is restricted to the liquid only. Saturated and subcooled liquid streams are instantaneously mixed while any superheat in the vapor entering the control volume is assumed to increase i
the enthalpy of the subcooled liquid.
2.1.3 FRAP-T Fuel Model The FRAP-T fuel model has been included in the RELAP4/ MOD 7 code to provide a best estimate capability for fuel rod gap conductance, deformation character-istics, and material properties. The fuel rod temperaturec are calculated with the normal RELAP4 conduction and heat transfer models.
2.1. 4 Improved Reactor Kinetics and Decay Heat Model The reactor kinetics model was reprogrammed with some minor errors corrected.
In addition, the old ANS decay heat standard was replaced by the new ANS standard. An option is available to calculate the initial decay heat based on reactor operating history.
li 2.1.5 CHF Correlation J
l A new best-estimate critical heat flux (CHF) correlation was developed and I
incorporated in the code.
It evolved from a statistical correlation of a large data base, as a function of local conditions.
...y..._.-...
g ;..:.
e
(
Harold R. Denton 3
2.1.6 Liquid Entrainment Correlation The Steen-Wallis correlation, used in the RELAP4/M006, has limitations because it is dependent on the core inlet flow.
The modified Heat Transfer Research Institute (HTRI) correlation incorporated in RELAP4/ MOD 7 removes this limitation.
It correlates the entrainment fraction with the local superficial vapor velocity through the Kutateladze number and, therefore, may give significant entrainment even if the core inlet mass flow rate is zero or negative.
The correlation parameters were calibrated against the boiloff portion of the FLECHT-SEASET Test 35557. Using these parameters, RELAP4/M007 gave good '
I predictions of the reflood stage of the same experiment (Ref. 8). However, subsequent calculations performed in connection with the International Standard Problem #10 (PXL Test K9A) showed that the modified HTRI correlation gave poor -.-
-- _ results and the Steen-Wallis correlation provided more accurate predictions --
(Ref. 9). Hence, benefits of the modified HTRI correlation may be limited, to very low reflood velocities and the input parameters may be facility dependent.
Reference 1 will provide guidance to the user for selection of input parameters.
2.1.7 Critical Flow Stacnation Pressure Model The critical flow models in RELAP4 are based on stagnation pressure rather than static pressure.
Previous versions of RELAP4 did contain an option for calculation of the stagnation pressure. Mcwever, when the ficw areas upstream of and at the choking plane were approximately equal, supersonic velocities in the upstream control volumes were sometimes calculated, especially when using the nonequilibrium-based critical flow models, yielding erroneous results and oscillations. To avoid this problem in RELAP4/M007, the upstream-volume flow area is updated to values large enough to prevent supersonic velocities while maintaining the correct friction losses.
The model breaks down at low pressures, such as 50 psia.
I 2.1.8 Other Model Improvements These improvements include:
more accurate steam tables and steam-table interpola-tion methods, minor corrections to the centrifugal pump model, an improved
[
junction specific volume calculation, addition of a fill model where mass flux i
can be specified as a function of the mixture level in an adjacent volume, heat transfer model improvements, and various minor model improvements and error corrections.
'f
(
2.2 RELAP4/M007 FEATURES i
Several new features have been developed for dELAP4/M007 to provide more convenience to users, to save computer time, and to improve deficiencies identified in References 6 and 7.
A description of the major new features lj developed for RELAP4/M007 is given below:
e i
2.2.1 PWR Self-Initialization c
i I
The initial conditions in a RELAP4/M007 input model are overspecified in order to maintain the generality and usefulness of the program. However, the a
I..
.m j
~
w I
Harold R. Denton 4
overspecification may lead to inconsistencies that could, in turn, lead to erroneous initial perturbations or incorrect values being used.
In previous versions of RELAP4, eliminating the input model inconsistencies (balancing) was accomplishe.d through time-consuming iterative computer calculations.
The PWR self-initialization feature automates this procedure, balancing the energy by requiring the energy produced in the core and pumps to be equal to the energy transferred out of the steam generators, and balancing the pressure distribution by requiring the volume-to-volume pressure difference to be equal to the connecting junction pressure drop due to friction, momentum flux, elevation effects, and pumps.
The restrictions on using the energy balance (which may be applied with or without the pressure balance) are:
1.
An active core, pump, and steam generator must be considered,
' 2.T ~ The cie inlet temperature mu's't be 'specified, and'
~
~~
~
~
~
3.
The steam generator secondary side must be modeled as a single control volume, with the mixture level and mixture quality specified.
The restrictions on the pressure belance are:
1.
At least one flowing icop with an active pump must be specified.
2.
The pressure in one control volume in the locp must be specified as a reference.
3.
The mass flow rate in any parallel flow path (for example, core and
- - ~ ~
bypass flows) is specified as a reference.
Pressure balance is accomplished by adjusting control volume pressures and path friction loss coeffi-cients.
The NR self-initialization feature does calculate consistent initial conditions and at a very reasonable cost (typically in less than 20 CPU seconds on a CDC i
Cyber 176 series computer).
l 2.2.2 - Automatic Renodalization In previous versions of RELAP4, analysis of a complete LOCA required two separate calculations - a blowdown and refill analysis, and a reflood analysis. This was one of the major deficiencies pointed out in References 6 and 7 Data ii from the blowdown and refill analysis had to be hand-transferred to the reflood analysis, which caused all transient effects to be lost. The automatic renodali-I; zation feature has been developed to automate the transfer of data between
! j different nodalizations of the same problem, while maintaining transient i
effects. This feature requires that two (or more) nodalizations be prepared t
initially and transfer or mapping functions be constructed for the problem.
i The primary quantities conserved in mapping are mass and spatial volume, l
junction flow and then form losses, internal energy of the fluid and the
~I l'
energy stored in heat slabs. The input required with the renodalization is complex and must be checked by running a test case to ensure all data are being correctly transferred and all required input is present.
Renodalization
'l
Harold R. Denton 5
may be activated at any time during the transient.
For a LOCA analysis, renodalization is usually activated either during the refill period or at the beginning of reflood.
2.2.3 Decoupled Heat Transfer Time Step The time constant for significant changes in heat transfer conditions is often larger than the time constant 'or significant changes in hydraulic conditions.
For this reason, the heat transfer time advancement has been decoupled from t.
the hydraulic time advancement. A heat transfer time step control, based on temperature changes, calculates a time step bounded by a user-specified maximum time step and the hydraulic time step size.
This feature may result in a large saving in computer time as fewer calls to the relatively time-consuming
_ heat transfer routines are made. _ _--
2.2.4 SI Units Ootion
~
An SI units option has been incorporated.
Either British (" Engineering") or SI units may be selected, for either input or output.
Internal calculations in the code and the plot-restart tape are still made with the British units.
2.2.5 Other RELAP4/MCD7 Features Several other features and error corrections have been included. These involve water packing, mass depletion and mixture level crossing models, the generalized input feature which allows nonconsecutive numbering of volumes, functions and heat slabs, conversion of the RELAP4/M007 source to CDC update format to ease updating the program and to maintain better traceability of changes, and the restructuring of subroutines to improve computer efficiency.
3.0 RESULTS 3.1 CHECKOUT (DEVELOPMENTAL ASSESSMENT) CALCULAlIONS Several check problems were analyzed to ensure that the models and features developed for RELAP4/M007 were programmed correctly and that the new models I
and features did not introduce errors in the existing models.
The check i
calculations listed below do not represent a thorough assessment of the program, F
and no major sensitivity studies were performed.
1.
Simulation of Refill-Creare Tests 2.
Calculation of the Semiscale Mod-1 System Blowdown Test 5-02-9 f
3.
Calculation of the Semiscale Mod-1 System Integral LOCA Test S-06-3 i
4.
Calculation of the LOFT Integral LOCA Test L2-3 i
I 5.
Best-Estimate Prediction of LOCA in the Zion I Plant The results of these calculations indicate that incorporation of the new models improved predictions of 1) pressure transient during the bicwdown stage
,E:
(
\\
Harold R. Denton 6
and 2) downcomer penetration during the refill stage.
The calculations also demonstrate that the code is able to perform an integral, one pass, calculation from blowdown through reflood using two different nodalizations. More details on the results of developmental assessment calculations are presented in Appendix A.
4.0 RECOMMENDATIONS i
4.1 APPLICATION TO REGULATORY NEEDS The RELAP4/M007 code has been designed as a best-estimate large-break LOCA analysis tool.
The code more accurately calculates the blowdown and refill stages than does the RELAP4/M006 code and provides for an improved quantitative appraisal of ECCS which had been requested in References 4 and 5.
The calcula.__
~ ~ ~ ~
-tional capability of the reflood stage is almost the same as in the RELAP4/M006
~ - ~
s code.
The RELAP4/M007 code is user convenient due to its self-initialization and automatic-renodalization features.
Running time improvements were made by decoupling the heat transfer tie:e step control. With the self-initialization feature, the user spends less time in preparation of the input decks. Automatic renodalization feature permits integral LOCA calculations, from beginning of bicwdown through end of reflood.
The code is suitable to perform a statistical LOCA study as requested in Reference 3.
The RELAP4/M007 code can be usad for the following best-estimate calculations ~~
~~
associated with a large-break PWR LOCA:
BLOWDOWN:
In general, c:ocd temperature calculations will be obtained in the lower and midcore regions, including the het spot. As with RELAP4/M006, clad temperature predictions will be high in the upper part of the core and low in the lower part of the core, Most other parameters are expected
[
to be calculated reasonably well.
REFILL: Addition of nonequilibrium and vertical slip models substan-r tially improved the calculational capability of the refill stage. However, I
this conclusion has been obtained from limited developmental assessment.
I REFLOOD: Acceptable results can be obtained if the user selects appropriate parameters and options. At present, this selection is an art and detailed j
guidelines will be provided (Ref. 1).
The RELAP4/M007 code contains many models important in small-break LOCA analysis.
However, the limited number of analyses performed thus far have indicated both
}
successes and failures. For example, the code could not correctly predict the core uncovery; that is, mass inventory in the vessel, in the Semiscale S-SB-P1 test. However, assumptions made to simulate heat losses may play an important role in predicting the inventory distribution. The invent ~ory dist.ribution is T
[
also dependent on horizontal-flov~ stratification in the hot legs.
In predicting g
'_ Semiscale S-SB-P1 test data, the flow conditions in the hot legs were assumed to be homogeneous.
Scoping studies were performed to assess the existing f
i
~
(
Harold R. Denton 7
horizontal flow stratification model.
These studies, described in Reference 10, show that the horizontal flow stratification model in the code is not adequate to model physical phenomena, such as reflux boiling and critical discharge through small breaks associated with large horizontal pipes described in Reference 11.
The code cannot account for the presence of a noncondensible The control system model considered in the code may be too simplistic in gas.
The reactor point kinetics and reactivity feedback models in some cases.
RELAP4/M007 may be too simplistic for accurate calculation of transients without early reactor scram.
These difficulties do not preclude using RELAP4/ MOD 7 for small-break LOCAs or transients, but they will require greater care when applying the program to areas for which it was not specifically designed.
4.2 FURTHER IMPROVEMENTS AND CHECXOUT, IN PROGRESS
- _ :. - _ :^.-- - ~=--
- - ~
-In Reference 12, NRR has requested that additional capabilities be incorporated for analyses of system transients and steam-line-break accicents.
These capabilities include:
.i
'a.
Boron injection and tracking capability and inclusion of boron reactivity effects b.
Self-initiali::stion ca; ability for the secondary side of PWR plants.
Self-initialization capability for BWR plants.
c.
Boron injection and tracking capability is incorporated in a new version of the code. The remaining capabilities will be incorporated in a future version.
The RES technical contact for this work is Fuat Odar.
t Robert B. Minogue, Director Office of Nuclear Regulatory Research i
m' i
i
. i ll
.-m--.
- ,, f..
.m
(
REFERENCES 1.
RELAP4/M007 - A Best Estimate Computer Program to Calculate Thermal and Hydraulic Phenomena in a Nuclear Reactor or Related System (to be issued, draft copy available in RES office).
2.
RELAP4/M006 - A Computer Program for Transient Thermal-Hydraulic Analysis of Nuclear Reactors and Related Systems, User's Manual, EG&G Idaho, Inc.
Technical Report CDAP-TR-003, January 1978.
3.
Case, E. G., "NRR Requirements for Loss-of-Coolant Acc* dent Analysis Computer Programs (RR-NRR-77-5)," NRC Memo to S. Levine, June 23, 1977.
_ _ 4.
Bender, M., " Status of Generic Items Relating to Light ' Water Reactors:--- -
~ ~ ~ ~
~ Report NoJ'5," Letter from Chairman, ACRS, to M. A. Rowden, Chairman, NRC, Feb. 24, 1977.
5.
"NRC Program for the Resolution of Generic Issues Related to Nuclear Power Plants," Report to Congress, NUREG-0410, Jan.1,1978.
t 6.
Research Information Letter, 39, RELAP4/M006, November 1978.
7.
Research Information Letter, 90, RELAP4/M006 Assessment, May 22, 1980.
8.
D. M. Ogden, Data Comparison of FLECHT-SEASET Boiloff Test 35557 Using RELAP4/M007, EGG-CAAP-5152, May 1980.
9.
M. Hsu, F sttest Analysis of International Standard Problem 10 Using RELAP4/M007, EGG-CAAP-5210, July 1980.
10.
C. Tsai, "RELAP4/M007 Horizontal Stratification Scoping Study," WR-CD-028 g
(to be issued, draft copy available in RES office).
11.
N. Zuber, " Problem in Modeling of Small Break LOCA," NUREG-0724, August 1980.
- 12. Memo from Z. R. Rosztoczy to S. Fabic,
Subject:
RELAP4/M007 Modeling Modifications for Licensing Needs, Feb. 22, 1980.
I i
j s
i I
L..
<emo
~
(
iI APPENDIX-A l
3.1.1 Simulation of Refill-Creare Tests ii These calculations involved simulation of the 1/15-scale Creare Countercurrent Flow (CCF) test system.
A specified subcooled (ATsub=110*F) liquid injection rate was used in conjunction with the nonequilibrium model for a series of calculations varying the injected steam flow.
The flow regime dependent vertical slip model was applied at vertical junctions. The following results were obtained:
1.
With a steam injection rate of 0.6 lbm/sec or larger (the nondimensional gas flow, J* > 0.24) ECC water did not reach the lower plenum.
2FWith a steam injection rate'of 0.24 lbm/sec-~or smaller (J*c 0.14) al1 ~-~ ~ ~--- '
~
of the injected ECC water was delivered to the lower plenum.
3.
With an injected steam flow rate between 0.24 lbm/sec and 0.6 lbm/sec (0.14 < J < 0.24) relatively stable flow rates were calculated and gc partial delivery of the injected ECC water to the icwer plenum occurred.
Two sets of calculations were performed.
In the first set, the new vertical slip correlation was used, while in the second set, the old churn-turbulent slip correlation available in previous versions of RELAP4 was employed.
All other inputs were identical. Both calculation sets utilized nonequilibrium models in the same volumes involving the cold leg, downcomer, and the lower plenum.
Results of these calculations are compared with the Creare measured data in i
Figure 1.
These results indicate that the flow regime-dependent vertical slip l!
correlations give a better representation of ECC bypass and delivery than the l
churn-turbulent correlation. Hence, it is expected that the RELAP4/M007 code will calculate the downcomer penetration and refill stage more accurately than the RELAP4/M006 code.
However, since the downcomer penetration involves also other processes (for example, multidimensional flow, detailed nonequilibrium mixing and condensation, turbulence, and flow regimes that are more complex than those the code can consider) the above check calculations do not examine
[
the validity of the slip treatment, per se.
i.
3.1.2 Calculation of Semiscale Mod-1 System Blowdown Test 5-02-9 ii In this case, the focus was on checking the applicability of new models and to ensure that the existing models were not accidentally and erroneously changed. _
t The results of two comparisons with test data are presented in Figures 2 and 3.
Figure 2 is a plot of the RELAP4/M007 calculated upper plenum pressure, compared with both the test data and with the RELAP4/M006 calculation results.
There is an excellent agreement until 25 seconds into the transient. The ll RELAP4/M006 calculation slightly underpredicts the pressure due to nonequilibrium l'
effects after 25 seconds when ECC injection begins.
RELAP4/M007 accurately predicts the pressure after 25 seconds, indicating the improvements introduced I
)
t.i
.. L
.(
A-2 with the nonequilibrium model.
Figure 3 shows comparisons of clad surface temperatures at 11 and 16 inch levels.
Both the RELAP4/M006 and RELAP4/ MOD 7 calculations are within the measured data spread. The somewhat higher temperature calculated by the RELAP4/M007 code is caused by a slightly different system mass distribution that produced a smaller core inlet flow.
In general, the RELAP4/M007 overpredicts the surface temperatures in the upper part of the core, underpredicts in the lower part of the core and there is a good aqreement in the central part of the core.
These are the same conclusions reached for the RELAP4/M006 code (Refs. 6 and 7).
This very limited information shows that the RELAP4/ MOD 7 code accurately calculates the blowdown stage of the S-02-9 test. Note that the calculations of the blowdown period have also been quite satisfactory using the RELAP4/M006 code. However, it can be seen that the accuracy of calculations of pressure.._ _ -
- and core inlet flow rate is improved by addition of the new vertical slip correlation and the nonequilibrium model in the RELAP4/M007 code.
This should lead to more accurate predictions of clad temperatures.
3.1.3 Semiscale Mod-1 System Intecral LOCA Test S-06-3 The primary purpose of this checkout was to demonstrate the capability of RELAP4/M007 to perform an integral calculation from blowdcwn through reflood.
The Semiscale S-06-3 Test was a simulated 200 percent double-ended cold leg break with ECC injection in the intact loop.
The input model for the calculation
~~~
described herein was created from two different nodalizations (one for blowdown /~
refill and one for reflood) used in the previous, RELAP4/M006 analysis of the same test. The calculations were switched from the blowdown / refill nodalization j
to the reflood nodalization during the refill.
For reflood calculations the g
Steen-Wallis correlation was used because the new modified HTRI correlation gave very unrealistic results.
Therefore, the Steen-Wallis correlation was used in all developmental assessment calculations.
I The RELAP4/ MOD 7 calculation was run to 111 seconds of transient time, at which time all low power rods were quenched. A comparison of the calculated and the measured upper plenum pressures is shown in Figure 4.
The agreement is fairly
't good during the blowdown period.
However, after the onset of ECC injection, I
at about 18 seconds, the RELAP4/M006 code predicted depressurizaticn due to
!l its thermal equilibrium assumption, while the RELAP4/H007 calculation, using
~
the nonequilibrium model, showed a much smaller deviation from the data.
i Figures 5 and 6 show comparisons of calculated vs. measured low power rod surface temperatures at low and upper-middle elevations cf the core, respectively.
These low powered rods (32) are representative of the core average condition.
For low elevations, the RELAP4/ MOD 7 results are, in general, in better agreement with the measured experimental data than the RELAP4/M006 results. The calculation g
of the quench time is also better. At the upper-middle elevation, the RELAP4/M007 l
calculation indicated the heat slab quenched at 110 seconds, while the measured
- 4 quench time is 160 seconds. For high elevations.(not shown), the predicted quenching times with RELAP4/M006, and RELAP4/ MOD 7 were 260 and 80 seconds, respectively, while the measured quench time was 125 seconds.
L
g....
(
A-3 For the four high powered rods which were represented by a hot channel, both RELAP4/M007 and M006 underpredicted the peak clad temperatures in the lower region of the core by about the same amount.
In the upper and upper-middle parts of the core, both codes overpredicted the data, the RELAP4/M007 predictions I
of the' peak clad temperature being higher in both parts of the core.
- Hence, the predictions of peak clad temperatures by the RELAP4/M007 in the hot channel for this test do not represent any improvement over those predicted by the RELAP4/M006.
RELAP4/M006 predicted a minor rewet for all elevations at about 10 seconds while RELAP4/M007 calculated this type of rewet only at the low elevation.
Test data do not indicate c rewet at this time. Mcwever, at the lower portion of the hot channel, the data indicate a sudden quench at about 95 seconds L_
which both codes failed to predict.. _. _. _ _ _
-~
3.1.4 LOFT Integral LOCA Test L2-3 This case was sek.ed to test incorporation of the FRAP-T fuel model and to check the automatic renodalization feature.
The calculation was continued t
until the final core quench which occurred approximately 55 seconds after experiment initiation.
The nodalizations for the blowdown / refill and reflood periods were originally t
developed by the LOFT Program for L2-3 experiment pretest predictions and did not contain the modeling of the hot channel. As with the cther developmental assessment tests, Steen-Wallis correlation was used for the reflood calculations.
Some comparisons of the calculated vs. measured results are shown in Figures 7 through 9.
Figure 7 pertains to the upper plenum pressure where the difference between the calculated and experimental results is traced to the difference in i
predicted vs. experimental break flows. The fuel rod cladding temperatures for the middle and upper 1/3 of the core are shown in Figures 8 and 9.
The calculated temperatures are compared with measurements for the peripheral bundles where the radial power peaking factor was about equal to one.
In general, the calculation underpredicts the blowdown peak clad temperatures, although the time for CHF is predicted reasonably well, and does not reproduce the total effect of core-wide rewet at about 8 seconds.
This calculation showed that the new models and features, particularly the automatic renodalization feature, the FRAP-T fuel model, and the nonequili-brium model were successfully incorporated to the code.
'I 3.1.5 Best Estimate Prediction of LOCA in 7 ion I Plant I
i The Zion I LOCA (large cold leg break) calculation was performed in the previous l
e BE/EM study (Ref. A-1) using RELAP4/M005. The system nodalization developed for that study was maintained, as nearly as possible, for the current investiga-tion; however, RELAP4/M007 analytical models and features were selected for
{,.
this calculation, where appropriate.
Steen-Wallis entrainment correlation was used in the calculations of the reflood stage.
1 e
d.
.-..a-..
v
(
A-4 t
The average core red surface temperatures from.both calculations are presented in Figure 10.
The large differences between the temperatures are due to the positive core inlet flow prediction by RELAP4/M007 between 4 and 11 seconds.
This difference in prediction is primarily due to different break flow models i
and to a lesser extent to new vertical slip and nonequilibrium models.
These calculations also indicate that the results are consistent with LOFT large-break loss-of-coolant experiment (LOCE) tests and that the code can perform integral LOCA calculations.
Further information on Zion I calculations is presented in Reference A-2.
References A-1. G. W. Johnsen, F. W. Childs, J. M. Broughton, A Cceparison of Best Estimate and Evaluation Model LOCA Calculatters: The BE/EM Study, EG&G Idaho,
--~ ~
-~
Inc., Technical Report PG-R76-009, December 1976.
A-2 T. L. DeYoung, "RELAP4/M007 Developmental checkout:
Zion I" (to be issued, draf t copy available in RES office).
l.
I I
I l
l I
f l
I i
- ~ ~ ~
h
=
L
(....;....
\\-
il l'
2 4
l-1
\\
f 1
i 3
04 o
I I
I J
k D CREARE DATA POINTS A RELAP4/ MOD 7 NEW BLIP M
X Retap4/MOoz n
stIP oto i
U j
i n
J i
i a
)
2 4
5 i
O 0.2 j
4 m
m x
g Z
A' i
It!
1 E
4 H
o X
4 Z
0 li Z
X i
1 I
I I
l 0
0.0 0.00 0.05 0.10 0.l15 0.20 t
- i.
JFD* (H0HDIHENSIDHAL LIQUID FLOH)
SLIP UELOCITY CORRELATION'!
RELAP4/H007 CHECK 0UT j
Figure 1 i
j' i
g**
6A..
W
(
n' Ce I
i l
l I
v i
)1L n
M g
'i l a
l one I
q
+oo M
2 3
>2X l
m ANN
/
m lwTT U
AL L
c.
H,J J 10 Ci g
<~J N
I
=
Q C: C::
of c
o 2
)
a n
s a
e 80 o> u) h u
q.*
y 8
o w
a N
w
_2 I
rc
- o i
P ~ ~~
I-(O cy 5
l n
~
o>
m b.
- u.
ii o
i 1
n I
l I
i o
O O
O c3 o
o o
g C
C O
C.
o c.
t]
O u~
C to o
to i
o4 m'
o r-u m
p IS (Yd) SEOS63Hd l
'~"
~
~'
)
I l
l 1000 L
1[
\\ l-1 t
900 i
4~
m :
l 3"
3-I' a
~
25 800 l4 j
g I-l 1
P l
M
]
W 700 2
5 l
i r
i O
I' l
600 l-Heasu' red. Data Til-F4-14G g
2-fleasured Data Til-C3-13U 3-fteasured Data Til-E8-14K i
4-RELAP 4/fl006 i
5-RELAP 4/ MOD 7 j
500 l.
j I,.
ji 400
- i A
a A
O.0 10.0 20.0 30.0 40.0 i
TIME AFTER RUPTURE (S)
Figure 3.
Average Channel Clad TemperatureJ.11-16 inch level i
1l
'~
I
- - - - ~ ~ -
,t i
1 I
i l:
li 12000 I
I l
l l
l 1
I I
l 10000 1
2 8000
.w I
- DATA (PV-UP+10) 6000
RELAP4/ MOD 7 0
RELAP4/ MOD 6 i
=
b 1
a.
fi 4
j 4000 1
I l
2000 i
,'53hM000 M '
I I
I I
I o
j D
5 10 15 20 25 30 35 I 40 45 50 4
i TIME (S)
I, SEMISCALE S-06-3 i
i t.
Figure 4 - Upper Plenum Pressure 9
w,..,,,,..,
..m.,
i & ?'
k l..E.... _.
L _ __ _.
(
l g
e f' W
i i
aa a
.g
.~
i 1.'.
oo e s I
eNNe 8'i' i
- u. u. c c h
e pfgg j
.J.
v co ss c
AA 3
--vt s
i
)
e q
mozz 4
l d
N to- -
G
( J i
e.
8 oc i
-Njo f
c-m 57 o
i c
t to 3.
~
n't o
=o
^ s nm mm e-y v
Ea o
w 25 4
e wa s
j 3
3 E C o5
.-* O a F-(
sm I
~.
-g_
r
=
c a m
h<,. (f w
e w
(
g i
i s
\\
l t
( J t
.t I
8 i
i I
J eM
.t
- o cu t-I i
I N
I i
.o
)
w t
i 4
00 SHA1YH3dH31 I
1 i
h
r:
I i
i i
'r, i
\\,
l I
i i
1400 g
g g
g i
g l
[
4 1200 j
,t
~,
/
is i
l k
o i
t I
s
+
i j
s
't l
I L
1000
(
.i j
d 3
's t
I g _
e"%. - a s
s t
l N
t v
l h
800 J
~
h-i y
4 i
to 4
1 i
600 i
DATA (TH-D8-28)
RELAP4/ MOD 7 I
O RELAP4/ MOD 6 i
n 400 8
i i
t
(
i I
I I
I I
I
.I I
200 k
6 20 40 60
'80 100 120 j140 160 180 TIME (S)
SENISCALE S-06-3 l
l Figure 6 - Low Power Rod Surface Temperature, l
(i 24-30 in. Elevation s
S
l
- ~
~ ~'
O.
I s
l I
17500 i '
I I
I I
I i
o 15000
'4 l
'. 1 12500 I
I l
2 10000 S
i
- DATA (PE-lUP-00l A)
~
u,
RELAP4/ MOD 7 0!
7500 e
t l
s 5000
's U
i s
s 2500 I
i I
I 0
-~~I O
10 20 30 40 50 60 I
t
- t..
TIME (S)
Figure 7-CalculatedandExperimenkal Upper Plenum Pressures
{;
er--
4
s.
i I
l e
1 I
I I
t.
o t "%,.
I 1 n
P e
\\
9
\\
_- 4 _'
l p* {
l l
~
JA
{
e E-x n_
._ p rs J,
g__
\\
=c N
E" t
)
J t
C 2 IE
\\
o s_
t fr3 g
"e "as o
58,
n
=
m xx~
w Bu i
l'
/
LAQ N
TH v
t
~~o n
J S en%
8 2
- 222, N
w ec g
(
'= w w w
--a
[o 8@g g
WAAW k
u o
-mu g
nn<
>a e--
I abMg G
uur emeW N
I aco>
7 o,.e e 4 m
r i
s-ee E
u.m m s n
m g
i me-o S
p a: eo g
wwwz t
P H t-N g
wwwT g
444k f
s.
H >= I-=
3 I
/
_f 444W QQQE a
,,y t-n
-Ne g~~
.e I
i 1
a 1
.O o
o o
o o
o k
i o
o o
o o
o r
o.
co N
to n
e g.
t F L P'
(N) 3E01YB3dH21 t
A 9., -,,
..__.__-.,.,_,._,,.9y._p,.. _
m g,,,.
,p.,,.,,,
__,_g..
L
(
e i
I I
4 n
,e s i
,b 1
l-n
'< J e
\\{-
e oE
- 8 u, _. _
g
-e i
e TR a
t v
t 8-i fb t
oa t
e NS m
w
^M v5 8
f JJ a.
m 5,
i
- m. c N
v e
c
-o n
J m'
g I
.CCNEl N
w' >
- 3
~W c
an g
W a. W u
-u vu o
~ o.
=o 4
H u c.
g
--em e
';;i B WNT N
u l-
- C TSO
\\ [N-0084 E4"-
muxP n
e i
mmaa
=
i eo en WWWE' l
FFHN
---T r
i
- 444k,
- e i
FFH p
444 COGE
't
,l n
f
-Ne e I g s.
e e
i'i 3
l I
h.
i I
o o
e o
o i
O O
G G
O r-ve n
v i !
i (X1 3HO1YE3dM31 i
i I
t
(
n 8
- l.. - - -.. _ _ _ _.. _ _ _ _,, - _ _ _ _..... _,. _ _ _ _., _, _ _ _ _, _, _,,,,,,,, _ _, _ _ _ _, _, _,, _ _. _ _ _ _ _ _ _ _ _ _ _ _,,
T' 'h
.7
-,.3.....
i l
l j !,
l i
t.
1100 I
I l
l l
g i
i 1
ZION / MOD 7 2
BE/EM see.
l n
^
A
- ow m
!9 w
cr 3,.
a H<
aw O-764 rw F
3 Gee s
I I
I I
i m
8 5
to 15 26 25 3e 35 f
TIME AFTER RUPTURE CSEC) l l
,i l
e.
Figure 10 - Average-core-rod surface temperatures - center of core.
i
_