View source for L-PI-04-017, Development and Qualification of a Gothic Containment Evaluation Model for the Prairie Island Nuclear Generating Plants.

{{Adams
| number = ML042570292
| issue date = 04/30/2004
| title = Development and Qualification of a Gothic Containment Evaluation Model for the Prairie Island Nuclear Generating Plants.
| author name = Ofstun R
| author affiliation = Westinghouse Electric Co, LLC
| addressee name = 
| addressee affiliation = NRC/NRR
| docket = 05000282, 05000306
| license number = DPR-042, DPR-060
| contact person = 
| case reference number = L-PI-04-017
| document report number = WCAP-16219-NP
| document type = Topical Report
| page count = 135
}}

=Text=
{{#Wiki_filter:Westinghouse Non-Proprietary Class 3 WCAP-16219-NP                                       April 2004 Development and Qualification of a GOTHIC Containment Evaluation Model for the Prairie Island Nuclear Generating Plants g    Westinghouse

WESTINGHOUSE NON-PROPRIETARY CLASS 3 WCAP-16219-NP Development and Qualification of a GOTHIC Containment Evaluation Model for the Prairie Island Nuclear Generating Plants R. Ofstun Systems and Safety Analysis April 2004 Westinghouse Electric Company LLC P.O. Box 355 Pittsburgh, PA 15230-0355 0 2004 Westinghouse Electric Company LLC All Rights Reserved 6403-NP.doc-042904

iii TABLE OF CONTENTS LIST OF TABLES ................
iii LIST OF FIGURES .................                                                                                                                iii LIST OF ACRONYMS .................................................................................... ..................................
iii I          INTRODUCTION ................................................................                                                       1-3 2         DESCRIPTION OF THE PRAIRIE ISLAND GOTHIC CONTAINMENT EVALUATION MODELS                                 ..                                                                                 2-3 2.1      NODING STRUCTURE ............................                                                                               2-3 2.2      VOLUME INPUT ............................                                                                                   2-3 2.3      INITIAL CONDITIONS ............................                                                                             2-3 2.4      FLOW PATHS ............................                                                                                    2-3 2.5      HEAT SINKS ............................                                                                                    2-3 2.6      MODELING SUMP RECIRCULATION ............................                                                                   2-3 2.7      BOUNDARY CONDITIONS ............................                                                                           2-3 3         COMPARISON OF MODELING ASSUMPTIONS FOR CONTAINMENT ANALYSES ........ 3-3 3.1      CODE COMPARISON ................................................................                                           3-3 3.2      INPUT COMPARISON ................................................................ 3-3 4         BENCHMARK COMPARISON ................................................................                                               4-3 4.1      LOCA CONTAINMENT RESPONSE BENCHMARK COMPARISON ....................... 4-3 4.2      MSLB CONTAINMENT RESPONSE BENCHMARK COMPARISON ...................... 4-3 5          SAMPLE TRANSIENT RESULTS USING THE PRAIRIE ISLAND GOTHIC CONTAINMENT EVALUATION MODEL.                                                    ..................................... 5-3 5.1      LOCA SAMPLE TRANSIENT RESULTS ......................................                                                       5-3 5.2      MSLB SAMPLE TRANSIENT RESULTS ......................................                                                       5-3 6         

==SUMMARY==
OF RESULTS AND CONCLUSIONS .                                                  ..................................... 6-3 7         REFERENCES .......................................                                                                                  7-3 APPENDIX A           GOTHIC PRAIRIE ISLAND MSLB CONTAINMENT RESPONSE MODEL INPUT TABLES .......................................                                                               A-3 APPENDIX B           GOTHIC PRAIRIE ISLAND LOCA CONTAINMENT RESPONSE MODEL INPUT TABLES........................................................................................... B-3 April 2004 WCAP-16219-NP 64033NNPdod0c20424 W6CA3NPd1629O                                                                                                                         April 2004

v LIST OF TABLES Table 2-1          Heat Sink Data for the Prairie Island Containment Model ....................................... 2-3 Table 2-2          Thermal Conductivity, Density, and Volumetric Heat Capacity (VHC) for Heat Sink Materials ............                                                                    2-3 Table 3-1          Input Comparison ............                                                                       3-3 April 2004 WCAP-l 6219-NP                                                                                                  April 2004 6403-NP.doc-042904

vii LIST OF FIGURES Figure 2-1         Prairie Island GOTHIC Containment Model for MSLB Events ...................................... 2-3 Figure 2-2         Prairie Island GOTHIC Containment Model for LOCA Events ...................................... 2-3 Figure 3-1         Fan Cooler Heat Removal Rate per Train .....................................................                3-3 Figure 4-1         Containment Pressure Comparison for LOCA Benchmark Case .................................... 4-3 Figure 4-2         Containment Temperature Comparison for LOCA Benchmark Case ................                  .............. 4-3 Figure 4-3         Containment Sump Temperature Comparison for Benchmark LOCA Case .................... 4-3 Figure 4-4         Containment Pressure Comparison for LOCA Benchmark Case with Input Change to Simulate the CONTEMPT Temperature Flash Model .4-3 Figure 4-5         Containment Temperature Comparison for LOCA Benchmark Case with Input Change to Simulate the CONTEMPT Temperature Flash Model .4-3 Figure 4-6         Containment Sump Temperature Comparison for LOCA Benchmark Case with Input Change to Simulate the CONTEMPT Temperature Flash Model .................. 4-3 Figure 4-7         Containment Pressure Comparison for MSLB Benchmark Case .................................... 4-3 Figure 4-8         Containment Temperature Comparison for MSLB Benchmark Case .............................. 4-3 Figure 4-9         Containment Sump Temperature for MSLB Benchmark Case ........................................ 4-3 Figure 4-10        Containment Pressure Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option .4-3 Figure 4-11        Containment Temperature Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option .4-3 Figure 4-12        Containment Sump Temperature Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option .4-3 Figure 4-13        Containment Pressure Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option and 100% Spray Efficiency .4-3 Figure 4-14        Containment Temperature Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option and 100% Spray Efficiency .4-3 Figure 4-15        Containment Sump Temperature Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option and 100% Spray Efficiency ............................................                                     4-3 Figure 5-1         Containment Pressure Transient Response for LOCA ............................................               5-3 Figure 5-2         Containment Temperature Transient Response for LOCA ............................................            5-3 Figure 5-3         Containment Sump Temperature Transient Response for LOCA .................................... 5-3 April 2004 WCAP-l 621 9-NP                                                                                                  April 2004 6403.NP.doc.042904

viii LIST OF FIGURES (cont.)
Figure 54         Containment Pressure Transient Response for MSLB ........................................        5-3 Figure 5-5        Containment Temperature Transient Response for MSLB ........................................ 5-3 Figure 5-6         Containment Pressure Transient Response for MSLB ........................................        5-3 WCAP-16219-NP                                                                                                April 2004 6403-NP.doc-042904

ix LIST OF ACRONYIS CCW                component cooling water CS                 containment spray DBA                design basis accident DLM                diffusion layer model EPRI               Electric Power Research Institute FRV                feedwater regulator valve HTA                heat transfer area HX                 heat exchanger LCO                limiting condition of operation LOCA               loss-of-coolant accident M&E                mass and energy MDLM               mist diffusion layer model MSLB               main steamline break NAI                Numerical Applications Incorporated NMC                Nuclear Management Company NRC                Nuclear Regulatory Commission NSP                Northern States Power RAI                Request for Additional Information RCS                reactor coolant system RH                 relative humidity RHR                residual heat removal USAR               Updated Safety Analysis Report VHC                volumetric heat capacity WCAP-16219-NP                                            April 2004 6403-NP.doc-042904

I-I 1        INTRODUCTION The GOTHIC program is rapidly becoming the industry standard for performing both inside and outside containment pressure and temperature design basis analyses. The program is being developed by Numerical Applications Incorporated (NAI) with funding by the Electric Power Research Institute (EPRI). The GOTHIC program consists of a pre-processor, which helps with input generation; a solver, which performs the calculations; and a post-processor, which produces the output tables and plots. The GOTHIC Technical Manual (Reference 1) provides a description of the governing equations, constitutive models, and solution methods in the solver. The GOTHIC Qualification Report (Reference 2) provides a comparison of the solver results with both analytical solutions and experimental data. The GOTHIC User Manual (Reference 3) provides information to help the user develop models for various applications.
Nuclear Management Company (NMC) approved a plan to develop a GOTHIC containment evaluation model for the Prairie Island Nuclear Generating Plants. The GOTHIC model would replace the current CONTEMPT containment evaluation model and be used to perform the containment peak pressure and temperature analyses, equipment qualification analyses, and peak liner temperature calculation.
The Prairie Island GOTHIC containment evaluation model was constructed based on the recently accepted Kewaunee GOTHIC containment evaluation model, including the limitations and conditions specified in the associated Nuclear Regulatory Commission (NRC) Safety Evaluation Report (Reference 4). Both the Kewaunee and Prairie Island containment evaluation models use the diffusion layer model (DLM) heat and mass transfer option in GOTHIC.
The Kewaunee containment evaluation model was developed with GOTHIC version 7.Op2 and submitted for NRC review. The NRC conditions for acceptance, primarily the change from the mist diffusion layer model (MDLM) to the DLM correlation, along with other improvements and error corrections, were incorporated into GOTHIC version 7.1pl and subsequent analyses for the Kewaunee Updated Final Safety Analysis Report were performed with that version. The differences between GOTHIC version 7.0 and 7.1 are documented in Appendix B of the GOTHIC User Manual (Reference 3).
The Prairie Island GOTHIC containment evaluation model was compared with the CONTEMPT containment evaluation model to identify modeling differences. Changes were then made to the GOTHIC containment evaluation model to create a model for benchmark comparison with the CONTEMPT model.
The GOTHIC benchmark model results were compared with the CONTEMPT model results for both a loss-of-coolant accident (LOCA) and main steamline break (MSLB) transient.
Finally, representative LOCA and MSLB cases were run using the Prairie Island containment evaluation model to demonstrate that the new evaluation model produces an acceptable (less than containment design pressure), conservative analysis result.
April 2004 WCAP-16219-NP                                                                                   April 2004 6403-NP.doc-042904

2-1 2         DESCRIPTION OF THE PRAIRIE ISLAND GOTHIC CONTAINMENT EVALUATION MODELS The modeling assumptions used in the Prairie Island containment evaluation models for LOCA and MSLB events are described below. The Prairie Island containment evaluation models are based on the previously accepted Kewaunee containment evaluation models (Reference 4) and use GOTHIC version 7.1patchl with the DLM heat and mass transfer correlation. The GOTHIC input tables for the LOCA and MSLB sample cases presented in Section 5 are provided in Appendices A and B, respectively.
2.1       NODING STRUCTURE The Prairie Island GOTHIC containment model for an MSLB event is shown in Figure 2-1. The containment (Volume 1) is modeled with a single lumped parameter node. Two boundary conditions (I F and 2F) are used to represent the sources of mass and energy (M&E) from the break. The spray injection system is modeled with boundary condition (3F). Flow paths connect the boundary conditions to the containment volume. Twenty-four heat sinks and a fan cooler component are also shown.
The Prairie Island GOTHIC containment model for a LOCA event is shown in Figure 2-2. Additional boundary conditions, volumes, flow paths, and components are used to model accumulator nitrogen release and sump recirculation. The recirculation system model uses GOTHIC component models for the residual heat removal (RHR) and component cooling water (CCW) heat exchangers and the CCW pumps.
Recirculation flow from the sump is modeled using a boundary condition.
2.2       VOLUME INPUT GOTHIC requires the volume, height, diameter, and elevation input values for each node. The containment free volume specified in the Prairie Island Updated Safety Analysis Report (USAR) was used as input for the volume of node 1. Standard values of 100 feet, 10 feet, and 0 foot were used for the GOTHIC height, diameter, and elevation input values; the model is not sensitive to these representative input values.
A conservatively calculated pool surface area is used to model interfacial heat and mass transfer to the liquid pools on the various floor surfaces inside containment. The conductors representing the floors are essentially insulated from the vapor afterthe pools develop; however, there can still be condensation or evaporation from the surface of the liquid pools. The pool area input value represents the sum of the three floor conductor surface areas (22, 23, and 24) from Table 2-1. Using this method to model the interfacial heat and mass transfer between the pools and the atmosphere was previously approved by the NRC for the Kewaunee containment design basis accident (DBA) and equipment qualification analyses (Reference 4).
The LOCA containment response model includes an additional volume to represent the nominal volume of pressurized nitrogen gas in the accumulators that is released to the containment after the accumulator water volume has been discharged. Standard values of 20 feet, 20 feet, and 0 foot were used for the GOTHIC height, diameter, and elevation input values for the accumulator volume, respectively; the model is not sensitive to these representative input values.
WCAP-16219-NP 6403-NP.doc-042904

J   a 2-2 The LOCA containment response model input values for the RHR and CCW systems' volume, height, diameter, and elevation are not important for modeling the sump temperature response after recirculation.
Standard values of 100 feet 3 , 10 feet, 10 feet, and 0 foot were used, respectively; the model is not sensitive to these representative input values.                                                                   l 2.3       INITIAL CONDITIONS IL The containment initial conditions are listed below:
* Pressure - 16.7 psia (biased high for peak pressure and liner temperature calculation)                  l
                    - 14.2 psia (biased low for peak vapor temperature calculation)
* Relative Humidity - 30% (biased low for peak pressure and liner temperature calculation)                I
                              - 100% (biased high for peak vapor temperature calculation)                         I
* Temperature - 120'F                                                                                     I The LOCA containment response model contains volumes representing the accumulators, the RHR system, and the CCW system. The initial pressure of the volume representing the accumulators (366 psia)            A was calculated assuming an ideal gas law expansion of the initial accumulator gas volume to the full accumulator volume. The gas temperature was set to a conservatively high value (2601F). The RHR                  ],
system volumes were initially filled with hot (200'F) water at the containment pressure. Most of the              1 CCW system volumes in the LOCA containment response model were also initially filled with hot water,              I, but at a higher pressure of 60 psia; the CCW surge tank was modeled as half full.                                  I 2.4       FLOW PATHS                                                                                               I, Flow paths connect the boundary conditions to the containment volume. The flowrate is specified by the boundary condition, so most of the flow path input is not important. Standard values are used for the area, hydraulic diameter, friction length, and inertia length of the flow path. The elevation of the break       \I flow paths is arbitrarily set to 30 feet and the elevation of the spray flow path is arbitrarily set to 90 feet.
Flow path 14 in the LOCA containment response model represents the accumulator injection lines to the reactor coolant system (RCS). The valve in this flow path is opened to model the release of the accumulator nitrogen gas to the containment after the accumulator water volume is depleted. This is different from the Kewaunee model, in which the accumulator gas was released at a specified rate over a fixed time period using a boundary condition. The flow path friction length, inertia length, hydraulic diameter, and flow area input values are based on actual line data.
Flow path 6 in the LOCA containment response model connects the containment volume to the RHR                       J system model. The valve in this flow path is opened at the start of recirculation. The upstream elevation of this flow path was set to the bottom of the containment (0.0 foot) and the downstream elevation was set to the middle of the RHR piping volume (-5.0 feet). The height, area, hydraulic diameter, friction, and inertial length input values for this flow path were arbitrarily set to 0.1 foot, 1.0 foot 2, 1.0 foot, 1.0 foot 'I and 1.0 foot, respectively.                                                                                         I WCAP-1 6219-NP 6403-NP.doc-042904

2-3 2.5       HEAT SINKS The CONTEMPT model documented in Reference 5 provides data for 24 containment heat sinks. The heat sink data for the Prairie Island GOTHIC containment model is based on conservatively low surface areas and nominal thicknesses and is summarized in Table 2-1. Since the heat sink area input provided in CONTEMPT represents the total exposed surface area of the heat sink, the actual thickness of the two-sided heat sinks was divided by 2 to preserve the heat sink volume, and the Volume B boundary was modeled as adiabatic for the GOTHIC model.
A thin air gap is assumed to exist between the steel and concrete for steel-jacketed heat sinks in the GOTHIC containment evaluation model; this gap is not modeled in the CONTEMPT containment evaluation model. A standard value of 100 Btu/hr-ft2~-F is assumed for the minimum gap conductance.
The gap width is determined by dividing the gap thermal conductivity by the gap conductance.
The volumetric heat capacity and thermal conductivity for most of the heat sink materials are taken from Reference 5 as summarized in Table 2-2 and taken from engineering publications. Values for the density of concrete and steel were taken from Reference 6 and the specific heat values were calculated based on the volumetric heat capacity values specified in Reference 5.
The density and specific heat input values for the paint, the zinc coating on galvanized steel, and the air gap (between the concrete and steel) were not specified in Reference 5. The specific heat for the paint was assumed to be 1.0 BtuAbm-0 F; the density was calculated based on the volumetric heat capacity given in Reference 5. The thermal properties for zinc and air were taken from Reference 6 at a temperature of approximately 200'F.
Heat and Mass Transfer Correlations GOTHIC has a number of heat transfer coefficient options that can be used for containment analyses.
These include the film, direct, Tagami, and user-specified heat transfer coefficient options. GOTHIC also has a number of condensation options that can be used for containment analyses. These include the Uchida, Gido-Koestel, and DLM options.
The direct heat transfer coefficient set is used, along with the DLM mass transfer correlation, for all of the heat sinks inside containment. This heat and mass transfer methodology was recently reviewed and approved for use in the Kewaunee containment DBA analyses (Reference 4). The DLM correlation does not require the user to specify a revaporization input value, as was done in previous analyses using the Uchida correlation.
The direct heat transfer coefficient set (with split option) is used for the heat sinks representing floors.
The split option will shut off heat transfer to the floor once the containment liquid volume fraction exceeds 0.0001.
Insulated surfaces are modeled with a constant (0.0 Btu/hr-ft2 -F) heat flux.
WCAP-1621 9-NP 6403-NP.doc-042904

24 2.6       MODELING SUMIP RECIRCULATION
[
                                                                  ]ac 2.7       BOUNDARY CONDITIONS The LOCA and MSLB transient M&E releases are calculated separately and input to the GOTHIC                    4 containment models via boundary conditions. The break mass and enthalpy are input to the containment         I model through forcing functions on flow boundary conditions IF and 2F.                                       It1 This report addresses the applicability of GOTHIC for Prairie Island and is not intended to address other methods used to determine GOTHIC inputs such as the M&E releases. The M&E release inputs used in the benchmarking cases (Sections 4.1 and 4.2) and sample cases (Sections 5.1 and 5.2) were calculated specifically for Prairie Island and are considered representative of Prairie Island and thus, sufficient for the purposes of this report.
JI The M&E input used in the GOTHIC and CONTEMPT LOCA benchmarking case is from the existing Prairie Island licensing basis LOCA containment response, as documented in Appendix K of the Prairie Island USAR. The M&E input used in the LOCA GOTHIC containment evaluation model sample case was calculated using the Westinghouse method for LOCA M&E, as described in WCAP-10325 WCAP-16219-NP 6403-NP.doc-042904

2-5 (Reference 7). While the M&E release input used in the containment evaluation model sample case is not presently the licensing basis LOCA release input, it is considered adequately conservative, yet representative of actual expected LOCA release input, for Prairie Island.
The M&E input used in the GOTHIC and CONTEMPT MSLB benchmarking case is from the existing Prairie Island licensing basis MSLB containment response, as documented in Section 14.5.5.3.1 of the Prairie Island USAR (Reference 8). The M&E used as input into the MSLB GOTHIC containment evaluation model sample case was calculated using the Westinghouse method for MSLB M&E, as described in WCAP-8822 (Reference 9). While the M&E release input used in the containment evaluation model sample case is not presently the licensing basis MSLB release input, it is considered adequately conservative, yet representative of actual expected MSLB release input, for Prairie Island.
The liquid portion of the break flow is released as drops with an assumed diameter of 100 microns (0.00394 inches). This is consistent with the methodology approved for Kewaunee (Reference 4) and is based on data presented in Reference 10.
The containment fan coolers are modeled with a cooler component. There are two fan coolers per train and two trains are normally available. The fan coolers are modeled to actuate on the containment HI-1 pressure setpoint (including biased high uncertainty) and to begin removing heat from the containment after an accumulation of conservatively long delays. The fan cooler heat removal rate per train is given as a function of containment temperature in Reference 5 and is shown in Figure 3-1. This fan cooler heat removal rate is conservatively low to account for the two-phase flow analysis that was completed to address NRC Generic Letter 96-06. The heat removal rate is read into a GOTHIC function and a multiplier, based on the number of fan cooler trains running, is used to calculate the heat removal rate from containment.
Containment spray is modeled with a flow boundary condition. There is one spray pump per train and two trains are normally available. The spray pumps are modeled to actuate on the containment HI-2 setpoint (including biased high uncertainty) and to begin injecting water from the spray header after an accumulation of conservatively long delays. The spray flowrate is biased low and entered into GOTHIC as a function of time. From Section 6.4.2.2.2 of Reference 8, the Prairie Island containment spray nozzles produce a mean drop diameter of 700 microns. Therefore, in the GOTHIC model, the spray is injected in the form of drops with a conservatively large diameter of 1,000 microns (0.0394 inches). A function multiplier is used to convert the input volumetric spray flowrate (in gpm) to units of cubic feet per second and factors in the number of trains running.
WCAP-16219-NP 6403-NP.doc-042904

2-6 Table 2-1        Heat Sink Data for the Prairie Island Containment Model Actual               Actual                      Actual   Model Area                               Paint       Primer          Zinc'       Stainless Steel       Carbon Steel           Air     Concrete  Total Conductor          (fty)            Sides             (in.)         (in.)         (in.)             (in.)               (in.)            (in.)      (in.)    (in.)
1          40,300                I             0.011                                                                1.5                               1.5110 2            17,300               1             0.011                                                              0.75                               0.7610 3            1,260               02             0.011                                                             0.1875           0.0021       12    12.2006 4            5,150                1                                                          0.1875                                0.0021       12    12.1896 5            1,010                1                                                            0.25                                0.0021       12    12.2521 6            5,780                1             0.011                                                              0.375                              0.3860 7            4,055               2              0.011                                                              0.336                               0.1790 8            16,925              2              0.011                                                               0.5                               0.2610 9           28,500               2              0.011                                                              0.75                               0.3860 10           2,000               2              0.011                                                                1.5                              0.7610 l    l        500                2              0.011                                                                2                                 1.0110 12            1,695               2             0.011                                                                0.5                               0.2610 13            12,400              2             0.011                                                              0.1875                              0.1048 14           14,000               1                                         0.004                                   0.1                               0.1040 15          27,100               2                                          0.004                                 0.1084                              0.0582 16            1,346                1            0.011                                                               2.75                               2.7610 17            1,256                1            0.011                                                                1.39                               1.4010 18             973                2             0.011                                                              0.625                               0.3235 19            19,320              2                                                                                         .                    12    6.0000 20           15,000                1                                                                                                             12    12.0000 21            3,810               2             0.028                                                                                            12    6.0280 22           3,322.5              i3            0.028                                                                                            12     12.056 23            5,830                1            0.028                                                                                            12    12.0280 24           7,687.5       _         _          0.028                                                                                             6     6.056 Notes:
: 1. The zinc plating thickness was not dcefincd in Refercnce 5, it was estimated to bc 4 mils.
: 2. Rcfcrencc 5 states that the reactor vessel liner (conductor 3) is not to be modeled as a heat sink, so both surfaces were modeled as insulated.
: 3. This heat sink represents a floor on one side and a ceiling on the other side. The area of one side is given.
WCAP-16219-NP                                                                                                                                                 April 2004 6403-NP.doc-042904 i                                         -e-w  f,                /e w         " /e          #         t,                           we                                    e

2-7 Table 2-2         Thermal Conductivity, Density, and Volumetric Heat Capacity (VHC) for Heat Sink Materials Conductivity           Density               VHC              Specific Heat
_                           (Btu/hr-ft-0 F)         (Ibm/cf)           (Btu/cf- 0 F)       (Btu/lbm-0 F)
Paint- Steel                     0.29                 28.0                 28                    1.0 Paint - Concrete                 0.29                 28.0                 28                    1.0 Concrete                         0.8                   144                28.8                   0.2 Carbon Steel                      26                  490                 56.4                 0.115 Stainless Steel                  9.4                  488                 60.12                0.123 Zinc                              64                  446                 N/A                  0.091 Air(gap)                       0.0174                 0.06                N/A                  0.241 WCAP-16219-NP                                                                                     April 2004 6403-NP.doc.042904

2-8 Figure 2-1 Prairie Island GOTHIC Containment Model for MSLB Events WCAP-16219-NP                                                                      April 2004 6403-NP.doc-042904

2-9 ac Figure 2-2 Prairie Island GOTHIC Containment Model for LOCA Events April 2004 WCAP-l    621 9-NP WCAP-16219-NP                                                                      April 2004 6403-NP.doc-042904

3-1 3         COMPARISON OF MODELING ASSUMPTIONS FOR CONTAINMENT ANALYSES The Prairie Island licensing basis containment response analyses had been performed previously with the CONTEMPT code. This section presents a comparison of the differences between the CONTEMPT and GOTHIC containment models.
3.1       CODE COMPARISON The GOTHIC and CONTEMPT code features and modeling techniques are compared in Reference 11.
The key differences are summarized below:
* GOTHIC was designed as a general purpose thermal-hydraulic transient analysis code.
CONTEMPT was specifically designed for containment analyses and contains four special hard-wired volume types for this purpose.
* GOTHIC performs M&E balances for the droplet phase; the spray efficiency is a user input in CONTEMPT. Similarly, GOTHIC performs M&E balances for steam/air bubbles rising through a pool; CONTEMPT assumes 100 percent condensation efficiency.
* For interface heat and mass transfer, CONTEMPT assumes the interface is always at saturation.
GOTHIC calculates the interface temperature from first principles. Spray droplets condense steam based on the difference in steam partial pressure at the drop surface and in the atmosphere.
The droplets can eventually come to thermal equilibrium with the containment atmosphere and fall into the sump.
* CONTEMPT uses a user-specified phase separation model for blowdown (pressure flash or temperature flash). GOTHIC determines the blowdown phase separation based on fundamental models for interface heat and mass transfer at the containment conditions. Liquid break flow released to the containment evaporates based on the difference in steam partial pressure at the drop surface and in the atmosphere. The droplets can eventually come to thermal equilibrium with the containment atmosphere and fall into the'sump.
* GOTHIC models momentum balanced, multi-phase flow between compartments and to the atmosphere. GOTHIC uses break flow tables or correlations for choked flow. CONTEMPT only models the vapor phase flow between compartments and to the atmosphere.
* GOTHIC calculates the tube and shell heat transfer coefficients for heat exchangers. The overall heat transfer coefficient is a user input for CONTEMPT heat exchangers.
* Each version of GOTHIC is qualified against a wide range of tests (both experimental and analytical) and the results are documented in the code Qualification Report (Reference 2).
Although CONTEMPT results from various versions have been compared with various experimental and analytical tests, there is no corresponding documentation to demonstrate the qualification of each version of CONTEMPT.
April 2004 WCAP-l 621 9-NP                                                                                     April 2004 6403-NP.doc-042904

3-2 3.2       INPUT COMPARISON                                                                                       L Table 3-1 presents a comparison of the key input values for the CONTEMPT and GOTHIC containment                  I models. The key differences between modeling input assumptions are summarized below:                             L
* The DLM is used to calculate condensation mass transfer between the heat sinks and the atmosphere in the Prairie Island GOTHIC containment evaluation model. The DLM model is                'I}
described in Reference I and the qualifications for its use in the containment design basis            I, analyses are described in References 2 and 3. In addition to the results presented in the GOTHIC      IL code Qualification Report, the DLM correlation was compared with additional test data in response to a Request for Additional Information (RAI) from the NRC on the Kewaunee submittal (Reference 12). CONTEMPT uses the Tagami and Uchida heat and mass transfer                   I correlations.
* GOTHIC uses something similar to the CONTEMPT pressure flash assumption to model phase separation for the break flow. The liquid portion of the break flow is released as 100 p droplets that come into temperature equilibrium with the atmosphere before falling into the sump.                i
* A larger pool surface area for interface heat and mass transfer is modeled in GOTHIC. The GOTHIC pool surface area input value is the sum of the horizontal floor heat sink areas.
Condensation occurs when the containment steam partial pressure is higher than the saturation pressure at the pool surface; evaporation occurs when the containment steam partial pressure is         t lower than the saturation pressure at the pool surface.                                                I
* The accumulator nitrogen release is modeled for the LOCA event in the GOTHIC containment              I model. This was not considered in the CONTEMPT model for the LOCA event.                               I.
* GOTHIC models the fan cooler and spray actuation based on the containment pressure setpoint.          ]
Both CONTEMPT and GOTHIC model delays for containment heat removal from the fan coolers and spray. The fan cooler heat removal delay is modeled to account for signal generation, diesel      i generator and fan startup, and other factors. The spray injection delay is modeled to account for       9, signal generation, diesel generator and pump startup and filling of the spray lines.                    L
* GOTHIC models the incoming spray as droplets with a 1,000 p diameter. GOTHIC models the                J interface heat and mass transfer between the drops and the atmosphere. CONTEMPT models a             .1' fixed spray efficiency rate.                                                                           1
* The long-term LOCA heat removal rate from the RHR and CCW heat exchangers is flow and                 IL temperature dependent. GOTHIC calculates the overall heat transfer coefficient for each type of      '
heat exchanger based on the flowrates, number of tubes and tube hydraulic diameter. The              J.
CONTEMPT model uses a fixed UA for the RHR heat exchanger and assumes a conservative                 i CCW flowrate and temperature.
Some of the differences between the modeling assumptions can be eliminated by input changes in the              A, GOTHIC model. The impact of these changes is shown in the benchmark tests described in Section 4.                 1, April 2004 WCAP- 16219-NP WCAP-1    6219-NP                                                                                   April 2004 6403-NP.doc-042904

3-3 Table 3-1         Input Comparison Description                      CONTEMPT                        GOTHIC Containment Data Noding Structure                      Single lumped                  Single lumped Volume (cf)                           1,320,000                      1,320,000 Height (fi)                           N/A                            100 Pool Area (sqft)                      8,659(')                       16,840 Heat Sink Geometry Area           Thickness     Area              Thickness Number              Description      (sqft)   Sides  (inches)    (sqft)      Sides    (inches) 1(2)          Containment Cylinder   40,300      1           1.5  40,300        1         1.511 2(2)          Containment Dome       17,300      1         0.75   17,300        1         0.761 3P)           Vessel Cavity Liner      1,260     0            12    1,260       0     12.20059 4(4)           Refueling Canal          5,150     1     12.1875     5,150        1     12.18959 5(4)          Refueling Canal          1,010     1        12.25     1,010       1     12.25209 6(2)          Exposed Piping           5,780     1       0.375     5,780        1         0.386 7(2).(S)       Steel Structures         4,055     2       0.336     4,055        2         0.179 8(2).(5)      Steel Structures       16,925      2          0.5   16,925        2         0.261 9(2),(5)       Steel Structures       28,500      2         0.75   28,500        2         0.386 lo(2).(5)     Steel Structures         2,000     2           1.5   2,000        2         0.761 Steel Structures           500     2          2.0       500       2         1.011 12   Handrails & Ladders      1,695     2          0.5     1,695       2         0.261 13(2)(5)      Grating                12,400      2      0.1875    12,400        2      0.10475 14(6)         Cable Trays & Conduit  14,000      1          0.1   14,000        1         0.104 15(5) (6)     Ductwork               27,100      2      0.1084    27,100        2        0.0582 16(2)         Accumulators (2)         1,346     1         2.75     1,346       1         2.761 17(2)         Accumulators (2)         1,256     1         1.39   - 1,256       1         1.401 I8(2),(5)     Accumulators (2) skirt     973     2       0.625        973       2        0.3235 19(5)         Heavy Walls            19,320      2            12  19,320        2            6.0 20            Heavy Walls            15,000      1            12  15,000        1           12.0 21(5)_(v)     Heavy Walls              3,810     2            12   3,810        2         6.028 22(51(7).(12)  Heavy Floors             6,645     2            12 3,322.5        1        12.056 23'           Heavy Floors             5,830     1            12   5,830        1        12.028 24(5),M,(12)   Light Floors           15,375      2             6 7,687.5        1         6.056 WCAP-1621 9-NP                                                                              April 2004 6403-NP.doc.042904

3-4 Table 3-1         Input Comparison (cont.)
Description                         CONTEMPT                                 GOTHIC Heat Sink Material Properties Thermal               Volumetric             Thermal            Volumetric Conductivity              Heat              Conductivity             Heat Material             (Btu/hr-ft-0 F)          Capacity           (Btu/hr-ft-0 F)         Capacity Paint - Steel                         0.29                    28                  0.29                  28 Paint - Concrete                      0.29                    28                  0.29                  28 Concrete                              0.8                   28.8                   0.8                 28.8 Carbon Steel                           26                    56.4                  26                  56.4 Stainless Steel                       9.4                   60.12                  9.4                60.12 Zinc                                  N/A                   N/A                    64                 40.59 Air (gap)                             N/A                   N/A                  0.0174              0.01446 Heat Transfer Coefficients  - LOCA      Tagami/Uchida                          Direct + DLM Heat Transfer Coefficients - MSLB       Uchida                                 Direct + DLM Revaporization Fraction                 0.08                                   N/A Initial Conditions Initial Pressure (psia)(8)              14.2, 16.85                            14.2, 16.70(9)
Initial Temperature   (0 F)              120                                   120 Initial Humidity- LOCA (% RH)           30                                     30 Initial Humidity- MSLB (% RH)           30                                     100("&deg;)
Boundary Conditions Break Flow Phase Separation             Liquid released as a stream, uses a    Liquid released as 100 [L diameter temperature flash assumption           drops, uses a pressure flash assumption Accumulator Nitrogen Release            Not Modeled                            Modeled for LOCA Fan Cooler Initiation                   60 seconds after event initiation(i)   19.7 psia with a 60-second delay Fan Cooler Heat Removal Rate            See Figure 3-1                         See Figure 3-1 Spray Flow Initiation - LOCA            75 seconds after event initiation      38.7 psia with a 72-second max.
l_                                      delay Spray Flow Initiation - MSLB            80 seconds after event initiation      38.7 psia with a 73-second max.
l_                                      delay Spray Flowrate (gpm/pump)               1,200                                  1,200 Spray Drop Modeling                     100% efficiency                        1,000 g drop diameter WCAP-16219-NP                                                                                               April 2004 6403-NP.doc-042904

3-5 Table 3-1           Input Comparison (cont.)
Description                                  CONTEMPT                                     GOTHIC LOCA Sump Recirculation Modeling Transfer to Recirculation after event             3,000                                      840/1,680 (max)( 13) initiation (seconds)                                                                         1,440 (min)
RHR Flowrate (Ibm/s)                              208                                        277.9/555.8 (max)( 13) 207.3 (min)
RHR Heat Exchanger UA                             689,842                                    Code calculated using flow area, (Bt/hr-F)                                                                                    Dh and HTA of HX CCW Flowrate (lbm/s)                              313.9                                      305.5 CCW Heat Exchanger UA                             N/A(t15)                                   Code calculated using flow area, l (Btu/hr-__)_                                                                                 Da and HTA of HX Other CCW Heat Loads (Btu/hr)                     N/A(' 5)                                   500,000 Service Water Flowrate (lbm/s)                    N/A(15)                                    275 Service Water Temperature (OF)                    N/A(15)                                    95(14)
Notes:
(1) The pool evaporation/condensation modeling is conservatively used for the CONTEMPT MSLB limiting temperature case but is not used for the CONTEMPT MSLB limiting pressure case or the LOCA containment response case.
(2) An 11 mil paint thickness is modeled for painted steel heat sinks in both GOTHIC and CONTEMPT, but not shown in the CONTEMPT thickness.
(3) Heat transfer to/from this heat sink is not modeled.
(4) A thin (2.09 mil) air gap is modeled between the steel plate and the concrete in GOTHIC. This simulates a gap conductance of 100 Btu/hr-sqft-0 F.
(5) The two-sided conductors are modeled with half thickness to conserve heat sink volume in both GOTHIC and CONTEMPT.
(6) A 4 mil zinc coating for the galvanized steel is modeled GOTHIC.
(7) A 28 mul paint thickness is modeled for painted concrete heat sinks in GOTHIC and CONTEMPT, but not shown in the CONTEMPT thickness.
(8) The MSLB limiting temperature case conservatively uses the containment minimum pressure, 14.2 psia.
(9) The value used in GOTHIC reflects the Technical Specification Limit on containment pressure of2 psig (LCO 3.6.4).
(10) Sensitivity cases have shown that GOTHIC calculates higher MSLB peak temperatures when an initial humidity of 100 percent was assumed.
(11) The CONTEMPT MSLB limiting temperature case conservatively starts the CS pump at 120 seconds after event initiation.
(12) This two-sided conductor represents a floor/ceiling inside containment and is modeled with the full thickness and one side area in GOTHIC.
(13) The LOCA M&E releases are calculated with minimum and maximum safety injection flowrates. The RHR flowrate varies with time in the maximum safety injection case.
(14) Service water temperature may vary for the long-term equipment qualification analysis.
(15) The CCW and Service Water systems were not explicitly modeled in CONTEMPT. A conservatively high CCW water temperature of 122.81F was used for the entire transient.
Note:
It is not the intent to restrict use of any of these inputs that are specific to plant equipment design or operation that may require changes in the future due to plant modifications conducted under 10 CFR 50.59 or due to any input calculations corrections or conservative enhancements.
WCAP-16219-NP                                                                                                                April 2004 6403-NP.doc-042904

3-6 1
I'l 20000 L
rr~                                                        I EI 15000                                                                       1 E-/
100           150                200                 250       300 Vapor Temperature (F)                           4.
I I
J.
J Figure 3-1 Fan Cooler Heat Removal Rate per Train              It WCAP-16219-NP                                                          April 2004 6403-NP.doc-042904

4-1 4         BENCHMARK COMPARISON A benchmark comparison between the CONTEMPT and GOTHIC code results was made to determine the impact of differences in the modeling methodology that cannot be completely eliminated with input changes (i.e., the modeling of liquid flashing and droplets). The containment response for both the LOCA and MSLB events was compared.
4.1       LOCA CONTAINMENT RESPONSE BENCHMARK COMPARISON The GOTHIC containment response model described in Section 2 was modified for the benchmark comparison to be like the CONTEMPT LOCA containment response model documented in Reference 13.
The following input changes were made to the GOTHIC model for the LOCA benchmark comparison case:
: 1.        The LOCA mass and enthalpy input data from Reference 13 was input to GOTHIC.
: 2.        Entrainment of the liquid break flow was limited to the blowdown portion of the event (22 seconds) by using a multiplier function on the break drop diameter. This was done to try to simulate the liquid temperature flash option in CONTEMPT.
: 3.        All the GOTHIC direct/DLM heat transfer coefficients were changed to Tagami/Uchida to match the CONTEMPT model. The Tagami peak was set to the end of blowdown (22 seconds) with an integrated energy release value of I .76E8 Btu, based on an integration of the LOCA energy release input data up to the end of blowdown.
: 4.        The GOTHIC revaporization fraction input value was set to 0.08 to match CONTEMPT.
: 5.        The GOTHIC fan cooler start time was set to 60 seconds to match CONTEMPT. Only one train (two fan coolers) was assumed to be operating, as in CONTEMPT.
: 6.        The GOTHIC spray start time was set to 75 seconds to match CONTEMPT. Only one pump was assumed to be running with a spray flowrate of 1,200 gpm, as in CONTEMPT.
: 7.        The GOTHIC containment volume liquid-vapor interface area was set to 8,659 square feet to match CONTEMPT.
: 8.        The GOTHIC initial containment pressure was set to 16.85 psia to match CONTEMPT.
The LOCA benchmark case was run until the start of recirculation. The containment pressure, temperature, and sump temperature comparisons are shown in Figures 4-1 through 4-3. The GOTHIC results are plotted as a solid line and the CONTEMPT results are represented by Xs. GOTHIC predicted a lower transient containment pressure and temperature response than CONTEMPT, but a higher transient sump temperature response. This indicated that GOTHIC was putting more energy into the water than the air. It was postulated that the difference might be due to the CONTEMPT temperature flash option for liquid break flow.
April 2004 WCAP-1     6219-NP WCAP-16219-NP                                                                                    April 2004 6403-NP.doc 042904

4-2 A modification was made to the GOTHIC input model to better simulate the temperature flash option in CONTEMPT. A circular flow path with a volumetric fan was added for this sensitivity case. The fan flowrate was set to 60,000,000 cfm and the flow path drop deposition was set to -1 (simulates 100 percent     I) de-entrainment). All of the liquid drops that enter the circular flow path are de-entrained and put into the sump. The volumetric fan is tripped off when the spray starts.
The containment pressure comparison with the simulated temperature flash model is shown in Figure 4-4.
GOTHIC predicted the same pressure response as CONTEMPT during the blowdown phase of the large                 l LOCA; however, GOTHIC still predicted a slight decrease in pressure between 20 and 60 seconds. The            I GOTHIC peak pressure is slightly lower than CONTEMPT (59.85 psia at 160.1 seconds versus 60.2 psia             l at 159 seconds) and the GOTHIC pressure is increasing slowly at the end of the transient.                      I The containment temperature comparison is shown in Figure 4-5. GOTHIC did not predict the increase in temperature that occurred just prior to the spray injection in CONTEMPT. The GOTHIC peak                    I temperature was 275.51F at 160.1 seconds versus a CONTEMPT peak temperature of 274.40 F at 75.4 seconds.
The containment sump temperature comparison is shown in Figure 4-6. The GOTHIC calculated sump                 I.
temperature is closer to, but still slightly higher than CONTEMPT.
Therefore, with these input changes, the GOTHIC model can be made to match the CONTEMPT results for the containment LOCA response reasonably well.                                                             l 4.2       MSLB CONTAINMENT RESPONSE BENCHMARK COMPARISON The GOTHIC containment response model described in Section 2 was modified for the benchmark                    l comparison to be like the CONTEMPT MSLB containment response model documented in Reference 5.
The following input changes were made to the GOTHIC model for the MSLB benchmark comparison case:
Id
: 1.       The faulted MSLB mass and enthalpy input data from Reference 5 was added to GOTHIC.
: 2.        The intact MSLB mass and enthalpy input data from Reference 5 was added to GOTHIC.                 ,
: 3.        A boundary condition and flow path were added to the GOTHIC model to incorporate the auxiliary feedwater mass and enthalpy input data from Reference 5.                                  1
: 4.        All the GOTHIC directlDLM heat transfer coefficients in GOTHIC were changed to direct/Uchida to match the CONTEMPT model.                                                            I
: 5.        The GOTHIC revaporization fraction input value was set to 0.08 to match CONTEMPT.
: 6.        The GOTHIC fan cooler start time was set to 60 seconds to match CONTEMPT. Only one train             I' (two fan coolers) were assumed to be operating, as in CONTEMPT.                                    -
WCAP-16219-NP                                                                                     Apil 2004 6403-NP.doc-042904

4-3
: 7.        The GOTHIC spray start time was set to 80 seconds to match CONTEMPT. Only one pump was assumed to be running with a spray flowrate of 1,200 gpm, as in CONTEMPT.
: 8.       A trip to terminate spray flow at 1,200 seconds was added to GOTHIC to match CONTEMPT.
: 9.        The GOTHIC containment volume liquid-vapor interface area was set to 0.0 square feet to eliminate heat and mass transfer to the break pool, as in CONTEMPT.
: 10.      The GOTHIC initial containment pressure was set to 16.85 psia to match CONTEMPT.
The MSLB benchmark case was run for 3,000 seconds. The containment pressure, temperature, and sump temperature comparisons are shown in Figures 4-7 through 4-9. The GOTHIC results are plotted as a solid line and the CONTEMPT results are represented by Xs. GOTHIC predicts a slower increase in pressure during the early part of the transient (0 to 80 seconds). GOTHIC does not predict the rapid increase in air temperature predicted by CONTEMPT at between 20 and 80 seconds. The GOTHIC calculated sump temperature is higher than predicted by CONTEMPT. This indicates that GOTHIC is storing more energy in the water than the air. The GOTHIC and CONTEMPT calculated pressure and temperature transient response converge again after spray starts. The GOTHIC calculated peak pressure is slightly lower than CONTEMPT (60.2 psia at 140.2 seconds versus 60.6 psia at 140.4 seconds).
A modification was made to the GOTHIC model to better simulate the temperature flash option in CONTEMPT. A circular flow path with a volumetric fan was added for this sensitivity case. The fan flowrate was set to 60,000,000 cfm and the flow path drop deposition was set to -1 (simulates 100 percent de-entrainment). All of the liquid drops that enter the circular flow path are de-entrained and put into the sump. The volumetric fan is tripped off when spray starts.
The containment pressure, temperature, and sump temperature comparisons are shown in Figures 4-10 through 4-12. Now, GOTHIC predicts the same initial increase in pressure and temperature prior to spray injection, but still does not predict the same response to the spray injection as CONTEMPT. The sump temperature comparison is much better. After spray injection occurs, the GOTHIC pressure continues to increase, but at a slower rate. The GOTHIC calculated peak pressure for this sensitivity case is about 1.3 psi higher than CONTEMPT (61.9 at 138.2 seconds).
A second modification was made to the GOTHIC model to simulate the 100 percent spray efficiency in CONTEMPT. The spray drop diameter input value was reduced from 1,000 microns to 100 microns.
The containment pressure, temperature, and sump temperature comparisons are shown in Figures 4-13 through 4-15. With the change in the spray drop diameter input value, GOTHIC predicts the same response to spray injection as CONTEMPT. The GOTHIC calculated peak pressure for this sensitivity case is about 0.2 psi higher than CONTEMPT (60.8 psia at 140.2 seconds).
Therefore, with these input changes, the GOTHIC model can be made to match the CONTEMPT results for the containment MSLB response reasonably well.
April 2004 WCAP-1 6219-NP                                                                                    April 2004 6403.NP.doc-042904

4-4
                                                                                  .1 AL 1                 ~Containment Pressure.l GOTHIC         CONTEMPT                                  l 110 'rme Fiue4tContainment Pressure
  -GOTHIC 7.1Patrhl(QA) OdW2=003 11:41:55 100 CoprsnfrIC (sec)1 1000 ecmr aea, 1e+00 s
4 WCAP-16219-NP                                                          April 2004 6403-NP.doc-042904

4-5 2                        Containment Temperature GOTHIC         CONTEMPT o                  X-- X 0t 0
04 EL 0
LO 0         I     I  IIIII ,,1 I I ,, ,,,     ,,     , 1,1,,,,1   , . , .,,   ,
1                       10                100           1000           le+004 Time (see)
GOTHIC 7.1Pat h1(QA) Oc=2O2003 l1A1i5 Figure 4-2 Containment Temperature Comparison for LOCA Benchmark Case April 2004 WCAP-1    6219-NP WCAP-16219-NP                                                                                        April 2004 6403-NP doc4042904

4-6 C1.1                                          .1 I
LOL C14.1 U.1 xL                                  I, CLL                                                       .1 E    0 Ln GOTHIC CONTEMPT                            1II   II 111
                                                                                                       .1
                                                                                                      'I, C)1 Time (sec).
GOTHIC 7.1 POW l(QA) OCVJ~2=DD 11.41.5.[
                                                                                                       .1
                                                                                                      .1~
14 It F igue4D    otimn              upTmeaueCmaio              o  ecmr    OACs WCAP-1 6219-NP                                                                             April 2004 6403-NP.doc-042904

4-7 Containment Pressure GOTHIC       CONTEMPT o            X-   -x LO In Iin ELt U)_
I I  H i l I I I I If l     I I I I 1I 1   II 1                       10            100             1000    Ie.004 Time (sec)
GOTHIC 7.1Pah1(QA) Octt2OZ03 11:52:0 Figure 4-4        Containment Pressure Comparison for LOCA Benchmark Case with Input Change to Simulate the CONTEMPT Temperature Flash Model WCAP-I 6219-NP                                                                             April 2004 6403-NP.doc-042904

                                                                                                -- U-4-8                                                                                               1' 2                   Containment Temperaturel GOTHIC      CONTEMPT.l o          X- -X LL 0L MX1L GOTHIC               1:20 CONTEMPTQA                                                       Ld I         Il I II C1.1           I I"-11 I I IIHl                        I1I 0                                                         .1 I,
Figure 4-5        Containment Temperature Comparison for LOCA Benchmark Case with Input            l Change to Simulate the CONTEMPT Temperature Flash Model April 2004 WCAP- 16219-NP WCAP-16219-NP                                                                        April 2004 6403-NP.doc-04290.4

4-9 3                          Sump Water Temperature GOTHIC    CONTEMPT x- -x C,.'
0                                 x 0x E                 i CD             II a ..1l l I II            *il I I, I.. 1   'I 'I '"
I  "
1                 10        100             1000          1e+004 Time (sec)
GOTHIC 7.lPakhl Q(ANow2O1ZXS3 12At3 Figure 4-6           Containment Sump Temperature Comparison for LOCA Benchmark Case with Input Change to Simulate the CONTEMPT Temperature Flash Model WCAP-16219-NP                                                                              April 2004 6403-NP.doc-042904

4-10                                                                                                    L I,
I
                                                                                                        .1 Containment Pressure                                                         .
GOTHIC     CONTEMPT                                                           .
C,,)
                                                                                                         .I, l
0.1 0              50      C..1 100            150       200       250 GOTHIC 7.11Parth1(OA) Ocf=01148:28 O
Time (see)                                      .1 U                                                        .1~
Figue47CotaimentPresureComarisn fr MSB, encharkCas April 2004 WCAP-162 19-NP WCAP-16219-NP                                                                               April 2004 6403-NP.doc-042904

4-11 2                         Containment Temperature GOTHIC      CONTEMPT X     -~X       -xXXt NI (N
U)
E 0               50       100           150 200        250 Tine (sec)
C0`HC7.1Paft1ct(QA) MdfSf=O3 12A828 Figure 4-8 Containment Temperature Comparison for MSLB Benchmark Case WCAP-1 6219-NP                                                                        April 2004 6403-NP.doc-042904

I 4-12                                                                                               L L
3                      Sump Water Temperature GOTHIC      CONTEMPT                                                          I o               X--X xxxxxxxxxxxxxxxxxxII
                                                                                                    .1 0                                                    .1 50          100           15D        200     250            1 La0 I,
GoT~l7.ln~hl(A~ocn5 o..1         14828Time (sec).
GO    7.49     FgrC      12.4828 CPotnn(A)     Scum   Tm5e2                 B     a
                                                                                                    .1 49CotaimentSum FigueTemeraure        or SLB encmarkCas WCAP-16219-NP                                                                            April 2004 6403-NP.doc-042904

4-13 Containment Pressure GOTHIC       CONTEMPT
        .i5,        c)mo                _XX                 XX~
2X    -,{
LO a.a 0,0               5010                 10      0025 0                 50   100             150    200          250 Time (sec)
GOTHIC 7.1Psth1(QA) Odt152003 11:1727 Figure 4-10        Containment Pressure Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option WCAP-16219-NP                                                                          April 2004 6403-NP.doc-042904

L 4-14                                                                                               1 xL 2Containment                           Temperature LOL GOTHIC           CONTEMPT o              ~X- -X                                                  .
CLL C1.1 1
0                    50           100            150 200 250 Time (sec).
GOT11C 7.1 Pa"I (QA) Od/15J203 I1A:T.27.
Figure 4-11          Containment Temperature Comparison for MSLB Benchmark Case with Input           1 Changes to Simulate the CONTEMPT Temperature Flash Option April 2004 WCAP-16219-NP                                                                           April 2004 6403.NP.doc-042904

4-15 3                       Sump Water Temperature GOTHIC    CONTEMPT 04 6,X                                   ->x
                                                       -X.-_  xxxxxxxxxxxxxxxx E
0             50        100            150    200     250 Time (sec)
Oc~ts2003 11727 GtWC7.1Palh1e0Aj Figure 4-12 Containment Sump Temperature Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option WCAP-16219-NP                                                                     April 2004 6403-NP.doc-042904

4-16                                                                                             '
AL
                                                                                                .L L
1                   Containment Pressure GOTHIC    CONTEMPT x- -x 0 ___100                               150      200          250 Time (sec)l Figure 4-13      Containment Pressure Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option and 100% Spray Effciency WCAP-1 6219-NP                                                                       April 2004 6403-NP.doc-442904

4-17 2                           Containment Temperature GOTHIC        CONTEMPT
                                  --       x- -x  X E                   10 EE 0                 50     100           150   200          250 Time (see)
GOThiC 7.1 Patchc(OA) Fe~btllb004 1428.1 Figure 4-14           Containment Temperature Comparison for MSLB Benchmark Case with Input Changes to Simulate the CONTEMPT Temperature Flash Option and 100% Spray Efficiency April 2004 WCAP-1 621 9-NP                                                                           April 2004 6403.NP.doc-042904

I 4-18 I
L L
3                           Sump Water Temperature                                                 I GOTHIC         CONTEMPT X--x I
CI I,
I E0 I,
in L
I C_
I, I-L1 I
o                 so     100           150   200          250             I Time (sec)                                    I, GOTIC 7.1Pldl t(Q  Febf1 t20G4 1425;35 I'
I, IL J,
I,
                                                                                                    .1~
Figure 4-15          Containment Sump Temperature Comparison for MSLB Benchmark Case with Input     -I, Changes to Simulate the CONTEMPT Temperature Flash Option and 100% Spray Efficiency                                                                       t WCAP-I 6219-NP                                                                           April 2004 6403-NP.doc-042904

5-1 5        SAMPLE TRANSIENT RESULTS USING THE PRAIRIE ISLAND GOTHIC CONTAINMENT EVALUATION MODEL This section presents containment response results for representative LOCA and MSLB transients using the Prairie Island GOTHIC containment evaluation model described in Section 2.0. The GOTHIC input tables for these sample cases are given in Appendices A and B and summarized in Section 3. The purpose of these sample cases is to demonstrate that the GOTHIC model has analysis margin to the containment design basis pressure limits.
5.1      LOCA SAMPLE TRANSIENT RESULTS The M&E input used in the LOCA GOTHIC containment evaluation model sample case was calculated using the Westinghouse method for LOCAM&E, as described in WCAP-10325 (Reference 7). The double-ended pump suction (DEPS) LOCA case was run through the transfer to sump recirculation.
While the M&E release input used in the containment evaluation model sample case is not presently the licensing basis LOCA release input, it is considered adequately conservative, yet representative of actual expected LOCA release input, for Prairie Island.
A DEPS LOCA is initiated from full power. A loss of offsite power and failure of an emergency diesel generator is assumed. The fan coolers start to remove heat at 60.3 seconds and the spray starts to inject at 75.6 seconds.
The containment pressure response is shown in Figure 5-1. The containment peak pressure is 40.7 psig at 14 seconds. This is less than the containment design pressure of 46 psig. The containment temperature response is shown in Figure 5-2. The containment temperature remains less than 260'F for the entire transient. The containment sump temperature response is shown in Figure 5-3. The containment sump temperature remains less than 2501F for the entire transient.
5.2      MSLB SAMPLE TRANSIENT RESULTS The M&E used as input into the MSLB GOTHIC containment evaluation model sample case was calculated using the Westinghouse method for MSLB M&E, as described in WCAP-8822 (Reference 9). A number of MSLB M&E cases, covering various break sizes and power levels, were generated to determine the limiting conditions for containment pressure. While the M&E release input used in the containment evaluation model sample case is not presently the licensing basis MSLB release input, it is considered adequately conservative, yet representative of actual expected MSLB release input, for Prairie Island.
The limiting MSLB event in terms of containment pressure represents a 4.6 ft2 break at 70-percent power.
A single failure of a feedwater regulator valve (FRV) is assumed. The fan coolers start to remove heat at 61 seconds and spray starts to inject at 73 seconds.
The containment pressure response for the peak pressure case is shown in Figure 5-4. The peak pressure for this case was 45.8 psig at 89 seconds. This is less than the containment design pressure of 46 psig.
The containment temperature response is shown in Figure 5-5. The containment temperature remains less than 2680 F for the entire transient. The containment sump temperature response is shown in Figure 5-6.
The containment sump temperature remains less than 2551F for the entire transient.
WCAP-1 621 9-NP                                                                                   April 2004 6403-NP.doc.042904

L 5-2 IL L-l I
I L
L.
I, Containment Pressure GOTHIC I
L 1,
                                                                                                      .1 Ul)
L I
IL U
L11 I,
                                                                                                     .1 U
1              10     100           1000   1le+004    1le+005 11 Time (sec)
GOTHIC 7.1Patch1(QA) Mar/19/204 13:32-20 L
I.
I I11 II
                                                                                                      .1~
I, I,
                                                                                                         .1 Figure 5-1 Containment Pressure Transient Response for LOCA                    I, I
I WCAP-1 6219-NP 6403-NP.doc-042904 April 2004  I

5-3 2                         Containment Temperature GOTHIC CO 0to 04 c'j M.
1               to0      100           1000 le+004 1 e+005 Time (sec)
GOTHIC 7.1 Patchl (OA) Mar/19/2004 13:3220 Figure 5-2 Containment Temperature Transient Response for LOCA WCAP-1 6219-NP                                                                            April 2004 6403-NP.doc-042904

s5..1                                                       t/
1~
I L
II 3                 ~Sump Water Temperature                                                >
0L 0U CD0.        0 GOTHIC                                      7.1                          arh a
E1 E            in i                 10Ql        100      1000  I e+004 p   e +005 0~1/
GOTHIC 7.1 Patchl1(OA) Mar11 9/2004 13:32:20 0.1          Tm sc
                                                                                                     .1 Figre5-SmpTepeatreTrnsen Cntinen   Rspns fr OC WCAP- 16219-NP                                                                           April 2004 6403-NP.doc-042904

5-5 Containment Pressure GOTHIC F~LO/
0 0             25      So      75      100       125     150 Time (sec)
GOTHIC 7.11Patchl (aA) Mar/19/2X0416:31:45 Figure 5-4 Containment Pressure Transient Response for MSLB April 2004 WCAP-1 621 9-NP                                                                         April 2004 6403-NP.doc-042904

L 5-6                                                                                            L L
L L
                                                                                              .1, 2                   Containment Temperature                                              I GOTHIC9 04 0+
c'J o.1~
Q..
lt12 Figure 5-5 Containment Temperature Transient Response for MSLB[
WCAP- 16219-NP                                                                     April 2004 6403-NP.doc-042904

5-7 3                      Sump Water Temperature GOTHIC 0
LO CNI E           _a E
CD cn 0            25         50      75    100      125      150 Time (sec)
GOTHIC 7.1Patch (A) Mar/19/2004 16:31:45 Figure 5-6 Containment Pressure Transient Response for MSLB WCAP-16219-NP                                                                           April 2004 6403-NP.doc-042904

6-1 6         

==SUMMARY==
OF RESULTS AND CONCLUSIONS The original Prairie Island containment evaluation model was based on the CONTEMPT code. A new Prairie Island containment evaluation model, based on the NRC-approved containment evaluation model for Kewaunee, was built using the GOTHIC code. Most of the input data for the Prairie Island GOTHIC containment evaluation model was taken from the CONTEMPT LOCA and MSLB containment model input decks.
The GOTHIC code modeling assumptions and input values were compared with the CONTEMPT code modeling assumptions and input values. Differences were identified in modeling condensation heat and mass transfer to the heat sinks, flashing of the liquid break flow, and condensation on the spray droplets.
The heat and mass transfer correlations can be changed in GOTHIC to match the CONTEMPT models; however, GOTHIC does not have the same flashing or spray condensation models as CONTEMPT. To determine the effect of these differences, the GOTHIC Prairie Island containment evaluation model was modified for benchmark comparisons with the original CONTEMPT LOCA and MSLB containment evaluation models. In addition to the various input changes required to add the M&E releases and change the heat transfer correlations, a circular flow path with drop de-entrainment was used to simulate the temperature flash option in CONTEMPT and the containment spray drop diameter input value was reduced by a factor of 10 to simulate the 100 percent spray efficiency in CONTEMPT. With these changes, the GOTHIC model results were reasonably close to those predicted by CONTEMPT. For the LOCA event, GOTHIC predicted a slightly lower (0.35 psi) peak pressure and a slightly higher (1.1 0F) peak temperature. For the MSLB event, GOTHIC predicted a slightly higher (0.2 psi) peak pressure and a slightly lower (10 0 F) peak temperature. GOTHIC predicted a slightly higher sump temperature forboth the LOCA and MSLB transients.
The GOTHIC Prairie Island containment evaluation model was used to produce sample results for the LOCA and MSLB transients. The GOTHIC model predicted that the containment would remain below design pressure for both cases.
In conclusion, as shown in the benchmark cases, the GOTHIC code maintains a similar degree of conservatism with respect to the present licensing basis CONTEMPT code. The modeling assumptions used in the proposed GOTHIC evaluation model have adequate engineering justification for being more mechanistic or conservative, as has been previously determined in the Kewaunee License Amendment No. 169, which was reviewed and approved by the NRC. Since the Prairie Island and Kewaunee plants are nearly identical Westinghouse two-loop pressurized water reactor plants in design and operation, it is reasonable to apply the same methodology for performing the containment design basis calculations.
Therefore, the GOTHIC evaluation model described in this report will provide conservative results to assure containment design basis limits are met as required by the Prairie Island licensing basis General Design Criterion 10 with consideration of General Design Criterion 41, as described in Section 1.5 of the Prairie Island USAR.
April 2004 WCAP-1 6219-NP                                                                                     April 2004 6403-NP.doc.042904

7-1 7         REFERENCES
: 1.       NAI 8907-06, Revision 13, "GOTHIC Containment Analysis Package Technical Manual,"
Version 7.1, January 2003.
: 2.        NAI 8907-09, Revision 7, "GOTHIC Containment Analysis Package Qualification Report,"
Version 7.1, January 2003.
: 3.        NAI 8907-02, Revision 14, "GOTHIC Containment Analysis Package User Manual," Version 7.1, January 2003.
: 4.        NRC Letter from Anthony C. McMurtray (NRC) to Thomas Coutu (NMC), Enclosure 2, Safety Evaluation, September 29, 2003.
: 5.        NMC folder CF.PX.00.RSE.016 "MSLB Cont Press Response, Updated," February 2001.
: 6.        Kreith, "Principles of Heat Transfer," Third Edition, 1973.
: 7.        WCAP-1 0325-P-A, "Westinghouse LOCA Mass and Energy Release Model for Containment Design March 1979 Version," May 1983.
: 8.        Prairie Island Updated Safety Analysis Report, Revision 25, May 2003.
: 9.        WCAP-8822, "Mass and Energy Releases Following a Steam Line Rupture," September 1976.
: 10.       Brown and York, "Sprays formed by Flashing Liquid Jets," AICHE Journal Volume 8, #2, May 1962.
: 11.      NAI 9301-05, "Comparison of Results from GOTHIC to Results from Other Codes,"
January 1996.
: 12.      NMC Letter NRC-03-077, "Response to Request for Additional Information Related to NMC Request for the use of GOTHIC 7 for the Kewanuee Nuclear Power Plant Containment Design Basis Accident Analyses," July 24, 2003.
: 13.      NMC CONTEMPT Case Number: /CASE/PX/00/CONTEQ/CNP/C010006.
April 2004 WCAP-l 6219-NP                                                                              April 2004 6403-NP.doc.042904

A-l APPENDIX A GOTHIC PRAIRIE ISLAND MSLB CONTAINMENT RESPONSE MODEL INPUT TABLES WCAP-16219-NP                                     April 2004 6403-NP.doc-042904

U A-2 Pl Containment DBA Model - MSLB Peak Pressure for OSG Apr/16/2004 13:08:12 GOTHIC Version 7.1 Patcthl(OA) - September 2003 L
File: D.\DATA\GOTHIC\Prairte Island\OSG MSLB Cases\FCFlp\peak-press-mslb PI Containment DBA Model - MSLB Peak Pressure for OSG                                                                                    IJ Control Volumes Vol                                     Vol           Elev          Ht        Hyd. D.       L/V IA     Burn Description                 (ft3)          (ft)         Cft)         (ft)         (ft2)      Opt 1    lContainment Vol           l1320000.       l        0. I       100. I       10. l    16840. lNONE Laminar Leakage Leak Rate           Ref            Ref          Ref        Sink/                                       Leak Vol       Factor          Press           Temp        Humid       Source      Model       Rep     Subvol       Area I       (%/hr)         (psia)            (F)          C%)          BC       Option      Wall    Option       (ft2) l1              0. I               I              I           I           I           I       I           IDEFAULT Turbulent Leakage Leak Rate          Ref             Ref          Ref       Sink/                                        Leak Vol        Factor          Press           Temp        Humid       Source      Model       Rep     Subvol       Area 1       (%/hr)         (psia)            (F)          (%)          BC       Option     Wall     Option       (ft2)       fL/D 1               0.                 I              I           I           I           I       IUNIFORM    IDEFAULT     I Fluid Boundary Conditions - Table 1 Press.              Temp.              Flow            S J    ON     OFF     Elev.
BC#        Description               (psia)        FF      (F)      FF   (lbn/s)      FF P 0 Trip        Trip       Mft) 1F     SG Side Break F                    60.                El     2            1    1 N Y         0       0       30.
2F     Condenser Side                     60.                El     2            0   1 N Y                          30.
3F     Spray from RWST                  14.7               120         V0.0044       4 N Y          2       0       90.
Fluid Boundary Conditions - Table 2 Liq. V.          Stm.             Drop D.          Cpld      Flow            Heat             Outlet BC#      Frac. FF     P.R.        FF      (in)     FF BC#        Frac.      FF   CBtu/s)      FF Quality      FF 1F          1.                  1         0.C039                                                    DEFAULT 2F          1. l1                         0.0039                                                    DEFAULT WCAP- 16219-NP                                                                                                                   April 2004 6403-NP.doc-042904

A-3 Pl Containment DBA Evaluation Model - MSLB Peak Pressure                                                         2 Feb/18/2004 11:44:57 GOTHIC Version 7.1Patch1 (OA) - September 2003 File: D:XDATA\GOTHICAPrairie Island\MSLB Cases\peak-pressrnsib Fluid Boundary Conditions             -  Table 2 (cont.)
Liq. V.           Stxn.          Drop D.        Cpld      Flow         Heat        Outlet BC#      Frac. FF      P.R.      FF     (in)     FF BC1       Frac. FF (Btu/s)   FF Quality  FF 3F            -1.l              1l     10.0394 l      l       1                           DEFAULT Fluid Boundary Conditions - Table 3 Gas Pressure Ratios Air BC#     Gas 1     FF     Gas 2      FF    Gas 3     FF   Gas 4     FF IF          1.
2F          1.
3F           1.
Fluid Boundary Conditions - Table 4 Gas Pressure Ratios BC#     Gas 5     FF     Gas 6      FF   Gas 7      FF   Gas 8     FF IF 2F 3F Flow Paths - Table 1 F.P.                                 Vol        Elev        Ht        Vol     Elev      Ht Description                A         (ft)      (ft)        B      (ft)     (ft) 1     SG Side Break F               1           30.                 1F        30.
2     Condenser Side                1           30.                 2F        30.
3     Spray Flow                    1           90.                 3F        90.
Flow Paths - Table 2 Flow        Flow           Hyd.        Inertia       Friction      Relative     Dep    Mom Strat.
Path        Area          Diam.         Length        Length        Rough-     Bend    Trn Flow (ft2)           (ft)           (ft)          (ft)         ness      (deg)   Opt   Opt I              1.I             1.             1            1                     0.      NONE WCAP-16219-NP                                                                                             April 2004 6403.NP.doc.042904

M A4 L
I,1 PI Containment DBA Model - MSLB Peak Pressure for OSG                                                                    3 Apr116/2004 13:08:12 GOTHIC Version 7.1 Patch I (QA) - September 2003 File: D NDATA\GOTHIC\Praide Islanc\OSG MSLB Cases\FCFIX\peak_press~mslb Flow Paths - Table 2              (cont.)
Flow       Flow            Hyd.         Inertia          Friction       Relative           Dep        Mom Strat Path       Area           Dian.          Length             Length        Rough-          Bend        Trn Flow
         #      (ft2)           (ft)            (ftl               (ft)           ness         (deg)       Opt  Opt 2               1.             1.               1.               1.                        0.       -  NONE 3               1.             1.               1.               1.                        0.       _  NONE Flow Paths - Table 3 Flow       Fwd.            Rev.                    Critical           Exit        Drop Path       Loss            Loss        Comp.           Flow           Loss     Breakup d     Coeff.          Coeff.        Opt.           Model        Coeff.       Model 1                                     OFF            OFF              0.       OFF 2                                     OFF            OFF              0.       OFF 3                                     OFF            OFF              0.       OFF Thermal Conductors -              Table 1 Cond                                  Vol       HT        Vol      HT Cond        S. A.       Init.
d       Description                 A        Co          B      Co Type         (ft2)       T.(F)       Or 1    Containment Cyl                 1       1           1        2     1      40300.           120. I 2    Containment Dom                 1        1          1        2     2      17300.           120. I 3    Reactor Vessel                 1        2           1        2     3        1260.          120. I 4    Refueling Canal                 1       1           1        2     4        5150.          120. I 5    Refueling Canal                 1       1           1        2     5        1010.          120. I 6    Exposed Piping                 1        1           1        2     6        5780.          120. I 7    Structural       Stee          1        1           1        2     7        4055.          120. I 8    Structural       Stee           1       1           1        2     8      16925.           120. I 9    Structural       Stee           1       1           1       2      9      28500.           120. I 10    Structural Stee                1        1           1       2    10         2000.          120. I 11    Structural Stee                I         1          1        2   11           500.         120. I 12    Handrails                       1       1           1        2   12         1695.          120. I 13    Grating                        1        1           1       2    13       12400.           120. I 14    Conduit       - Cable          1        1           1        2   14       14000.           120. I 15    Ductwork                       1        1           1        2   15       27100.           120. I 16    Accumulators 1                 1        1           1        2   16         1346.          120. I 17    Accumulators 2                 1        1           1        2   17         1256.          120. I 18    Accumulators 3                 1        1           1       2    18           973.         120. I 19    Heavy Walls                    1        1           1       2    19       19320.           120. I 20    Heavy Walls                    1        1           1       2    20       15000.           120. I 21    Heavy Walls w/P                1        1           1       2    22         3810.          120. I 22    Heavy Floors                   1        3           1        1   22       3322.5           120. I WCAP-16219-NP                                                                                                        April 2004 6403-NP.doc-042904

A-5 PI Containment DBA Model -MSLB Peak Pressure for OSG                                               4 Apr/16/2004 13:08:12 GOTHIC Version 7.1 Patchl (OA) - September 2003 File: D:.DATA\GOTHICXPrairie Island\OSG MSLB CasesXFCFl)Cpeak-pressjimsib Thermal Conductors - Table 1 (cont.)
Cond                                 Vol      HT     Vol   HT Cond        S. A. Init.
        #        Description               A       Co      B    Co Type        (ft2)  T.(F)  Or 23     Heavy Floors                 1        3       1    2      21      5830. 120. I 24     Light Floors                 1        3       1    1      23     7687.5   120. I Thermal Conductors - Table 2 Cond     Therm. Rad.              Emiss.         Therm. Rad.          Emiss.
        #         Side A               Side A            Side B            Side B 1            No                                     No 2            No                                     No 3            No                                     No 4            No                                     No 5            No                                     No 6            No                                     No 7            No                                     No 8            No                                     No 9            No                                     No 10            No                                     No 11            No                                     No 12            No                                     No 13            No                                     No 14            No                                     No 15            No                                     No 16            No                                     No 17            No                                     No 18            No                                     No 19            No                                     No 20            No                                     No 21            No                                     No 22            No                                     No 23            No                                     No 24            No                                     No WCAP-16219-NP                                                                                   April 2004 6403-NP.doc-042904

I A-6                                                                                                                        I PI Containment DBA Evaluation Model - MSLB Peak Pressure Feb/18/2004 11:44:57 5  I GOTHIC Version 7.1 Patch1 (QA)- September 2003 File: DADATA\GOTHIC\Prairie IslandWMSLB Cases'peakpress.mslb Heat Transfer Coefficient Types - Table 1 (cont.)
Heat                            Cnd/            Sp       Nat              For Type    Transfer        Nominal            Cnv Cnd        Cnv       Cnv              Cnv        Rad
       #     Option           Value        FF Opt Opt         HTC       Opt              Opt        Opt 31Direct         I              I    IADD    DLM   l      IVERT SURF     IPIPE FLOW       ION Heat Transfer Coefficient Types - Table 2 Min          Max       Convection         Condensation Type    Phase       Liq          Liq         Bulk Temp           Bulk Temp 4  Opt       Fract        Fract       Model       FF      Model      FF 1   YAP                                Tg-Tf               Tb-Tw 2
3 SPLIT             0. le-004       Tg-Tf               Tb-Tw
                                                                                                                            .1 Heat Transfer Coefficient Types - Table 3 Char.       Nat Conv          For Conv        Nom              Minimum          Char.
Type    Length      Coef       Exp   Coef      Exp    Vel       Vel   Conv HTC         Height 4     (ft)         FF       FF      FF      FF   I ft/s)     FF   (B/h-f2-F)          Ift)
DEFAULT             1.
DEFAULT HTC Types - Table 4 Total           Peak                   Initial         BD        Post-BD    Post-BD Type     Ccnst            Heat           Time       Exp          Value         Exp          EXP      Direct
       #      CT            IBtu)           Isec)       XT       (B/h-f2-F)        yt           xt         FF 1    72.5                                      0.62                          1.       0.025 2
3    72.5                                      0.62                          1.       0.025 WCAP-16219-NP                                                                                                   April 2004 6403-NP.doc-042904

A-7 PI Containment DBA Model - MSLB Peak Pressure for OSG                                                 6 Apr/16/2004 13:08:12 GOTHIC Version 7.1Patchl (OA) - September 2003 File: D:tDATA\GOTHIC\Prairie lslandOSG MSLB Cases\FCFIX\peakpress_mstb Thermal Conductor Types Type                                       Thick.         O.D.                  Heat     Heat
        #       Description            Geom       (in)         (in)      Regions (Btu/ft3-s)   FF 1     Painted CS               WALL        1.511            0.       2             0.
2     Painted CS               WALL        0.761            0.       2             0.
3     CS Lined Concre          WALL     12.2006             0.      4              0.
4     SS Lined Concre          WALL     12.1896             0.       3             0.
5     SS Lined Concre          WALL     12.2521             0.       3             0.
6     Painted CS               WALL        0.386            0.      2              0.
7     Painted CS               WALL        0.179            0.       2             0.
8     Painted CS               WALL        0.261            0.      2              0.
9     Painted CS               WALL        0.386            0.       2             0.
10     Painted CS               WALL        0.761            0.       2             0.
11     Painted CS               WALL        1.011            0.      2              0.
12     Painted CS               WALL        0.261            0.      2              0.
13     Painted CS               WALL     0.10475             0.       2             0.
14     Galvanized               WALL        0.104            0.       2             0.
15    Galvanized                WALL       0.0582            0.      2              0.
16     Painted CS               WALL        2.761            0.      2              0.
17     Painted CS               WALL        1.401            0.      2              0.
18     Painted CS               WALL       0.3235            0.      2              0.
19     6-in Concrete            WALL             6.          0.      1              0.
20     12-in Concrete           WALL           12.           0.      1              0.
21     12-in Painted C          WALL       12.028            0.      2              0.
22     12-in Painted C          WALL       12.056            0.      3              0.
23     6-in Painted Co          WALL        6.056            0.      3              0.
Thermal Conductor Type 1
Painted CS Mat.         Bdry.         Thick       Sub-     Heat Region         #           (in)         (in)       regs. Factor 1         1                 0.       0.011         1       0.
2         2            0.011            1.5       15       0.
Thermal Conductor Type 2
Painted CS Mat.         Bdry.         Thick       Sub-     Heat Region         #           (in)         (in)      regs. Factor 1      1                    0.      0.011          1       0.
WCAP-16219-NP                                                                                      April 2004 6403-NP.doc.042904

A-8 Pl Containment DBA Evaluation Model - MSLB Peak Pressure                        7 Febl8/2004 11:44:57 GOTHIC Version 7.1 Patctl(OA) - September2003 File: D:\DATA\GOTHIC\Prairie Island\MSLB Cases\peakpressmsIb Thermal     Conductor Type          (cont.)
2 Painted CS Mat.         Bdry.         Thick        Sub-     Heat Region                     (in)          (in)       regs. Factor 2         2             0.011            0.75         8        0.
Thermnal Conductor         Type 3
CS Lined Concrete Mat.         Bdry.         Thick        Sub-     Heat Regicn          N           (in)          (in)       regs. Factor 1         1                  0.       0.011          1        0.
2        2              0.011      0.1875            2         0.
3         5          0.1985        0.0021            1        0.
4        4           0.2006               12.      24         0.
Thermal Conductor Type 4
SS Lined Concrete Mat.         Bdry.         Thick         Sub-     Heat Region          x           (in)          (in)       regs. Factor 1        3                   0. 0.1875            2        0.
2         5           0.1875        0.0021            1        0.
3        4           0.1896               12.      24         0.
Thermal Conductor Type 5
SS Lined Concrete Mat.         Bdry.         Thick         Sub-     Heat Region         #           (in)           (in)       regs. Factor 1        3                   0.         0.25         3        0.
2         5                0.25     0.0021            1        0.
3         4           0.2521               12.      24         0.
WCAP-16219-NP                                                              April 2004 6403-NP.doc-042904

A-9 PI Containment DBA Evaluation Model - MSLB Peak Pressure                        8 Feb/18/2004 11:44:57 GOTHIC Version 7.1 Patchl (OA) - September 2003 File: D:ADATA\GOTHIC\PraIrie Island\MSLB Cases\peakpressnmslb Thermal Conductor Type 6
Painted CS Mat.         Bdry.         Thick       Sub-     Heat Region          y          (in)           (in)     regs. Factor 1         1                  0.       0.011       1         0.
2         2             0.011         0.375       4         0.
Thermal Conductor Type 7
Painted CS Mat.         Bdry.         Thick       Sub-     Heat Region         #           (in)           (in)     regs. Factor 1         1                  0.       0.011         1       0.
2         2             0.011         0.168         2       0.
Thermal Conductor Type 8
Painted CS I.at.        Bdry.         Thick       Sub-     Heat Region         #           (in)           (in)     regs. Factor 1         1                  0.       0.011          1      0.
2         2            0.011           0.25         3       0.
Thermal Conductor Type 9
Painted CS Mat.         Bdry.         Thick       Sub-     Heat Region         #           (in)           (in)     regs. Factor 1      1                  0.       0.011         1       0.
2         2            0.011          0.375         4       0.
WCAP-16219-NP                                                           - April 2004 6403-NP.doc-042904

U I
A-10 L
                                                                                    'L Pt Containment DBA Evaluation Model - MSLB Peak Pressure                       9 1L Feb/1 8/2004 11:44:57 GOTHIC Version 7.1 Patchl(OA) - September 2003 File: D:IDATA\GOTHIC\Prairie Island\MtSLB Cases~peakpressmslD Thermal Conductor Type 10 Painted CS Mat.         Bdry.          Thick      Sub-     Heat Region           I          (in)           (in)     regs. Factor 1          1                 0.        0.011        1        0.
2          2            0.011           0.75        8        0.
Thermal Conductor Type 11 Painted CS Mat.         Bdry.          Thick      Sub-     Heat Region           #          (in)           (in)     regs. Factor 1          1                 0.        0.011        1       0.
2          2           0.011               1. 10        0.
Thermal Conductor Type 12 Painted CS Mat.         Bdry.          Thick      Sub-     Heat Region           #          (in)           (in)     regs. Factor 1          1                 0.       0.011         1       0.
2         2            0.011            0.25        3       0.
Thermal Conductor Type 13 Painted CS Mat.         Bdry.         Thick      Sub-      Heat Region          #           (in)           (in)     regs. Factor 1         1                  0.       0.011         1       0.
2          2            0.011       0.09375          1       0.
WCAP-16219-NP                                                            April 2004 6403-NP.doc-04290.4

A-1l Pi Containment DBA Evaluation Model - MSLB Peak Pressure                      10 Feb/1812004 11:44:57 GOTHIC Version 7.1 Patchl (OA) - September 2003 File: D:tDATA\GOTHICPrairie Island\MSLB Casestpeak-press-nmsIb Thermal Conductor Type 14 Galvanized Mat.        Bdry.          Thick        Sub-     Heat Region         #          (in)            (in)      regs. Factor 1        6                  0.        0.004          1       0.
2         2            0.004             0.1         1       0.
Thermal Conductor Type 15 Galvanized Mat.        Bdry.          Thick        Sub-     Heat Region         #          (in)            (in)      regs. Factor 1         6                 0.        0.004          1       0.
2         2            0.004        0.0542           1       0.
Thermal Conductor Type 16 Painted CS Mat.        Bdry.          Thick       Sub-      Heat Region          #          (in)            (in)      regs. Factor 1         1                 0.        0.011       1          0.
2         2            0.011           2.75      28          0.
Thermal Conductor Type 17 Painted CS Mat.        Bdry.          Thick       Sub-      Heat Region          #          (in)            (in)      regs. Factor 1         1                 0.        0.011       1          0.
2         2            0.011           1.39      14          0.
WCAP-16219-NP                                                            April 2004 6403-NP.doc.042904

A-12 PI Containment DBA Evaluation Model - MSLB Peak Pressure                      11 Feb/i8/2004 11:44:57                                                              .11 GOTHIC Version 7.1lPatchl(OA) -September 2003 File: D:\DATA\GOTHIC\Prairie )sland\MSLBCases\poak-press-mslb Thermal Conductor Type 18 Painted CS I.
                                                                                    .1 Mat.         Bdry.        Thick        Sub-      Heat Region          #          (in)         (in)       regs. Factor 1          1                0.      0.011         1          0.             I 2         2            0.011       0.3125         3          0.
                                                                                   .1 Thermal Conductor Type 19 6-in Concrete Mat.         Bdry.        Thick        Sub-      Heat Region         K           (in)         (in)       regs. Factor 1     1   4     1           0.1          6.         12       0.
Thermal Conductor Type 20 12-in Concrete Mat.         Bdry.        Thick        Sub-      Heat Region         If          (in)         (in)       regs. Factor 1         4                 0.          12. 1 24             0.
Thermal Conductor Type 21 12-in Painted Concrete Mat.         Bdry.        Thick        Sub-      Heat Region         K           (in)         (in)       regs. Factor 1         1                 0.      0.028             1      0.
2         4            0.028           12.          24       0.
WCAP-I 6219-NP                                                          April 2004 6403-NP.doc-042904

A-13 PI Containment DBA Model - MSLB Peak Pressure for OSG                                                                      12 Apr/1612004 13:08:12 GOTHIC Version 7.1 Patchl (GA) - September2003 File: D:-DATA\GOTHIC\Prairie Islarfd\OSG MSLB Cases\FCFIX\peak.pressmsfb Thermal Conductor Type 22 12-in Painted Concrete Mat.         Bdry.           Thick        Sub-      Heat Region          *          (in)            (in)       regs. Factor 1          1                  0.       0.028            1        0.
2          4            0.028              12.        24         0.
3          1          12.028           0.028            1        0.
Thermal Conductor Type 23 6-in Painted Concrete Mat.         Bdry.           Thick        Sub-      Heat Region         #           (in)            (in)       regs. Factor 1         1                   0.       0.028        1            0.
2         4            0.028                6. 12            0.
3         1             6.028          0.028        1            0.
Cooler/Heater Heater                                            On     Off      Flow     Flow      Heat      Heat Cooler                                 Vol.      Trip Trip        Rate     Rate     Rate       Rate    Phs Ctrlr I       Description                #         I       #      (CFM)     FF     (Btu/s)      FF     Opt     Loc IC    IFan Cooler Tr           I    1     l   1            le+005 I         l        2.1    3    lVTS  I1 Volume Initial Conditions Vapor       Liquid       Relative      Liquid       Ice          Ice Vol         Pressure         Temp.         Temp.      Humidity      Volume     Volume     Surf.A.
           #           (psia)            (F)          F            (%)       Fract. Fract.        (ft2) def              16.7           120.         120.            30.          0.         0.           0.
Initial Gas Pressure Ratios Vol          Air I        Gas 1       Gas 2        Gas 3      Gas 4      Gas 5      Gas 6     Gas 7      Gas 8 def              1.           0. I       0. I        0.         0.        0.         0.         0.
WCAP-162 19-NP                                                                                                        April 2004 6403-NP.doc-042904

L A-14 L
L 1L PI Containment DBA Evaluation Model - MSLB Peak Pressure                                                         13 Feb/18/2004 11:44:57 GOTHIC Version 7.1 PatchI (OA) - September 2003 File: D:\DATANGOTHIC\Prairie lstand\MSLB Cases'peak press-mslb Noncondensing Gases Gas         Description             Symbol     Type       Mol.       Lennard-Jones       Parameters No.                                                     Weight           Diameter              e/K (Ang)               (K) 1   Air                      I   Air     IPOLY 1      28.97   l             3.617 1            97.
Noncondensing Gases - Cp/Visc.               Equations Gas               Cp   Equation        (Required)                Visc.      Equation     (Optional)
No.         Tmin          Tmax                  Cp              Tmin           Tmnax        Viscosity (R)          (R)           (Btu/lbm-R)              (R)           (R)        (lbm/ft-hr) 11         360.       2280.      0.238534-6.2006 Materials Type #                  Description                   Gap 1       Paint                                    NO 2       Carbon Steel                             NO 3       Stainless Steel                          NO 4       Concrete                                 NO 5       Gap                                      NO 6       Zinc                                     NO Material Type 1
Paint Temp.          Density              Cond.            Sp. Heat (F)         (lb-n/ft3)     (Btu/hr-ft-F)        (Btu/lbm-F)
: 32.             28.                  0.29                  1.
500.             28.                  0.29                  1.
WCAP-1621 9-NP                                                                                             April 2004 6403-NP.doc-042904

A-15 Pi Containment DBA Evaluation Model - MSLB Peak Pressure                     14 Feb/18/2004 11:44:57 GOTHIC Version 7.1 Patchl (OA) - September 2003 File: DO:\DATA\GOTHIC\Prairie Island\MSLB Cases\peak-press-mslb Material Type 2
Carbon Steel Temp.        Density                Cond.             Sp. Heat (F)       (lbm/ft3)         (Btu/hr-ft-F)        (Btu/lbm-F)
: 32.           490.                    26.            0.115 500.           490.                    26.            0.115 Material Type 3
Stainless Steel Temp.        Density                Cond.             Sp. Heat (F)       (lbm/ft3)         (Btu/hr-ft-F)        (Btu/lbm-F)
: 32.           488.                    9.4            0.123 500.           488.                    9.4            0.123 Material Type 4
Concrete Temp.        Density                Cond.            Sp. Heat (F)       (lbm/ft3)         (Btu/hr-ft-F)        (Btu/lbm-F)
: 32.           144.                    0.8              0.2 500.            144.                    0.8              0.2 Material Type 5
Gap Temp.        Density                Cond.            Sp. Heat (F)       (lbm/ft3)         (Btu/hr-ft-F)        (Btu/lbm-F)
: 32.           0.06                0.0174             0.241 500.           0.06                0.0174             0.241 WCAP-16219-NP                                                          April 2004 6403-NP.doc-042904

L A-16 L
L L
Pi Containment DBA Evaluation Model - MSLB Peak Pressure                                                  15 Feb/11 12004 11:44:57 GOTHIC Version 7.1 Patchl (OA) - September 2003 File: D:\DATA\GOTHICPrairie IsIand\MSLB Cases\peak_pressmsib Material Type 6
Zinc Temp.          Density             Cond.            Sp. Heat (F)        (lbm/ft3)      (Btu/hr-ft-F)        (Btu/lbm-F)
: 32.           446.                    64.           0.091 500.           446.                    64.           0.091 Ice Condenser Parameters Initial           Bulk        Surface Area           Heat Temp.         Density         Multiplier       Transfer (F)       (lbm/ft3)         Function          Option
: 15.         33.43                          UCHIDA Component Trips Trip                         Sense     Sensor    Sensor     Var.      Set   Delay  Rset Cond Cond I     Description        Var.      1 Loc. 2 Loc. Limit   Point    Time  Trip Trip Type 1     Fan Cooler         PRESS         1                UPPER      19.7    60.
2     Spray Injec        PRESS         1                UPPER      38.7    72.           AND Functions FF#         Description            Ind. Var.       Dep. Var. Points 0     Constant                                        -           0 1     Framatome Mass           Time    (sec)    Flow Rate        76 2     Framatome Entha          Time    (sec)    Enthalpy (       76 3     Fan Cooler Heat          Vapor    Sat. Heat Remov        14 4     Spray Flow Rate          Time    (sec)    Flow Rate         5 WCAP-16219-NP                                                                                       April 2004 6403-NP.doc-042904

A-17 PI Containment DBA Evaluation Model - MSLB Peak Pressure                    16 Feb/18004 11:44:57 GOTHIC Version 7.1 Patchl (OA) - September 2003 File:DA\DATA\GOTHIC\Prairie Island\MSLB Cases'peakpress_mslb Function 1
Framatome Mass Limiting Pressure Case (tabl Ind. Var.: Time (sec)
Dep. Var.: Flow Rate (lbm/s)
Ind. Var.        Dep. Var.          Ind. Var.        Dep. Var.
: 0.          1253.3                  1.       1253.3
: 2.          1220.1                  3.       1226.2
: 4.          1595.1                  5.       2136.1
: 6.         2364.3                   7.       2589.9
: 8.           2702.                  9.        2567.
: 10.            2349.                12.        2249.5
: 14.          2042.5                 16.        1577.5
: 18.          1377.5                 20.         1308.
: 22.           1366.5                24.        1084.5
: 26.          1145.5                 28.         1087.
: 30.          1109.5                 32.        1013.5
: 34.            1266.                36.         891.5
: 38.              831.               40.          741.
: 42.              689.               44.           696.
: 46.            757.5                48.           650.
: 50.              629.               55.           607.
: 60.            578.4                65.         553.2
: 70.            529.8                75.         508.8
: 80.            490.2                85.         472.8
: 90.            457.2                95.         442.6 100.             429.6               105.         420.2 110.             411.4               115.         404.6 120.              399.               125.         393.8 130.             389.2               135.         385.2 140.               382.              145.         378.8 150.             376.8               160.         375.9 170.             374.9               180.         373.3 190.             371.5               200.         369.2 220.             365.2               240.          360.
260.             354.3               280.          348.
300.             340.3               320.         329.2 340.             312.7               360.         291.6 380.             267.2               400.         230.7 420.             201.8               440.         114.9 460.               25.3              480.           12.2 500.                 12.             520.            12.
540.              12.1               560.            12.
580.               12.1              600.           12.1 WCAP-1 621 9-NP                                                      April 2004 6403-NP.doc-042904

aa A-18
                                                                                ,L Pi Containment DBA Evaluation Model - MSLB Peak Pressure L
17 Fob/18/2004 11:44:57 GOTHIC Version 7.IPatchi(OA) - September 2003 File: D:\DATA\GOTHIC\Praine Islanr\MSLB Cases'peak pressrnslb Function 2
Framatome Enthalpy Limiting Pressure Case                                   L Ind. Var.: Time (sec)
Dep. Var.: Enthalpy (BTU/lbm)
L Ind. Var.         Dep. Var.         Ind. Var.       Dep. Var.              F
: 0.        1190.04                  1.       1190.04
: 2.        1188.31                  3.       1167.82
: 4.         959.28                  5.        792.65
: 6.         739.69                  7.        699.14
: 8.         679.68                  9.        692.79
: 10.         719.88                12.         726.65
: 14.         751.46                16.         851.32
: 18.         897.68                20.         907.15
: 22.         871.72                24.         996.36
: 26.         944.96                28.         964.77
: 30.         939.25                32.         983.23
: 34.         843.84                36.       1044.03
: 38.        1078.04                40.       1153.17
: 42.        1195.72                44.       1168.97
: 46.        1082.38                48.       1190.92
: 50.        1201.19                55.       1201.52
: 60.        1202.52                65.       1202.35
: 70.         1202.57                75.       1203.26
: 80.        1202.65                85.       1203.05
: 90.        1203.06                95.       1202.89 100.         1203.26              105.        1202.57 110.         1203.31              115.        1203.51 120.         1203.51              125.        1202.84 130.         1203.49              135.        1203.53 140.         1203.14              145.        1203.27 150.        1203.29               160.        1203.78 170.        1203.25               180.        1203.32 190.        1203.23               200.        1202.87 220.        1203.31               240.        1203.06 260.        1202.96               280.        1202.87 300.        1203.23               320.        1204.92 340.        1208.38               360.        1211.93 380.        1214.49               400.        1211.53 420.        1224.54               440.          1193.3 460.        1193.68               480.        1255.14 500.          1262.5              520.          1262.5 540.        1257.26               560.        1266.67 580.        1261.41               600.        1257.26 WCAP- 16219-NP                                                       April 2004 6403-NP.doc-042904

A-19 PI Containment DBA Evaluation Model - MSLB Peak Pressure                                                         18 Feb/12004 11:44:57 GOTHIC Version 7.lPatchl (OA) - September 2003 File: D:\DATA\GOTHIC\Prairie Island\MSLB Cases\peak-press_...slb Function 3
Fan Cooler Heat Removal vs. Vapor Temp.
Ind. Var.: Vapor Sat. Temp. (F)
Dep. Var.: Heat Removal (BTU/s)
Ind. Var.        Dep. Var.          Ind. Var.       Dep. Var.
: 0.                0.           100.                500.
120.            2000.             140.             4500.
160.            6250.             180.             8250.
200.          11000.              220.            14000.
240.          17750.            240.1              9861.
260.          12500.              265.            13194.
270.          13611.            270.1                   0.
300.                 0.           500.                  0.
Function 4
Spray Flow Rate Ind. Var.: Time (sec)
Dep. Var.: Flow Rate (gpm)
Ind. Var.        Dep. Var.          Ind. Var.       Dep. Var.
: 0.                0.            0.1             1200.
3600.             1200.           3601.                   0.
1000000.                  0.
Control Variables CV                              Func. Initial       Coeff.      Coeff.                 Upd. Int.
      #      Description             Form       Value             G       aO       Min    Max    Mult.
1    Drop Mass                    ault          0.      le+006          0. -le+03 le+032         0.
2    Liquid Mass                  mult          0.      le+006          0. -le+03 le+032         0.
3    Heat Rate 1                  mult          0.             1.       0. -le+03 le+032         0.
4    Heat Rate 2                  mult          0.             1.       0. -le+03 le+032         0.
5    Heat Rate 3                  mult          0.             1.       0. -le+03 le+032         0.
6    Heat Rate 4                  mult          0.             1.       0. -le+03 le+032         0.
7    Heat Rate 5                  mult          0.             1.      0.  -le+03 le+032         0.
8    Heat Rate 6                  mult          0.             1.      0.  -le+03 le+032         0.
9    Heat Rate 7                  mult          0.             1.      0.  -le+03 le+032         0.
10     Heat Rate 8                  mult          0.             1.      0.  -le+03 le+032         0.
11     Heat Rate 9                  mult          0.             1.      0.  -le+03 le+032         0.
12     Heat Rate 10                 mult          0.             1.      0.  -le+03 le+032         0.
13     Heat Rate 11                 mult          0.             1.      0.  -le+03 le+032         0.
WCAP-16219-NP                                                                                             April 2004 6403-NP.doc-042904

L A-20 L
PI Containment DBA Evaluation Model - MSLB Peak Pressure Feb/1182004 11:44:57 GOTHIC Version 7.1Patchl (A) - September 2003 19 L
File: DA\DATA\GOTHIC\Prairie lsland\MSLB Casestpeakpressmsib                                                          L Control Variables       (cont.)                                            L CV                              Func. Initial       Coeff.      Coeff.                   Upd. Int.
LI
* Description             Form        Value        G            aO        Min    Max    Mult.               L 14     Heat  Rate   12            mult           0.          1.            0. -le+03 le+032         0.
I 15     Heat  Rate   13            mult           0.          1.            0. -le+03 le+032         0.
16     Heat  Rate   14            mult           0.          1.            0. -le+03 le+032         0.
17     Heat  Rate   15            rnult          0.          1.            0. -le+03 le+032         0.
18     Heat  Rate   16            mult           0.          1.            0. -le+03 le+032         0.
19     Heat  Rate   17            mult           0.          1.            0. -le+03 le+032         0.
20     Heat  Rate   18            mult           0.          1.            0. -le+03 le+032         0.            L.
21 22 Heat Heat Rate Rate 19 20 mult mult 0.
0.
1.
1.
0.
0.
                                                                                   -le+03
                                                                                   -le+03 le+032 le+032 0.
0.
I 23 24 Heat Heat Rate Rate 21 22 mult mult 0.
0.
1.
1.
0.
0.
                                                                                   -le+03
                                                                                   -le+03 le+032 le+032 0.
0.
I, 25     Heat  Rate   23            mult           0.          1.            0. -le+03 le+032         0.            L 26     Heat  Rate   24             rmult         0.          1.            0. -le+03 le+032         0.
27 28 Drop Drop Enth Q
Dif             sum mult 0.
0.
1.
1.
0.
0.
                                                                                   -le+03
                                                                                   -le+03 le+032 le+032 0.
0.
I, I-Function Components Control Variable        1 Drop Mass mult Y=G* (alXlla2X2*. ... *anXn)
Gothic-s                   Variable          Coef.
0                 Name                    location             a l                                Ad             cVl               1.
2                                Rd             cVl                1.
Function Components Control Variable 2 Liquid Mass mult Y=G* (alXl'a2X2. .      . anXnl Gothic-s                   Variable          Coef.
Name                    location             a 1                               Al              cVl                .
2                                Rl             cVl               1.
WCAP- 16219-NP                                                                                              April 2004 6403-NP.doc-042904

A-21 PI Containment DBA Evaluation Model - MSLB Peak Pressure                         20 Feb/18/2004 11:44:57 GOTHIC Version 7.1Patchl (OA) - September 2003 File: DADATA\GOTHIC\Prairie IsIand\MSLB Cases\peak-press-mslb Function Components Control Variable 3 Heat Rate 1 mult Y=G*IalXl*a2X2*... *anXn)
Gothics                   Variable         Coef.
       #                 Name                   location            a 1                         Qwv l)                ccl          1.
2                      Rmuls(l)                  ccl          1.
Function Components Control Variable 4 Heat Rate 2 mult Y=G*(alXl*a2X2*... *anXn)
Gothic_s                  Variable         Coef.
Name                   location           a 1                          Qwv(l)                cC2          1.
2                      Rmuls(l)                  cC2 l.
Function Components Control Variable 5 Heat Rate 3 mult Y=G*(alXl*a2X2*... *anXn)
Gothic-s                  Variable         Coef.
Name                   location            a 11                         Qwv(l)   I           cC3           1.
2                      Rmuls(l)                  cC3 l.
WCAP-16219-NP                                                             April 2004 6403-NP.doc-042904

A-22 L
Pi Containment DBA Evaluation Model - MSLB Peak Pressure                     21 Fet/18/2004 11:44:57 GOTHIC Version 7.1 Patchl(OA) - September 2003 File: D:\DATA\GOTH IC\Praine IsIand\MSLB Cases'peakpressmsIb L
L Function Components Control Variable 6                                    L Heat Rate 4 mult U
Y=G*(alXl*a2X2*... *anXn)                                   I Gothic-s                  Variable      Coef.
I Name                  location         a L
12 2
Qwv(l)
Rmuls(l) cC5 cC5 1.
1.
L Function Components Control Variable 7 Heat Rate 5 mult Y=C*(alXl*a2X2*... *anXn)
Gothics                   Variable      Coef.
       #                 Name                   location         a 1                       Qwv(l)               cC4        1.
2                      Rmuls(l)                cC4        1.
Function Components Control Variable 8 Heat  Rate 6 mult Y=G*(alXl*a2X2 ... *anXn)
Gothics                   Variable      Coef.
       #                 Name                   location         a 1                         Qwv(l)               cC7        1.
2                      Rmuls(l)                cC7        1.
WCAP-16219-NP                                                          April 2004 6403-NP.doc-042904

A-23 Pi Containment DBA Evaluation Model - MSLB Peak Pressure                       22 Feb/18/2004 11:44:57 GOTHIC Version 7.1 Patch1 (OA) - Septembor 2003 File: D:\DATA\GOTHIC\Prairie Island\MSLB Cases\peakjpress.mstb Function Components Control Variable 9 Heat Rate 7 mult Y=G* (alXl*a2X2*.. . *anXn)
Gothics                    Variable       Coef.
       #                 Name                    location         a 1                         Qwv(l)               cC8         1.
2                     Rmuls(l)                 cC8         1.
Function Components Control Variable 10 Heat Rate 8 mult Y=G* (alXl*a2X2* ... *anXn)
Gothics                    Variable       Coef.
4                 Name                    location         a 1                         Qwv(l)               cC9         1.
2                     Rmuls(l)                 cC9         l.
Function Components Control Variable 11 Heat Rate 9 mult Y=G* (alXl*a2X2* ... *anXn)
Gothics                    Variable       Coef.
       #                 Name                    location         a 1                         QV(l)              cclo          1.
2                     Rmuls(l)               cClO          1.
WCAP-1 6219-NP                                                           April 2004 6403-NP.doc-042904

L A-24 I
IL IL Pl Containment DBA Evaluation Model - MSLB Peak Pressure                        23 Feb/18/2004 11:44:57 GOTHIC Version 7.1 Patchl (OA) - September 2003 File: D:ADATA\GOTHIC\Prairie Island\MSLB Cases'peak-press-mslb Function Components Control Variable 12 Heat Rate 10 mult Y=G* (alXl*a2X2* ... .anXn)
Gothics                    Variable       Coef.
If                Name                    location          a 1                          Qwv(1               cCl           .
2                       Rmuls(l)               cC11          .
Function Components Control Variable 13 Heat Rate 11 mult Y=G*(alXl*a2X2* ... *anXn)
Gothics                    Variable       Coef.
       #f                Name                    location         a 1                          Qwv (1)            cC12          1.
2                      Rmuls(l)               cCl2          1.
Function Components Control Variable 14 Heat Rate 12 mult Y=G* (alXl*a2X2* .. . *anXn)
Gothics                    Variable       Coef.
f                 Name                    1ccation         a 1                          Qwv(1)             cC13          1.
2                      Rmuls(l)               cC13          1.
WCAP-16219-NP                                                            April 2004 6403-NP.doc-042904

A-25 Pi Containment DBA Evaluation Model - MSLB Peak Pressure                       24 Feb/t8/2004 11:44:57 GOTHIC Version 7.1 Patch1 (OA) - September 2003 File: D:%DATA\GOTHIC\Prairie Island\MSLB Cases\peak-press-mslb Function Components Control Variable 15 Heat Rate 13 mult Y=G*(alXl*a2X2*... *anXn)
Gothics                    Variable       Coef.
Name                   location          a 1                          Qwv l)             cC14          1.
2                       Rmuls(l)              cC14          1.
Function Components Control Variable 16 Heat Rate 14 mult Y=G*(alXl*a2X2*... *anXn)
Gothics                    Variable       Coef.
       #                 Name                    location          a 1                          QWV(1)             cC15          1.
2                       Rmuls(l)              cC15          1.
Function Components Control Variable 17 Heat Rate 15 mult Y=G*(alXl*a2X2*... *anXn)
Gothic-s                   Variable       Coef.
       #                 Name                    location          a 1                          QWV(l)             cC21          1.
2                       Rmuls(l)              cC21          1.
WCAP-16219-NP                                                            April 2004 6403-NP.doc.042904

I A-26                                                                                IF I
PI Containment DBA Evaluation Model - MSLB Peak Pressure                        25 Feb11812004 11:44:57 GOTHIC Version 7.1 Patchl (OA) - September 2003 File: DAOATAWGOTHIC\Prairie lsland\MSLB Cases'peak_ press-msIb Function Components Control Variable 18 Heat Rate 16 mult Y=G* (alXl*a2X2* ... *anXn)
Gothics                    Variable       Coef.
Name                    location         a 1                         Qwv(l)              cC23          l.
2                      Rmuls(l)               cC23          1.
Function Components Control Variable 19 Heat Rate 17 mult Y=G*(alXl*a2X2* ... -anXn)
Gothic-s                   Variable       Coef.
Name                    location         a 1                          Qwv(l)              cCl9          1.
2                      Rmuls(l)                cCl9          1.
Function Components Control Variable 20 Heat Rate 18 mult Y=G* (alXl*a2X2*. . . *anXn)
Gothics                    Variable       Coef.
Name                    location         a 1                          Qwv(l)              cC22          1.
2                      Rmuls l)                cC22          1.
WCAP-16219-NP                                                             April 2004 6403-NP.doc-042904

A-27 Pi Containment DBA Evaluation Model -MSLB Peak Pressure                       26 Feb18004 11:44:57 GOTHIC Version 7.1 Patch 1(OA) - September 2003 File: D:\DATA\GOTHIC\Prairie IsIand\MSLB Cases\peakpress-mstb Function Components Control Variable 21 Heat Rate 19 mult Y=G*(alXl*a2X2*... *anXn)
Gothics                    Variable      Coef.
* Name                   location         a I                         Qwv(1)             cC20         1.
2                     Rmuls(1)               cC20         1.
Function Components Control Variable 22 Heat Rate 20 mult Y=G*(alXl1a2X2*... *anXn)
Gothics                    Variable      Coef.
Name                    location        a 1                         QWV(l)             cC24         1.
2                     Rmuls(l)               cC24         1.
Function Components Control Variable 23 Heat Rate 21 mult Y=G*(alXl*a2X2... .*anXn)
Gothics                     Variable      Coef.
Name                    location         a 1                         QWV(l)             cC16         1.
2                     Rmuls(l)               cC16         1.
WCAP-1 621 9-NP                                                         April 2004 6403-NP.doc-042904

I A-28 U
I Pi Containment OBA Evaluation Model - MSLB Peak Pressure Feb/18/2004 11:44:57 GOTHIC Version 7.1 Patchl(OA) - Septomber 2003 27 I
File: DADATA\GOTHIC\Prairie Island\MSLB Cases\peak._pressmsIb L
Function Components                                      I, Control Variable 24 Heat Rate 22 mult I
Y=G*(alXl*a2X2*... *anXn)                                      L Gothics                   Variable       Coef.
Name                    location         a 1                         Qwv(l)              cC17          1.             L 2                     Rmuls l)               cC17          1.              I, Function Components Control Variable 25 Heat Rate 23 mult Y=G (alXlla2X2*...*anXn)
Gothics                   Variable       Coef.
       #                 Name                   location          a 1                        QWV(l)              cC18          1.
2                     Rmuls(l)               cC18          1.
Function Components Control Variable 26 Heat Rate 24 mult Y=G*(alXl*a2X2-... anXn)
Gothics                   Variable       Coef.
Name                   location         a 1                        Qwv(l)                cC6         1.
2                     Rmuls(l)                 cC6         1.
WCAP-16219-NP                                                            April 2004 6403-NP.doc-042904

A-29 Pi Containment DBA Evaluation Model - MSLB Peak Pressure                                                              28 Feb/18/2004 11:44:57 GOTHIC Version 7.1 Patch1 (QA) - September 2003 File: D:\DATA\GOTHIC\Prairie Island\MSLB Cases\peak.press..msIb Function Components Control Variable 27 Drop Enth Diff sum Y=G*(aO+alXl+a2X2+...+anXn)
Gothic_s                   Variable          Coef.
Name                    location              a 1                               Hd              cVl l.
2                           Bc_hd              cB3F             -1.
Function Components Control Variable 28 Drop Q mult Y=G*(alXl*a2X2*... *anXn)
Gothics                    Variable          Coef.
Name                    location             a 1                          Cvval             cv27                 .
2                          Wdjnc               cJ3               1.
Run Control Parameters         (Seconds)
Time         DT       DT           DT       End       Print     Graph       Max     Dump     Phs Chng Int        Min      Max         Ratio      Time        Int      Int        CPU      Int    Time Scale 1       0.001      0.05            1.      600.          5.        1. 1200.       0. DEFAULT Solution Options Time      Solution    Imp Conv        Imp Iter      Pres Sol       Pres Conv     Pres Iter   Differ    Burn Dom       Method         Limit         Limit         Method          Limit        Limit     Scheme   Sharp 1    lSEMI-IMP                 0.            1l    DIRECT     I           0.I           1l FOUP  I     0.
Run Options Parameter                              Value Restart Time (sec)                                      I          0 WCAP-16219-NP                                                                                                  April 2004 6403-NP.doc-042904

I A-30                                                                                           I, I1 Pi Containment DBA Evaluation Model - MSLB Peak Pressure                                  29  I, Feb/18/2004 11:44:57 GOTHIC Version 7.lPatchl(OA) - September 2003 File: D:\DATA\GOTHIC\Prairie Island\MSLB Cases~peak-press-mslb L
Run Options        (cont.)
I, Parameter                              Value                            LI Restart Time Step 4                                                0                        L Restart Time Control                                            NEW Revaporization Fracticn Fog Model DEFAULT OFF I
Maximum Mist Density(lbm/ft3)                              DEFAULT Drop Diam. From Mist(in)
Minimum HT Coeff.(B/h-ft2-F)
DEFAULT 0
I, Reference Pressure(psia)                                    IGNORE                           L Forced Ent.      Drop Diam.(in)                           0.09996 Vapor Phase Head Correction                                INCLUDE                           I Kinetic Energy Vapor Phase IGNORE INCLUDE I-Liquid Phase                                              INCLUDE                           I Drop Phase Force Equilibrium INCLUDE IGNORE I,
Drop-Liq. Conversion                                      INCLUDE                            I.
QA Logging                                                       OFF Debug Output Level                                                  0                        'I Restart Dump on CPU Interval               (sec)             3600.0 T1 1/
Graphs Graph                                                      Curve    Number Title              Mon        1            2          3    4  5 0     M & E Imbalance                     EM           EE 1     Containment Pre                    PRl 2     Containment Tem                    TVl 3     Sump Water Temp                    TLE 4     Sump Level                         LL1 5     Fan Cooler Heat                   CQlC 6     Spray Flow Rate                    FD3 7     Metal and Concr                    HAl          HA4 8     Drop and Liquid                    cvl          cv2 9     Heat Rate for C                    cv3         cvl0        cv15 cv20 10                                       cv28         cv27 11                                        ED1 12                                        TA1          TA2 WCAP-1 621 9-NP                                                                     April 2004 6403-NP.doc-042904

A-31 PI Containment DBA Evaluation Model- MSLB Peak Pressure                     30 Feb/i 8/2004 11:44:57 GOTHIC Version 7.1 Patchl (OA) - September 2003 Fie: D:ADATA\GOTHIC\Prairie IslandMSLB Cases\peak_pressnmsib Envelope Sets Set                                            Set       No.
No.               Description                 Type      Items Data Files File                                               Inter-  Output II               Name                Type       polate   Files 1      datal_case3rpo_               TIME          YES  SINGLE 2      data2_case3rpo_               TIME          YES  SINGLE 3      pressurecase3r                TIME          YES  SINGLE WVCAP-16219-NP                                                       April 2004 6403-NP.doc-042904

B-I APPENDIX B GOTHIC PRAIRIE ISLAND LOCA CONTAINMENT RESPONSE MODEL INPUT TABLES WCAP-16219-NP                                     April 2004 6403-NP.doc-042904

                              ---- I I
B-2 a, I I
I
                                    'I IL I,
I I.
I.
I, I
I, I,
I, I,
L IL I,
                                     .t I
I, U.
I 1I I,
I I
1U WCAP-16219-NP      April 2004        ',t 6403-NP.doc-042904
                                    -I'

B-3 ac WCAP-1 621 9-NP    April 2004 6403-NP.doc-042904

I B4 I
a,c     L t.1 L.
I I.
I.
I-I, I.1 I
It LI
                                      'I
                                      .1, I
1~
I
                                      .I1L
                                     .I It'
                                     'IL I-
                                      .1~
I I.
I
                                     .1I II I -
WCAP-16219-NP      April 2004, 6403-NP.doc-042904 1t

B-5 a,c April 2004 WCAP-1 6219-NP     April 2004 6403-NP.doc-042904

B-6                                 1.
a,c.
                                    .1~
L 1-
                                    ..1
                                    .1
                                    .1
                                    .1 I,
                                    .1~
                                    .1,
                                    .1
                                     .1~
                                    .1 I
I, I,
1 WCAP-162 19-NP      April 2004       l, 6403-NP.doc-04290 4

B-7 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

I B-8                              "L I
                                  'I I-I I-I.
L I,
I I
                                  .1~
I II I1, I,
I
                                  ,I, I,
I, I,
I, I
I, Lj I
                                  "I
                                  'I I,
I '
WCAP-16219-NP        April 2004    I 6403-NP.doc-042904
                                   'I

B-9 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-10 a,c I I,
                                  'I-,
I 1I 1,
L
                                  'L 1I I,
L It L
L1
                                  'L.
IL
                                  'L I,
I1 I,
                                  -11 I
                                   'L "L
I, I.
L, I
I I
I WCAP-1 6219-NP     April 2004      1, 6403-NP.doc-042904
                                  'L

B-Il wS        a,c K'
K' K>
i_'
K' K'
k-K' K,'A-61-N          pl20 KY3N~oc020 V>

I I,
B-12 a,c L I-I-
II I,
I I,
Id I,
I, I
I,
                                   'I
                                   'I I
I, I)
I,
                                   .1, I-It I
I I
I
                                   'I
                                   'I, I,
I, WCAP-1 6219-NP     Apr i112004     -1' 6403-NP.doc-042904 L

B-13 a,c April 2004 WCAP-1 6219-NP          April 2004 6403-NP.doc-042904

I, B-14                             1L a,c L L
F L
I, I,
I, L
I; I-F
                                  'I I
1.
I I,
I-I I,
                                  .1 I,
I, I
                                   'I, I
I, I,
                                  'I I
I WCAP-1 6219-NP     Apnil 2004       1 6403-NP.doc-042904

B-15 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-16                                   I a,c     L 1
I-I'
                                     .1 I)
I.
11 1~
                                    .1 I.
                                    'I I.1 I
I.
I I,
I, I,
                                    'I I,
I.t I,
I
                                    .1I, I,
                                      .1 I,
I,
__ I I
WCAPI16219-NP 6403-NP.doc-042904 April 2004 J
ItId

B-17 ac WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-18 axc WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-19 a,c WCAP-16219-NP      April ZUU4 6403-NP.doc-042904

IF B-20                                F I
I I
I
                                   'It ah,c I I
It I
I I
I I
I)
IL I;
I, IL It, It I
                                   'I I
WCAP-16219-NP      April 2004        J 6403-NP.doc-042904

B-21 ac WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-22                              I 1~
a,c  i
                                  'IL F
                                  'It IL IL J)
JF F
F IL
                                  .I April 2004 WCAP-16219-NP 6403-NP.doc-042904

B-23 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-24                              FL a, IL FL FL I
IL
                                  'IL I)
FL I)
IL IL I
Id I
                                  'IL I,
IL
                                  'I)
IL
                                  'I I,L I,
WCAP-1 6219-NP      April 2004 6403-NP.doc-042904I

v            B-25 v               a,c v-vY vY vY vY vY vY KY vY vY vY vY vY vY vY KY vY WCP1619NPArl20 vI60-Pdc020 v'

B-26 axc WCAP-16219-NP      April 2004 6403-NP.doc-042904

B3-27 v
a,c v
v v~
v-v-
v-v-
v-vY
'v/
vY vY vY  WCAP-1 6219-NP       Arl20 pl20 6403-NP.doc-042904 VY

B-28                                   F a,c 1
I I
J.I l
IL It
_              _ .1  IL WCAP-16219-NP        April 2004il 6403-NP.doc-042904                       1

B-29 a,c WCAP-1 621 9-NP    April 2004 6403-NP.doc-042904

B-30 ac April 2004 WCAP-16219-NP      April 2004 6403-NP.doc4042904

B-31 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-32     'IL a,t I)
It It IL IL I)
IL
         -1I, It I)
IL IL I'
It It I~t It IL I1
         .1 IL It'I, 11
         'It I
         "IL It

KYV B-33 a,c KY1 KYv KYl kY/
NKY KYV KY~
WCAP-16219-NP      April 2004 6403-NP.doc-042904

IL B-34                             I a,c I,
                                     .1~
I, "ti
                                  'I I-)
                                  "t I
1t I
IL
                                  '.t
                                  \.t I) 11 IL
                                  'Ii It I
                                  'I it It IJ I
IL It
                                  .,L it
                               -j :,L It WCAPd16219-NP      April 2004 6403-NP.doc-042904

B-35 axc WCAP-16219-NP      April 2004 6403-NP.doc-042904

It B-36I a,c    I 1)
I, I
                                      'IL I
IL I
I, It
                                      'I-I, It IL I-IL I)
                                      'I, I)
It IL
__  IL WCAP-1 6219-NP      April 2004 6403-NP.doc-042904I

B-37 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

lb NIL B-38                             I st
                                 'I I,
                                 'I 1'
I' 4L I/
I' MI I,
It 1)
                                  .11 VI
                                 'It It sIt I,
NI I'
It 1/
                                 't
                                 'IL It
                              -i .1 I-)
It J1.
NI 11 NI1 WCAPI16219-NP      April 2004     4,L 6403-NP.doc-042904 1t

B-39 a,c Api .ZUU4 WCAP-16219-NP      April 2Uo4 6403-NP.doc-042904

B40 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-41 a,c WCAP-1 6219-NP     April 2004 6403-NP.doc-042904

B-42 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

B43 ac I
K1.'
k-1 kl,~
KY-MY.1 WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-44 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

B45 axc I
WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-46                              I a,c
                                  -IL I,
It
                                  -I
                                  .I I,
VJo IitI, It J,
I' 4'
Ii It I.
                                   'I i4 4
A\}
utI
                                  'It I
I
                                   .J
                                  'IL' It A
                                  ,II I
I, 11
                                  ,I JiA J
WCAP-I 6219-NP     April 2004     4 6403-NP.doc-042904 Jit

B-47 a,c A - -! %%flrIA WCAP-16219-NP      Apnl suuq 6403-NP.doc-042904

B-48 a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

0 B-49
_            a,c C;
k-kI' k-k-
k_/
k-;
   -t K-'
K_
      -WCAP-1           6219-NP    April 2004 6403-NP.doc-042904

y B-50
                                     ;.1 a,c I,
                                     't 1~
It Is IL 1
                                      -I I'
IJ I
It 4,
I1 Is II.
I
                                     '1, I~i4 3'
I+
It I
J
                                      'I
                                      .,II It
                                      >1 3
_J.J WCAP-16219-NP      April 2004         j 6403-NP.doc-042904 J

B-51I a,c WCAP-16219-NP      April 2004 6403-NP.doc-042904

B-52 a,c April 2004 WCAP-1 6219-NP     April 2004 6403-NP.doc-042904

B-53 a,c April 2004 WCAP-1 621 9-NP WCAP-16219-NP      April 2004 6403-NP.doc 042904

B-54 ac WCAP-1 6219-NP     April 2004 6403-NP.doc-042904}}