ML19063B345

From kanterella
Revision as of 23:37, 4 May 2019 by StriderTol (talk | contribs) (Created page by program invented by StriderTol)
Jump to navigation Jump to search
NUREG/IA-511 - International Agreement Report - Simulation of ROSA-2 Test-2 Experiment: Application to Nuclear Power Plant.
ML19063B345
Person / Time
Issue date: 03/31/2019
From: Gallardo S, Lorduy M, Querol A, Kirk Tien, Verdu G
Office of Nuclear Regulatory Research, Universitat Politecnica de Valencia
To:
Meyd, Donald
References
NUREG/IA-0511
Download: ML19063B345 (65)


Text

NUREG/IA-05 11 International Agreement Report Simulation of ROSA-2 Test-2 Experiment

Application to Nuclear Power Plant Prepared by:

S. Gallardo, A. Querol, M. Lorduy and G. Verdú Universitat Politcnica de Valncia Instituto Universitario de Seguridad Industrial, Radiofísica y Medioambiental Camí de Vera s/n 46022 Valencia, SPAIN K. Tien, NRC Project Manager Division of Systems Analysis Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, DC 20555

-0001 Manuscript Completed: April 2018 Date Published

March 2019 Prepared as part of The Agreement on Research Participation and Technical Exchange Under the Thermal-Hydraulic Code Applications and Maintenance Program (CAMP)

Published by U.S. Nuclear Regulatory Commission

AVAILABILITY OF REFERENCE MATERIALSIN NRC PUBLICATIONS NRC Reference Material As of November 1999, you may electronically access NUREG-series publications and other NRC records at

NRC's Library at www.nrc.gov/reading-rm.html. Publicly released records include, to name a few, NUREG-series publications; Federal Register notices; applicant, licensee, and vendor documents and correspondence; NRC correspondence and internal memoranda; bulletins

and information notices; inspection and investigative

reports; licensee event reports; and Commission papers

and their attachments.

NRC publications in the NUREG series, NRC regulations, and Title 10, "Energy," in the Code of Federal Regulations may also be purchased from one of these two sources. The Superintendent of Documents Mail Stop IDCC

Washington, DC 20402-0001

Internet:

bookstore.gpo.gov

T elephone: (202) 512-1800

Fax: (202) 512-2104 The National Technical Information Service 5301 Shawnee R dAlexandria, VA 22312-0002 www.ntis.gov800-553-6847 or, locally, (703) 605-6000A single copy of each NRC draft report for comment isavailable free, to the extent of supply, upon written

request as follows:

Address: U.S. Nuclear Regulatory Commission Washington, DC 20555-0001 E-mail: distribution.resource@nrc.gov

Facsimile: (301) 415-2289 Some publications in the NUREG series that are posted at NRC's Web site address www.nrc.gov/reading-rm/

doc-collections/nuregs are updated periodically and may differ from the last printed version. Although references to material found on a Web site bear the date the material

was accessed, the material available on the date cited

may subsequently be removed from the site.

Non-NRC Reference Material Documents available from public and special technical libraries include all open literature items, such as books, journal articles, transactions, Federal Register notices, Federal and State legislation, and congressional reports.

Such documents as theses, dissertations, foreign reports

and translations, and non-NRC conference proceedings

may be purchased from their sponsoring organization.

Copies of industry codes and standards used in a substantive manner in the NRC regulatory process are

maintained at-The NRC Technical Library Two White Flint North 11545 Rockville Pike

Rockville, MD 20852-2738 These standards are available in the library for reference

use by the public. Codes and standards are usually

copyrighted and may be purchased from the originating organization or, if they are American National Standards, from-American National Standards Institute 11 West 42nd Street New York, NY 10036-8002 www.ansi.org (212)642-4900 Legally binding regulatory requirements are stated only in laws; NRC regulations; licenses, including technical

speci views expressed in contractor prepared publications in

this series are not necessarily those of the NRC.

The NUREG series comprises (1) technical and adminis-trative reports and books prepared by the staff (NUREG-

XXXX)or agency contractors (NUREG/CR-XXXX), (2) proceedings of conferences (NUREG/CP-XXXX), (3) reports

resulting from international agreements (NUREG/IA-XXXX),

(4)brochures (NUREG/BR-XXXX), and (5) compilations of legal decisions and orders of the Commission and Atomic and Safety Licensing Boards and of Directors' decisions under Section 2.206 of NRC's regulations (NUREG-0750).

DISCLAIMER:

This report was prepared under an international cooperative agreement for the exchange of

technical information. Neither the U.S. Government nor any

agency thereof, nor any employee, makes any warranty, expressed or implied, or assumes any legal liability or

responsibility for any third party's use, or the results of such

use, of any information, apparatus, product or process

disclosed in this publication, or represents that its use by

such third party would not infringe privately owned rights.

NUREG/IA-051 1 International Agreement Report Simulation of ROSA-2 Test-2 Experiment

Application to Nuclear Power Plant Prepared by:

S. Gallardo, A. Querol, M. Lorduy and G. Verdú Universitat Politcnica de Valncia Instituto Universitario de Seguridad Industrial, Radiofísica y Medioambiental Camí de Vera s/n 46022 Valencia, SPAIN K. Tien, NRC Project Manager Division of Systems Analysis Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, DC 20555

-0001 Manuscript Completed: April 2018 Date Published

March 2019 Prepared as part of The Agreement on Research Participation and Technical Exchange Under the Thermal

-Hydraulic Code Applications and Maintenance Program (CAMP)

Published by U.S. Nuclear Regulatory Commission

iii ABSTRACT The purpose of this work is to overview the results provided by a first approach to the simulation of an Intermediate Break Loss

-Of-Coolant Accident (IBLOCA) in a 3

-loop PWR Nuclear Power Plant (NPP) using the thermal

-hydraulic code TRACE5 patch 2 and the Symbolic Nuclear Analysis Packages software (SNAP) version 2.1.2.

The IBLOCA transient applied to the standard PWR TRACE5 model is the Test 2 (IB

-CL-03) handled at the Large Scale Test Facility (LSTF) in the frame of the OECD/NEA ROSA

-2 Project. Test 2 simulates a 17% cold leg IBLOCA under the assumption of the single

-failure of High Pressure Injection and Low Pressure Injection systems and total failure of the Auxiliary Feedwater.

The LSTF is a Full Height Full Pressure (FHFP) facility designed to simulate a 4-loop W-type PWR (Tsuruga unit II NPP). The volumetric scaling factor is 1/48. The four primary loops of the reference PWR are scaled in LSTF by two equal

-volume loops. The core power used to simulate the decay power is 10 MW, corresponding to 14% of the 1/48 volumetrically scaled reference PWR rated power.

The simulation results are provided throughout several graphs, where the main system variables, such as pressures, pressure vessel liquid levels and temperatures are shown. These results represent a contribution to assess the predictability of computer codes such as TRACE5.

v FOREWORD Thermalhydraulic studies play a key role in nuclear safety. Important areas where the significance and relevance of TH knowledge, data bases, methods and tools maintain an essential prominence are among others:

assessment of plant modifications (e.g., Technical Specifications, power uprates, etc.);

analysis of actual transients, incidents and/or start

-up tests; development and verification of Emergency Operating Procedures; providing some elements for the Probabilistic Safety Assessments (e.g., success criteria and available time for manual actions, and sequence delineation) and its applications within the risk informed regulation framework; training personnel (e.g., full scope and engineering simulators); and/or assessment of new designs.

For that reason, the history of the involvement in Thermalhydraulics of CSN, nuclear Spanish Industry as well as Spanish universities, is long. It dates back to mid 80's when the first serious talks about Spain participation in LOFT

-OCDE and ICAP Programs took place. Since then, CSN has paved a long way through several periods of CAMP programs, promoting coordinated joint efforts with Spanish organizations within different periods of associated national programs (i.e., CAMP-Espana). From the CSN perspective, we have largely achieved the objectives. Models of our plants are in place, and an infrastructure of national TH experts, models, complementary tools, as well as an ample set of applications, have been created. The main task now is to maintain the expertise, to consolidate it and to update the experience. We at the CSN are aware on the need of maintaining key infrastructures and expertise, and see CAMP program as a good and well consolidated example of international collaborative action implementing recommendations on this issue.

Many experimental facilities have contributed to the today's availability of a large thermal

-hydraulic database (both separated and integral effect tests). However, there is a continuous need for additional experimental work and code development and verification, in areas where no emphasis have been made along the past. On the basis of the SESAR/FAP 1 reports "Nuclear Safety Research in OECD Countries: Major Facilities and Programmes at Risk" (SESAR/FAP, 2001) and its 2007 updated version "Support Facilities for Existing and Advanced Reactors (SFEAR) NEA/CSNI/R(2007)6", CSNI is promoting since the beginning of this century several collaborative international actions in the area of experimental TH research. These reports presented some findings and recommendations to the CSNI, to sustain an adequate level of research, identifying a number of experimental facilities and programmes of potential interest for present or future international collaboration within the nuclear safety community during the coming decade. The different series of PKL, ROSA and ATLAS projects are under these premises. CSN, as Spanish representative in CSNI, is involved in some of these research activities, helping in this international support of facilities and in the establishment of a large network of international collaborations. In the TH framework, most of these actions are either covering not

1 SESAR/FAP is the Senior Group of Experts on Nuclear Safety Research Facilities and Programmes of NEA Committee on the Safety of Nuclear Installations (CSNI

vi enough investigated safety issues and phenomena (e.g., boron dilution, low power and shutdown conditions, beyond design accidents), or enlarging code validation and qualification data bases incorporating new information (e.g., multi

-dimensional aspects, non

-condensable gas effects, passive components).

This NUREG/IA report is part of the Spanish contribution to CAMP focused on:

Analysis, simulation an d investigation of specific safety aspects of PKL2/OECD and ROSA2/OECD experiments.

Analysis of applicability and/or extension of the results and knowledge acquired in these projects to the safety, operation or availability of the Spanish nuclear power plants. Both objectives are carried out by simulating the experiments and conducting the plant application with the last available versions of NRC TH codes (RELAP5 and/or TRACE).

On the whole, CSN is seeking to assure and to maintain the capability of the national groups with experience in the thermalhydraulics analysis of accidents in the Spanish nuclear power plants. Nuclear safety needs have not decreased as the nuclear share of the nations grid is expected to be maintained if not increased during next years, with new plants in some countries, but also with older plants of higher power in most of the countries. This is the challenge that will require new ideas and a continued effort.

________________________________

Rosario Velasco García, CSN Vice-president Nuclear Safety Council (CSN) of Spain

vii TABLE OF CONTENTS ABSTRACT ................................

................................................................

...............................

iii FOREWORD ................................................................................................

..............................

v LIST OF FIGURES

................................................................................................

....................

ix LIST OF TABLES

................................

................................................................

......................

xi EXECUTIVE

SUMMARY

................................

................................................................

........ .xiii ACKNOWLEDGMENTS

................................................................................................

...........

x v ABBREVIATIONS AND ACRONYMS

................................

....................................................

xv ii 1 INTRODUCTION

................................................................................................

...............

1-1 2 LSTF DESCRIPTION

................................................................................................

........ 2-1 3 LSTF TRANSIENT

................................

................................................................

............

3-1 4 THE STANDARD 3

-LOOP PWR MODEL

................................

......................................... 4-1 5 SCALING CONSIDERATIONS

................................

.........................................................

5-1 6 RESULTS AND DISCUSSION

................................

..........................................................

6-1 6.1 Transient Conditions ...............................................................................................

6-1 6.2 Steady-State ................................

................................................................

...........

6-1 6.3 Transient ................................

................................................................

.................

6-2 6.4 System Pressures ................................

................................................................

... 6-2 6.5 Break Mass Flow Rate

................................

............................................................

6-3 6.6 Primary Loop Mass Flow Rate s ..............................................................................

6-4 6.7 Pressure Vessel Collapsed Liquid Levels

................................................................

6-5 6.8 Temperatures

................................

................................................................

..........

6-6 6.9 Hot and Cold Legs Liquid Levels

................................

.............................................

6-9 6.10 U-Tubes Collapsed Liquid Level

................................

...........................................

6-10 6.11 Study of the ECCS Actuation ................................

................................................

6-11 6.12 Void Fraction................................

................................

................................

......... 6-15 6.13 LSTF Model vs Standard PWR M odel ................................

.................................. 6-20 7 CONCLUSIONS

................................

................................................................

................

7-1 8 REFERENCES

................................

................................................................

..................

8-1

ix LIST OF FIGURES Figure 1 Schematic View of the LSTF Facility ..................................................................... 2-1 Figure 2 LSTF Experimental Data

Primary a nd Secondary Pressures and PCT ................

3-2 Figure 3 Standard 3-Loop PWR Model Nodalization

................................

...........................

4-2 Figure 4 Primary a nd Secondary Pressures

. ................................

....................................... 6-3 Figure 5 Break Mass Flow Rate

................................

..........................................................

6-4 Figure 6 Primary Loop Mass Flow Rates

................................

............................................

6-5 Figure 7 Core a nd Downcomer Collapsed Liquid Levels

.....................................................

6-6 Figure 8 Maximum Fuel Rod Surface a nd Core Exit Temperatures

................................

.... 6-7 Figure 9 Fuel Rod Surface Temperature a t Different Axial Positions

................................

... 6-8 Figure 10 Fuel Rod Surface Temperature for Different HTSTR .............................................

6-8 Figure 11 Hot Legs Collapsed Liquid Level

...........................................................................

6-9 Figure 12 Cold Legs Collapsed Liquid Level

................................

....................................... 6-10 Figure 13 Steam Generator U

-Tube Up-Flow Side Collapsed Liquid Levels

.......................

6-11 Figure 14 Primary a nd Secondary Pressures

................................

...................................... 6-12 Figure 15 Core Collapsed Liquid Levels

................................

.............................................

6-12 Figure 16 Downcomer Collapsed Liquid Levels

................................

................................

.. 6-13 Figure 17 Maximum Fuel Rod Surface and Core Exit Temperatures

................................

.. 6-13 Figure 18 Accumulator Injection System Mass Flow Rate

................................

...................

6-14 Figure 19 High Pressure Injection System Mass Flow Rate

................................

................

6-14 Figure 20 Void Fraction i n 3-Loop PWR Plant at 0 s

................................

...........................

6-16 Figure 21 Void Fraction i n 3-Loop PWR Plant at 25 s ................................

.........................

6-17 Figure 22 Void Fraction i n 3-Loop PWR Plant at 200 s .......................................................

6-18 Figure 23 Void Fraction i n 3-Loop PWR Plant at the End o f the Transient

..........................

6-19 Figure 24 Primary and Secondary Pressures

...................................................................... 6-20 Figure 25 Break Mass Flow Rate Obtained i n LSTF and i n t he 3-Loop PWR Plant Scaled

. 6-21 Figure 26 Core Collapsed Liquid Levels

.............................................................................

6-22 Figure 27 Downcomer Collapsed Liquid Levels

................................

.................................. 6-22 Figure 28 Maximum Fuel Rod Surface and Core Exit Temperatures

................................

.. 6-23 Figure 29 Hot Legs Collapsed Liquid Levels

................................

....................................... 6-24 Figure 30 Cold Legs Collapsed Liquid Levels

..................................................................... 6-24

xi LIST OF TABLES Table 1 Control Logic and Sequence of Major Events in the LSTF Experiment

................. 3-1 Table 2 Main Characteristics of LSTF, 4-Loop NPP and 3-Loop NPP ................................ 5-3 Table 3 Steady- tate Conditions i n the TRACE5 M odel ....................................................6-1 Table 4 Chronological Sequence of Events ....................................................................... 6-2

xiii EXECUTIVE

SUMMARY

The purpose of this work is to test the capability of the thermal

-hydraulic code TRACE5 patch 2 in a first approach to the simulation of a cold leg Intermediate Break LOCA (IBLOCA) in a standard PWR plant of three loops. For this goal, a TRACE5 model of a standard 3

-loop PWR has been adapted to reproduce the Test 2 transient, in the frame of the OECD/NEA ROSA

-2 Project. Test 2 (IB-CL-03) was conducted in the Large Scale Test Facility (LSTF) of the Japan Atomic Energy Agency (JAEA). LSTF is a Full Height Full Pressure (FHFP) facility reproducing Tsuruga unit II Nuclear Power Plant (W-type 4-loop PWR). It allows to preserve time, power and mass inventory as in the LSTF, because the fluid exhibits the same properties at full pressure.

Test 2 simulates a 17% cold leg IBLOCA using an upward break nozzle mounted with the cold leg inner surface. Single

-failure of both High Pressure Injection and Low Pressure Injection systems and total failure of the Auxiliary Feedwater are assumed. Results of the simulation using TRACE5 are shown in several graphs, including primary and secondary pressures, discharged inventory, primary mass flow rates, and collapsed liquid levels (in the pressure vessel, hot legs, steam generators U

-tubes, etc.).

xv ACKNOWLEDGMENTS This paper contains findings that were produced within the OECD

-NEA ROSA-2 Project. The authors are grateful to the Management Board of the ROSA

-2 Project for their consent to this publication, and thank the Spanish Nuclear Regulatory Body (CSN) and the Asociación Espanola d e la Industria Eléctrica (UNESA) for the technical and financial support.

xvii ABBREVIATIONS AND ACRONYMS A FW Auxiliary Feedwater AIS Accumulator Injection System AM Accident Management BE Best Estimate CAMP Code Assessment and Management Program CET Core Exit Temperature CPU Central Processing Unit CRGT Control Rod Guide Tubes CSN Nuclear Safety Council, Spain DBE Design Basis Event ECCS Emergency Core Cooling System HPI High Pressure Injection IBLOCA Intermediate Break Loss

-Of-Coolant Accident JAEA Japan Atomic Energy Agency JAERI Japan Atomic Energy Research Institute JC Jet Condenser LOCA Loss-Of-Coolant Accident LPI Low Pressure Injection LSTF Large Scale Test Facility MFW Main Feedwater MSIV Main Steam Isolation Valve NPP Nuclear Power Plant NRC U.S. Nuclear Regulatory Commission NV Normalized to the Steady State Value PA Auxiliary Feedwater Pump PCT Peak Cladding Temperature PF Feedwater Pump PGIT Primary Gravity Injection Tank PJ High Pressure Charging Pump PL High Pressure Injection Pump PORV Pilot Operated Relief Valve PV Pressure Vessel PWR Pressurized Water Reactor PZR Pressurizer RHR Residual Heat Removal RV Relief Valve SBLOCA Small Break Loss

-Of-Coolant Accident SG Steam Generator SI Safety Injection SNAP Symbolic Nuclear Analysis Package SRV Safety Relief Valve ST Storage Tank

1-1 1 INTRODUCTION Since the accident a t Three Mile Island Unit

-2 (TMI-2) i n 1979, many test facilities hav e produced experimental d ata on S mall B reak Los s-O f-C oolant A ccidents (SBLOCA) for c ode assessment and dev elopment. H owever, the ex perimental d ata for Intermediat e Break Los s-O f-C oolant A cc idents (IBLOCA) i s quite l imited although t he t hermal-h ydraulic r esponses c ould differ bet ween SBLOCA and IBLOCA.

B ecaus e full-s cale testing is us ually i mpossibl e to perform, there i s a need t o test sca l ed models o f prototype systems.

With this aim, t he Japan Atomic Energy Agency (JAEA) i n the frame o f the OECD/NEA ROSA

-2 Project prepared detailed thermal-h y draulic dat a. The pur pose o f t his w ork i s t o t est t he capability of t he thermal-h y draulic c ode TRACE5 patch 2 [1, 2] i n the s imulation of a cold leg IBLOCA i n a standard 3-l oop PWR pl ant (representative of Spanish NPPs).

Fo r this goal, a generic 3-l oop PWR TRACE5 model has been adapted t o reproduce t he Test 2 , c onducted i n the Lar ge S cale Test Facility (LSTF) [4] i n t he frame o f t he OECD/NEA RO SA-2 Project.

Test 2 (IB-C L-0 3) [3] s imulates a PWR 17% c old l eg IBLOCA under t he as sumption o f a s ingle-f ailure of t he High Pres sure Injection (HPI) and th e Low P ress ure Injection (LPI) systems and the total f ailur e of t he Auxiliary Feedwater (AFW) system. LSTF is a Full Height Ful l Pressure (FHFP) facility designed t o simulat e Tsurug a unit I I NPP, which is a 4-l oop W-type PWR of 3423 MWt. The volumetric scaling factor i s 1/48. T he four primary l oops o f t he reference PWR a r e represented by two loops in LSTF. The initial core power used to simulate the decay core pow er i s 10 MW, corresponding t o 14% o f the 1/48 volumetrically scaled reference PWR rated power.

A generic 3-loop PWR TRACE5 model has been adapted applying the power

-to-volume scaling strategy t o simulate t he transient conducted i n the LSTF. Th e simulation results are shown in different graphs, including primary and secondary pressures, discharged inventory, primary mass flow rates, liquid levels (in the pressure vessel, hot leg, steam generators U-tubes, etc.) and temperatures.

2-12 LSTF DESCRIPTION The LSTF simulates a W-type PWR reactor, of four loops and 3423 MW of thermal power. It is characterized using prototypical

-scaled components with full

-height, 1/48 scaled volume and full-pressure conditions to its reference PWR plant (Tsuruga unit II NPP). The four primary loops of the reference PWR are represented by two loops with a volume factor of 1

/24. Figure 1 shows a schematic view of the LSTF facility.

Figur e 1 Schematic View of the LSTF Facility The main features of LSTF are the following:

Elevations: preserved.Volumes: scaled by 1/48 to the reference PWR.Flow area: scaled by 1/48 in the pressure vessel and by 1/24 in the steam generators.The flow area in the hot and cold legs is scaled to conserve the ratio of the length (L) toCore power: scaled by 1/48 and limited to 14% of the scaled core power of the reference PWR.Fuel assembly: The total number of fuel rods is scaled by 1/48. There are 1008 heated rods.The pr imary c oolant system o f LS TF c onsists o f the primary l oop A w ith the pressurizer (PZR) and the s ymmetrical p rimary l oop B. B oth l oops i nclude a primary c oolant pum p (PC) a nd 141 U-t ubes f or each steam g enerator. T he secondary-co olant s ys tem i ncludes a jet c ondenser (J C), a feed w ater pu m p (PF), t he aux iliary feedwater pumps (PA) and t he pi ping system.The Emergency Core Cooling System (ECCS) consists of the f ollowing sub-s ystems: the high-p ressure charging pump (PJ), the high-p ressure injection pump (PH), the accumulators (ACC), the low-p ressure injection pump (PL), the Residual Heat Removal (RHR) sys tem and the 2-2 primary gravity injection tank (PGIT). The primary coolant discharged from the primary system is stored in the break flow Storage Tank (ST).

The pressure vessel has five regions: the upper head located above the upper core support plate; the upper plenum situated between the upper core support plate and the upper core plate; the active core; the lower plenum and the downcomer annulus region that surrounds the core and the upper plenum. The LSTF vessel includes 8 spray nozzles (of 3.4 mm inner

-diameter) in the upper head, and 8 Control Rod Guide Tubes (CRGTs), which lead the flow path between the upper head and the upper plenum.

3-1 3 LSTF TRANSIENT Test 2 conducted in the LSTF simulates a 17% cold leg IBLOCA under the assumption of single-failure of the ECCS (HPI, AIS and LPI systems) and the total failure of the AFW. The complete control logic of the transient is listed in Table 1. The system pressures and the PCT evolution are shown in Figure 2. The break unit is connected to the cold leg and the orifice flow area corresponds to 17% of the volumetrically

-scaled cross

-sectional area of the reference PWR cold leg. The transient starts at time 0 with the break valve opening in the cold leg (in the loop without pressurizer) and increasing the rotational speed of the reactor coolant pumps.

F ew seconds later, the scram signal is generated. This signal produces the initiation of the core power decay curve, the initiation of the primary coolant pumps coastdown, the turbine trip, the closure of MSIV and the termination of the MFW.

When the primary pressure is lower than 12.27 MPa, the Safety Injection (SI) signal is generated. The HPI system is activated few seconds after the SI signal generation, only in the loop with pressurizer. The accumulators actuate when the primary pressure is lower than 4.5 MPa. The core power is automatically decreased by the LSTF Core Protection System when the maximum fuel rod surface temperature reaches 958 K. When the primary pressure is lower than 1.24 MPa, the LPI system actuates in the loop with pressurizer. Test 2 finishes with the break valve closure when primary and secondary pressures are stabilized.

Table 1 Control Logic and Sequence of Major Events i n the LSTF Experiment Break. Time zero Generation of scram signal.

Primary pressure = 12.97 MPa.

Initiation of core power decay curve simulation.

Generation of scram signal.

Initiation of Primary Coolant Pump coastdown. Generation of scram signal.

Turbine trip (closure of steam generators Main Steam Stop Valves).

Generation of scram signal.

Closure of Steam Generator Main Steam Isolation Valves (MSIVs).

Generation of scram signal.

Termination of Steam Generator Main Feedwater.

Generation of scram signal.

Generation of Safety Injection (SI) signal Primary pressure = 12.27 MPa. Initiation of High Pressure Injection (HPI) in loop w/pressurizer.

12 seconds after SI signal.

Initiation of the Accumulator Injection in loop w/pressurizer.

Primary pressure = 4.51 MPa.

Initiation of Low Pressure Injection (LPI) in loop w/pressurizer.

Primary pressure = 1.24 MPa.

3-2 Figure 2 LSTF Experimental Data: Primary and Secondary Pressures and PCT 050100150200250300350400450500 0 2 4 6 8 10 12 14 16 Pressure (Pa) 10 6 Primary Pressure Secondary Pressure050100150200250300350400450500 Time (s)400 500 600 700 800 900 1000 Temperature (K)

PCT 4-1 4 THE STANDARD 3

-LOOP PWR MODEL The pl ant o f reference used in this st udy i s a standard 3-loop PWR with a thermal power o f 2785 MWt. The TRACE5 model i ncludes t he pressure vessel, ho t and cold legs, pressurizer (in loop 3), reactor coolant pumps, loop seals, HPI, A I and LP I systems, stea m generators, etc. The pressure vessel has been modeled using a 3-D VESSEL component available in TRACE5.

The VESSEL component contains 16 axial levels, 4 radial rings and 4 azimuthal sectors. Three rings ar e use d t o simulate the upper head, upper plenum, active core and lower pl enum o f t he pressure vessel. The fourth ring i s use d t o simulate the downcomer annulus region. Low er plenum i s located between the axial levels 1 and 2. Axial levels 3 to 10 model the active core.

Upper pl enum i s simulated between the axial levels 11 and 14.

The levels 15 and 16 model the upper head. T he control rod guide tubes (CRGTs) ar e modeled using 6 PI PEs: 3 PIPEs allow t he flow pat h between upper pl enum and upper head (connect levels 13 to 16), while the other 3 connect t he core exit (level 11) with th e upper head (level 16). The pow er i s supplied to the active cor e using 12 Heat Structur e components (HTSTR), which simulat e t he fuel rods present i n the standard 3-loop PWR plant. The core includes 157 fuel assemblies with 17x17 lattice design.

A POWER component has been used to manage t he power from these HTSTRs.

I n this case, the power has been characterized by m eans o f a point kinetics model with constant initial reactivity and table-input reactivity as a function of time after a trip. The hot legs hav e been modeled with 3 PIPE components connect ed t o the pressure vessel and the U-tubes. I n l oop 3, the hot leg ha s a cross flow connection with the surge line of t he pressurizer.

A PIPE component i s u se d to simulate the pressurizer.

The pressurizer Relief Valves (RVs) and the Pilot Operated Relief Valve (PORV) are also considered. The cold legs hav e been simulated with 3 PIPEs link ed t o the pressure vessel, t he reactor coolant pumps and the ECCS. I n t he broken l oop (loop 1), the cold leg presents a cross flow connection with a VALVE component joined with a BREAK component t o simulat e the atmospheric conditions. In this valve, t he chok ed flow coefficients hav e be en fixed to default values (1.0 for subcooled and two-phase coefficients). Th e HPI and LP I systems have been simulated with FILL components. The AIS i s modeled using PIPE components with "accumulator" option. The discharge lines ar e modeled with check VALVE components connected to the col d legs. The reactor coolant pumps ar e simulated with PUMP components, considering specific head and torque homologous curves. The reactor coolant pumps a r e joi ned t o the loop seal.

The U-tubes o f each steam generator hav e been modeled using a PIPE component. The secondary side consists of a boiler, a downcomer, a steam separator and the steam lines. The boiler and the downcomer hav e been simulated with PIPE components, while a SEPD component has been used to m odel the steam separator.

T he relief and isolation valves have been modeled using VALVE components. The Relief Valve is linked t o a BREAK component, while the Main Steam Isolation Valve is connected to a PIPE component, which simulates the steam l ine header.

4-2standard 3-l oop PWR TRACE5 model us ing t he Symbolic N uclear A nalysis P ackag e (SNAP) [5] software version 2.1.2. Figure 3 Standard 3-Loop PWR Model Nodalization

5-1 5 SCALING CONSIDERATIO NS In LOCA scenarios, the most important consideration of the scaling criteria is to preserve both the power and the coolant mass inventory during the transient. The power

-to-volume scaling criterion is frequently used to preserve time, power and mass inventory in FHFP facilities regarding reference NPPs because they have the same fluid properties at full pressure. To perform a power

-to-volume analysis of the NPP model, scaling considerations should be made. In this frame, scaling factors between the LSTF and 3-loop NPP must be evaluated to assess the viability of the scaling analyses and can be applied to define the boundary conditions of the scaled model.

As LSTF is a FHFP facility of an actual NPP (Tsuruga unit II) and the scenario is an IBLOCA in the cold leg, the power

-to-volume scaling criterion has been chosen to develop a scale

-up TRACE5 model of LSTF. This criterion results from the application of conservation equations (1) to (4) under some requirements and implications:

Continuity equation:

+()=0 (1) Momentum equation:

+= 1 1()(2) Energy equation:

+=()+++ (3) State equation:

=( ,) (4) where is the density, u is the velocity, x is the coordinate, t is the time, p is the pressure, F is the friction coefficient, h is the enthalpy, thermal expansion coefficient and is the power density.

Substituting the next dimensionless parameters (denoted by an asterisk) in equations (1) to (4): = , = ,= , = , = , = ,= , = ,= gives a set of nondimensionalized equations (5) to (8):

+()=0 (5) 5-2+= 1 1()(6) +=()+++(7) =(,) (8) These equations contain other dimensionless parameters such as , and. The first two parameters are Froude and Euler numbers, respectively. The third is known as heat source number following Ishii and Kataoka terminology. Power-to-volume scaling method requires that all dimensionless parameters in equations (5) to (8)have to be equal in the LSTF model and the scale

-up LSTF model.====1 (9) Considering that s imilarity bet ween both s ystems has been achieved, t he following pow er-t o-v olume relations a re obtained: ==== (10) being power, Q mass flow rate, V volume, A area and the volumetric scaling factor.

Furthermore, in the scale

-up model, the Froude number is conserved in horizontal components. It implies varying the diameter and length of these components. Trying to conserve the Froude number, from the scale

-up mass flow rate calculated as Eq. (11), the scale

-up diameter, D, can be obtained as Eq. (12):

  • 4*=* 4**(11) = */ (12) where U is velocity in the scale

-up model, is the coolant density and d is LSTF diameter. Furthermore, from the volume equation Eq. (13) and trying to conserve the Froude number, lengths of the scale

-up piping system are obtained as Eq. (14):

  • 4*=* 4**(13) =*/ (14) being l and L, LSTF and scale-up length, respectively.

Considering the nominal power and volume of LSTF, Tsuruga unit II and the 3-loop NPP, the scaling ratios between them can be obtained. Table 2 summarizes the main characteristics of LSTF, 4-loop and 3

-loop NPP models.

5-3 I n terms of pressure and heights, the relation between LSTF and its reference NPP is 1

/1, so the initial value for pressures and temperatures in the 3

-loop PWR model is also the same than in the experiment performed in LSTF. The pressure of accumulators is the same in all cases.

The break size considered in Test 2 has been scaled applying the power

-to-volume criterion between LSTF, Tsuruga II and the 3-loop PWR NPP.

Table 2 Main Characteristics of LSTF, 4

-Loop NPP and 3-Loop NPP Parameter LSTF Reference 4

-loop NPP 3-loop NPP Kv (3-loop/LSTF)

Primary pressure (MPa).

15.5 15.5 15.5 1 Core power (MW).

10 3423 2785 39.0 (14% of 2785 MW). Number of loops. 2 4 3 Number of fuel assemblies

. 193 157 Fuel assembly array

. 17x 17 17x17 Core height (m)

. 3.66 3.66 3.66 1 Vessel volume (m 3). 2.75 131.7 106 38.54 Number of fuel rods.

1008 50952 41447 41.11 Volume of pressurizer (m

3) 1.2 51.0 39.7 33.08 Pressurizer heaters power (kW). 124 1800 1400 Average length U

-tubes (m).

20.2 20.2 20.2 1 Hot leg inner diameter (m).

0.207 0.737 0.737 8.11 8.11 8.11 1 Cold leg inner diameter (m).

0.207 0.698 0.698 Break area: 17% cold leg flow area reference NPP (m 2). 1.36 10-3 0.065 0.048 35.2 Accumulator pressure (MPa).

4.6 4.6 4.6 1

6-1 6 RESULTS AND DISCUSSION 6.1 Transient Conditions The following assumptions are imposed to reproduce Test 2 in the NPP model:

1.Break size is 17% of the cold leg flow area.

2.An upward long break nozzle is located on top of the cold leg in the broken loop withoutpressurizer (loop 1).

3.Loss of off

-site power concurrent with the scram.

4.HPI, AI and LPI systems are activated in intact loops.

5.Non-condensable gas inflow from accumulator tank may take place.

6.Total failure of AFW.

6.2 Steady- tate Table 3 s hows t he s teady-s tate conditions achiev ed using TRACE5.

Table 3 Steady-St ate C onditions i n t he T RACE5 Mo del Item 3-loop NPPTRACE5 model Core Power (MW t). 2785 Hot leg Fluid Temperature (K).

600 Cold leg Fluid Temperature (K).

569 Mass Flow Rate (kg/s).

4416 Pressurizer Pressure (MPa).

15.5 Pressurizer Liquid Level (m).

3.65 Accumulator System Pressure (MPa).

4.6 SG Secondary

-side Pressure (MPa).

6.8 SG Secondary

-side Liquid Level (m).

5.7 Steam Flow Rate (kg/s).

505.6 Main Feedwater Flow Rate (kg/s).

505.6 Main Feedwater Temperature (K).

490 Auxiliary Feedwater Temperature (K).

300 6-26.3 Transient Table 4 lists the chronological sequence of events during the transient performed using the PWR TRACE5 model. In addition, in this section several graphs are shown with the main system variables obtained with the standard PWR TRACE5 model.

Table 4 Chronological Sequence of Events Event 3-loop TRACE 5 m odelTime (s)Break valve open.

0 Scram signal.

5 Turbine trip and closure of steam generator MSIVs.

5 Initiation of a decrease in liquid level in steam generator U

-tube. 5 Initiation of coastdown of primary coolant pumps.

5 Termination of steam generator Main Feedwater.

5 Open of steam generator Relief Valves.

5 Initiation of core power decay.

5 Initiation of HPI system in intact loops.

25 Primary pressure became lower than steam generator secondary pressure.

About 40 Initiation of AIS in intact loops.

About 100 Maximum fuel rod surface temperature.

130 End of the test.

300 6.4 System Pressures In Figure 4, the simulated primary and secondary pressures are shown. As it can be seen, the primary pressure starts to decrease at time zero due to the loss of coolant through the break. When the primary pressure decreases below 15.92 MPa, the scram signal is produced. The MSIVs are closed, producing the increasing of the secondary pressure. Due to this increas e, the steam generator RVs are opened to maintain the pressure almost constant. When the secondary pressure is lower than 6.8 MPa, the RVs are closed, and they are maintained closed till the end of the transient.

6-3 Figure 4 Primary and Secondary Pressures The primary pressure soon becomes lower than the steam generator secondary pressure (at about 50 s). From this moment on, the steam generators no longer serve as a heat sink. The Accumulator Injection system is initiated at about 100 s, when the primary pressure is lower than 4.6 MPa. At this moment, the depressurization becomes effective due to steam condensation caused by the coolant injection in the cold legs. When the primary pressure is about 1 MPa, the Low Pressure Injection system is activated.

6.5 Break Mass Flow Rate Figure 5 shows the mass flow rate through the break. As it can be seen, the change from liquid single-phase to two

-phase flow is produced in a very short time after the break. The two

-phase fluid regime is maintained until the hot legs are emptied. At this moment, the fluid regime changes to one

-phase vapor. The change from two

-phase to one phase vapor is produced at around 60 s. In order to adjust the break mass flow rate with TRACE5, the discharge coefficients for single

-phase liquid and two

-phase [2, 3] have been fixed to 1.0 (default value). These discharge coefficients and the application of the power

-to-volume strategy to the break size provide similar results than in the experiment (Figure 25).

0 50 100 150 200 250 300 Time (s)0 2 4 6 8 10 12 14 16 Pressure (Pa) 10 6 Typical PWR TRACE5 Primary Pressure Typical PWR TRACE5 Secondary PressureHPI activation in loop 2 and 3 AIS activation in loop 2 and 3 SI activation Closure of steam generators RVs

6-4 Figure 5 Break Mass Flow Rate 6.6 Primary Loop Mass Flow Rates Figure 6 shows the primary mass flow rates obtained with the model. As it can be seen, the mass flow rates start to decrease simultaneously with the pumps coastdown. After the pumps coastdown only natural circulation occurs. There are no significant differences among loops, even considering the loop with the break.

0 50 100 150 200 250 300 Time (s)0 500 1000 1500 2000 2500 3000 3500 4000 4500 Mass flow rate (kg/s)

Standard PWR TRACE5ACC activation in loop 2 an d 3 6-5 Figure 6 Primary Loop Mass Flow Rates 6.7 Pressure Vessel Collapsed Liquid Levels The simulated pressure vessel collapsed liquid levels are shown in Figure 7. As it can be seen, the liquid level in the core sharply decreases after the break due to flashing of fluid because of a fast-primary depressurization, while in the downcomer, the minimum liquid level is achieved at 120 s. From this moment on, the refill of the active core is produced. The entrance of cold water injected in two of the cold legs produces the vapor condensation in the active core and its collapsed liquid level rapidly increases.

0 50 100 150 200 250 300 Time (s)0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 Mass flow rate (kg/s)

Standard PWR TRACE5 Primary Loop 1 Standard PWR TRACE5 Primary Loop 2 Standard PWR TRACE5 Primary loop 3 6-6Figure 7 Cor e and Downcomer Collapsed Liquid Levels 6.8 Temperatures Figure 8 shows t he Maximum Peak Cladding (PCT) and t he Core Exit Temperature (CET). A s i t can be seen, bot h t emperatures s tart to i ncrease when the core c learance t akes pl ac e, at 40 s approximately.

T he maximum v alue of the PCT is r eached at 130 s and t he value is about 940 K.The CET excursion shows a similar behav ior t han the P CT, reaching i ts maximum v alue at asimilar time.0 50 100 150 200 250 300 Time (s)0 1 2 3 4 5 6 7 8 9 Liquid level (m)

Core Collapsed Liquid Level Downcomer Collapsed Liquid Level 6-7 Figure 8 Maximum Fuel Rod Surface and Core Exit Temperatures Figure 9 shows the cladding temperature at eight axial positions of a given fuel rod. In this case, a HTSTR component belonging to the second ring has been chosen due to the radial power distribution is maximum in it. The axial positions correspond to the eight axial levels in which is divided the active core. Furthermore, fuel rod temperatures for different HTSTRs (two of each ring) measured in the fifth axial level of the pressure vessel are shown in Figure 10. As it can be seen, the maximum fuel rod temperatures correspond to the HTSTR 51 and 83, which are in the second ring.

0 50 100 150 200 250 300 Time (s)300 400 500 600 700 800 900 1000 Temperature (K)

Standard PWR TRACE5 PCT Standard PWR TRACE5 CETACC activation in loop 2 and 3

6-8 Figure 9 Fuel Rod Surface Temperature at Different Axial Positions Figure 10 Fuel Rod Surface Temperature for Different HTSTR 0 50 100 150 200 250 300 Time (s)400 500 600 700 800 900 1000 Temperature (K)

HTSTR 51 pos 1 HTSTR 51 pos 2 HTSTR 51 pos 3 HTSTR 51 pos 4 HTSTR 51 pos 5 HTSTR 51 pos 6 HTSTR 51 pos 7 HTSTR 51 pos 8 0 50 100 150 200 250 300 Time (s)450 500 550 600 650 700 750 800 850 900 Temperature (K)

HTSTR 41 HTSTR 73 HTSTR 51 HTSTR 83 HTSTR 61 HTSTR 93 6-96.9 Hot and Cold Legs Liquid Levels Figure 11 and 12 show the collapsed liquid levels in the hot and cold legs, respectively. As it can be seen, liquid levels rapidly decrease due to the flashing produced by the fast primary depressurization just after the break. Regarding the cold leg liquid levels, it is observed that the three cold legs are emptied at around 70 s and at 120 s. When the activation of the AIS is produced, the water injected refills the three cold legs.

Figure 11 Hot Legs Collapsed Liquid Level 0 50 100 150 200 250 300 Time (s)0 0.5 1 1.5 2 2.5 Liquid level (m)

Standard PWR TRACE5 Loop 1 Standard PWR TRACE5 Loop 2 Standard PWR TRACE5 Loop 3Broken loop

6-10Figure 12 Col d Legs Collapsed Liquid Level 6.10 U-Tu bes Collapsed Liquid Level The c ollapsed liquid levels o f the U-t ubes a re shown in Figure 13.

A s i t c an be seen, the U-t ubes a re emptied v ery ear ly i n t he transient.

A t 7 0 s app roximately, t he U-t ubes of l oop 1 are completely em pty , av oiding the heat t ransfer bet ween the pr imary and secondary s ides. However, in loops 2 and 3 the U-tubes emptying is delayed until 150 s, approximately.

0 50 100 150 200 250 300 Time (s)0 0.5 1 1.5 2 2.5 Liquid level (m)

Standard PWR TRACE5 Loop 1 Standard PWR TRACE5 Loop 2 Standard PWR TRACE5 Loop 3 6-11 Figure 13 Steam G enerator U-T ube Up-F low Side Collapsed Liquid L evels 6.11 Study o f the E CCS Ac tuation I n the LSTF experiment, t he ECCS (HPI, AIS and LPI) are actuated only i n the intact l oop. However, in the 3-loop PWR TRACE5 model, i t i s modeled actuating in both intact l oops and results obtained ar e similar t o the experiment.

I n this section, t he effect o f the ECCS actuation in two or one intact l oops i s studied. Figure 14 shows t he pressures obtained i n both cases (ECCS actuation i n one or two intact loop s). A s i t ca n be seen, there ar e no important differences between them.

Cor e and downcomer collapsed liquid levels are shown in Figures 15 and 16, respectively.

A s i t can be seen, when the ECCS i s actuated i n two loops t he core and t he downcomer ar e refilled before. Figure 17 shows t he maximum PCT and the CET. A s i t i s expected, t he maximum values o f bot h temperatures ar e reached when the ECCS i s actuated only in one loop.

Consequently, the drop of these temperatures i s delayed in comparison to t he ca se with t he ECCS activated in two loops. Figures 18 and 19 sh ow the accumulator injection flow rate and t he Hig h Pressure Injecti on flow rate, respectively in the PWR model. 0 50 100 150 200 250 300 Time (s)0 2 4

6 8 10 12 Liquid level (m)

Standard PWR TRACE5 SG 1 Standard PWR TRACE5 SG 2 Standard PWR TRACE5 SG 3 6-12 Figure 14 Primary and Secondary Pressures Figure 15 Core Collapsed Liquid Levels 0 50 100 150 200 250 300 Time (s)0 2 4 6 8 10 12 14 16 Pressure (Pa) 10 6 Primary Pressure 2 ACC + 2 HPI Secondary Pressure 2 ACC + 2 HPI Primary Pressure 1 ACC + 1 HPI Secondary Pressure 1 ACC + 1 HPI 0 50 100 150 200 250 300 Time (s)0 0.5 1 1.5 2 2.5 3 3.5 Liquid level (m) 2 ACC + 2 HPI 1 ACC + 1 HPI 6-13 Figure 16 Downcomer Collapsed Liquid Levels Figure 17 Maximum Fuel Rod Surface and Core Exit Temperatures 0 50 100 150 200 250 300 Time (s)4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 Liquid level (m) 2 ACC + 2 HPI 1 ACC + 1 HPI 0 50 100 150 200 250 300 Time (s)400 500 600 700 800 900 1000 1100 1200 Temperature (K)

CET 2 ACC + 2 HPI PCT 2 ACC + 2 HPI CET 1 ACC + 1 HPI PCT 1 ACC + 1 HPI 6-14 Figure 18 Accumulator Injection System Mass Flow Rate Figure 19 High Pressure Injection System Mass Flow Rate 0 50 100 150 200 250 300 Time (s)0 50 100 150 200 250 300 350 400 450 Mass flow rate (kg/s)

Standard PWR TRACE5 Loop 1 Standard PWR TRACE5 Loop 2 Standard PWR TRACE5 Loop 3 0 50 100 150 200 250 300 Time (s)0 5 10 15 Mass flow rate (kg/s)

Standard PWR TRACE5 Loop 1 Standard PWR TRACE5 Loop 2 Standard PWR TRACE5 Loop 3 6-156.12 Void Fraction Figures 20, 21, 22 and 23 show the void fraction obtained at different moments during this transient. Firstly, Figure 22 shows the void fraction at the initiation of the test. As it can be seen, the primary system is full of liquid at this time. Then, the situation at 25 s is shown in Figure 21, when the pressurizer is empty. At this time, the liquid is located in the loop seals and in the bottom part of the pressure vessel. Figure 22 shows the situation when the accumulators of loops 2 and 3 are empty. As it can be seen, the loop seals are emptied, while the liquid is located in the core and in the pressure vessel bottom. Figure 23 shows the void fraction at the end of the transient.

6-16 Figure 20 Void Fraction i n 3-Loop PWR Plant at 0 s 6-17 Figure 21 Void Fraction in 3-Loop PWR Plant at 25 s

6-18 Figure 22 Void Fraction in 3

-Loop PWR Plant at 200 s

6-19 Figure 23 Void Fraction in 3

-Loop PWR Plant at the End of the Transient 6-206.13 LSTF Model vs Standard PWR Model In this section, a qualitative comparison between the experimental values (LSTF) and the simulation of both LSTF and Standard 3

-loop PWR is shown.

There is a good agreement in timing between the experiment (LSTF) and the simulations of LSTF and the standard PWR. Nevertheless, some differences are found. Primary pressure (Figure 24) in the 3

-loop model drops slightly earlier than in LSTF (experiment and simulation). This different behavior can be attributed to the different pump coast down used in the 3

-loop and LSTF models. Pressure discrepancies can be also attributed to the different modelization of the break valve used in each case. In the LSTF model, the break has been simulated following the experiment conditions, while in the standard PWR model the break has been simulated as an orifice.

The accumulator system is initiated at about 120 s in LSTF (experiment and simulation) and about 100 s in the 3

-loop PWR model, when the primary pressure is lower than 4.51 MPa. At this moment depressurization becomes effective due to steam condensation caused by the coolant injection in the cold legs. Finally, when primary pressure is about 1 MPa, Low Injection System is activated.

Figure 24 Primary and Secondary Pressures Figure 25 shows the comparison between the mass flow rates through the break obtained in the LSTF and in the standard PWR divided by the volumetric factor scaling. As it can be seen, the break mass flow rate obtained in the standard PWR divided by the volumetric scaling factor shows a good agreement with the mass flow rate obtained in the LSTF.

0 50 100 150 200 250 300 0 2 4 6 8 10 12 14 16x 10 6Time (s)Pressure (Pa)

LSTF Experiment Primary PressureLSTF Experiment Secondary PressureLSTF TRACE5 Primary PressureLSTF TRACE5 Secondary PressureStandard PWR TRACE5 Primary PressureStandard PWR TRACE5 Secondary PressureSI activation HPI activation Closure of SG RVs ACC activation 6-21 Figure 25 Break Mass Flow Rate Obtained In LSTF a nd i n the 3-Loop PWR Plant Scaled Figures 26 and 27 show the core and downcomer collapse liquid levels of LSTF (experiment and model) and the standard PWR model, respectively. As it can be seen, the general trend is similar in three cases, although in the standard PWR plant model the liquid level is recovered slightly rather than the LSTF (experiment and model). These discrepancies can be due to the mass flow rate injected by the AIS system in the standard PWR model.

0 50 100 150 200 250 300 0 20 40 60 80 100 120Time (s)Mass flow rate (kg/s)

LSTF TRACE5Standard PWR TRACE5 6-22 Figure 26 Core Collapsed Liquid Levels Figure 27 Downcomer Collapsed Liquid Levels 0 50 100 150 200 250 300 00.5 11.5 22.5 33.5 4Time (s)Liquid level (m)

LSTF ExperimentLSTF TRACE5Standard PWR TRACE5 0 50 100 150 200 250 300 3 4 5 6 7 8 9 10 11Time (s)Liquid level (m)

LSTF ExperimentLSTF TRACE5Standard PWR TRACE5 6-23 Regarding the temperatures of the system, in Figure 28 the maximum peak cladding (PCT) and the core exit temperatures (CET) are shown for the three cases considered (experiment, LSTF and standard PWR simulations). Both temperatures start to increase when the core clearance takes place, at 40 s approximately, in the experiment and in the standard PWR TRACE5 model.

In the LSTF TRACE5 model, this increase is delayed some seconds. The maximum value of the PCT in the standard PWR simulation is reached at 130 s and the value is about 940 K. CET shows a similar behavior than the PCT, reaching its maximum value at a similar time. As it can be seen, standard PWR temperatures obtained using TRACE5 are very similar to the experimental data. Figure 28 Maximum Fuel Rod Surface and Core Exit Temperatures Figures 29 and 30 show the comparison between hot and cold legs collapsed liquid levels achieved in the experiment and in both simulations (LSTF and 3

-loop PWR). As it can be seen, the refill in hot legs of the standard PWR is very similar to the LSTF. Regarding the cold leg, differences are more important, due to the discrepancies in the mass flow rate of the accumulator injection system.

0 50 100 150 200 250 300 300 400 500 600 700 800 900 1000Time (s)Temperature (K)

LSTF Experiment PCTLSTF TRACE5 PCTLSTF TRACE5 CETStandard PWR TRACE5 PCTStandard PWR TRACE5 CET 6-24 Figure 29 Hot Legs Collapsed Liquid Levels Fi gure 30 Cold Legs Collapsed Liquid Levels 0 50 100 150 200 250 300 350 400 450 500 00.20.40.60.8 11.21.4Time (s)Normalized Values LSTF Experiment ALSTF Experiment BLSTF TRACE5 ALSTF TRACE5 BStandard PWR TRACE5 1Standard PWR TRACE5 2Standard PWR TRACE5 3 0 50 100 150 200 250 300 350 400 450 500 00.10.20.30.40.50.60.70.80.9 1Time (s)Normalized Values LSTF Experiment ALSTF Experiment BLSTF TRACE5 ALSTF TRACE5 BStandard PWR TRACE5 1Standard PWR TRACE5 2Standard PWR TRACE5 3 7-1 7 CONCLUSIONS This w ork p resents results c orresponding t o t he simulation of an Intermediate Break Los s-O f-C oolant A ccident (IBLOCA) i n the c old leg a pplied to a standard 3-l oop P WR Nuclear P ower Plant (NPP). The c apability of t he thermal-h ydraulic c ode TRACE5 patch 2 i n the s imulation o f a Cold Leg IBLOCA i n a NPP has been t es ted. The ex periment has been ad apt ed t o this c ase us ing t he power-t o-v olum e scaling m ethodology. In general, TRACE5 can reproduce the significant events during an IBLOCA using the 3

-loop PWR model. The simulation results are consistent with the experimental data obtained in the IBLOCA transient reproduced in the LSTF.

Although differences app ear in design and scaling among the LSTF and the standard PWR TRACE5 model, i t has been test ed that using a volumetric scaling criterion the important physical phenomena are common t o t he LSTF and the 3-l oop PWR TRACE 5 model. T he volumetric s caling criterion has been us ed i n the break s ize and the ECCS i njection r ates (HPI, AI S and LPI). It has bee n observed that t he m ass flow r ate t hroug h the break obt ained i n the 3-l oop PWR di vided by t he volumetric s cali ng factor fits v ery w ell t he break mass f low rate obtained in the LSTF.

8-1 8 REFERENCES 1.Nuclear Regulatory Commission, Division of Risk Assessment and Special Projects.Office of Nuclear Regulatory Research. U. S Nuclear Regulatory Commission, TRACEV5.0. "Theory Manual. Field Equations, Solution Methods and Physical Models", NuclearRegulatory Commission, U.S.

(2007).2.Nuclear Regulatory Commission, Division of Risk Assessment and Special Projects,Office of Nuclear Regulatory Research, TRACE V5.0. User's Manual. Volume 1: InputSpecification, Nuclear Regulatory Commission, U. S (2007).

3.Thermohydraulic Safety Research Group, Nuclear Safety Research Center, Quick

-lookData Report of ROSA

-2/LSTF Test 2 (Cold Leg Intermediate Break LOCA IB

-CL-03 inJAEA), Japan Atomic Energy Agency, JAEA (2010).

4.The ROSA-V Group, ROSA

-V Large Scale Test Facility (LSTF) System Description forthe 3rd and 4th Simulated Fuel Assemblies, JAERI-Tech, Japan (2003).

5.Nuclear Regulatory Commission and Applied Programming Technology, SymbolicNuclear Analysis Package Version 2.1.2. (SNAP) (2012).

Simulation of ROSA

-2 Test-2 Experiment: Application to Nuclear Power Plant NUREG/IA-0511 S.Gallardo, A. Querol, M. Lorduy and G. Verdú Universitat P olitcnica de Valncia Instituto Universitario de Seguridad I ndustrial, R adiofís ica y M edioambiental Camí de Vera s/n 46022 Valencia, S PAINK.Tien, Project M anagerThe purpose of this work is to overview the results provided by a first approach to the simulation of an Intermediate Break L oss-Of-Coolant Accident (IBLOCA) in a 3-loop PWR Nuclear Power Plant (NPP) using the thermal-hydraulic code TRACE5 patch 2 and the Symbolic Nuclear Analysis Packages software (SNAP) version 2.1.2. The IBLOCA transient applied to the standard PWR TRACE5 model is the Test 2 (IB-C L-03) handled at the Large Scale Test Facility (LSTF) in the frame of the OECD/NEA ROSA

-2 Project. Test 2 simulates a 17% cold leg IBLOCA under the assumption of the single-f ailure of High Pressure In jection and Low Pressure In jec tion systems and total failure of the Auxiliary Feedwater. The LSTF is a Full Height Full Pressure (FHFP) facility designed to simulate a 4-loop W-type PWR (Tsuruga unit II NPP). The volumetric scaling factor is 1/48. The four primary loops of the reference PWR are scaled in LSTF by two equal-volume loops.

The core power used to simulate the decay power is 10 MW, corresponding to 14% of the 1/48 volumetrically scaled reference PWR rated power.

The simulation results are prov ided throughout several graphs, where the main system variables , such as press ures, pressure vessel liquid levels and temperatures are shown. These results represent a contribution to assess the predictability of computer codes such as TRACE5 Accident M anagement (AM) Committee on th e Safety of N uclear I nstallations (CSNI)

Consensus D E Segurid ad Nuclear (Nuclear Safety C ouncil,Spain, (C SN) Lar ge Scale Test F acility (LSTF) High P ressure Charging P ump (PJ)

High P ressure Injection Pump (PL)

Intermediate Break Loss-O f-C oolant-A ccident (IBLOCA) March 2019 Technical

NUREG/IA-05 11 International Agreement Report Simulation of ROSA-2 Test-2 Experiment

Application to Nuclear Power Plant Prepared by:

S. Gallardo, A. Querol, M. Lorduy and G. Verdú Universitat Politcnica de Valncia Instituto Universitario de Seguridad Industrial, Radiofísica y Medioambiental Camí de Vera s/n 46022 Valencia, SPAIN K. Tien, NRC Project Manager Division of Systems Analysis Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, DC 20555

-0001 Manuscript Completed: April 2018 Date Published

March 2019 Prepared as part of The Agreement on Research Participation and Technical Exchange Under the Thermal-Hydraulic Code Applications and Maintenance Program (CAMP)

Published by U.S. Nuclear Regulatory Commission

AVAILABILITY OF REFERENCE MATERIALSIN NRC PUBLICATIONS NRC Reference Material As of November 1999, you may electronically access NUREG-series publications and other NRC records at

NRC's Library at www.nrc.gov/reading-rm.html. Publicly released records include, to name a few, NUREG-series publications; Federal Register notices; applicant, licensee, and vendor documents and correspondence; NRC correspondence and internal memoranda; bulletins

and information notices; inspection and investigative

reports; licensee event reports; and Commission papers

and their attachments.

NRC publications in the NUREG series, NRC regulations, and Title 10, "Energy," in the Code of Federal Regulations may also be purchased from one of these two sources. The Superintendent of Documents Mail Stop IDCC

Washington, DC 20402-0001

Internet:

bookstore.gpo.gov

T elephone: (202) 512-1800

Fax: (202) 512-2104 The National Technical Information Service 5301 Shawnee R dAlexandria, VA 22312-0002 www.ntis.gov800-553-6847 or, locally, (703) 605-6000A single copy of each NRC draft report for comment isavailable free, to the extent of supply, upon written

request as follows:

Address: U.S. Nuclear Regulatory Commission Washington, DC 20555-0001 E-mail: distribution.resource@nrc.gov

Facsimile: (301) 415-2289 Some publications in the NUREG series that are posted at NRC's Web site address www.nrc.gov/reading-rm/

doc-collections/nuregs are updated periodically and may differ from the last printed version. Although references to material found on a Web site bear the date the material

was accessed, the material available on the date cited

may subsequently be removed from the site.

Non-NRC Reference Material Documents available from public and special technical libraries include all open literature items, such as books, journal articles, transactions, Federal Register notices, Federal and State legislation, and congressional reports.

Such documents as theses, dissertations, foreign reports

and translations, and non-NRC conference proceedings

may be purchased from their sponsoring organization.

Copies of industry codes and standards used in a substantive manner in the NRC regulatory process are

maintained at-The NRC Technical Library Two White Flint North 11545 Rockville Pike

Rockville, MD 20852-2738 These standards are available in the library for reference

use by the public. Codes and standards are usually

copyrighted and may be purchased from the originating organization or, if they are American National Standards, from-American National Standards Institute 11 West 42nd Street New York, NY 10036-8002 www.ansi.org (212)642-4900 Legally binding regulatory requirements are stated only in laws; NRC regulations; licenses, including technical

speci views expressed in contractor prepared publications in

this series are not necessarily those of the NRC.

The NUREG series comprises (1) technical and adminis-trative reports and books prepared by the staff (NUREG-

XXXX)or agency contractors (NUREG/CR-XXXX), (2) proceedings of conferences (NUREG/CP-XXXX), (3) reports

resulting from international agreements (NUREG/IA-XXXX),

(4)brochures (NUREG/BR-XXXX), and (5) compilations of legal decisions and orders of the Commission and Atomic and Safety Licensing Boards and of Directors' decisions under Section 2.206 of NRC's regulations (NUREG-0750).

DISCLAIMER:

This report was prepared under an international cooperative agreement for the exchange of

technical information. Neither the U.S. Government nor any

agency thereof, nor any employee, makes any warranty, expressed or implied, or assumes any legal liability or

responsibility for any third party's use, or the results of such

use, of any information, apparatus, product or process

disclosed in this publication, or represents that its use by

such third party would not infringe privately owned rights.

NUREG/IA-051 1 International Agreement Report Simulation of ROSA-2 Test-2 Experiment

Application to Nuclear Power Plant Prepared by:

S. Gallardo, A. Querol, M. Lorduy and G. Verdú Universitat Politcnica de Valncia Instituto Universitario de Seguridad Industrial, Radiofísica y Medioambiental Camí de Vera s/n 46022 Valencia, SPAIN K. Tien, NRC Project Manager Division of Systems Analysis Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, DC 20555

-0001 Manuscript Completed: April 2018 Date Published

March 2019 Prepared as part of The Agreement on Research Participation and Technical Exchange Under the Thermal

-Hydraulic Code Applications and Maintenance Program (CAMP)

Published by U.S. Nuclear Regulatory Commission

iii ABSTRACT The purpose of this work is to overview the results provided by a first approach to the simulation of an Intermediate Break Loss

-Of-Coolant Accident (IBLOCA) in a 3

-loop PWR Nuclear Power Plant (NPP) using the thermal

-hydraulic code TRACE5 patch 2 and the Symbolic Nuclear Analysis Packages software (SNAP) version 2.1.2.

The IBLOCA transient applied to the standard PWR TRACE5 model is the Test 2 (IB

-CL-03) handled at the Large Scale Test Facility (LSTF) in the frame of the OECD/NEA ROSA

-2 Project. Test 2 simulates a 17% cold leg IBLOCA under the assumption of the single

-failure of High Pressure Injection and Low Pressure Injection systems and total failure of the Auxiliary Feedwater.

The LSTF is a Full Height Full Pressure (FHFP) facility designed to simulate a 4-loop W-type PWR (Tsuruga unit II NPP). The volumetric scaling factor is 1/48. The four primary loops of the reference PWR are scaled in LSTF by two equal

-volume loops. The core power used to simulate the decay power is 10 MW, corresponding to 14% of the 1/48 volumetrically scaled reference PWR rated power.

The simulation results are provided throughout several graphs, where the main system variables, such as pressures, pressure vessel liquid levels and temperatures are shown. These results represent a contribution to assess the predictability of computer codes such as TRACE5.

v FOREWORD Thermalhydraulic studies play a key role in nuclear safety. Important areas where the significance and relevance of TH knowledge, data bases, methods and tools maintain an essential prominence are among others:

assessment of plant modifications (e.g., Technical Specifications, power uprates, etc.);

analysis of actual transients, incidents and/or start

-up tests; development and verification of Emergency Operating Procedures; providing some elements for the Probabilistic Safety Assessments (e.g., success criteria and available time for manual actions, and sequence delineation) and its applications within the risk informed regulation framework; training personnel (e.g., full scope and engineering simulators); and/or assessment of new designs.

For that reason, the history of the involvement in Thermalhydraulics of CSN, nuclear Spanish Industry as well as Spanish universities, is long. It dates back to mid 80's when the first serious talks about Spain participation in LOFT

-OCDE and ICAP Programs took place. Since then, CSN has paved a long way through several periods of CAMP programs, promoting coordinated joint efforts with Spanish organizations within different periods of associated national programs (i.e., CAMP-Espana). From the CSN perspective, we have largely achieved the objectives. Models of our plants are in place, and an infrastructure of national TH experts, models, complementary tools, as well as an ample set of applications, have been created. The main task now is to maintain the expertise, to consolidate it and to update the experience. We at the CSN are aware on the need of maintaining key infrastructures and expertise, and see CAMP program as a good and well consolidated example of international collaborative action implementing recommendations on this issue.

Many experimental facilities have contributed to the today's availability of a large thermal

-hydraulic database (both separated and integral effect tests). However, there is a continuous need for additional experimental work and code development and verification, in areas where no emphasis have been made along the past. On the basis of the SESAR/FAP 1 reports "Nuclear Safety Research in OECD Countries: Major Facilities and Programmes at Risk" (SESAR/FAP, 2001) and its 2007 updated version "Support Facilities for Existing and Advanced Reactors (SFEAR) NEA/CSNI/R(2007)6", CSNI is promoting since the beginning of this century several collaborative international actions in the area of experimental TH research. These reports presented some findings and recommendations to the CSNI, to sustain an adequate level of research, identifying a number of experimental facilities and programmes of potential interest for present or future international collaboration within the nuclear safety community during the coming decade. The different series of PKL, ROSA and ATLAS projects are under these premises. CSN, as Spanish representative in CSNI, is involved in some of these research activities, helping in this international support of facilities and in the establishment of a large network of international collaborations. In the TH framework, most of these actions are either covering not

1 SESAR/FAP is the Senior Group of Experts on Nuclear Safety Research Facilities and Programmes of NEA Committee on the Safety of Nuclear Installations (CSNI

vi enough investigated safety issues and phenomena (e.g., boron dilution, low power and shutdown conditions, beyond design accidents), or enlarging code validation and qualification data bases incorporating new information (e.g., multi

-dimensional aspects, non

-condensable gas effects, passive components).

This NUREG/IA report is part of the Spanish contribution to CAMP focused on:

Analysis, simulation an d investigation of specific safety aspects of PKL2/OECD and ROSA2/OECD experiments.

Analysis of applicability and/or extension of the results and knowledge acquired in these projects to the safety, operation or availability of the Spanish nuclear power plants. Both objectives are carried out by simulating the experiments and conducting the plant application with the last available versions of NRC TH codes (RELAP5 and/or TRACE).

On the whole, CSN is seeking to assure and to maintain the capability of the national groups with experience in the thermalhydraulics analysis of accidents in the Spanish nuclear power plants. Nuclear safety needs have not decreased as the nuclear share of the nations grid is expected to be maintained if not increased during next years, with new plants in some countries, but also with older plants of higher power in most of the countries. This is the challenge that will require new ideas and a continued effort.

________________________________

Rosario Velasco García, CSN Vice-president Nuclear Safety Council (CSN) of Spain

vii TABLE OF CONTENTS ABSTRACT ................................

................................................................

...............................

iii FOREWORD ................................................................................................

..............................

v LIST OF FIGURES

................................................................................................

....................

ix LIST OF TABLES

................................

................................................................

......................

xi EXECUTIVE

SUMMARY

................................

................................................................

........ .xiii ACKNOWLEDGMENTS

................................................................................................

...........

x v ABBREVIATIONS AND ACRONYMS

................................

....................................................

xv ii 1 INTRODUCTION

................................................................................................

...............

1-1 2 LSTF DESCRIPTION

................................................................................................

........ 2-1 3 LSTF TRANSIENT

................................

................................................................

............

3-1 4 THE STANDARD 3

-LOOP PWR MODEL

................................

......................................... 4-1 5 SCALING CONSIDERATIONS

................................

.........................................................

5-1 6 RESULTS AND DISCUSSION

................................

..........................................................

6-1 6.1 Transient Conditions ...............................................................................................

6-1 6.2 Steady-State ................................

................................................................

...........

6-1 6.3 Transient ................................

................................................................

.................

6-2 6.4 System Pressures ................................

................................................................

... 6-2 6.5 Break Mass Flow Rate

................................

............................................................

6-3 6.6 Primary Loop Mass Flow Rate s ..............................................................................

6-4 6.7 Pressure Vessel Collapsed Liquid Levels

................................................................

6-5 6.8 Temperatures

................................

................................................................

..........

6-6 6.9 Hot and Cold Legs Liquid Levels

................................

.............................................

6-9 6.10 U-Tubes Collapsed Liquid Level

................................

...........................................

6-10 6.11 Study of the ECCS Actuation ................................

................................................

6-11 6.12 Void Fraction................................

................................

................................

......... 6-15 6.13 LSTF Model vs Standard PWR M odel ................................

.................................. 6-20 7 CONCLUSIONS

................................

................................................................

................

7-1 8 REFERENCES

................................

................................................................

..................

8-1

ix LIST OF FIGURES Figure 1 Schematic View of the LSTF Facility ..................................................................... 2-1 Figure 2 LSTF Experimental Data

Primary a nd Secondary Pressures and PCT ................

3-2 Figure 3 Standard 3-Loop PWR Model Nodalization

................................

...........................

4-2 Figure 4 Primary a nd Secondary Pressures

. ................................

....................................... 6-3 Figure 5 Break Mass Flow Rate

................................

..........................................................

6-4 Figure 6 Primary Loop Mass Flow Rates

................................

............................................

6-5 Figure 7 Core a nd Downcomer Collapsed Liquid Levels

.....................................................

6-6 Figure 8 Maximum Fuel Rod Surface a nd Core Exit Temperatures

................................

.... 6-7 Figure 9 Fuel Rod Surface Temperature a t Different Axial Positions

................................

... 6-8 Figure 10 Fuel Rod Surface Temperature for Different HTSTR .............................................

6-8 Figure 11 Hot Legs Collapsed Liquid Level

...........................................................................

6-9 Figure 12 Cold Legs Collapsed Liquid Level

................................

....................................... 6-10 Figure 13 Steam Generator U

-Tube Up-Flow Side Collapsed Liquid Levels

.......................

6-11 Figure 14 Primary a nd Secondary Pressures

................................

...................................... 6-12 Figure 15 Core Collapsed Liquid Levels

................................

.............................................

6-12 Figure 16 Downcomer Collapsed Liquid Levels

................................

................................

.. 6-13 Figure 17 Maximum Fuel Rod Surface and Core Exit Temperatures

................................

.. 6-13 Figure 18 Accumulator Injection System Mass Flow Rate

................................

...................

6-14 Figure 19 High Pressure Injection System Mass Flow Rate

................................

................

6-14 Figure 20 Void Fraction i n 3-Loop PWR Plant at 0 s

................................

...........................

6-16 Figure 21 Void Fraction i n 3-Loop PWR Plant at 25 s ................................

.........................

6-17 Figure 22 Void Fraction i n 3-Loop PWR Plant at 200 s .......................................................

6-18 Figure 23 Void Fraction i n 3-Loop PWR Plant at the End o f the Transient

..........................

6-19 Figure 24 Primary and Secondary Pressures

...................................................................... 6-20 Figure 25 Break Mass Flow Rate Obtained i n LSTF and i n t he 3-Loop PWR Plant Scaled

. 6-21 Figure 26 Core Collapsed Liquid Levels

.............................................................................

6-22 Figure 27 Downcomer Collapsed Liquid Levels

................................

.................................. 6-22 Figure 28 Maximum Fuel Rod Surface and Core Exit Temperatures

................................

.. 6-23 Figure 29 Hot Legs Collapsed Liquid Levels

................................

....................................... 6-24 Figure 30 Cold Legs Collapsed Liquid Levels

..................................................................... 6-24

xi LIST OF TABLES Table 1 Control Logic and Sequence of Major Events in the LSTF Experiment

................. 3-1 Table 2 Main Characteristics of LSTF, 4-Loop NPP and 3-Loop NPP ................................ 5-3 Table 3 Steady- tate Conditions i n the TRACE5 M odel ....................................................6-1 Table 4 Chronological Sequence of Events ....................................................................... 6-2

xiii EXECUTIVE

SUMMARY

The purpose of this work is to test the capability of the thermal

-hydraulic code TRACE5 patch 2 in a first approach to the simulation of a cold leg Intermediate Break LOCA (IBLOCA) in a standard PWR plant of three loops. For this goal, a TRACE5 model of a standard 3

-loop PWR has been adapted to reproduce the Test 2 transient, in the frame of the OECD/NEA ROSA

-2 Project. Test 2 (IB-CL-03) was conducted in the Large Scale Test Facility (LSTF) of the Japan Atomic Energy Agency (JAEA). LSTF is a Full Height Full Pressure (FHFP) facility reproducing Tsuruga unit II Nuclear Power Plant (W-type 4-loop PWR). It allows to preserve time, power and mass inventory as in the LSTF, because the fluid exhibits the same properties at full pressure.

Test 2 simulates a 17% cold leg IBLOCA using an upward break nozzle mounted with the cold leg inner surface. Single

-failure of both High Pressure Injection and Low Pressure Injection systems and total failure of the Auxiliary Feedwater are assumed. Results of the simulation using TRACE5 are shown in several graphs, including primary and secondary pressures, discharged inventory, primary mass flow rates, and collapsed liquid levels (in the pressure vessel, hot legs, steam generators U

-tubes, etc.).

xv ACKNOWLEDGMENTS This paper contains findings that were produced within the OECD

-NEA ROSA-2 Project. The authors are grateful to the Management Board of the ROSA

-2 Project for their consent to this publication, and thank the Spanish Nuclear Regulatory Body (CSN) and the Asociación Espanola d e la Industria Eléctrica (UNESA) for the technical and financial support.

xvii ABBREVIATIONS AND ACRONYMS A FW Auxiliary Feedwater AIS Accumulator Injection System AM Accident Management BE Best Estimate CAMP Code Assessment and Management Program CET Core Exit Temperature CPU Central Processing Unit CRGT Control Rod Guide Tubes CSN Nuclear Safety Council, Spain DBE Design Basis Event ECCS Emergency Core Cooling System HPI High Pressure Injection IBLOCA Intermediate Break Loss

-Of-Coolant Accident JAEA Japan Atomic Energy Agency JAERI Japan Atomic Energy Research Institute JC Jet Condenser LOCA Loss-Of-Coolant Accident LPI Low Pressure Injection LSTF Large Scale Test Facility MFW Main Feedwater MSIV Main Steam Isolation Valve NPP Nuclear Power Plant NRC U.S. Nuclear Regulatory Commission NV Normalized to the Steady State Value PA Auxiliary Feedwater Pump PCT Peak Cladding Temperature PF Feedwater Pump PGIT Primary Gravity Injection Tank PJ High Pressure Charging Pump PL High Pressure Injection Pump PORV Pilot Operated Relief Valve PV Pressure Vessel PWR Pressurized Water Reactor PZR Pressurizer RHR Residual Heat Removal RV Relief Valve SBLOCA Small Break Loss

-Of-Coolant Accident SG Steam Generator SI Safety Injection SNAP Symbolic Nuclear Analysis Package SRV Safety Relief Valve ST Storage Tank

1-1 1 INTRODUCTION Since the accident a t Three Mile Island Unit

-2 (TMI-2) i n 1979, many test facilities hav e produced experimental d ata on S mall B reak Los s-O f-C oolant A ccidents (SBLOCA) for c ode assessment and dev elopment. H owever, the ex perimental d ata for Intermediat e Break Los s-O f-C oolant A cc idents (IBLOCA) i s quite l imited although t he t hermal-h ydraulic r esponses c ould differ bet ween SBLOCA and IBLOCA.

B ecaus e full-s cale testing is us ually i mpossibl e to perform, there i s a need t o test sca l ed models o f prototype systems.

With this aim, t he Japan Atomic Energy Agency (JAEA) i n the frame o f the OECD/NEA ROSA

-2 Project prepared detailed thermal-h y draulic dat a. The pur pose o f t his w ork i s t o t est t he capability of t he thermal-h y draulic c ode TRACE5 patch 2 [1, 2] i n the s imulation of a cold leg IBLOCA i n a standard 3-l oop PWR pl ant (representative of Spanish NPPs).

Fo r this goal, a generic 3-l oop PWR TRACE5 model has been adapted t o reproduce t he Test 2 , c onducted i n the Lar ge S cale Test Facility (LSTF) [4] i n t he frame o f t he OECD/NEA RO SA-2 Project.

Test 2 (IB-C L-0 3) [3] s imulates a PWR 17% c old l eg IBLOCA under t he as sumption o f a s ingle-f ailure of t he High Pres sure Injection (HPI) and th e Low P ress ure Injection (LPI) systems and the total f ailur e of t he Auxiliary Feedwater (AFW) system. LSTF is a Full Height Ful l Pressure (FHFP) facility designed t o simulat e Tsurug a unit I I NPP, which is a 4-l oop W-type PWR of 3423 MWt. The volumetric scaling factor i s 1/48. T he four primary l oops o f t he reference PWR a r e represented by two loops in LSTF. The initial core power used to simulate the decay core pow er i s 10 MW, corresponding t o 14% o f the 1/48 volumetrically scaled reference PWR rated power.

A generic 3-loop PWR TRACE5 model has been adapted applying the power

-to-volume scaling strategy t o simulate t he transient conducted i n the LSTF. Th e simulation results are shown in different graphs, including primary and secondary pressures, discharged inventory, primary mass flow rates, liquid levels (in the pressure vessel, hot leg, steam generators U-tubes, etc.) and temperatures.

2-12 LSTF DESCRIPTION The LSTF simulates a W-type PWR reactor, of four loops and 3423 MW of thermal power. It is characterized using prototypical

-scaled components with full

-height, 1/48 scaled volume and full-pressure conditions to its reference PWR plant (Tsuruga unit II NPP). The four primary loops of the reference PWR are represented by two loops with a volume factor of 1

/24. Figure 1 shows a schematic view of the LSTF facility.

Figur e 1 Schematic View of the LSTF Facility The main features of LSTF are the following:

Elevations: preserved.Volumes: scaled by 1/48 to the reference PWR.Flow area: scaled by 1/48 in the pressure vessel and by 1/24 in the steam generators.The flow area in the hot and cold legs is scaled to conserve the ratio of the length (L) toCore power: scaled by 1/48 and limited to 14% of the scaled core power of the reference PWR.Fuel assembly: The total number of fuel rods is scaled by 1/48. There are 1008 heated rods.The pr imary c oolant system o f LS TF c onsists o f the primary l oop A w ith the pressurizer (PZR) and the s ymmetrical p rimary l oop B. B oth l oops i nclude a primary c oolant pum p (PC) a nd 141 U-t ubes f or each steam g enerator. T he secondary-co olant s ys tem i ncludes a jet c ondenser (J C), a feed w ater pu m p (PF), t he aux iliary feedwater pumps (PA) and t he pi ping system.The Emergency Core Cooling System (ECCS) consists of the f ollowing sub-s ystems: the high-p ressure charging pump (PJ), the high-p ressure injection pump (PH), the accumulators (ACC), the low-p ressure injection pump (PL), the Residual Heat Removal (RHR) sys tem and the 2-2 primary gravity injection tank (PGIT). The primary coolant discharged from the primary system is stored in the break flow Storage Tank (ST).

The pressure vessel has five regions: the upper head located above the upper core support plate; the upper plenum situated between the upper core support plate and the upper core plate; the active core; the lower plenum and the downcomer annulus region that surrounds the core and the upper plenum. The LSTF vessel includes 8 spray nozzles (of 3.4 mm inner

-diameter) in the upper head, and 8 Control Rod Guide Tubes (CRGTs), which lead the flow path between the upper head and the upper plenum.

3-1 3 LSTF TRANSIENT Test 2 conducted in the LSTF simulates a 17% cold leg IBLOCA under the assumption of single-failure of the ECCS (HPI, AIS and LPI systems) and the total failure of the AFW. The complete control logic of the transient is listed in Table 1. The system pressures and the PCT evolution are shown in Figure 2. The break unit is connected to the cold leg and the orifice flow area corresponds to 17% of the volumetrically

-scaled cross

-sectional area of the reference PWR cold leg. The transient starts at time 0 with the break valve opening in the cold leg (in the loop without pressurizer) and increasing the rotational speed of the reactor coolant pumps.

F ew seconds later, the scram signal is generated. This signal produces the initiation of the core power decay curve, the initiation of the primary coolant pumps coastdown, the turbine trip, the closure of MSIV and the termination of the MFW.

When the primary pressure is lower than 12.27 MPa, the Safety Injection (SI) signal is generated. The HPI system is activated few seconds after the SI signal generation, only in the loop with pressurizer. The accumulators actuate when the primary pressure is lower than 4.5 MPa. The core power is automatically decreased by the LSTF Core Protection System when the maximum fuel rod surface temperature reaches 958 K. When the primary pressure is lower than 1.24 MPa, the LPI system actuates in the loop with pressurizer. Test 2 finishes with the break valve closure when primary and secondary pressures are stabilized.

Table 1 Control Logic and Sequence of Major Events i n the LSTF Experiment Break. Time zero Generation of scram signal.

Primary pressure = 12.97 MPa.

Initiation of core power decay curve simulation.

Generation of scram signal.

Initiation of Primary Coolant Pump coastdown. Generation of scram signal.

Turbine trip (closure of steam generators Main Steam Stop Valves).

Generation of scram signal.

Closure of Steam Generator Main Steam Isolation Valves (MSIVs).

Generation of scram signal.

Termination of Steam Generator Main Feedwater.

Generation of scram signal.

Generation of Safety Injection (SI) signal Primary pressure = 12.27 MPa. Initiation of High Pressure Injection (HPI) in loop w/pressurizer.

12 seconds after SI signal.

Initiation of the Accumulator Injection in loop w/pressurizer.

Primary pressure = 4.51 MPa.

Initiation of Low Pressure Injection (LPI) in loop w/pressurizer.

Primary pressure = 1.24 MPa.

3-2 Figure 2 LSTF Experimental Data: Primary and Secondary Pressures and PCT 050100150200250300350400450500 0 2 4 6 8 10 12 14 16 Pressure (Pa) 10 6 Primary Pressure Secondary Pressure050100150200250300350400450500 Time (s)400 500 600 700 800 900 1000 Temperature (K)

PCT 4-1 4 THE STANDARD 3

-LOOP PWR MODEL The pl ant o f reference used in this st udy i s a standard 3-loop PWR with a thermal power o f 2785 MWt. The TRACE5 model i ncludes t he pressure vessel, ho t and cold legs, pressurizer (in loop 3), reactor coolant pumps, loop seals, HPI, A I and LP I systems, stea m generators, etc. The pressure vessel has been modeled using a 3-D VESSEL component available in TRACE5.

The VESSEL component contains 16 axial levels, 4 radial rings and 4 azimuthal sectors. Three rings ar e use d t o simulate the upper head, upper plenum, active core and lower pl enum o f t he pressure vessel. The fourth ring i s use d t o simulate the downcomer annulus region. Low er plenum i s located between the axial levels 1 and 2. Axial levels 3 to 10 model the active core.

Upper pl enum i s simulated between the axial levels 11 and 14.

The levels 15 and 16 model the upper head. T he control rod guide tubes (CRGTs) ar e modeled using 6 PI PEs: 3 PIPEs allow t he flow pat h between upper pl enum and upper head (connect levels 13 to 16), while the other 3 connect t he core exit (level 11) with th e upper head (level 16). The pow er i s supplied to the active cor e using 12 Heat Structur e components (HTSTR), which simulat e t he fuel rods present i n the standard 3-loop PWR plant. The core includes 157 fuel assemblies with 17x17 lattice design.

A POWER component has been used to manage t he power from these HTSTRs.

I n this case, the power has been characterized by m eans o f a point kinetics model with constant initial reactivity and table-input reactivity as a function of time after a trip. The hot legs hav e been modeled with 3 PIPE components connect ed t o the pressure vessel and the U-tubes. I n l oop 3, the hot leg ha s a cross flow connection with the surge line of t he pressurizer.

A PIPE component i s u se d to simulate the pressurizer.

The pressurizer Relief Valves (RVs) and the Pilot Operated Relief Valve (PORV) are also considered. The cold legs hav e been simulated with 3 PIPEs link ed t o the pressure vessel, t he reactor coolant pumps and the ECCS. I n t he broken l oop (loop 1), the cold leg presents a cross flow connection with a VALVE component joined with a BREAK component t o simulat e the atmospheric conditions. In this valve, t he chok ed flow coefficients hav e be en fixed to default values (1.0 for subcooled and two-phase coefficients). Th e HPI and LP I systems have been simulated with FILL components. The AIS i s modeled using PIPE components with "accumulator" option. The discharge lines ar e modeled with check VALVE components connected to the col d legs. The reactor coolant pumps ar e simulated with PUMP components, considering specific head and torque homologous curves. The reactor coolant pumps a r e joi ned t o the loop seal.

The U-tubes o f each steam generator hav e been modeled using a PIPE component. The secondary side consists of a boiler, a downcomer, a steam separator and the steam lines. The boiler and the downcomer hav e been simulated with PIPE components, while a SEPD component has been used to m odel the steam separator.

T he relief and isolation valves have been modeled using VALVE components. The Relief Valve is linked t o a BREAK component, while the Main Steam Isolation Valve is connected to a PIPE component, which simulates the steam l ine header.

4-2standard 3-l oop PWR TRACE5 model us ing t he Symbolic N uclear A nalysis P ackag e (SNAP) [5] software version 2.1.2. Figure 3 Standard 3-Loop PWR Model Nodalization

5-1 5 SCALING CONSIDERATIO NS In LOCA scenarios, the most important consideration of the scaling criteria is to preserve both the power and the coolant mass inventory during the transient. The power

-to-volume scaling criterion is frequently used to preserve time, power and mass inventory in FHFP facilities regarding reference NPPs because they have the same fluid properties at full pressure. To perform a power

-to-volume analysis of the NPP model, scaling considerations should be made. In this frame, scaling factors between the LSTF and 3-loop NPP must be evaluated to assess the viability of the scaling analyses and can be applied to define the boundary conditions of the scaled model.

As LSTF is a FHFP facility of an actual NPP (Tsuruga unit II) and the scenario is an IBLOCA in the cold leg, the power

-to-volume scaling criterion has been chosen to develop a scale

-up TRACE5 model of LSTF. This criterion results from the application of conservation equations (1) to (4) under some requirements and implications:

Continuity equation:

+()=0 (1) Momentum equation:

+= 1 1()(2) Energy equation:

+=()+++ (3) State equation:

=( ,) (4) where is the density, u is the velocity, x is the coordinate, t is the time, p is the pressure, F is the friction coefficient, h is the enthalpy, thermal expansion coefficient and is the power density.

Substituting the next dimensionless parameters (denoted by an asterisk) in equations (1) to (4): = , = ,= , = , = , = ,= , = ,= gives a set of nondimensionalized equations (5) to (8):

+()=0 (5) 5-2+= 1 1()(6) +=()+++(7) =(,) (8) These equations contain other dimensionless parameters such as , and. The first two parameters are Froude and Euler numbers, respectively. The third is known as heat source number following Ishii and Kataoka terminology. Power-to-volume scaling method requires that all dimensionless parameters in equations (5) to (8)have to be equal in the LSTF model and the scale

-up LSTF model.====1 (9) Considering that s imilarity bet ween both s ystems has been achieved, t he following pow er-t o-v olume relations a re obtained: ==== (10) being power, Q mass flow rate, V volume, A area and the volumetric scaling factor.

Furthermore, in the scale

-up model, the Froude number is conserved in horizontal components. It implies varying the diameter and length of these components. Trying to conserve the Froude number, from the scale

-up mass flow rate calculated as Eq. (11), the scale

-up diameter, D, can be obtained as Eq. (12):

  • 4*=* 4**(11) = */ (12) where U is velocity in the scale

-up model, is the coolant density and d is LSTF diameter. Furthermore, from the volume equation Eq. (13) and trying to conserve the Froude number, lengths of the scale

-up piping system are obtained as Eq. (14):

  • 4*=* 4**(13) =*/ (14) being l and L, LSTF and scale-up length, respectively.

Considering the nominal power and volume of LSTF, Tsuruga unit II and the 3-loop NPP, the scaling ratios between them can be obtained. Table 2 summarizes the main characteristics of LSTF, 4-loop and 3

-loop NPP models.

5-3 I n terms of pressure and heights, the relation between LSTF and its reference NPP is 1

/1, so the initial value for pressures and temperatures in the 3

-loop PWR model is also the same than in the experiment performed in LSTF. The pressure of accumulators is the same in all cases.

The break size considered in Test 2 has been scaled applying the power

-to-volume criterion between LSTF, Tsuruga II and the 3-loop PWR NPP.

Table 2 Main Characteristics of LSTF, 4

-Loop NPP and 3-Loop NPP Parameter LSTF Reference 4

-loop NPP 3-loop NPP Kv (3-loop/LSTF)

Primary pressure (MPa).

15.5 15.5 15.5 1 Core power (MW).

10 3423 2785 39.0 (14% of 2785 MW). Number of loops. 2 4 3 Number of fuel assemblies

. 193 157 Fuel assembly array

. 17x 17 17x17 Core height (m)

. 3.66 3.66 3.66 1 Vessel volume (m 3). 2.75 131.7 106 38.54 Number of fuel rods.

1008 50952 41447 41.11 Volume of pressurizer (m

3) 1.2 51.0 39.7 33.08 Pressurizer heaters power (kW). 124 1800 1400 Average length U

-tubes (m).

20.2 20.2 20.2 1 Hot leg inner diameter (m).

0.207 0.737 0.737 8.11 8.11 8.11 1 Cold leg inner diameter (m).

0.207 0.698 0.698 Break area: 17% cold leg flow area reference NPP (m 2). 1.36 10-3 0.065 0.048 35.2 Accumulator pressure (MPa).

4.6 4.6 4.6 1

6-1 6 RESULTS AND DISCUSSION 6.1 Transient Conditions The following assumptions are imposed to reproduce Test 2 in the NPP model:

1.Break size is 17% of the cold leg flow area.

2.An upward long break nozzle is located on top of the cold leg in the broken loop withoutpressurizer (loop 1).

3.Loss of off

-site power concurrent with the scram.

4.HPI, AI and LPI systems are activated in intact loops.

5.Non-condensable gas inflow from accumulator tank may take place.

6.Total failure of AFW.

6.2 Steady- tate Table 3 s hows t he s teady-s tate conditions achiev ed using TRACE5.

Table 3 Steady-St ate C onditions i n t he T RACE5 Mo del Item 3-loop NPPTRACE5 model Core Power (MW t). 2785 Hot leg Fluid Temperature (K).

600 Cold leg Fluid Temperature (K).

569 Mass Flow Rate (kg/s).

4416 Pressurizer Pressure (MPa).

15.5 Pressurizer Liquid Level (m).

3.65 Accumulator System Pressure (MPa).

4.6 SG Secondary

-side Pressure (MPa).

6.8 SG Secondary

-side Liquid Level (m).

5.7 Steam Flow Rate (kg/s).

505.6 Main Feedwater Flow Rate (kg/s).

505.6 Main Feedwater Temperature (K).

490 Auxiliary Feedwater Temperature (K).

300 6-26.3 Transient Table 4 lists the chronological sequence of events during the transient performed using the PWR TRACE5 model. In addition, in this section several graphs are shown with the main system variables obtained with the standard PWR TRACE5 model.

Table 4 Chronological Sequence of Events Event 3-loop TRACE 5 m odelTime (s)Break valve open.

0 Scram signal.

5 Turbine trip and closure of steam generator MSIVs.

5 Initiation of a decrease in liquid level in steam generator U

-tube. 5 Initiation of coastdown of primary coolant pumps.

5 Termination of steam generator Main Feedwater.

5 Open of steam generator Relief Valves.

5 Initiation of core power decay.

5 Initiation of HPI system in intact loops.

25 Primary pressure became lower than steam generator secondary pressure.

About 40 Initiation of AIS in intact loops.

About 100 Maximum fuel rod surface temperature.

130 End of the test.

300 6.4 System Pressures In Figure 4, the simulated primary and secondary pressures are shown. As it can be seen, the primary pressure starts to decrease at time zero due to the loss of coolant through the break. When the primary pressure decreases below 15.92 MPa, the scram signal is produced. The MSIVs are closed, producing the increasing of the secondary pressure. Due to this increas e, the steam generator RVs are opened to maintain the pressure almost constant. When the secondary pressure is lower than 6.8 MPa, the RVs are closed, and they are maintained closed till the end of the transient.

6-3 Figure 4 Primary and Secondary Pressures The primary pressure soon becomes lower than the steam generator secondary pressure (at about 50 s). From this moment on, the steam generators no longer serve as a heat sink. The Accumulator Injection system is initiated at about 100 s, when the primary pressure is lower than 4.6 MPa. At this moment, the depressurization becomes effective due to steam condensation caused by the coolant injection in the cold legs. When the primary pressure is about 1 MPa, the Low Pressure Injection system is activated.

6.5 Break Mass Flow Rate Figure 5 shows the mass flow rate through the break. As it can be seen, the change from liquid single-phase to two

-phase flow is produced in a very short time after the break. The two

-phase fluid regime is maintained until the hot legs are emptied. At this moment, the fluid regime changes to one

-phase vapor. The change from two

-phase to one phase vapor is produced at around 60 s. In order to adjust the break mass flow rate with TRACE5, the discharge coefficients for single

-phase liquid and two

-phase [2, 3] have been fixed to 1.0 (default value). These discharge coefficients and the application of the power

-to-volume strategy to the break size provide similar results than in the experiment (Figure 25).

0 50 100 150 200 250 300 Time (s)0 2 4 6 8 10 12 14 16 Pressure (Pa) 10 6 Typical PWR TRACE5 Primary Pressure Typical PWR TRACE5 Secondary PressureHPI activation in loop 2 and 3 AIS activation in loop 2 and 3 SI activation Closure of steam generators RVs

6-4 Figure 5 Break Mass Flow Rate 6.6 Primary Loop Mass Flow Rates Figure 6 shows the primary mass flow rates obtained with the model. As it can be seen, the mass flow rates start to decrease simultaneously with the pumps coastdown. After the pumps coastdown only natural circulation occurs. There are no significant differences among loops, even considering the loop with the break.

0 50 100 150 200 250 300 Time (s)0 500 1000 1500 2000 2500 3000 3500 4000 4500 Mass flow rate (kg/s)

Standard PWR TRACE5ACC activation in loop 2 an d 3 6-5 Figure 6 Primary Loop Mass Flow Rates 6.7 Pressure Vessel Collapsed Liquid Levels The simulated pressure vessel collapsed liquid levels are shown in Figure 7. As it can be seen, the liquid level in the core sharply decreases after the break due to flashing of fluid because of a fast-primary depressurization, while in the downcomer, the minimum liquid level is achieved at 120 s. From this moment on, the refill of the active core is produced. The entrance of cold water injected in two of the cold legs produces the vapor condensation in the active core and its collapsed liquid level rapidly increases.

0 50 100 150 200 250 300 Time (s)0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 Mass flow rate (kg/s)

Standard PWR TRACE5 Primary Loop 1 Standard PWR TRACE5 Primary Loop 2 Standard PWR TRACE5 Primary loop 3 6-6Figure 7 Cor e and Downcomer Collapsed Liquid Levels 6.8 Temperatures Figure 8 shows t he Maximum Peak Cladding (PCT) and t he Core Exit Temperature (CET). A s i t can be seen, bot h t emperatures s tart to i ncrease when the core c learance t akes pl ac e, at 40 s approximately.

T he maximum v alue of the PCT is r eached at 130 s and t he value is about 940 K.The CET excursion shows a similar behav ior t han the P CT, reaching i ts maximum v alue at asimilar time.0 50 100 150 200 250 300 Time (s)0 1 2 3 4 5 6 7 8 9 Liquid level (m)

Core Collapsed Liquid Level Downcomer Collapsed Liquid Level 6-7 Figure 8 Maximum Fuel Rod Surface and Core Exit Temperatures Figure 9 shows the cladding temperature at eight axial positions of a given fuel rod. In this case, a HTSTR component belonging to the second ring has been chosen due to the radial power distribution is maximum in it. The axial positions correspond to the eight axial levels in which is divided the active core. Furthermore, fuel rod temperatures for different HTSTRs (two of each ring) measured in the fifth axial level of the pressure vessel are shown in Figure 10. As it can be seen, the maximum fuel rod temperatures correspond to the HTSTR 51 and 83, which are in the second ring.

0 50 100 150 200 250 300 Time (s)300 400 500 600 700 800 900 1000 Temperature (K)

Standard PWR TRACE5 PCT Standard PWR TRACE5 CETACC activation in loop 2 and 3

6-8 Figure 9 Fuel Rod Surface Temperature at Different Axial Positions Figure 10 Fuel Rod Surface Temperature for Different HTSTR 0 50 100 150 200 250 300 Time (s)400 500 600 700 800 900 1000 Temperature (K)

HTSTR 51 pos 1 HTSTR 51 pos 2 HTSTR 51 pos 3 HTSTR 51 pos 4 HTSTR 51 pos 5 HTSTR 51 pos 6 HTSTR 51 pos 7 HTSTR 51 pos 8 0 50 100 150 200 250 300 Time (s)450 500 550 600 650 700 750 800 850 900 Temperature (K)

HTSTR 41 HTSTR 73 HTSTR 51 HTSTR 83 HTSTR 61 HTSTR 93 6-96.9 Hot and Cold Legs Liquid Levels Figure 11 and 12 show the collapsed liquid levels in the hot and cold legs, respectively. As it can be seen, liquid levels rapidly decrease due to the flashing produced by the fast primary depressurization just after the break. Regarding the cold leg liquid levels, it is observed that the three cold legs are emptied at around 70 s and at 120 s. When the activation of the AIS is produced, the water injected refills the three cold legs.

Figure 11 Hot Legs Collapsed Liquid Level 0 50 100 150 200 250 300 Time (s)0 0.5 1 1.5 2 2.5 Liquid level (m)

Standard PWR TRACE5 Loop 1 Standard PWR TRACE5 Loop 2 Standard PWR TRACE5 Loop 3Broken loop

6-10Figure 12 Col d Legs Collapsed Liquid Level 6.10 U-Tu bes Collapsed Liquid Level The c ollapsed liquid levels o f the U-t ubes a re shown in Figure 13.

A s i t c an be seen, the U-t ubes a re emptied v ery ear ly i n t he transient.

A t 7 0 s app roximately, t he U-t ubes of l oop 1 are completely em pty , av oiding the heat t ransfer bet ween the pr imary and secondary s ides. However, in loops 2 and 3 the U-tubes emptying is delayed until 150 s, approximately.

0 50 100 150 200 250 300 Time (s)0 0.5 1 1.5 2 2.5 Liquid level (m)

Standard PWR TRACE5 Loop 1 Standard PWR TRACE5 Loop 2 Standard PWR TRACE5 Loop 3 6-11 Figure 13 Steam G enerator U-T ube Up-F low Side Collapsed Liquid L evels 6.11 Study o f the E CCS Ac tuation I n the LSTF experiment, t he ECCS (HPI, AIS and LPI) are actuated only i n the intact l oop. However, in the 3-loop PWR TRACE5 model, i t i s modeled actuating in both intact l oops and results obtained ar e similar t o the experiment.

I n this section, t he effect o f the ECCS actuation in two or one intact l oops i s studied. Figure 14 shows t he pressures obtained i n both cases (ECCS actuation i n one or two intact loop s). A s i t ca n be seen, there ar e no important differences between them.

Cor e and downcomer collapsed liquid levels are shown in Figures 15 and 16, respectively.

A s i t can be seen, when the ECCS i s actuated i n two loops t he core and t he downcomer ar e refilled before. Figure 17 shows t he maximum PCT and the CET. A s i t i s expected, t he maximum values o f bot h temperatures ar e reached when the ECCS i s actuated only in one loop.

Consequently, the drop of these temperatures i s delayed in comparison to t he ca se with t he ECCS activated in two loops. Figures 18 and 19 sh ow the accumulator injection flow rate and t he Hig h Pressure Injecti on flow rate, respectively in the PWR model. 0 50 100 150 200 250 300 Time (s)0 2 4

6 8 10 12 Liquid level (m)

Standard PWR TRACE5 SG 1 Standard PWR TRACE5 SG 2 Standard PWR TRACE5 SG 3 6-12 Figure 14 Primary and Secondary Pressures Figure 15 Core Collapsed Liquid Levels 0 50 100 150 200 250 300 Time (s)0 2 4 6 8 10 12 14 16 Pressure (Pa) 10 6 Primary Pressure 2 ACC + 2 HPI Secondary Pressure 2 ACC + 2 HPI Primary Pressure 1 ACC + 1 HPI Secondary Pressure 1 ACC + 1 HPI 0 50 100 150 200 250 300 Time (s)0 0.5 1 1.5 2 2.5 3 3.5 Liquid level (m) 2 ACC + 2 HPI 1 ACC + 1 HPI 6-13 Figure 16 Downcomer Collapsed Liquid Levels Figure 17 Maximum Fuel Rod Surface and Core Exit Temperatures 0 50 100 150 200 250 300 Time (s)4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 Liquid level (m) 2 ACC + 2 HPI 1 ACC + 1 HPI 0 50 100 150 200 250 300 Time (s)400 500 600 700 800 900 1000 1100 1200 Temperature (K)

CET 2 ACC + 2 HPI PCT 2 ACC + 2 HPI CET 1 ACC + 1 HPI PCT 1 ACC + 1 HPI 6-14 Figure 18 Accumulator Injection System Mass Flow Rate Figure 19 High Pressure Injection System Mass Flow Rate 0 50 100 150 200 250 300 Time (s)0 50 100 150 200 250 300 350 400 450 Mass flow rate (kg/s)

Standard PWR TRACE5 Loop 1 Standard PWR TRACE5 Loop 2 Standard PWR TRACE5 Loop 3 0 50 100 150 200 250 300 Time (s)0 5 10 15 Mass flow rate (kg/s)

Standard PWR TRACE5 Loop 1 Standard PWR TRACE5 Loop 2 Standard PWR TRACE5 Loop 3 6-156.12 Void Fraction Figures 20, 21, 22 and 23 show the void fraction obtained at different moments during this transient. Firstly, Figure 22 shows the void fraction at the initiation of the test. As it can be seen, the primary system is full of liquid at this time. Then, the situation at 25 s is shown in Figure 21, when the pressurizer is empty. At this time, the liquid is located in the loop seals and in the bottom part of the pressure vessel. Figure 22 shows the situation when the accumulators of loops 2 and 3 are empty. As it can be seen, the loop seals are emptied, while the liquid is located in the core and in the pressure vessel bottom. Figure 23 shows the void fraction at the end of the transient.

6-16 Figure 20 Void Fraction i n 3-Loop PWR Plant at 0 s 6-17 Figure 21 Void Fraction in 3-Loop PWR Plant at 25 s

6-18 Figure 22 Void Fraction in 3

-Loop PWR Plant at 200 s

6-19 Figure 23 Void Fraction in 3

-Loop PWR Plant at the End of the Transient 6-206.13 LSTF Model vs Standard PWR Model In this section, a qualitative comparison between the experimental values (LSTF) and the simulation of both LSTF and Standard 3

-loop PWR is shown.

There is a good agreement in timing between the experiment (LSTF) and the simulations of LSTF and the standard PWR. Nevertheless, some differences are found. Primary pressure (Figure 24) in the 3

-loop model drops slightly earlier than in LSTF (experiment and simulation). This different behavior can be attributed to the different pump coast down used in the 3

-loop and LSTF models. Pressure discrepancies can be also attributed to the different modelization of the break valve used in each case. In the LSTF model, the break has been simulated following the experiment conditions, while in the standard PWR model the break has been simulated as an orifice.

The accumulator system is initiated at about 120 s in LSTF (experiment and simulation) and about 100 s in the 3

-loop PWR model, when the primary pressure is lower than 4.51 MPa. At this moment depressurization becomes effective due to steam condensation caused by the coolant injection in the cold legs. Finally, when primary pressure is about 1 MPa, Low Injection System is activated.

Figure 24 Primary and Secondary Pressures Figure 25 shows the comparison between the mass flow rates through the break obtained in the LSTF and in the standard PWR divided by the volumetric factor scaling. As it can be seen, the break mass flow rate obtained in the standard PWR divided by the volumetric scaling factor shows a good agreement with the mass flow rate obtained in the LSTF.

0 50 100 150 200 250 300 0 2 4 6 8 10 12 14 16x 10 6Time (s)Pressure (Pa)

LSTF Experiment Primary PressureLSTF Experiment Secondary PressureLSTF TRACE5 Primary PressureLSTF TRACE5 Secondary PressureStandard PWR TRACE5 Primary PressureStandard PWR TRACE5 Secondary PressureSI activation HPI activation Closure of SG RVs ACC activation 6-21 Figure 25 Break Mass Flow Rate Obtained In LSTF a nd i n the 3-Loop PWR Plant Scaled Figures 26 and 27 show the core and downcomer collapse liquid levels of LSTF (experiment and model) and the standard PWR model, respectively. As it can be seen, the general trend is similar in three cases, although in the standard PWR plant model the liquid level is recovered slightly rather than the LSTF (experiment and model). These discrepancies can be due to the mass flow rate injected by the AIS system in the standard PWR model.

0 50 100 150 200 250 300 0 20 40 60 80 100 120Time (s)Mass flow rate (kg/s)

LSTF TRACE5Standard PWR TRACE5 6-22 Figure 26 Core Collapsed Liquid Levels Figure 27 Downcomer Collapsed Liquid Levels 0 50 100 150 200 250 300 00.5 11.5 22.5 33.5 4Time (s)Liquid level (m)

LSTF ExperimentLSTF TRACE5Standard PWR TRACE5 0 50 100 150 200 250 300 3 4 5 6 7 8 9 10 11Time (s)Liquid level (m)

LSTF ExperimentLSTF TRACE5Standard PWR TRACE5 6-23 Regarding the temperatures of the system, in Figure 28 the maximum peak cladding (PCT) and the core exit temperatures (CET) are shown for the three cases considered (experiment, LSTF and standard PWR simulations). Both temperatures start to increase when the core clearance takes place, at 40 s approximately, in the experiment and in the standard PWR TRACE5 model.

In the LSTF TRACE5 model, this increase is delayed some seconds. The maximum value of the PCT in the standard PWR simulation is reached at 130 s and the value is about 940 K. CET shows a similar behavior than the PCT, reaching its maximum value at a similar time. As it can be seen, standard PWR temperatures obtained using TRACE5 are very similar to the experimental data. Figure 28 Maximum Fuel Rod Surface and Core Exit Temperatures Figures 29 and 30 show the comparison between hot and cold legs collapsed liquid levels achieved in the experiment and in both simulations (LSTF and 3

-loop PWR). As it can be seen, the refill in hot legs of the standard PWR is very similar to the LSTF. Regarding the cold leg, differences are more important, due to the discrepancies in the mass flow rate of the accumulator injection system.

0 50 100 150 200 250 300 300 400 500 600 700 800 900 1000Time (s)Temperature (K)

LSTF Experiment PCTLSTF TRACE5 PCTLSTF TRACE5 CETStandard PWR TRACE5 PCTStandard PWR TRACE5 CET 6-24 Figure 29 Hot Legs Collapsed Liquid Levels Fi gure 30 Cold Legs Collapsed Liquid Levels 0 50 100 150 200 250 300 350 400 450 500 00.20.40.60.8 11.21.4Time (s)Normalized Values LSTF Experiment ALSTF Experiment BLSTF TRACE5 ALSTF TRACE5 BStandard PWR TRACE5 1Standard PWR TRACE5 2Standard PWR TRACE5 3 0 50 100 150 200 250 300 350 400 450 500 00.10.20.30.40.50.60.70.80.9 1Time (s)Normalized Values LSTF Experiment ALSTF Experiment BLSTF TRACE5 ALSTF TRACE5 BStandard PWR TRACE5 1Standard PWR TRACE5 2Standard PWR TRACE5 3 7-1 7 CONCLUSIONS This w ork p resents results c orresponding t o t he simulation of an Intermediate Break Los s-O f-C oolant A ccident (IBLOCA) i n the c old leg a pplied to a standard 3-l oop P WR Nuclear P ower Plant (NPP). The c apability of t he thermal-h ydraulic c ode TRACE5 patch 2 i n the s imulation o f a Cold Leg IBLOCA i n a NPP has been t es ted. The ex periment has been ad apt ed t o this c ase us ing t he power-t o-v olum e scaling m ethodology. In general, TRACE5 can reproduce the significant events during an IBLOCA using the 3

-loop PWR model. The simulation results are consistent with the experimental data obtained in the IBLOCA transient reproduced in the LSTF.

Although differences app ear in design and scaling among the LSTF and the standard PWR TRACE5 model, i t has been test ed that using a volumetric scaling criterion the important physical phenomena are common t o t he LSTF and the 3-l oop PWR TRACE 5 model. T he volumetric s caling criterion has been us ed i n the break s ize and the ECCS i njection r ates (HPI, AI S and LPI). It has bee n observed that t he m ass flow r ate t hroug h the break obt ained i n the 3-l oop PWR di vided by t he volumetric s cali ng factor fits v ery w ell t he break mass f low rate obtained in the LSTF.

8-1 8 REFERENCES 1.Nuclear Regulatory Commission, Division of Risk Assessment and Special Projects.Office of Nuclear Regulatory Research. U. S Nuclear Regulatory Commission, TRACEV5.0. "Theory Manual. Field Equations, Solution Methods and Physical Models", NuclearRegulatory Commission, U.S.

(2007).2.Nuclear Regulatory Commission, Division of Risk Assessment and Special Projects,Office of Nuclear Regulatory Research, TRACE V5.0. User's Manual. Volume 1: InputSpecification, Nuclear Regulatory Commission, U. S (2007).

3.Thermohydraulic Safety Research Group, Nuclear Safety Research Center, Quick

-lookData Report of ROSA

-2/LSTF Test 2 (Cold Leg Intermediate Break LOCA IB

-CL-03 inJAEA), Japan Atomic Energy Agency, JAEA (2010).

4.The ROSA-V Group, ROSA

-V Large Scale Test Facility (LSTF) System Description forthe 3rd and 4th Simulated Fuel Assemblies, JAERI-Tech, Japan (2003).

5.Nuclear Regulatory Commission and Applied Programming Technology, SymbolicNuclear Analysis Package Version 2.1.2. (SNAP) (2012).

Simulation of ROSA

-2 Test-2 Experiment: Application to Nuclear Power Plant NUREG/IA-0511 S.Gallardo, A. Querol, M. Lorduy and G. Verdú Universitat P olitcnica de Valncia Instituto Universitario de Seguridad I ndustrial, R adiofís ica y M edioambiental Camí de Vera s/n 46022 Valencia, S PAINK.Tien, Project M anagerThe purpose of this work is to overview the results provided by a first approach to the simulation of an Intermediate Break L oss-Of-Coolant Accident (IBLOCA) in a 3-loop PWR Nuclear Power Plant (NPP) using the thermal-hydraulic code TRACE5 patch 2 and the Symbolic Nuclear Analysis Packages software (SNAP) version 2.1.2. The IBLOCA transient applied to the standard PWR TRACE5 model is the Test 2 (IB-C L-03) handled at the Large Scale Test Facility (LSTF) in the frame of the OECD/NEA ROSA

-2 Project. Test 2 simulates a 17% cold leg IBLOCA under the assumption of the single-f ailure of High Pressure In jection and Low Pressure In jec tion systems and total failure of the Auxiliary Feedwater. The LSTF is a Full Height Full Pressure (FHFP) facility designed to simulate a 4-loop W-type PWR (Tsuruga unit II NPP). The volumetric scaling factor is 1/48. The four primary loops of the reference PWR are scaled in LSTF by two equal-volume loops.

The core power used to simulate the decay power is 10 MW, corresponding to 14% of the 1/48 volumetrically scaled reference PWR rated power.

The simulation results are prov ided throughout several graphs, where the main system variables , such as press ures, pressure vessel liquid levels and temperatures are shown. These results represent a contribution to assess the predictability of computer codes such as TRACE5 Accident M anagement (AM) Committee on th e Safety of N uclear I nstallations (CSNI)

Consensus D E Segurid ad Nuclear (Nuclear Safety C ouncil,Spain, (C SN) Lar ge Scale Test F acility (LSTF) High P ressure Charging P ump (PJ)

High P ressure Injection Pump (PL)

Intermediate Break Loss-O f-C oolant-A ccident (IBLOCA) March 2019 Technical