ML20137W130
ML20137W130 | |
Person / Time | |
---|---|
Site: | Vermont Yankee File:NorthStar Vermont Yankee icon.png |
Issue date: | 03/31/1997 |
From: | Fujita N, Kapitz D, Lefrancois M YANKEE ATOMIC ELECTRIC CO. |
To: | |
Shared Package | |
ML20137W107 | List: |
References | |
YAEC-1339, NUDOCS 9704180039 | |
Download: ML20137W130 (68) | |
Text
. - - _ _ _ _ - _ _ _ - _---
b.
[:
a 0
[: 6 i
Application of the FIBWR2 Core Hydraulics Code
[.[
to BWR Reload Analysis
, -1 E
l h January 31,1997 Major Contributors: M. P. LeFrancois N. Fujita D. M. Kapitz E
E F
I u
h: 9704190039 970414 N -
PDR ADOCK 05000271 ~ % i P pop
)
( Prepared By: N6- ~
/YMCOso //J/ 97 .
'M. P. LeFranco'i's ~ '(Dite)
Transient Analysis Group
\
\,
Prepared By: bA ' F/
/ 7/ / 7 i D. M. Kapitz (Date) '
Transient Analysis Group Prepared By: k N. Fujita / (Date)
{ Transient Analysis Group
(
Approved By: 8/#I A+t"N M I d&#M / 3/J97 P. A. Bergeron / 6' '(Datt)
Transient . Analysis Group Approved By: tb -
/ 8/ '
2 i
. R. Chapma% / (Dde)
Nuclear Engineering Depanment Yankee Atomic Electric Company Nuclear Services Division 580 Main Street k Bolton, Massachusetts 01740 b
[
ii
4: DISCLAIMER OF RESPONSIBILITY h .
-- This document was prepared by Yankee ' Atomic Electric Company (" Yankee"). The use of 1
i 1 information contained in this document by anyone other than Yankee, or the Organization for which l g(-
L this document was prepared under contract, is not authorized, and with respect to any unauthorized )
ya neither Yankee nor its officers, directors, agents, or employees assume any obligation, ]
- . responsibility, or liability or make any warranty or representation as to the accuracy or completeness of the material contained in this document.
b 1
[
[ 4 b
I I
i
[
[i c;
[
L g
E F _ _ - - - _ - - _
ABSTRACT The benchmark and validation effort for the transient critical power ratio methods using FIBWR2 computer code, written specifically for Boiling Water Reactor steady state and transient thermal hydraulic analysis, is summarized in this report. FIBWR2 is an upgraded and enhanced version of the FIBWR code, enhancing the steady state portions to accommodate new fuel designs with part length fuel rods and enhanced water rod designs and upgrading to include transient capability. The benchmarks for the transient simulations result in excellent agreement to accepted YAEC methodology, GE ATLAS test facility critical heat flux data, and exact analytic solutions.
~
b t
[ !
[ iv
( .
ACKNOWLEDGEMENTS
, The authors express appreciation to Mr. Richard Weader (Yankee Atomic), Mr. Steven Ganthner (Onsite Engineering), Mr. Ben Gitnick (Scientech Inc.) and Mr. Nick Newman (Scientech Inc.).
i Each contributed to the significant engineering effort for the validation of the FIBWR2 computer code and analysis methodology summarized in this report.
l V
? _
TABLE OF CONTENTS fD.81 H AB STRACT . . . . .. . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iv ACKNOWLEDGEMENTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Y h TABLE ' OF CONTENTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . vi LIST OF TAB LES ' . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . viii E
. LI ST . OF FIG UR ES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ix 1.0 INTR ODUCTION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
- 1.1 Overview of Transient CPR Analysis Methodology . . . . . . . . . . . . . . . . . . . . . . . . . 3 1.2 - FIBWR2 Background .......................................... 4 1.2.2 FIBWR2 D,evelguent and Verification Objectives . . . . . . . . . . . . . . . . . 4 1.2.3 ' Approval Status of FIBWR and RETRAN Hot Channel Methodology . . . . . 6 1.3 . Expected Uses of FIBWR2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 1.4 Organization of Remainder of. Report . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 ,
1 2.0 DESCRIITION OF FIBWR2 CODE ENHANCEMENTS . . . . . . . . . . . . . . . . . . . . . . . . 12 2.1. FIBWR2 Conservation Equations . . . . . . . . . . . -. . . . . . . . . . . . . . . . . . . . . . . . . 12
{.
2.2 FIBWR2 Core Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 2.3 FIBWR2 Solution Scheme . . . . . . . . . . . . . . .......................... 14 i
{ ' 3.0 COMPARISON OF FIBWR2 WITH ANALYTIC SOLUTIONS ................... 20 3.1 Initial Conditions of Transient Two Phase Problems . . . . . . . . . . . . . . . . . . . . . . . 20 3.2 Comparison of the Steady-State Void Profile . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 3.3 Description of Analytical Solution Technique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 3.4 Flow Decay Tran sient . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 3.5 Oscillatory Flow Transient . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 3.6 Summary of Results for the Analytic Comparisons . . . . . . . . . . . . . . . . . . . . . . . . 28 a
4.0 COMPARISON OF FIBWR2 TO CRITICAL HEAT FLUX DATA . . . . . . . . . . . . . . . . . 35 4.1 Simulation of Steady State and Transiat P.i Rod Tests . . . . . . . . . . . . . . . . . . . . . 35 4.2 Sensitivity Studies . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37 4.3 Summary of 16 Rod Test Simulations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 5.0 FIBWR2 TRANSIENT HOT CHANNEL PERFORMANCE . . . . . . . . . . . . . . . . . . . . . . 42
. 5.1 Hot Channel Methodology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42 5.2 Hot Channel Re sults . . . . . . . . . . . . . .' . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43 I
L
TABLE OF CONTENTS (Continued)
[.
z.
5.2.1 : Comparison of Thermal Hydraulic Parameters . . . . . . . . . . . . . . . . . . . . 44 5.2.2 Comparison of Calculated CPR Values . . . . . . . . . . . . . . . . . . . . . . . . . . 46
. 6.0 CONCLU S ION S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
( - -
7.0 REFERENCES
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
[
(
(-
(
r 4-
{:
1
(- 1 1
Vii
..m _ __ _ _ _ _ _ _ . _ _ _ _ _ _ _ - _ - - - - - _ - . _ - _ _ -
- q. '
LIST OF TABLES I
Numhcr Title fags.
1.1- - Enhancements to FIBWR to create FIBWR2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 l' 3.1 Initial Conditions Assumed for Transient Two Phase Problems . . . . . . . . . . . . . . . . . . . 29
- h. 5.1 Calculated CPR Values' for Pressurization Events . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47 5.2. ' ' Calculated CPR Values for Subcooling Events . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48 h
1 i
l
')
o L
[
[ viii e - _ _ - _ _ _ _ - - - - _ _ _ - - - - _ - - - _ .- - - - - - - -
g, I ' LIST OF FIGURES
! )
Number Title Eagg 1
~1-1.. Flow. Diagram of Existing Reload Methodology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 1 1-2 Flow Diagram of Reload Methodology with FIBWR2 . . . . . . . . . . . . . . . . . . . . . . . . . 10-1 Flow Paths in a BWR Core . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 l 201- Nodalization Diagram - . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18- !
2 FIBWR2 Compound Lower Unheated Node . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 3-1 Steady State Void Fraction Comparison . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 j'- l3-2 Flow Decay Transient - Void Fraction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31 a
- 3-3) . Flow Decay Transient ~- Exit Mode Flow Rates . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32 j 3-4 Flow Oscillation Transient - Void Fraction - 25 Nodes . -. . . . . . . . . . . . . . . . . . . . . . . 33' i 3-5 Flow Oscillation Transient - Void Fraction - 50 Nodes . . . . . . . . . . . . . . . . . . . . . . . . 34 l L4-1l
- FIBWR2 Steady State' 16 Rod CPR Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39 4-2: FIBWR2 Transient 16 Rod Boiling Transition Onset Results . . . . . . . . . . . . . . . . . . . . 40 43 ? ATLAS 16 Rod Test CPR Predictions - Sensitivity Study Results . . . . . . . . . . . . . . . . . 41 5-1 GLRWOBP Boundary Conditions Channel Power . . . . . '. . . . . . . . . . . . . . . . . . . . . . . 49-52- GLRWOBP Boundary Conditions - Upper Plenum Pressure . . . . . . . . . . . . . . . . . . . . . 50 5-3 GLRWOBP Results - Active Channel flow Rate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 5-4' ' GLWOBP Results - Channel Exit Void . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52 55 LOSC Boundary Conditions - Channel Power . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53 5 LOSC Boundary Conditions - Upper Plenum Pressure . . . . . . . . . . . . . . . . . . . . . . . . . 54 5 LOSC Results - Active Channel Flow Rate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55 L5 8- LOSC Results - Channel Exit Void Fraction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56 b l L
1 i
[ ix
1.0 INTRODUCTION
The calculation of Critical Power Ratios (CPR) reached during transient conditions for Boiling Water Reactors (BWRs) requires the application of thermal hydraulic models to evaluate the key parameters used in the applicable CPR correlation for the most limiting fuel assemblies (hot channels). ' The existing YAEC methodology, at approved for application to BWR reload analysis,
-utilizes the steady state FIBWR (Flow In Boiling Water Reactors) code to establish initial pressure drops and flow distributions for the transient calculation of thermal hydraulic conditions. The systems level transient thermal hydraulic conditions for the hot channel are then calculated using the l RETRAN code. The resultant CPR is calculated using a vendor proprietary correlation with the TCPYA01 code. In this methodology, the FIBWR code is used to accurately model the complex flow geometry of the BWR core under steady state conditions with its multiple parallel flow channels and significant water tube and leakage bypass flow paths.
L The FIBWR2 code is a new version of the FIBWR code designed to extend the BWR specific
)
t models'used in the original FIBWR code to the analysis of transient thermal hydraulic conditions.
L The structure of the FIBWR2 program also accommodates the installation of subroutines for the
' calculation of CPR within the FIBWR2 coding. These features allow the FIBWR2 code to replace two separate codes (RETRAN and TCPYA01) used for hot channel analysis. FIBWR2 also includes features which allow explicit modeling of advanced fuel designs with axially variable active flow r .
E areas due to past length fuel rods and alternative water tube designs.
l The purpose of this report is to describe the key models and demonstrate the applicability of i the FIBWR2 program to the analysis of transient thermal hydraulic conditions representative of BWR fuel and core geometric including the implementation of CPR correlation within the FIBWR2
[ program coding. ' Key features added to the original FIBWR code to enable evaluation of transient t
thermal hydraulic conditions are identified and the nodalization and solutions scheme summarized.
[ The performance of these FIBWR2 models is demonstrated by comparison of the results obtained using FIBWR2 for two problems with analytical solutions. The results of simulations of experimental test data obtained from a 16 rod test assembly configuration under both steady state and i transient conditions using FIBWR2_(with the internal calculation of CPR) is presented and compared with both the test results and a similar simulation performed with the RETRAN/rCPYA01 models.
f e
The performance of the FIBWR2 is then compared with qualified RETRAN base models for representative transient analysis performed for recent actual operating cycle conditions. The comparison of FIBWR2 and RETRAN based results for the reload transient analysis include minimum CPR values. For both the FIBWR2 and RETRAN based evaluations, the CPR is cuculated using the same vendor proprietary CPR correlation. Application of the correlation for the RETRAN based evaluations is achieved by use of the TCPYA01 code while for the FIBWR2 results, the correlation was incorporated within the FIBWR2 program.
l t
i I
2 I
I
t g -1.1 Overview of Transient CPR Analysis Methodology
- A flow diagram of the current YAEC transient analysis methodology is presented in -
Figure 1-1. The CPR methodology' includes three components involving analysis of thermal' f '
- j. hydraulic conditions, they are:
' Steady state thermal hydraulic conditions for the specific core loading design using the FIBWR thermal hydraulics code. FIBWR provides core flow and pressure drop-information'to the lattice physics code, the fuel performance analysis and initial conditions for the RETRAN systems level calculations.
l l
[.
Simulation of the overall plant system response to an Abnormal Operational Transient l
'(AOT) using the RETRAN system level model, including the simulation of neutronic
- feedback using one dimensional kinetics. This calculation provides time dependent normalized power and pressure values as input boundary conditions for the CPR
]
. methodology.
1 i
h -
Calculation of the transient change in CPR for the limiting fuel bundle using the p RETRAN/TCPYA01 code combination. FIBWR provides the initial power, active flow and bypass flow based on specific calculations for the fuel type analyzed.
l I-h The transient CPR calculation is referred to as the " hot channel" analysis since it simulates the response of the single highest power fuel assembly within the core loading to a plant transient.
The hot channel analysis provides a CPR operating limit for the transient, such that the CPR safety
[ _ limit for fuel cladding integrity (FCISL), the lowest allowable CPR defined by the fuel vendor, is not I
achievable should the event occur. The most restrictive operating limit is determined by analysis of j L a BWR's limiting AOT for a given reactor core design's fuel cycle.
L j Figure 1-2 shows the revised approach to the reload analysis using the FIBWR2 methodology I
to explicitly address the advanced fuel design. FIBWR2 replaces the FIBWR code in setting the
- initial steady state thermal hydraulic conditions for the lattice physics, fuel performance, systems level and hot channel calculations. This application of FIBWR2 is identical to that of FIBWR.
3 i
I Additionally, since FIBWR2 has transient thermal hydraulic analysis capability and contains the vendor CPR correlation it provides an alternate to the RETRAN/I'CPYA01 combination. It should be noted that a separate report [13] documents a summary of the qualifications carried out for the steady state portion of the use of FIBWR2.
1.2 FIBWR2 Backcround The FIBWR code was developed as an accurate, steady state multichannel core hydraulic simulator for BWRs. The FIBWR methodology was designed specifically for the complex flow geometry in BWR cores, with many parallel fuel channels and other small but significant flow paths, such as water tubes and leakage flows to the bypass. The two-phase flow and pressure drop models selected for inclusion in FIBWR were specifically chosen for their excellent agreement under BWR non-LOCA transient conditions. FIBWR has been distributed by Yankee Atomic and the Electric t
l Power Research Institute and widely used in the industry. FIBWR2 is a new version of FIBWR developed specifically to allow simulation of advanced fuel types with the same accuracy as FIBWR.
l It is also intended to enhance the preser.t hot channel analysis of CPR methodology during postulated l traasients through the application of the more detailed representation of the fuel and core flow paths found in the FIBWR, l
1.2.2 FIBWR2 Development and Verification Objectives l
In the years since the original development of the original FIBWR, BWR fuel designs have become more complex, and computer speed and memory have greatly increased. In 1990, work began on FIBWR2, with the following six main objectives:
l l -
To accommodate the new, more complex fuel designs and water tube types currently offered by BWR fuel vendors, To add transient capability, including simulation of transient variable axial power l shape, 4
y - -
To provide a more detailed representation of the transient core hydraulics for CPR calculations,' -l t 1 To simplify the transient CPR methodology by structuring the program to facilitate-implementation of CPR correlations in the code,.
To solve larger problems with more flow channels and more water tubes by improving the speed of the solution technique, I
L To improve code portability and maintainability by conversion to FORTRAN 77.
These modifications provide the capability to analyze new fuel designs and perform full core multi-channel thermal hydraulic analysis. FIBWR2 has capability to represent each fuel assembly in a BWR core discretely. A summary of the new features of FIBWR2 which achieve these objectives L are presented in Table 1.1.
I L . Although new capabilities were added and improvements were made to the original FIBWR coding and solution techniques, steady-state core flow, void and pressure drop distribution results L have not changed. The FIBWR2 development effort was required to use the same pressure drop and j Lvoid distribution models as the FIBWR code. These models were considered to be well validated I and were previously documented in the original FIBWR report [1,2]. Only minimal changes and additions to original FIBWR models were made to support transient calculation capability. In L
I Section 2.0, the original FIBWR models are briefly reviewed and the modifications which comprise J
FIBWR2 changes are documented.
[
)
i FIBWR2 to FIBWR cross-code verification calculations were performed. Small differences '
were noted in the models of the original FIBWR code, which were revised during the development p of FIBWR2. Additional' differences were the result of the enhanced solution technique.
, _ Modifications to enable the modeling of new fuel designs and water tube types were verified L' '
against steady-state analytical solutions. In addition, pressure drop measurements for several
- operating states of Vermont Yankee operation were used to validate steady-state FIBWR2 results. A i: 5 O-
validation effort carried out for the steady state portions of the FIBWR2 code has been documented sepasate of this report [13]. It should be noted that the steady state exercise of FIBWR2 is consistent j Lwith the current method using FIBWR.
Transient solutions were verified against exact analytical solutions. . Cross code comparisons y to transient the hot channel RETRAN model were used to verify the fuel rod temperatum calculation I and transient pressure drop / flow distribution solution technique. Cross code void profile comparisons i
rgainst RETRAN were also used to verify the implementation of the drift flux model in FIBWR2. -!
(
l:
The FIBWR2 code was validated for hot channel analysis by comparison to the published 16
(_.'
rod CHF Test Facility data [3,4,5,6].' These comparisons show that FIBWR2 with the vendor CPR correlation installed can adequately predict the onset of boiling transition. l 1.2.3 Approval Status of FIBWR and RETRAN Hot Channel Methodology L ~ The FIBWR[1,2] code was qualified by validation against analytical solutions and by comparison to measured plant data and other codes. Yankee received NRC approval of the FIBWR application submittal [7], and has used FIBWR as part of the NRC approved licensing methodology package since 1981. The RETRAN and TCPYA01 code combination was validated against Critical Heat Flux test facility results [3,4,5] and approved for the use [8] in establishing transient Minimum Critical Power Ratio based operating limits. RETRAN/TCPYA01 has been used for the past eleven
[ reload cycles for Vermont Yankee, i
1.3 Expected Uses of FIBWR2 g FIBWR2 will be used for thermal-hydraulic evaluations of BWR cores during steady state operation and abnormal operational transients (AOTs). The steady state evaluations will provide b input to lattice physics codes, fuel performance codes and transient codes. The steady state
~
application'of FIBWR2 is identical to the FIBWR application. The steady state application is r- 3 1
' docu[mented in'[13] The code will be used for transient CPR evaluations as described in this report.
FIBWR2 will, therefore, be used for steady state flow distribution analysis and for transient CPR-L evaluations.
( 6' r
[ FIBWR2 will not be used to analyze: .!
i
' Acoustic and/or Sonic Flow Phenomena, such as Choking or Water Hammer
{-
Horizontal Flow ,
h -
Significant Vapor Superheat (Exit Qualities Greater Than 1.0)
LOCA Blowdown /Reflood, ADS actuation, CCFL
[ -
Fuel Rod Heat Transfer Beyond Critical Power I Fuel / Clad Deformation or Melt .
[.
1.4 Organization of Remainder of Report i
Section 2 discusses the modifications to the FIBWR code to include transient capability.
Section 3 presents the results of the comparison of FIBWR2 to two analytical solutions, an exponential flow decay and a sinusoidal flow oscillation. Benchmarks of FIBWR2 to the General
-l
[ Electric ATLAS test facility critical heat flux experiments are presented in Section 4. Comparison to )
i the' approved YAEC transient CPR methodology is summarized in Section 5. A summary of the l r l conclusions for each of the benchmarks and qualification studies are reported in Sections 3-5 are l presented in Section 6. 1
[-
F L
I L
r L
f u
g.
L
[. 7 e
r
i TABLE 1.1 Enhancements to FIBWR to create FIBWR2
.)
1 Improved Geometric Flexibility -
- - ' Axially varying flow geometry- l Advanced water tubes / water cross features !
- Multiple fuel rod geometries / fuel rod propeo.ies per flow channel J Improved Numerics
- For convergence at low-flow / low-power conditions For transient multi-channel Critical Power Ratio (CPR)' calculations needed to determine compatibility of new fuel designs in a mixed core For transient heat flux and Linear Heat Generation Rate (LHGR) edits !
l Automated Linkage to Physics Codes
- . 3 D power readings from nuclear simuletor files j) t Automated Linkage to RETRAN
- l. -
Boundary conditions for transients e
single (average) channel edits l.
All FORTRAN 77 for Ease of Maintenance l
-l 1
l- l 1
l L
[.
8
' m.
M
SWTICII Intual
[ Simulate to 3-D Nwer Dist. IAPURX l ---> STABIIIIY FIBWR Power / : I.C.'s (Distribution T/H ir MBWR ,r SIMUIATE 3r Kinetics, Fuel Rod inlet Steady Bypass 3-D Nodal Power Dist.
'Ihermal < State g, g y,g, - IAPURW Effects T-H Initial 3-D SUCKir STABIIIIY Point EXCLUSION
> Kinetics REGION Generation ANALYSIS Point N ir Bypass Mow Core Avg. - ASME hgap S l
M
PRESSURE l ANALYSIS Bourdary Condnions RETRAN B.C.'s liot Channel ir MECHANICAL Hot Channel NI N
- Hot "ITIERMAL bgap Channel OVERPOWER Model ANALYSIS ir RE'ITIAN to TCPYA Fonnat
! h's, v s. p*s
'ICPYA ' ' Transient i Transient Critical l Cl R via ; Power i OExL Rauo Analysis l
l l FIGURE l-1 l
Flow Diacram of Existine Reload Methodology a 1
i 9 >
l i
SWTICII Initial
'C 8 1APURX Simulate to I .C. s --+ STAB M FIBWR2 Distribution Powe)r 4 T/H n
,, {***y*g*gj Qr****** ****j SIMULATE 3r
- Kinetics.
Fuel Rod ,Idlet Steady ! Bypass . 3.D Nodal Power Dist.
IAPURW Thermal ' . State .Mw PhWes Effects : T.H :
............... ..........) Initial 3-D Power Dist. 1r SUCK 1r e Point EXCLUSION
- Kinetics REGION Generation ANALYSIS Point l RETRANir i
OIE"*'
- ASME Core Avg -
Systern > TRANSIENT hgap Model FRESSURE ANALYSIS j
.f......
l FIBWR2 B.C.'s :
i . ir : MECHANICAL llot Channel Hot l 4 "I1IERMAL x
> hgap : Channel : OVERPOWER
- Model i ANALYSIS
- ir l Transient :
i criucal i l ! Power i i ....................... M 0........)
FIGURE 1-2 Flow Diacram of Reload Methodolony with FIBWR2 l
1 l
l l
10
Top of Core .
7 1
lZCH1 zou t 9
1 r -
NOTE: Bottom' entry Peripheral Fuel ZHET ff j supports are welded into the core support plate. For these bundles, path numbero 1, 2, Channel 5 and 7 do not exist.
n Spring Plug ZUHB 8
Core Support t puuniinn Lower H Tie ZCEO Plate R g1 34 L.
y Bottom of core ZSTU I
NN P
" Shroud Fuel Support ' Control ~
lb Rod Guide
- Core Length = ZCHI + ^ Tube l
Fuel Length + ZCEO Control Rod Guide Tube - Fuel Support 1.
- 2. Control Rod Guide Tube - Core Support Fuel Length = ZUHA +
ZHET + ZUHB o 5' Plate 7
- 3. Core' Support Plate - Incore Guide Control Rod Tube Drive Housing 4. Core Support Plate - Shroud 5 Control Rod Guide Tube - Drive Housing
- 6. Fuci Support - Lower Tie Plate
- 7. Control Rod Drive Cooling Water
- 8. Channel - Lower Tie Plate
- 9. Lower Tie Place Holes
- 10. Spring Plug - Core Support FIGURE l-3 1
Flow Paths in a BWR Core 11
0
, : 2.0. DESCRWTION OF RBWR2 CODE ENHANCEMENTS This section describes the key features of FIBWR2 that allow detailed modeling of fuel and bypass paths, the capability to model enhanced fuel types, and addition of transient capability.
Several of these key features are compared against' the methodology employed in the current YAEC RETRAN hot channel model which is used to evaluate the response of the limiting fuel assembly to
' AOTs. Changes to the conservation equations used in the original RBWR program and the solution
. closure relationship to enable transient capability are summarized in Section 2.1. The RBWR2 core nodalization scheme, including the basic division of control cells, the specification of the active core zone, and the use of compound nodes to model the unheated portions of a fuel channel are described in Section 2.2. Section 2.3 describes the FIBWR2 solution ' scheme, which is similar to the original FIBWR code, but now accommoda'tes a larger problem size. A separate and significant discussion in j Section 2.4 focuses on the void and pressure drop models in FIBWR2 and RETRAN which forms the basis for the two codes differences in results.
1 I
(
2.1 RBWR2 Conservation Eauations E 1 1
The mass and energy equations of the original FIBWR program are modified in the FIBWR2 f
. program to allow thermal non-equilibrium between the liquid and vapor phases. Separate liquid and '
L vapor mass conservation equations are used and the energy equation modified to assume that the
[ vapor phase remains saturated while the liquid phase may be subcooled or superheated. The assumption that the vapor phase remains saturateci is reasonable for BWR vessel conditions for those 1
[. transients and accidents which maintain high pressure (>800 psia) and which do not involve acoustic or sonic flow phenomena such as might occur during the postulated Loss of Coolant Accident -
( -(LOCA). The FIBWR2 code therefore has four conservation equations compared with three conservation equations used in the FIBWR code.
(;
The treatment ofliquid and vapor conservation required the introduction of two terms in the conservation equations for the transfer of mass and heat between the two phases. These terms for the interphase mass and heat transfer are denoted GAMMA (ry) and H,,*Ag, respectively, with P, I L assumed proportional to H,,*A,r. ' At steady state conditions, the net interphasic heat transfer approaches its equilibrium value and the FIBWR2 equations are equivalent to the RBWR equations.
'12 p
r
- a b-fL L
I f The FIBWR2 program allows the specification of several different models for H,*A,. The recommended (and default) model for H,*Ag is based on the model included in the TRAC-BF1 program and described in Section 4.1 of Reference 2. Its basis is discussed in Section 4.1 of this reference. An alternate model is also available which transfers the proper amount of energy between I phases each time step to maintain thermal equilibrium. This model is useful for cross code-comparisons to codes such as RETRAN02 and VIPRE01, which assume thermal equilibrium is .
maintained between phases.
1 2.2 FIBWR2 Core Model i
i The FIBWR and FIBWR2 programs are specifically designed to model the complex flow -
paths of the BWR geometry. The FIBWR2 code also has enhance modeling capability relative to the original FIBWR code to explicitly model features of advanced fuel designs such as axially variable flow areas due to part length fuel rods and varied water tubes designs (atriums, watercrosses) which may have flow entrances and exits above or below the upper and lower tie plates. The FIBWR and-FIBWR2 programs model all significant flow paths between the lower and upper flow plena.
Based on the specific input for the modeling application, the FIBWR2 code methodology-allows the representation of the BWR core with a large number of active channel, water tube and l bypass flow paths between the lower.and upper plena. Internal to the FIBWR2 program, the flow paths are treated as cells with nodes and junctions with each cell a control volume for mass and energy.
(- The lower and upper plena are cach modeled as single cells. The number of cells for the
- active channel flow path is specified by the user and is usually specified such that the unheated j
7
. regicns in the active channels above and below the active fuel zone with heated fuel rods are modeled with a single cell. These upper and lower unheated zones are modeled as compound nodes.
The features of these nodes include flow area chages, flow restrictions and other geometric
[ characteristics associated with the specific representadon of the flow and fuel bypass paths leading to
) exiting from the fuel. The representation of the active chmuel lower unheated node in the FIBWR2 methodology.is shown in Figure 2-2. It includes the flow paths from inlet orifices to the top of the fuel lower tie plate. The upper unheated zone is relatively simple compared to the lower unheated 13 L
l region since only one flow path is needed. This upper node simulates the flow from the upper unheated region of the fuel to the flow expansion above the fuel channel.
The number of uniform cells used to model the active fuel zone is controlled by the user within limits specified during compilation of the program. The bypass flow paths are modeled using the same number of cells as the active channels. A typical value for the number of active and bypass flow path cells is 27, with 25 cells used to model the active fuel zone for consistency with reactor physics models. The number of cells used to model the water tube flow paths is specified during compilation of the program. A typical value for the number of water tube cells is 5. Cells may be connected by a " leakage flow path" which is assumed not to store either mass nor energy. A typical application of this nodalization scheme is illustrated in Figure 2-1.
~
The capability to specify leakage flow paths between tells is significant since it allows more detailed modeling of the pressure drops and flow raus of bypass flow paths through the core support s plate and the water tube flow path than is possible with general purpose thermal hydraulic codes such as RETRAN. As a result, the effects of changes in pressure drop distributions which occur as a result of system wide perturbations can be evaluated with FIBWR2 and allow the calculation of flow distributions into and leakage paths out of an individual fuel assembly during a transient.
2.3 FIBWR2 Solution Scheme The FIBWR2 code is intended specifically for use in BWR steady state distribution analysis and transient CPR determination. The BWR specific application uses several simplifying assumptions, delineated below, which allows a code structure requiring only a simple explicit Euler solution technique. Though this technique was chosen in part because it was accommodated by the simplifying assumptions, it was also chosen since it minimizes numerical dissipation which is especially important in boundary conditions which fluctuate such as fast pressurization events or ' 5' oscillations associated with an instability.
14 L_ - ________ _ --__ -__ _ __ __- - _ - _ - _ - - _ - _ _ - _ -
- r The RBWR2 solution of the conservation equations uses the following simplifying assumptions:
' All the core flow channels experience the same overall pmssure drop.
~ -
The flow is one-dimensional everywhem, although the liquid and vapor may travel at different velocities, and even in different directions. l I
Liquid / vapor density and temperature are at saturation at a uniform core pressure. !
l
- The void and pressure drop relationships used in the steady state analysis are valid for the .l
- transient analysis.
All of the default void and pressure loss relationships in RBWR have been retained in RBWR2, i including those for:
Initiation of subcooled boiling, 1
Relation of flow quality to equilibrium quality, Relation of void fraction to flow quality, Single-phase friction factor.
Two phase multiplier for frictional losses, and, l
Two phase multiplier for local (form) losses.
l Although flow paths are treated one-dimensional by MBWR2, the thermal-hydraulic effects of complex flow patterns have been incorporated into the models and input coefficients for predicting flow quality, void fraction, and pressure drop. As discussed previously, an additional model for the vapor generation from flashing or condensation is required for transient simulations. RBWR2 has u- 15 l
- incorporated several options, the default model (used in the CPR methodology) is derived from TRAC BF-I [12]." These assumptions are reasonable for analysis AOT and allow the fluid dynamics to be' simplified the a four equation model. . Separate mass conservation equations are used for the
. liquid and vapor phases, while the momentum and energy conservation equations have been combinsd. The numerical integration of the four conservation equations formed by the code insumptions is performed with the explicit Euler integration scheme for the parallel path flow r- network of nodes defined by MBWR2 input.
\
L The FIBWR2 numerical techniques utilize an implicit solution for the calculation of fuel rod temperatures which always results in a tri-diagonal matrix which can be solved rapidly. These conservation equations and fuel rod conduction models in the RBWR2 solution scheme enable a FIBWR2 based simulation containing several channels with about two hundred nodes and heat I structures to run acceptably on a DOS based personal computer. While this would be a relatively small FIBWR2 model, much larger models containing thousands of nodes can be readily L accommodated on a workstation.
3.4 Differences in the RETRAN and FIBWR2 Thermal Hydraulic Models
. In the process of establishing the FIBWR2 CPR methodology, the thermal hydraulic models
.used in RETRAN and FIBWR2 were reviewed to establish the potential differences and possible limitations in FIBWR2. The review found that the FIBWR2 thermal hydraulic models do differ from the models used by RETRAN02, primarily due to the BWR specific nature of the models in
{
FIBWR2. In RETRAN02, subcooled voids are not modeled in the hydraulic solution for channel quality and a generalized model, appropriate for both horizontal and vertical flow, for the two phase form loss multiplier is used. FIBWR2 includes a subcooled void model. Other models available in r
' FIBWR2 were developed specifically for vertical flow in rodded bundles. FIBWR2 uses flowing quclity in the Baroczy friction multiplier model, while the Yankee RETRAN model uses the static l quality. RETRAN does not explicitly calculate bypass and water tube leakage, nor calculates the quality change above the part length rods in new fuel designs. These differences will not generally
, affect system transient calculations, where the core is lumped into a single equivalent channel.
m However, these differences may impact the MCPR response, which is determined by the flow
( response of upper nodes the most limiting assembly.
[ 16
l>
L As part of the validation effort for the RBWR2 CPR methodology, the model differences noted above were assessed to determine their impact on CPR performance and whether the HBWR2 code provided the expected results. As will be discussed in more detail in the MBWR2 and h RETRAN CPR comparisons in Chapter 5, the majority of the difference in the transient performance sf the two codes was found to be in the increased accuracy of the two phase multiplier in RBWR2.
This single difference changed the fuel channel active flow and the magnitude of the transient CPRs.
i The change in magnitude of CPR was consistent for the two limiting types of events at Vermont Yankee, the increase in pressure and decrease in feedwater temperature AOTs. The other codes' l model differences did not identifiably impact the magnitude of CPR a given transient. ;
i i
i l
I l
l I
l 1
1 r.
17
M M M M M 'M M M M M E E E E E
, i
! i i
@ $ aEb NE ga [3 s' SM 3 e eO a e
a o a
- - e -
E 7 zf us 3< +n b m g- >
?
pl -
[ Channel i = 1 .j (
E.-.
- I k y '
- E I I 1 a k = KRL+1 Water Tube i = IRL+1 ,
Z E ;; r -
- 8. : u 3 =M + 4 >
[ o 1
=. c : - I E 8 y Channel i = 2 n .p l ; L
- 9. e
~
E i I
. , , . 7 -.I E E ! I Water Tube i = IRL+2 * ""
i B
3 E :
y -
i o
(
i S i
_+
g .
1 i :
l Channel i = IRL-1 -
i I
-li k q i i i i
- L I Water Tube i = IMAX k = KMAX i k g - x s-gb c . ." -
- t Bypass i = lRL .
d, I . .
Heated, Rodded i l ,
s .
Zone !
p 3 .-
Unheated, m . h m Rodded Zone S AuHB g g g
- Y- - gl.
13 -l -f-l- l'- lt E + Leak 8 '
Lower Unheated Lower + Leak 9 Node i Tie # !
Plate , Leak 6 k
j[
Core Plate l
,Jc....n,.-,... ., .
l
~_
l 4 -l i
! J
!< A + \ l INLET Fuel g l f i : Support me%H ImiraniNuo El_ _ ___ L, I I A tyg lm rl Proportion of lower plenum flow area for this channel l
FIGURE 2-2 FIBWR2 Compound Lower Unheated Node 1
i 1
19
- __ - _A
3.0 COMPARISON OF FIBWR2 WITH ANALYTIC SOLUTIONS
- A valid test of a numerical procedure is to evaluate its prediction for a simplified problem for.
which an analytical solution may be obtained.- Comparison to two analytic solutions are contained in this section to demonstrate that the enhance FIBWR conservation equations discussed in Section 2 have been correctly installed in FIBWR2. The first problem is an exponential flow decay transient.
For the second problem, an oscillatory flow transient.is simulated. The analytical solution is presented in Section 3.1, a verification of the steady state analytical solution in Section 3.2, and the results of the evaluation of each problem with FIBWR2 are compared to the results of the analytical solution in Sections 3.3 and 3.4.
3.1 Initial Conditions of Transient Two Phase Problems To define a transient two phase problem with a general analytical solution for void fraction and mass flow, a number of simplifying assumptions were made. Specifically, the following assumptions were made:
L ' All water physical properties are evaluated at an assumed constant system pressure of 1000 psia.
The liquid and vapor phases are in thermal equilibrium.
Slip between the phases is evaluated using the drift flux model with the values of C, and Vg
! held constant.
The power distribution 'along the length of the vertical heated channel is uniform and constant, l
l -
The fluid entering the bottom of the channel is at saturated conditions.
l; -
The system is at steady state conditions prior to the imposition of the inlet flow function.
1' 1
1 \
i j 20
- The flow into the channel is always positive. ,
l l
l Values of the parameters used to completely specify the initial conditions for the problems l consistent with the above assumptions are presented in Table 3.1.
3.2 Comparison of the Steady-State Void Profile :
The void profile of the initial steady state conditions can be determined analytically for each ,
axial location by first noting that the total power delivered to the coolant, Q, at an axial location, z, for a uniform and constant power distribution is Q=z,Q z, (1) where Q, and z, are the total channel power and length, respectively.
All the power deposited in the coolant is used to change the phase of coolant from liquid to vapor since the coolant at the entrance to the channel is assumed to be at saturated conditicas. At any axial location, the mass flow rate of the vapor is therefore the power that has been deposited divided by the heat of vaporization:
I o w* - h,-h, m I :
I 21 I - _- - - - --
I For the initial steady state conditions, the mass flow is constant and the mass flow rate of the liquid portion of the two phase flow is:
1 !
W,= Wu - W, (3y
'1 The drift velocity is then calculated:
I ' ~i.s f='P / E>
r (4) f i
l and the vapor velocity determined from the drift flux relationship:
v,=CpVg (5) l Using the vapor velocity and vapor mass flow rate, the void fraction is:
I w (6)
Av,p, I The void fraction profile calculated using the analytical solution presented above is compared to the initial void fraction profile evaluated with FIBWR2 in Figure 3-1. Inspection of Figure 3-1 indicates very good agreement between the analytical solution to the void fraction profile at the l initial steady state conditions and the results of the FIBWR2 evaluation. )
I 22 I
3.3. Description of Analytical Solution Techniaue
- The general analytical sch:ba for the void fraction profile under the assumptions specified in Section 3.1 is baEd on conservation of vapor mass:
a a g(A8P,)+p(Aap,v,) = AP (7)
When the expression for the vapor drift velocity of equation 5 is substituted into the above equation, the terms' can be rearranged to yield the following expression by recognizing that the area,
~ A, tie drift-flux parameters, C, and V ,g and liquid and vapor densities, p, and p,, are assumed constant:
4
" as da af T+(9+Vs )p = P -aCo (8)
The terms of the left side of equation 8 can be interpreted as a material derivative: ,
1 Da ,_P-aCy-0 Dt p, az (9) :
-with a corresponding velocity of:
dZ j " 9+Vs (io) b 23 .
^1 I
o j
=
i l
When the initial conditions are known, equation 9 can be integrated to yield the void fraction ct an axial location with that location determined by the relationship obtained by the integration of l
equation 10.
To perform these integrations, an expression for the drift velocity is required.
h=T
- P
= k, say (11) r Ph where k is constant. The value of j at the inlet is determined by the iniet mass flow rate, which is:
1 I
f(z=0,t) = I (W o+Wi e'8#; = fo+f e-8', say (12) 3 l
I when the inlet mass flow rate is assumed to be of the general form, W = Wo + W ie+. When l equation 11 is integrated and combined with equation 12 the result is: )
[I l
j(z,t) = fo+f ic'8'+kz (13)
I 1 which upon substitution into equation 10 yields the following expression for the axial location, z:
I l &
7-C kz o = Cgo+C/3e -#'+ V, (14) lI l
I ;
I l 4
l l
l l
This expression for the axial location, z, has the general solution:
(e@+qe++0 z (15)
- I where
l y-5 0+Cok (16)
I Vg+ Colo 0= (17)
Ck o
I and the value of ( is determined by the initial conditions.
.I Similarly, substitution of equation 13 into equation 9 yields:
Il" +aCo k = 1 (18)
Dt p, l
l which has the general solution:
I 1
l a = Ac'@+p (I9)
I i
I 25 '
.I l where:
P (20)
I P,Cok l
l l
and the value f A is determined by the initial conditions.
.E
- m i Values for the parameters and A must be determined if the void fraction is to be determined for a selected axial location and time. The value of ( can be determined by solving equation 15 directly for ( using the selected values of axial location and time. The value of A can be determined if the value of the void fraction can be determined for either any time or any axial location. Since the initial void fraction as a function of axial location is given by equation 6, if the value of the axial location, z is 20 as determined using equation 15, at time t=0, the value of the void fraction is known and equation 19 can be used to determine the value of A. If the value of the axial location, z is <0 at time t=0 as determined using equation 15, equation 15 can be used to iteratively solve for the time at which the value of z=0. Using this time and the assumption that water entering the inlet "l channel is at saturated conditions,i.e., the void fraction is zero, the value of A can be determined using equation 19. With the values of and A determined in this manner, the value of the void l fraction for the selected axial location and time are known by equations 15 and 19.
- I The above solution for the void fraction profile for an inlet flow function of the form, W = Wo + W ie+ is valid for both real and complex values of the parameters Wo, W i and p. A I
FORTRAN program was written to evaluate results based on this analytical solution.
3.4 Flow Decay Transient The FORTRAN program to evaluate a flow decay transient. The inlet flow function parameters selected to yield a flow decay condition consisted of all real values and were:
i W(t) = 111500 c-5' lb/hr (21) lI- 26 lI
1
{ '
The results of the analytical solution to the flow decay transient as evaluated using the FORTRAN code are presented along with the results of a HBWR2 evaluation of the same conditions in Figure 3-2 for the exit void fraction. The difference between the analytical and MBWR2 based veld fraction values were insignicant. Values of the analytical and FIBWR2 based exit mass flow rates are presented in Figure 3-3. The agreement between HBWR2 and the results based on the analytical solution are very good for the flow decay transient.
{ - 3.5 Oscillatory Flow Transieg The FORTRAN program was used to evaluate an oscillatory flow transient. The inlet flow
{
function parameters selected to yield the oscillatory flow function include both real and imaginary values. The real part of this flow function is equivalent to:
{-
i b
W(t) = 111500+108000 sin (5t) (22) l The results of the analytical solution to the oscillatory flow transient for the exit void fraction l are presented in Figure 3-4. The results obtained from evaluation of the oscillatory flow with a i
[ nominal FIBWR2 model with 24 axial nodes are also presented in Figure 3-4. The agreement between the results of the analytical solution and nominal HBWR2 based evaluations is very gocd l early in the transient. However, during the latter portion of the transient period evaluated, the
[
nominal FIBWR2 results do not predict the rapid changes in exit void fraction calculated with the
[ Cnalytical solution. By inspection of these results it can be determined that the void fraction is - l changing more rapidly than the relatively large control cell node volumes of the nominal FIBWR2 h model can resolve. - Therefore the inability of the nominal FIBWR2 model with 25 nodes to predict the rapid changes in exit void fraction is not unexpected. To demonstrate the impact of control cell b ; node volume, a HBWR2 model of the oscillatory flow transient with 50 nodes was preparea. The exit void fraction results obtained from this 50 node FIBWR2 model are also presented in Figure
- 3-4. Inspection of Figure 3-4, indicates that the agreement between the results based on the analytical solution and the 50 node FIBWR2 model are nearly exact for the duration of the transient modeled. l h l w
k 27.
n V -- -_ _ - -
3.6 Summary of Results for the Analvtic Comparisons The results of FIBWR2 evaluations were compared to results based on the analytical solution to two transient flow problems initiated from the same initial conditions. The void fraction profile calculated with FIBWR2 for the initial steady state conditions was essentially the same as the results obtained from the evaluation of the analytical solution. Very good agreement between the FIBWR2 l and analytical solution based results for the exit void fraction and exit flow were obtained for the I flow decay transient conditions evaluated. Agreement between the FIBWR2 and analytical solution exit void fraction results for the oscillatory flow transient conditions was shown to be dependent upon the number of nodes selected for use in the FIBWR2 model. Tivis dependence of the number of nodes is to be expected for thermal hydraulics codes employing a control cell method for this type of transient. These comparisons and the results obtained demonstrate the appropriateness and correct implementation of the conservation equations, thermal hydraulic correlations, and numerical ;
- procedures of FIBWR2. It should also be noted that for transients involving oscillatory behavior, the number of nodes chosen is particularly important. For flow decay transients, the choice of 25 nodes was appropriate as demonstrated by the excellent agreement with the analytical solution.
I
- I 28 I
l TABLE 3.1 I Initial Conditions Assumed for Transient Two Phase Problems Parameter / Property at 1000 nsia Value Units Water inlet enthalpy 542.6 Btu /lb Initial mass flow rate at inlet 111500 lb/hr Channel power 1.00 MW Channel length 12 ft Channel flow area 16.056 in:
Drift flux coefficient C, 1.13 --
Drift flux coefficient V g i ft/s Saturated liquid density 46.3156 lb/ft' Saturated vapor density 2.2423 lb/ft' Saturated liquid enthalpy 542.6 I Saturated vapor enthalpy 1192.8422 Btu /lb Btu /lb l
l l
t I
l l
1 29 I
0.45 8 r , , , 8 , ,
0.4 -
m
- a=s'O'O - ~
0.35 - 0 e.,e-c.
w E 0.3 -
a*c'0, a
<m e o FIBWR2
,#,a, -
k C3 $ 025 - "D ,c' ta M e
- a o :;o u.
$ $ 02 - ,#'g O, 1, > p' o , , -
9 0.15 - "
g' B p a
u ,
3- 0.1 -
- g O
, o ,0 _
e 0.05 - ,D o
,A - . ,
0.0 -' , ,- .
80 1h 120 140 160 Elevation (in)
W W W W M M M + M M M M M m m m7 0.75 ' ' ,
o o ,-
o,,*'
O.7 - o, '
o,,** -
T!
O,**
F o, ' .
c N 0.65 - c,.S -
8 __ Q'
% e o FIBWR2 9' ,
,- I 3 .g 06- -- Analybcal Solution ,0 * * ~
@ O M '
P' ;
w a C it ,o-
- a x ,
R 'T
> 0.55 -
. w ,o' -
< U V k 0.5 - cr
,G m
8 ,#,
E- F O.45 -
y -
,G '
c..O' O.4 i 8 e ,
0.0 0.5 1.0 1.5 2.0 2.5 3.0 Time (sec)
C E
t i I t t 3 m o ei b d, d,
?
- - 9 - y
,9
. ,9
,P
- ,P -
o cJ x ?
G S F P L oe' l e
,P
~
,P -
a ~f p
m -
0 r
- h i
~
,b 6
L
/
l -
- i l
,10 m
0 d a
p .
..e' F ,a**,0**
. o "I i i i i i d
- 8 8 N 8 S L .(s/qdmold
~ ,
! . i L
F FIGURE 3-3 Flow Decay Transient - Exit Mode Flow Rates I
- 32 L.
W
1 F t i 1 J' I I I I-. 1 I { } . 'l I ( .I { -l I I l . .. l I I E E E E E gg i i . e i i 3 . -~.
o o s
g G
o ,, ' o
. 9, 9 0.42 - FIBWR2 ,'
I oj,' o S.,i c ,,
- - * *
- 0 ,'
, o ; O ',
8-3
'n l ; o l o, ;
A 's. / . /
2 e o.4 - '. 0
" / l / -
< 3 S 4
,/ o 5' O G . ,/
", C E '.
o/ O
,/
O x m
a 2
< o '.
/
P. l j- % > 0.38 - 6 j j -
,n . o., .
g . ,
s ,
n '. / '
/
i'
= '.
,/ '.- ,l
. 0.36 - 4 o,f -
9 Vs '
, l
- 8. b~D 2
0.34 i i i i i 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Time (sec)
0.44
,d ', g' C J o
- 4 p' O o ,
O b
0 25 Node RBWR2 h' , O ,,, ' ',
0.42 - * ~ i
$ 0 50 Node RBWR2 O I, O 9' ,
E r4' g
-- Analytical Solution h'/ .
O q'
{ .
O/
4 's O/ 0; *
/
a ',
- o
- O O*
m e 0.4 - -
a
~
Q o
u D* / O e
. ,/
w A
E C
8 s
Of D'
=
=
0 a ,/
M '.
o w u
a .
' l B
- O /
l E a > c.x - 6 -
m l . ,
B o
t 8 I
I e o 3
g n * % a s
o 4 8
's,, l 3
0.36 - Q l -
2: '. :
g_ 6.g*
a 0.34 4
, I i 0.0 0.5 1.0 - 2.0 2.5 3.0 Time (sec)
4.0 COMPARISON OF FIBWR2 TO CRITICAL HEAT FLUX DATA Qualification of FIBWR2 was performed against selected Critical Heat Flux (CHF) experimental test data for the 16 rod electrically heated test section identified as 16R-2 [3,4,5]. The tests' selected for simulation were those tests with flow rates within the validated range of the -
General Electric proprietary Critical Power Ratio (CPR) GEXL correlation. These tests were conducted with values of mass flux in the range from 0.25 to 1.0 Mlbm/hrft ,2 pressures in the range from 800 to 1,000 psia, and inlet subcooling in the range from 22 to 165 Btu /lbm. The 16 rod tests are representative of the fuel assemblies used in boiling water reactors and were specifically designed and conducted to produce conditions which result in transition boiling as indicated by a rapid rise in the temperature of the heated surface. A base FIBWR2 model was developed and used in combination with test specific data to simulate 50 steady state and 27 transient tests.
4.1 Simulation of Steady State and Transient 16 Rod Tests
[' ,
Information describing the general configuration of the experimental test program [3,4,5] was used to construct a base FIBWR2 model. The 16 rod test bundle was constructed with direct heaters with variable cladding wall thickness along the length to provide a chopped cosine heat flux distribution with a peak to average value of 1.387. The heat flux distribution was simulated by specifying an axial power profile. The 16 rods were modeled using one hot rod and one lumped rod for the remaining 15 rods with 24 axial geometry nodes. The varying cladding thickness was
- simulated by an appropriate number of lattices and rod types to model the hot and lumped rods.
The onset of transition boiling and the associated changes in the heat transfer regime is
- predicted by empirical correlations such as the proprietary General Electric GEXL and GEXL-PLUS correlations. FIBWR2 is structured to allow implementation of fuel vendor specific correlations within the program coding. The calculated parameters required for boiling transition correlations can I
- therefore be provided either for post processing or used within the FIBWR2 program execution to calculate parameters such as CPR to assess the approach to the onset of transition boiling. !
Simulations of the steady state conditions under which the onset of boiling transition was observed during the 16 rod tests were performed using the FIBWR2 base model. The CPR values 35 L
calculated using FIBWR2 for each of the 50 steady state tests are plotted as a function of the !
l measured critical power for each test in Figure 5-1. Inspection of Figure 5-1 indicates that for all i 1
steady state tests, the onset of boiling transition, (i.e., CPR <l.00) was predicted within about 5%
and that the variation of the calculated CPR fmm a nominal value of 1.0 decreases slightly as the measured critica! power increases. For these 50 steady state tests, the mean of the FIBWR2 based CPR values was 1.007 with a standard deviation of 0.0186. These results indicate that the 16 rod test configuration was appropriately represented by the FIBWR2 base model and that the CPR correlation was properly implemented within the FIBWR2 program coding structure. For comparison, the mean and standard deviation of the CPR value for these same 50 steady state cases obtained from simulations using the RETRAN and TCPYA01 models were 1.006 and 0.0182 [9]. .
Transient 16 rod test conditions were simulated to demonstrate the capability of FIBWR2 in :
conjunction with a CPR correlation to adequately predict the onset of boiling transition under l l
changing thermal hydraulics conditions in a fuel assembly. This demonstration tests FIBWR2's ability to appropriately evaluate the changing thermal and hydraulic parameters at locations within the assembly and the fidelity of the CPR correlation. Therefore, the capability of a model to adequately predict the timing, and even the occurrence of the onset of boiling transition under test conditions where boiling transiticn was observed, is dependent on both the analysis of the transient i 1
thermal hydraulic conditions and the empirical CPR correlation. In the following paragraphs, the results of the FIBWR2 based simulations are compared with the test results and the results of similar RETRAN based simulations for a total 27 transient test cases. The 27 transient test cases which were simulated with both FIBWR2 and RETRAN consisted of 13 flow decay with constant power cases,12 flow and power decay cases and 2 flow and power increase cases. For each of these 27
~
transient cases, the time of onset of boiling transition was taken to be the time at which a rapid rise in surface temperature was indicated in the reported test data. BWR MCPR methodology uses a "best fit" CPR correlation, therefore, these simulations did not always predict CPR values equal to or less than 1.0. The time at which the CPR value was equal to 1.0 was taken as the calculated time of onset of boiling transition for the cases where the CPR went below 1.0.
The onset of boiling transition was calculated to occur during 20 of the 27 transient simulations. These times are plotted against the measured times of the onset of boiling transition in Figure 5-2. The variance (root mean square) of the differences between the measured times and the 36
{ f times calculated based on the FIBWR2 simulations for these 20 transient cases is 0.458 seconds.
The same 27 transient cases were simulated using the RETRAN and TCPYA01 models [10] with the -
[ onset of boiling transition calculated to occur during 19 of the 27 transient cases. The variance of
-the , differences between the measured and RETRAN and TCPYA01 calculated times of the onset of h boiling transition is 0.593. The small variance for both the FIBWR2 and RETRAN simulations indicate good agreement with the experimental results with slightly better agreement for the FIBWR2 b ' simulations relative to the RETRAN simulations.
b- Comparison of the CPR values and times of boiling tiansition calculated based on FIBWR2 simulations of experimental test data demonstrate the capability of FIBWR2 to accurately calculate
- thermal hydraulic parameters necessary for the prediction of the onset of boiling transition under both steady state and transient conditions. In addition, the good agreement between the results of the b FIBWR2 based simulations and the results of the RETRAN and TCPYA01 model simulations,
. previou ly approved for hot channel analysis, further demonstrates the applicability of FIBWR2 to steady state and transient thermal hydraulic analysis.
I L 4.2 Sensitivity Studies g
Sensitivity studies were performed to document the impact of three representative alternative choices that may be selected for simulations using FIBWR2. The three alternative modeling specifications considered were axial nodalization, void model and time step size. The 16 rod test cases for the 13 flow decay with constant power transients were selected for the sensitivity studies based on the relative severity of the change in thermal hydraulic conditions experienced in these p transients. For each case, the results of the FIBWR2 sensitivity studies were compared to a base case FIBWR2 simulation. Th'e base case FIBWR2 simulations used 24 axial nodes, the EPRI void F model and minimum and maximum time step sizes of 0.001 and 0.005 seconds, respectively.
E
{ The calculated times of the onset of boiling transition for the axial nodalization sensitivity study which reduced the base case number axial nodes from 24 to 12, the void model sensitivity F
study which specified the use of the homogeneous model and the time step size sensitivity study which increased the minimum and maximum time step sizes to 0.01 and 0.05 seconds, respectively, y ' are compared against the results for the base case simulations in Figure 5-3 for all cases for which F~ - 37 a
.. . _ _ _ _ _ _ _ _ _ - _ - - _ - - _ - - - - - - - - - - - - - - _---------__----------------------A
the onset of boiling transition was calculated. Inspection of Figure 5-3 indicates that specification of each of the alternate simulation options produces a negligible effect on the calculated time of boiling transition onset. Based on these studies it is concluded that the base case time steps and nodalization (0.01 seconds and 24 nodes) are adequate for a FIBWR2 hot channel model. As indicated in Section 2, the nodalization may be dependent on the type of transient to be simulated; thus; a sensitivity study may be required to ensure convergence of results. Convergence would be obtained if the parameter changed did not significantly change the results.
4.3 Summary of 16 Rod Test Simulations t
The FIBWR2 program was used to simulate selected steady state and transient test cases of 16 rod critical heat flux experiments. The results of these simulations were compared with both the test data and the results of RETRAN simulations of the same test data. The FIBWR2 based evaluations of the CPR for the steady state cases were in good agreement with the test data results and similar RETRAN simulations. The results of the FIBWR2 based evaluations of the time of onset of boiling transition for the transieni cases were also in good agreement with reported test data results. For those transient cases for which the onset of boiling transition was calculated, the variance between the calculated and measured times based on the FIBWR2 simulations was slightly less than for similar RETRAN simulations. The good agreement between the FIBWR2 results and both the test data and the results of the RETRAN simulations demonstrate the capability of the FIBWR2 models in conjunction with critical power correlations to adequately predict the onset of boiling transition. Furthermore, the results of sensitivity studies which specified different modeling options for axial nodalization, void model and time step size demonstrate the adequacy and l
robustness of FIBWR2 over a range of values and available alternative specifications. To obtain convergence of the FIBWR2 simulations a sensitivity study on time step and nodalization may be l required.
l l
=
38
(
l 8 8 8 e i e o f e' l
. g
- q. .
f' k.
- n' . o r %l ~d{
l e* !,,
- 5
- I
[
' - *' . S d .B
. . .'
- e
- O
. *i
]!=
r
- g
'
- O 2 l
~
e d L
l e l e
.e.
L ee
- - o!. q o gy l
L .
.' o i i
3 3 3 I i i g g H
,_ h $ h o enleA UdOW elels ApeelS EWMGld i
I FIGURE 4-1 FIBWR2 Steady State 16 Rod CPR Results r-s 39
[ ... 1....f. .Ii...f... 1 .. l....f.... g
- -N b _ _.-
L O O
~
O i
_ _m O
u
~
A O e F - -t .9 L o .
g O E O F u-1 E
_ _n=
j
[- :
N o 4
_- 1 _mg e
~
e _
EEE _,
i 004
~
P ^
l
" O
,,,,,,,.g,, ,,,,.,,,, y, ,,.. ,y....
C N O to T to N e O F (spuooes) sewillosuo uo!ysueJ16u!!!o8 EUM81.-l L
E E
m FIGURE 4-2 FIBWR2 Transient 16 Rod Boiline Transition Onset Results F
u a
L
~..
.aw, . . - . . . m_ .
3 E E 5 5 T p y 3 I I I E I 5 5 5 5 E 3 I E E I I E E I I T T 6
D
~
.a s
.O -
.g
- ~ "
w
- O -m f-g 3
O
- -> g
- n
! 3 a S 1-O -
c
~ - JC y
.O c
5 [
~ ne f 3
- emy* o (
= -mz
-s j Z E hC
.d. h
- OO> -
O
'''''''''''''''I*
DC N )0 tU A JC N 6- O
)sdnoces( semiT tesnO noitisnarT gnilioB ydutS ytivitilsneS
- I 4LDnMB V-C y11 ys [9 uop lasi 3dM dlap!apous - saus!pplA sinpA gasnys
!I I
5.0 FIBWR2 TRANSIENT HOT CHANNEL PERFORMANCE The applicability of the FIBWR2 code to the evaluation of transient hot channel thermal
{. hydraulic conditions, including' the calculation of CPR, is demonstrated by comparison to similar evaluations performed using the RETRAN and TCPYA01 hot channel ' methodology. The RETRAN i
[; and TCPYA01 models are currently used in the YAEC BWR reload analysis methodology [7,8,10].
Comparison of FIBWR2 and RETRAN hot channel hydraulic parameters for two characteristic
[ anticipated' transients is presented along with the change in CPR values for several limiting and near limiting postulated events. Differences in the thermal hydraulic parameters and in the change in-( CPR calculated for these events are identified and determined to be primarily the result of differences in the modeling of two phase multipliers for local form losses in the respective models; these will be
[- discussed further in Section 5.2.
h 5.1: Hot Channel Methodolony -
b The objective of the hot channel methodology is to evaluate the thermal hydraulic conditions of the fuel assembly that most nearly approaches the onset of boiling transition during postulated b events.' The approach of the hot channel to the onset of boiling transition is determined by the postulated event and the initial conditions of the hot channel and core.
L-In the current YAEC BWR reload analysis methodology a RETRAN model of the plant and E
its controls is utilized to establish the system wide response to the postulated event. This system wide response is used to set boundary conditions of upper and lower plena pressures and normalized power response for a second RETRAN model of the hot channel. The RETRAN hot channel model 4
consists of a single fuel assembly, a core bypass region and upper and lower plena. The pressure
' loss coefficients of the RETRAN hot chann:1 model are adjusted as necessary to make the pressure 7 drop and flow distributions agree with the results of a FIBWR analysis which models the various core flow paths in detall.
p The transient thermal hydraulic conditions for the single hot channel fuel assembly is then p evaluated using the RETRAN hot channel m(xlel and the boundary conditions set by the system wide RETRAN model transient response. The calculated hot channel thermal hydraulic performance o
42 I.
D 4
[1 . . . . . . . . , ____-- _ _ _ _ _ _ _ _ - - _ _ _ - - - _ - - - - _ _ _ - - - - - - - - - - -- -- -
)
parameters are subsequently evaluated using the TCPYA01 model to calculate the initial and minimum CPR values for the fuel assembly.
{ In an iterative process, the initial power level of the fuel assembly in the RETRAN hot channel'model is then adjusted until the change in CPR caused by the event ofinterest approaches
[ the predicted ' onset of boiling transition or safety limit minimum CPR value. The safety limit i minimum CPR value is the lowest allowable CPR and is set based on the combined uncertainties associated with the experimental data used to construct the CPR correlation, instrument, power distribution, core power and core flow uncenainties. The value of the safety limit minimum CPR is L tyPi cally on the order of 1.08. ~
( The capability of the FIBWR2 program to calculate both transient thermal hydraulic conditions and initial and minimum CPR values permits the RETRAN hot channel and TCPYA01 models to be replaced with a single code. Furthermore, the greater detail in modeling the BWR -
specific fuel and core geometry and the capability to use thermal hydraulic models specifically b developed for application to BWR conditions available with the FIBWR2 code allows greater fidelity g in modeling the pressure drop and flow distributions than can be obtained with the RETRAN hot t channel model, r.
U 5.2 Hot Channel Results F-L Comparisons of FIBWR2 based calculations associated with the evaluation cf hot channel g response are presented for both key thermal hydraulic parameters and for the resulting change in L
CPR for selected postulated transient events. Key thermal hydraulic parameters values are presented r for FIBWR2 and RETRAN transient evaluations of two postulated events with the same initial, and therefore the same transient, hot channel power levels. As expected, the two codes calculated
- r. - thermal hydraulic parameters are different due to the differences in the thermal hydraulic models.
This causes a difference in the minimum CPR calculated for each event. In addition to the comparison of thermal hydraulic parameters, the change in CPR values calculated using the FIBWR2 '
and RETRAN based hot channel models are also presented for several post"?ated events.
E L
L 43 L..
I
L5.2.1 Comparison of Thermal Hydraulic Parameters
{ Key values of the transient thermal hydraulic conditions for the evaluation of CPR are the active coolant flow and void fraction. Comparison of these parameters calculated with the FIBWR2 and RETRAN transient hot channel models are presented for a pressurization event and an increased
{
subcooling event. These events are the Generator Load Rejection Without Bypass (GLRWOBP) pressurization event and the Inadvertent High Pressure Core Injection (IHPCI), Loss of Feedwater Heater (LOFWH) and Loss of Stator Cooling (LOSC) events resulting in an increase in core
.subcooling. All cases were analyzed for Vermont Yankee.
{ The GLWOBP event was selected as a typical pressurization event. The power level and the upper plenum pressure during the event, determined from the system wide RETRAN evaluation, and used as boundary conditions for both the FIBWR2 and RETRAN hot channel models, are presented in Figures 5-1 and 5-2, respectively. The active coolant flow for the event calculated by the FIBWR2 and RETRAN transient hot channel models is presented in Figure 5-3. The initial active coolant flow rates are nearly identical. The exit void fraction is similarly presented for the two e
[ models in Figure 5-4. The initial void fractions are different due to the water tubes models in each
. code. In the RETRAN simulation of the water tube, the water tube flow is modelled to exit above L the channel active flow. In the FIBWR2 model, the water tube flow recombines with the coolant flow in the fuel assembly at the actual height associated with the fuel design, which is six inches i below the top of the active fuel. This difference in the water tube flow exit heights results in different channel exit void fractions observed in Figure 5-4, but does not impact the minimum CPR k which occurs below the exit of the water tubes.
As the GLRWOBP event progresses, the active coolant flow rates calculated by the two 7 models remain nearly the same until the pressure starts to increase at about 0.5 seconds. At this time, differences in the two phase pressure losses in the fuel assembly begin to occur. As expected,
.the models also indicate decreasing exit void fraction due to the pressurization starting at about 0.5 seconds. With this decrease in void fraction, the two-phase pressure loss of the hot channel begins to decrease.-
L F
s O
F
___ _=
Due to the decrease in two phase pressure loss'at this time, an increase in hot channel active coolant flow is observed for both the FIBWR2 and RETRAN models. Howevrr, the increase in active coolant flow calculated with the FIBWR2 model is larger than that calculated with the RETRAN model. - This difference in calculated active coolant flow is due to the modeling of the two phase form loss multiplier between the two codes. As described in Section 2.4, the RETRAN local
_ form loss two phase multiplier is less sensitive to changes in the void fraction than the corresponding
. FIBWR2 model. As a result, RETRAN calculates a smaller change in pressure loss and a subsequent smaller change in flow rate than FIBWR2. As previously discussed the FIBWR2 local form loss multiplier is a BWR fuel specific model and is a more accurate representation than RETRAN for this application. As will be discussed in Section 5.2.2, while the RETRAN code
_ provides a conservative prediction of transient CPR, FIBWR2 provides a more accurate result..
j . The Loss of Stator Cooling (LOSC) event was selected as a representative increase in
[ subcooling event since it generally results in the largest change in transient CPR. The power level l
and the upper plenum pressure during the event, determined from the RETRAN core wide evaluation and used as boundary conditions for both the FIBWR2 and RETRAN hot channel models, are presented in Figures 5 5 and 5-6, respectively. The active coolant flow for the event calculated by the FIBWR2 and RETRAN transient hot channel models is presented in Figure 5-7. As wr.s seen for the GLRWOBP event, the initial active coolant flow rates are nearly identical. The exit void fraction for the two models are presented in Figure 5-8. The difference in the exit void fraction is similar to that observed for the GLRWOBP pressurization event and as noted in the discussion of the results for that event is the result of differences between the two models in the modeling of the water tube flow.
l For the LOSC event the exit void fraction increases slowly as expected due to the core power rise. The greater sensitivity of the FIBWR2 local form loss two phase multiplier to void fraction 1
l 1
causes a slightly greater pressure increase than the RETRAN model. As a result, the FIBWR2 '
calculated hot channel active flow is slightly less than the active coolant flow calculated with the '
RETRAN model. . The FIBWR2 simulation thus predicts a larger change in transient CPR. These l . results will be discussed further in the following section.
l l
I L 45 L
p
E Due to the decrease in two phase pressure loss at this time, an increase in hot channel active I
coolant flow is observed for both the FIBWR2 and RETRAN models. However, the increase in active coolant flow calculated with the FIBWR2 model is larger than that calculated with the l
RETRAN model. This variance in calculated active coolant flow is due to the differences in the j modeling of the two phase form loss multipliers. As previously discussed in Section 2.4, the RETRAN local form loss two phase multiplier is less sensitive to changes in the void fraction than j the corresponding FIBWR2 model. As a resc!t, RETRAN calculates a smaller change in pressure loss and a subsequent smaller change in flow rate than FIBWR2. The FIBWR2 local form loss
} multiplier is a BWR fuel specific model and is a more accurate representation than RETRAN for this application. As will be discussed in Section 5.2.2, while the RETRAN code provides a conservative l
prediction of transient CPR, FIBWR2 provides a more accurate result.
l The Loss of Stator Cooling (LOSC) event was selected as a representative increase in subcooling event since it generally results in the largest change in transient CPR. The power level L and the upper plenum pressure during the event, determined from the RETRAN core wide evaluation and used as boundary conditions for both the FIBWR2 and RETRAN hot channel models, are
~
l presented in Figures 5-5 and 5-6, respectively. The active coolant flow for the event calculated by the FIBWR2 and RETRAN transient hot chancel modelt is presented in Figure 5-7 As was seen for
( the GLRWOBP event, the initial active coolant thw rates are nearly identical. The exit void fraction for the two models are presented in Figure 5-8. The d/ference in the exit void fraction is similar to that observed for the GLRWOBP pressurization event and as noted in the discussion of the results
- for that event is the result of differences between the two models in the modeling of the water tube flow.
For the LOSC event the exit void fraction increases slowly as expected due to the core power rise. The greater sensitivity of the FIDWR2 local form loss two phase multiplier to void fraction causes a slightly greater pressure increase than the RETRAN model. As a result, the FIBWR2 calculated hot channel active flow is slightly less than the active coolant flow calculated with the 4 RETRAN model. The FIBWR2 simulation thus predicts a larger change in transient CPR. These !
results will be discussed further in the following section.
M P
- 45 W
. 'na-.
a _"aa....\,a,. -
. _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ - _ _ _ _______.___-_____._____________________________,______________._____.____.__________a
5.2.2 Comparison of Calculated CPR Values Evaluations of the change in CPR using the transient capability of FIBWR2 as an alternative to the use of the RETRAN transient hot channel model and TCPYA01 were conducted for several postulated events. The results for selected pressurization and increased subcooling events for both the RETRAN and FIBWR2 based models are presented in Tables 5.1 and 5.2, respectively.
The largest difference in the change in CPR calculated using the two models was for the Turbine Trip without Bypass (TTWOBP) event at the end of full power life (EOFPL) reported in Table 5.1. The change in CPR calculated for this event was the most limiting for both the RETRAN and FIDWR2 based hot channel models. The reduction of about 0.051 in the calculated change in the CPR for this event using the FIBWR2 hot channel model is due to the greater sensitivity of the local form loss two phase multiplier of FIBWR2 relative to RETRAN. These results are consistent with the changes in hot channel active flow calculated with each code previously discussed in l
Section 2.4 and as shown in Section 5.2.1. While the RETRAN code predicts a lower active flow and more conservative transient CPR for the pressurization event, the local form loss two phase I multiplier model in FIBWR2 code is more applicable to the conditions found under BWR conditions l
than the more general model available with the YAEC RETRAN hot model.
In comparison, the FIBWR2 hot channel results for the increased subcooling even's result in slightly larger changes in CPR (about +0.01) than RETRAN. The results obtained for the subcooling transients are consistent with the results for the LOSC discussed in Section 5.2.1. The FIBWR2 model predicts a smaller active flow rate for all cases and subsequently an increase in the transient CPR. The FIBWR2 results provide the most accurate representation of the subcooling events.
46
TABLE 5.1 Calculated CPR Values for Pressurization Events Event / Notes
- RETRAN/TCPYA01 FIBWR2 Difference ACPR ACPR RETRAN-FIBWR2 TTWOBP/MST/ EOF"L 0.239 0.21 0.029 TTWOBP/MSTEOFPL-1 0.176 0.157 0.019 TTWOBP/MSTEOFPL-2 0.099 0.077 0.022 TTWOBP/MSTELLA 0.088 0.055 0.033 1TWOBP/67BEOFPL 0.307 0.256 0.051 TTWOBP/67BEOFPL-1 0.216 0.197 0.019
~ TTWOBP/67BEOFPL-2 0.141 0.12 0.021 TTWOBP/67BELLA 0.134 0.087 0.047
.GLRWOBP/MSTEOFPL 0.22 0.192 0.028 7
GLRWOBP/MSTEOFPL-1 0.152 0.127 0.025 P GLRWOBP/MSTEOFPL-2 0.071 0.038 0.033 GLRWOBP/MSTELLA 0.068 0.028 0.040 l GLRWOBP/67B/EOFPL 0.281 0.236 0.045 GLRWOBP/67BEOFPL-1 0.212 0.18 0.032 c GLRWOBP/67B/EOFPL-2 0.122 0.108 0.014 s GLRWOBP/67B/ELLA 0.127 0.068 0.059
- Notes:
J All events initiated from nominal full power and flow conditions unless otherwise noted GLRWOBP - Generator Load Rejection without Bypass 1TWOBP - Turbine Trip without Bypass F MST - Measured Scram Time ;
67B - Table 67B Scram Times EOFPL - End of Full Power Life; -1, minus 1 GWD/STU; -2, minus 2 GWD/STU
- ELLA - Extended Load Line Limit Analysis initial conditions of 100%
Power /87% Flow
?
47 L
?
C_-
TABLE 5.2 I
l Calculated CPR Values for Subcooling Events EventMotes.
RETRAN/TCPYA01 FIBWR2 Difference ACPR ACPR RETRAN-FTBWR2 LOFWH/EOFPL 0.135 0.141 -0.006 LOFWH/EOFPL-1 0.108 0.I15 -0.007 LOFWH/EOFPL-2 0.107 0.I14 -0.0 LOFWH/BOL 0.109 0.114 -0.005 IHPCI/EOFPL 0.112 0.118 -0.006 IHPCI/EOFPL-1 0.113 0.119 -0.006 IHPCI/EOFPL-2 0.118 0.119 -0.001 I IHPCI/BOL 0.I13 0.119 -0.006 LOSC/EOFPL 0.135 0.139 -0.004 LOSC/EOFPL-1 0.135 0.14 -0.005 LOSC/EOFPL-2 0.134 0.14 -0.006 l LOSC/BOL 0.135 0.14 -0.005 I
- Notes:
All events initiated from nominal full power and flow conditions unless otherwise noted l ;
LOFWH - Loss of Feedwater Heating j IHPCI - Inadvertent High Pressure Coolant Injection actuation :
LOSC - Loss of Stator Cooling EOFPL - End of Full Power Life; -1, minus 1 GWD/STU; -2, minus 2 GWD/STU l BOL - Beginning of Life ,
I B l h
48
l A i 1 a l 1 1 a a l e a 1 a l a e a a l n a 1 a O i v
m
- g e
o n
m d
n
. 3h
. 9 o
- 9 o
l o
o e o El 8 m,- o,-
(4Wm)Jomod leuusyo i..
I MGURE 5-1 GLRWOBP Boundary Conditions Channel Power 49
1 I
4 L
i' >
i
,,, l .,,,t . ,,t ,,,,t ,,,,
[-- .
r.
",I c
i r
~
f i
. -y*
L ^
m
.- o
-2 l
I o
, . . ,..,,,,,,,,,...,,,,. 6 i
g g 8 8 @
- - = = - -
(elsd) eJnsseld wnueld Jeddn FIGURE 5-2 GLRWOBP Boundary Conditions - Upper Plenum Pressure 6
s b
RETRAN Channel Flow Rate (Ib/sec) m e cn g-
. k g . S. E g .
8 E 7 h C i e i i ct f,
- l -
[ !
,! -N
(
ok
.l ;
-g
- b. : -
^
3:
l ~
< m l a m 's E a: ',
[ l; l -
1
' --g $ )
[
s ;
4
'k' -0
's -
i
_ _g
~
o i i I d
[ q I
q m e m o O d d (UH/87W) eleu mold leuueyo eA!py guMaid
[.
(. FIGURE 5-3 GLRWOBP Results - Active Channel flow Rate 51 b
i
([
9 l
l
[- -
/ -
_9 l n o'
4 -
a' _o n 9
~
WE .
l
~
-N Ex ,'
^
ll ,/
-
- a-n 0
- k
(.' l
-l _9-g o ',
t %
'. _o
[.. .
l '
. 1
( . l
/
~
L _ m d
( l I o
' i s e d .
9 m m e s l o a
-o d h o d 8 i d ;
I uolpeJ3 pl0A lix3 leuuey0 FIGURE 5-4
(
GLWOBP Results - Omanal Exit Void 52 f
L-
',- i , ,
g ,,
~
I
~
8 I'.
. o l
E e
E
_ o v
_ O N
M a i e o S
- e v m o (41MW) JeMod leuueyo FIGURE 5-5 LOSC Boundary Conditim - Channel Power 53
1
. . . . i,..,i....i.,, i,,,. g
~
~
(.
(
-8
{. -
_ _g 1
i:
f -
-?
~
[ .
L _ _g I
, , , . ',, ,,.,,. ,,,,,.,,,,. o b
e b
b
- S S b (ejsd) eJnsseJd wnueld Jeddn FIGURE 5-6 k' LOSC Boundary Conditions - Upoer Plenum Pressure 54 I _ _ _ -_ _ _ -_-- _ - - _ - -- _
l W31WYN V09Ae gysuuot dlom sele )18/S( j I8 N 8
N t
it d'
o i i i i 8 l
l l
.l l
l -8 l
! M l
bb 5g m=
i:
I f
_ _g
_ _g t t f I o
b m a* o b o o 8 o' 1 8 9
)RH/BLM( etaR wolf lennahC evitcA 2RWBIF diDnHH S-4 1oga yasnI:s - yaps.o aynuuol dlom 8810 95
,8-pe -
e f- - _g
! !g@
1-l W>
~
-8 l Em n
l -
~
l i:
l o
l l
l l
l l -R l l L :
l l
l i i . , o e E e m s o 6 h 6 E.
O o o uogoey ploA Ux3 leuueyo l.
FIGURE 5-8 LOSC Results - Channel Exit Void Fraction 56
\
6.0 CONCLUSION
S The FIBWR2 code extends the BWR specinc models used in the FIBWR code for use in the evaluation of transient thermal hydraulic conditions. FIBWR2 includes the capability to explicitly model features of advanced BWR fuel designs such as axially variable flow areas and alternative j
- water tube designs and the FIBWR2 program structure allows the implementation of CPR correlations within the code. The applicability of the FIBWR2 models and equations and application of FIBWR2 within the YAEC BWR transient analysis methodology for evaluation of transient hot channel CPR performance has been demonstrated by:
Comparison of FIBWR2 solutions to two thermal hydraulic problems (exponential flow decay and sinusoidal flow variation) and their analytical solutions to demonstrate the basic transient hydraulic equations and solution techniques have been correctly applied.
I Comparison to steady state and transient experimental test results for a 16 rod test assembly configuration and to the results of similar simulations performed using the RETRAN/TCPYA01 models. These comparisons demonstrate both the applicability of l the FIBWR2 thermal hydraulic models to be used for the evaluation of transient thermal hydraulic conditions and the proper implementation of the vendor proprietary CPR correlation within the FIBWR program structure. Furthermore, sensitivity studies l performed for selected transient tests demonstrate the applicability of alternate ;
modeling choices available with FIBWR2 relative to the specification of nodalization, )
void model and time step sizes.
l Comparison with RETRAN/TCPYA01 results for several characteristic limiting BWR l reload transients. Relative to the RETRAN/TCPYA01 results, the FIBWR2 results produce a smaller change (0.02 to 0.05) in the CPR for pressurization events (e.g.
turbine trip) and a larger change (0.01) in CPR for increased subcooling events.
l These differences in the CPR results are due to the more accurate modeling of the two phase form loss multiplier in FIBWR2 and the resulting changes in active coolant !
flow. l 57 L
These comparisons demonstrate the applicability of the FIBWR2 code to the' analysis of transient thermal hydraulic conditions during Abnormal Operational Transients. In addition, the' -
comparison of FIBWR2 results with RETRANfrCPYA01 for both the Vermont Yankee limiting i reload trensients and the 16 rod test assembly data specifically demonstrates the applicability of FIBWR2 based models as a revised approach in the YAEC BWR reload analysis methodology for hot channel CPR evaluation, l
i 58
l
~
. 7.0 . REFERENCES
- 1. ;EPRI NP-1923, FIBWR: "A Steady State Core Flow Distribution Code for Boiling Water Reactors, Code Verification and Qualification Report," July 1981.
- 2. EPRI NP-1924-CCM, "FIBWR: A Steady State Core Flow Distribution Code for Boiling Water Reactors, Computer Code _ User's Manual," July 1981.
- 3. GEAP-1022-II, " Deficient Cooling, lith Quarterly Progress Report, January 1 - March 31, 1972," General Electric Co., April 1972.
- 4. . GEAP-1022-12, " Deficient Cooling,12th Quarterly Progress Report, April 1 - June 30, 1972," General Electric Co., July 1972.
- 5. GEAP-13295, " Transient Critical Heat Flux - Experimental Results," General Electric Co.,
September 1972.
- 6. GEAP-21207, "BWR 8X8 Fuel Rod Simulation Using Electrical Heaters," General Electric, March 1976.
'7. NVY 82-157, USNRC to Vermont Yankee Nuclear Power Corporation, " Acceptance of YAEC BWR Reload Methods," September 15,1982.
- 8. Letter, USNRC to Vermont Yankee Nuclear Power Corporation, " Acceptance for Referencing of Topical Report YAEC-1693, ' Application of One-Dimensional Kinetics to Vermont Yankee Transient Analysis Methods,'" December 18,1989.
- 9. YAEC-1935, " Vermont Yankee Cycle 19 Core Performance Analysis Report," October 1996.
- 10. YAEC-1299P, " Methods for the Analysis of Boiling Water Reactors Transient Critical Power Ratio Analysis RETRAN - TCPYA01," Yankee Atomic Electric Company, March 30,1982. ,
- 11. EPRI NP-1580, "RETRAN - A Program for One-Dimensional Transient Thermal Hydraulic Analysis of Complex Fluid Flow Systems," December 1978.
- 12. NU. REG CR/4391, EGG-2680, ' TRAC-BFI/ MODI Models and Correlations," August 1992.
- 13. YAEC-1340," Application of the FIBWR2 Core Hydraulics Code to BWR Steady State Thermal Hydraulic Analysis," To Be issued.
59 l
l .. .
_ _ _ _ _ _ - _ - _ _ _ _ _ _ _