ML19056A075
ML19056A075 | |
Person / Time | |
---|---|
Issue date: | 02/28/2019 |
From: | Martinez-Quiroga V, Pretel C, Reventos F, Kirk Tien Office of Nuclear Regulatory Research, Technical University of Catalonia, Barcelona, Spain |
To: | |
Meyd, Donald | |
References | |
NUREG/IA-0499 | |
Download: ML19056A075 (41) | |
Text
NUREG/IA-0499 International Agreement Report Post-Test Calculation of the PKL-2 Test G7.1 Using RELAP5/MOD3.3 Prepared by:
V Martinez-Quiroga, F. Reventos, and C. Pretel Institute of Energy Technologies Technical University of Catalonia ETSEIB, Av. Diagonal 647, Pav. C 08028 Barcelona Spain Kirk Tien, NRC Project Manager Division of Systems Analysis Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, DC 20555-0001 Manuscript Completed: November 2017 Date Published: February 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 MATERIALS IN NRC PUBLICATIONS NRC Reference Material Non-NRC Reference Material As of November 1999, you may electronically access Documents available from public and special technical NUREG-series publications and other NRC records at libraries include all open literature items, such as books, NRCs Library at www.nrc.gov/reading-rm.html. Publicly journal articles, transactions, Federal Register notices, released records include, to name a few, NUREG-series Federal and State legislation, and congressional reports.
publications; Federal Register notices; applicant, Such documents as theses, dissertations, foreign reports licensee, and vendor documents and correspondence; and translations, and non-NRC conference proceedings NRC correspondence and internal memoranda; bulletins may be purchased from their sponsoring organization.
and information notices; inspection and investigative reports; licensee event reports; and Commission papers Copies of industry codes and standards used in a and their attachments. substantive manner in the NRC regulatory process are maintained at NRC publications in the NUREG series, NRC regulations, The NRC Technical Library and Title 10, Energy, in the Code of Federal Regulations Two White Flint North may also be purchased from one of these two sources. 11545 Rockville Pike Rockville, MD 20852-2738
- 1. The Superintendent of Documents U.S. Government Publishing Office These standards are available in the library for reference Mail Stop IDCC use by the public. Codes and standards are usually Washington, DC 20402-0001 copyrighted and may be purchased from the originating Internet: bookstore.gpo.gov organization or, if they are American National Standards, Telephone: (202) 512-1800 from Fax: (202) 512-2104 American National Standards Institute 11 West 42nd Street
- 2. The National Technical Information Service New York, NY 10036-8002 5301 Shawnee Road www.ansi.org Alexandria, VA 22312-0002 (212) 642-4900 www.ntis.gov 1-800-553-6847 or, locally, (703) 605-6000 Legally binding regulatory requirements are stated only in A single copy of each NRC draft report for comment is laws; NRC regulations; licenses, including technical available free, to the extent of supply, upon written specifications; or orders, not in NUREG-series publications.
request as follows: The views expressed in contractor prepared publications in this series are not necessarily those of the NRC.
Address: U.S. Nuclear Regulatory Commission The NUREG series comprises (1) technical and adminis-Office of Administration trative reports and books prepared by the staff (NUREG-XXXX) or agency contractors (NUREG/CR-XXXX), (2)
Multimedia, Graphics, and Storage &
proceedings of conferences (NUREG/CP-XXXX), (3) reports Distribution Branch resulting from international agreements (NUREG/IA-XXXX),
Washington, DC 20555-0001 (4) brochures (NUREG/BR-XXXX), and (5) compilations of E-mail: distribution.resource@nrc.gov legal decisions and orders of the Commission and Atomic Facsimile: (301) 415-2289 and Safety Licensing Boards and of Directors decisions under Section 2.206 of NRCs regulations (NUREG-0750).
Some publications in the NUREG series that are posted at NRCs Web site address www.nrc.gov/reading-rm/ DISCLAIMER: This report was prepared under an doc-collections/nuregs are updated periodically and may international cooperative agreement for the exchange of technical information. Neither the U.S. Government nor any differ from the last printed version. Although references to agency thereof, nor any employee, makes any warranty, material found on a Web site bear the date the material expressed or implied, or assumes any legal liability or was accessed, the material available on the date cited responsibility for any third partys use, or the results of such may subsequently be removed from the site. 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-0499 International Agreement Report Post-Test Calculation of the PKL-2 Test G7.1 Using RELAP5/MOD3.3 Prepared by:
V Martinez-Quiroga, F. Reventos, and C. Pretel Institute of Energy Technologies Technical University of Catalonia ETSEIB, Av. Diagonal 647, Pav. C 08028 Barcelona Spain Kirk Tien, NRC Project Manager Division of Systems Analysis Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, DC 20555-0001 Manuscript Completed: November 2017 Date Published: February 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
ABSTRACT The Thermalhydraulic Studies Group of Technical University of Catalonia (UPC) holds a large background in nuclear safety studies in the field of Nuclear Power Plant (NPP) code simulators.
This report analyzes with RELAP5mod3.3 the G7.1 PKL Test. This experiment is part of a Counterpart Test performed in LSTF and PKL Test Facilities within the framework of the OECD/NEA ROSA-2 and PKL-2 projects. Detailed core nodalizations and Pseudo 3D modeling have been object of study as well as the capabilities of the code for reproducing the correlation between the Core Exit Temperature (CET) and the Peak Cladding Temperature (PCT).
iii
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 80s 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 todays 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/FAP1 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 enough investigated safety issues and phenomena (e.g., boron dilution, low power and 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).
v
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 and investigation of specific safety aspects of PKL/OECD ROSA/OECD and ATLAS/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 vi
TABLE OF CONTENTS ABSTRACT ............................................................................................................................... iii FOREWORD .............................................................................................................................. v LIST OF FIGURESix LIST OF TABLES..ix EXECUTIVE
SUMMARY
.........xi ACKNOWLEDGMENTS .......................................................................................................... xiii ABBREVIATIONS AND ACRONYMS ..................................................................................... xv 1 INTRODUCTION ...............................................................................................................1 1.1 PKL-2 and ROSA-2 Counterpart Test .......................................................................1 1.2 The OECD PKL 2 Project ..........................................................................................2 2 FACILITY AND TEST DESCRIPTION ...............................................................................3 2.1 PKL III Test Facility ...................................................................................................3 2.2 Experimental Conditions ...........................................................................................4 2.3 Initial Conditions .......................................................................................................4 2.4 Test Phase ...............................................................................................................4 3 CODE INPUT MODEL DESCRIPTION ..............................................................................7 4 RESULTS ..........................................................................................................................9 5 RUN STATISTICS ........................................................................................................... 13 6 CONCLUSIONS .............................................................................................................. 15 7 REFERENCES ................................................................................................................ 17 vii
LIST OF FIGURES Figure 1 Diagram with the Different Counterpart Conditions (Courtesy of the OECD/NEA ROSA-2 Group) ....................................................................................................... 1 Figure 2 OECD PKL 2 Project Experiments (Courtesy of the OECD/NEA PKL-2 Group)....... 2 Figure 3 PKL Test Facility (Courtesy of the OECD/NEA PKL Group)..................................... 3 Figure 4 Nodalization of the Main Steam System................................................................... 7 Figure 5 Nodalization of the Bypass Heaters ......................................................................... 7 Figure 6 UPC PKL Pseudo 3D Nodalization Core Channels .................................................. 8 Figure 7 Primary Pressure ..................................................................................................... 9 Figure 8 Secondary Pressure ...............................................................................................10 Figure 9 Mass Flow Comparison ..........................................................................................10 Figure 10 Core Exit Temperature............................................................................................11 Figure 11 CET vs PCT Correlation .........................................................................................11 Figure 12 CET and PCT Versus Time.....................................................................................12 LIST OF TABLES Table 1 Chronology of the Main Events of Test G7.1 ............................................................ 5 Table 2 Initial Conditions of Test G7.1 .................................................................................. 9 Table 3 Run Statistics ..........................................................................................................13 ix
EXECUTIVE
SUMMARY
Experimental research activities are being performed in Germany by the OECD PKL 2 project with the aim of addressing thermal-hydraulic safety issues for current PWR and new PWR design concepts. These experiments are carried out at the PKL III test facility.
This report analyses the experiment G7.1 of the LSTF and PKL Counterpart Test, which was carried out as a part of the OECD/NEA ROSA-2 and PKL-2 projects. The aim of this international synergy was to analyze the effectiveness of Core Exit Temperature measurement in Accident Management strategies as well as the scaling effects that appear between counterpart transients performed at different sizes and designs.
UPC PKL Relap5mod33 nodalization, created for the simulation of the PKL boron dilution experiments, has been used and improved for this test. Two different nodalizations were prepared in order to check the capabilities of the code for simulating properly the correlation between the PCT and CET: one, the UPC PKL 1D nodalization, with just one core channel; and a second, the UPC Pseudo 3D nodalization, with concentric core channels and activating the transversal momentum equations.
Many other aspects related to the nodalization were adjusted and verified in order to improve results.
In general, the results of the simulation demonstrated that 1D nodalizations are good enough for describing the general behavior of the transient as well as the main events and phenomena. On the other hand, for simulating accurately the correlation between the CET and PCT, Pseudo 3D modeling is necessary.
xi
ACKNOWLEDGMENTS This paper contains findings that were produced within the OECD-NEA PKL-2 Project. The authors are grateful to the Management Board of the PKL-2 Project for their consent to this publication.
xiii
ABBREVIATIONS AND ACRONYMS ACC accumulator AM accident management CET core exit temperature ECCS emergency core cooling system EOP emergency operational procedure HPIS high pressure injection system HS heat structure LPIS low pressure injection system LSTF large scale test facility MS main steam MSIV main steam isolation valve MSLB main steam line break NEA Nuclear Energy Agency NPP nuclear power plant OECD Organization for Economic Cooperation and Development PCT peak cladding temperature PKL Primrkreislauf PZR pressurizer PWR pressurized water reactor RELAP reactor excursion and leak analysis program ROSA rig of safety assessment SBLOCA small break loss of coolant accident SG steam generator UP upper plenum UPC Universitat Politcnica de Catalunya (Technical University of Catalonia)
UT u-tubes xv
1 INTRODUCTION Several safety activities have been performed during the last decades under the auspices of the OECD to develop and improve computer codes. They include several experiments at integral test facilities like the Test G7.1, which forms part of the OECD PKL-2 project.
1.1 PKL-2 and ROSA-2 Counterpart Test In 2011 a Counterpart Test was performed in LSTF and PKL Test Facilities as a part of the OECD/NEA ROSA-2 and PKL-2 projects (Test 3 and test G7.1 respectively). The objective of both tests was devoted to analyze two aspects:
Core Exit Temperature measurement effectiveness in Accident Management of NPPs Scaling effects between PKL and LSTF Test Facilities In Accident Management strategies, core exit temperature measurement becomes critical for detecting core dryout and for avoiding that PCT rises until safety limits. In the Counterpart Test, the relationship between CET and PCT was object of study in order to analyze EOP set points.
The selected scenario was a hot leg SBLOCA. System failures as no high pressure safety injection and no automatic secondary-side safety cooldown were imposed. The particular test conditions for the PKL-2 Test G7.1 are described in section 2.2. The main phenomena of interest were:
Core boil-off with steam generation Steam flow toward hot-leg break Realistic pressure during core boil-off Relationship between PCT and CET Figure 1 Diagram with the Different Counterpart Conditions (Courtesy of the OECD/NEA ROSA-2 Group) 1
About the scaling issue, core dryout and AM actions were simulated under counterpart conditions at different sizes (PKL -1:145- and LSTF -1:48- ) and pressures (LSTF high and low pressure transient phases) in order to check how the scaling affects the CET and its related phenomena, and in order to sound out possible strategies for extrapolating the results to a commercial NPP (see figure one). The direct comparison of both transients was established between the whole transient of the PKL Test G7.1, and the low pressure transient phase of the LSTF Test 3.
1.2 The OECD PKL 2 Project The OECD PKL 2 is an international project carried out between 2008 and 2011 that includes several types of integral experiments (see figure two) on PKL III, ROCOM and PMK test facilities. The aim of the project is to address thermal-hydraulic safety issues for current PWR and new PWR design concepts. It also provides a wide database for the validation of computer codes and models for system integral analyses coupled with detailed analyses of local phenomena.
Begin of Project 1. PRG/MB 2. PRG/MB 3. PRG/MB 4. PRG/MB 5. PRG/MB 6. PRG/MB 7. PRG/MB 8. PRG/MB 2008 2009 2010 2011 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12 G1 (2 Tests)
PMK (2 Tests)
G2 (1 Test)
G3 (1 Test)
ROCOM Tests G4 (1 Test)
G5 (1 Test)
G6 (1 Test)
G7 (1 Test)
Figure 2 OECD PKL 2 Project Experiments (Courtesy of the OECD/NEA PKL-2 Group)
The topics of the experiments are:
G1 Systematic investigation of the heat transfer mechanisms in the SGs in presence of nitrogen, steam and water G2 Cooldown procedures with SGs isolated and emptied on the secondary side G3 Fast cooldown transients (main steam line break)
G4 Systematic study on heat transfer in the SG under reflux condenser conditions G5 Boron precipitation processes after LB - LOCA G6 RCS cooldown with void formation in RPV upper head G7 Counterpart Test with the OECD/NEA ROSA-2 Test 3 (hot leg SBLOCA) 2
2 FACILITY AND TEST DESCRIPTION 2.1 PKL III Test Facility PKL is an experimental power facility plant designed to simulate pressurized water reactors (PWR) under accidental conditions. The plant, situated in Erlangen (Germany), has been the scenario of several experiments in the last 25 years. The PKL facility replicates the entire primary system and most of the secondary system (except for the turbine and condenser) of a 1300-MW PWR plant, with elevations scaled 1:1 and diameters reduced by a factor of 12.
Therefore volumes and power have been reduced by a factor of 145. The number of rods in the core and the U-tubes in the steam generator has been divided by 145 too. The core has been modeled by 314 electrical heater rods. Unlike many experimental facilities with only two available loops (one for the break loop and one to simulate the other three loops) PKL simulates all four loops separately. This is very important in order to analyze asymmetrical transients, e.g.
with injection in two out of four loops. A diagram of the PKL test facility is shown in Figure 3.
Figure 3 PKL Test Facility (Courtesy of the OECD/NEA PKL Group)
The operating pressure of the PKL facility is limited to 45 bars on the primary side and to 56 bars on the secondary side. This allows simulation over a wide temperature range (322K to 522K) that is particularly applicable to the cooldown procedures investigated.
3
All emergency systems are represented and have a wide versatility referred to their functions and positions. There are 8 accumulators (one for each cold and hot leg). The pump injection system is available in all the hot and cold legs. Many break locations are available too.
PKL test facility has about 1500 measurement points that permit an exhaustive analysis of the tests. There are measurement devices for cladding, wall and fluid temperature, absolute and differential pressure, one and two phase mass flow, density and boron concentration.
Sixty of the measurement devices are identical to those that are used in a commercial plant to simulate what an operator would control in case of accident.
2.2 Experimental Conditions The hardware configuration of PKL is described in references [1] and [2]. Some important points are the following:
- Break: upward oriented SBLOCA (1.5 %) at Hot Leg without PZR.
- ECCs: HPIS full failure.
ACCs and LPIS set points fixed at 26 bars and 8 bars respectively.
ACC level and LPIS mass flow adjusted according to LSTF final conditions.
- Core power curve: constant at 1.8 % of the nominal power. Additional heat was added in order to compensate the differences between the environment losses of LSTF and PKL.
- SG depressurization: fully opening of 2 main steam relief valves when CET achieves 350
ºC. The depressurization of the 4 SGs was simulated connecting them via main steam header.
- Main steam relief valves: modified for having same ratio of LSTF valves.
2.3 Initial Conditions The PKL initial conditions were adjusted at its maximum pressures in order to reproduce as realistic as it can PWR SBLOCA reflux and condensation, core dry-out and Accident Management phases.
In relation to scaling, the mass inventory in the secondary side was adjusted using Kv factor in order to have the same ratio between liquid and energy storage in the SGs. In the primary system, mass inventory was tuned for having the same reflux condenser and break discharging initial conditions than the low pressure transient phase of the LSTF Test 3.
2.4 Test Phase The transient starts at t=0 seconds with the opening of the break valve. During the first 940 seconds, core is cooled under saturated conditions and reflux and condensation occurs in the UTs of the SGs. In this phase of the transient, primary pressure keeps constant over the 44 bars of the isolated secondary side. After that, core uncovery starts, and few seconds later (at 1020 seconds), primary pressure drops below secondary pressure. In this second phase, there is vapour superheating in the core, and the CET and the PCT rise above the temperature of saturation (with a delay of 270 seconds between both).
When the CET temperature achieves the AM temperature set point (1360 seconds), SGs are depressurized by MS relief valves opening. System pressures drop rapidly without a complete 4
core quenching, so CET and PCT temperatures do not drop to saturated conditions until accumulator injection system is activated (1500 seconds).
Finally, when primary pressure drops below 8 bars (2060 seconds), LPIS starts to inject water, compensating break losses and keeping constant plant parameters. At 5685 seconds, break valve is closed and the end of the transient is declared.
The main events are described in Table 1:
Table 1 Chronology of the Main Events of Test G7.1 Experimental UPC 1D UPC Pseudo-3D (s) nodalization (s) nodalization (s)
Start of the transient 0 0 0 Begin of core uncovery 940 800 940 Primary pressure below 1020 920 1010 secondary pressure Secondary side 1360 1190 1295 depressurization Start of accum. Injection 1500 1304 1450 ACC injection finished 1860 1712 1752 LPIS started 2060 1966 1993 End of the test 5685 5685 5685 5
3 CODE INPUT MODEL DESCRIPTION Two different nodalizations were prepared by making use of the UPC PKL RELAP5mod3.3 nodalization. This nodalization had been previously qualified with the simulation of the OECD/NEA PKL I tests F1.1 and F1.2, and the OECD/SETH tests E2.1, E2.2 and E2.3 (see references [3], [4] and [5] and [6]). Particular Counterpart nodalizations were:
UPC PKL 1D nodalization UPC PKL Pseudo-3D nodalization Several improvements were done for both nodalizations taking advantage of the particular specifications of the transient. In this way, main steam system was modified following the sketch shown in Figure 4. Moreover, realistic bypass heaters were added to the PZR and the SGs (as in Figure 5) in order to have a proper power balance between the environment heat losses and the power supplied by the exchangers.
Figure 4 Nodalization of the Main Steam System Figure 5 Nodalization of the Bypass Heaters 7
The differences between both nodalizations were exclusively related to the core and upper plenum (UP) modeling. UPC PKL 1D nodalization simulated them with one channel (in addition to the core bypass), having fuel and all passive heat structures (core barrel and unheated rods) linked to the same volumes. The fuel was modeled with three HSs, with the same power ratio and divided in 7 axial levels.
UPC PKL Pseudo-3D nodalization had the core and UP (until the CET thermocouple level) divided in three radial channels (see Figure 6), with one fuel HS for each channel. The HSs for the passive internal metal structures were split for each channel proportionally to the flow path of each one, and the core barrel was linked to the outer zone. The radial flow paths between cells were modeled and transversal momentum equations were activated following the recommendations of reference [7].The total number of core axial meshes was increased to 14 and the UP cell heights were adjusted so that the center of the node coincided with the elevation of the thermocouples in the test facility.
Figure 6 UPC PKL Pseudo 3D Nodalization Core Channels 8
4 RESULTS The results obtained for both nodalizations showed a quite close agreement reproducing the initial conditions (see Table 2) as well as the general behavior of the transient (see Figures 7 and 8). The Pseudo 3D nodalization provided closer results for the main events because it reduced the delay in the core uncovery (see Table 1). This was seen to be a consequence because of a higher vapor generation in the 1D nodalization during the phases of reflux and condensation and vapor superheating. It implied that, for similar break mass losses, liquid mass inventory decreased faster and core uncovery started before. In Figure 9, the differences between vapor generation and break mass flows are compared (the differences are calculated by substracting the values of the Pseudo 3D nodalization to the 1D nodalization).
Table 2 Initial Conditions of Test G7.1 Experimental UPC 1D UPC Pseudo 3D data nodalization nodalization Core power (Norm.) 1 0,996 0,996 Pressurizer pressure (Norm.) 1 1 1 Pressurizer liquid level (Norm.) 1 0,7 0,7 Secondary-side pressure (Norm.) 1 1 0.998 Secondary-side liquid level (Norm.) 1 0.998 1 Main feedwater temperature (Norm.) 1 1 1 Accumulators pressure (Norm.) 1 1 1 Accumulators temperature (Norm.) 1 1 1 LPIS pressure (initiation of system)
(Norm.) 1 1 1 LPIS temperature (Norm.) 1 1 1 Figure 7 Primary Pressure 9
Figure 8 Secondary Pressure Figure 9 Mass Flow Comparison Finally, the UPC PKL Pseudo 3D nodalization was qualified for reproducing the relationship between the CET and the PCT. This nodalization solved instabilities in the simulation of the overheated CET (Fig. 10), obtaining close results in the CETvsPCT curve (Fig. 11). The Pseudo 3D nodalization reproduced the same slope of the experimental data as well as the initial increase of the PCT. This agreement in the initial increase of the PCT is as consequence of reproducing qualitatively the reported delay between both temperatures (Fig. 12).
10
Figure 10 Core Exit Temperature Figure 11 CET vs PCT Correlation On the other hand, the Pseudo 3D nodalization was not qualified for reproducing closely the core quenching after SG depressurization action. Despite core refilling was simulated, in the calculation quench front achieved the top level of the active core before accumulators injection, showing a discrepancy with experimental results. In Fig. 10, each CET curve is associated with a vertical line that indicates the time in which the accumulators injection starts. The comparison shows that for both simulations, the temperatures dropped before accumulators injection, unlike experimental data, in which it did not occur.
11
Figure 12 CET and PCT versus time 12
5 RUN STATISTICS The calculations were performed on a Personal Computer with 3.0 GHz Intel Core Duo processor, 1.97 GB of RAM and Windows XP Service Pack 3 OS.
Table 3 shows main run statistics for all calculations performed in this report:
Table 3 Run Statistics Transient time CPU time Mass error ratio (s) (s) (emass/tmass)
UPC PKL 1D nodalization 15450.0 2325.66 7.5097*10-4 UPC PKL Pseudo-3D nodalization 15450.0 2920.2 5.6616*10-4 13
6 CONCLUSIONS The UPC PKL Relap5mod33, that was qualified for the PKL I boron dilution experiments, has been adjusted to the PKL-2 Test G7.1 proving its suitability to simulate the behavior of this transient. Results showed code and nodalization capabilities to reproduce main phenomena of the transient. Otherwise, some limitations were detected for reproducing condensation and core quenching during the fast depressurization induced by the steam generators valve opening.
Regarding core modeling, Pseudo 3D vessel modelling has shown to be a good tool for simulating core dryout and CET vs PCT correlation. On the other hand, some discrepancies have been detected with 1D modelling, as higher vapour generation and instabilities in the core exit temperatures during vapour superheating.
Closer results on CET vs PCT curve are a good starting point for later scaled plant applications.
15
7 REFERENCES
- 1. H. Kremin, H. limprecht, R. Güneysu, and K. Umminger. Description of the PKL III Test Facility. Technical report, Framatome ANP report, July 2001
- 2. S. Bernhard, and K. Umminger. Pkl II G7.1: SB-LOCA with Total Failure of HPSI (Counterpart Testing with ROSA/LSTF)-Quick Look Report-. Restricted Test Report, Areva, March 2012
- 3. Freixa, J. [et al.]. SBLOCA with Boron Dilution in Pressurized Water Reactors. Impact on Operation and Safety. "Nuclear Engineering and Design", April 2009, vol. 293, Num. 4, p.
749-760.
- 4. Freixa, J. SBLOCA with Boron Dilution in Pressurized Water Reactors. Impact to the Operation and Safety. Reventos, F. (director).Phd Thesis, Universitat Politcnica de Catalunya, DFEN, 2007.
- 5. Freixa, J. [et al.]. Boron Transport Model with Physical Diffusion for RELAP5. "Nuclear Technology", November 2007, vol. 160, Num. 2, p. 205-215
- 6. F. Reventós, J. Freixa, L. Batet, C. Pretel, D. Luebbesmeyer, D. Spaziani, J. Macek, F.
Lahovsky, F. Kasahara, K. Umminger, R. Wegner, An Analytical Comparative Exercise on the OECD-SETH PKL E2.2 Experiment, Nuclear Engineering and Design, Volume 238, Issue 4, April 2008, Pages 1146-1154,
- 7. RELAP5/MOD3 Code Manual. Volume IV: Models and Correlations. June, 1999.
SCIENTECH, Inc.
- 8. Martínez, V; Reventós, F.; Pretel, C.; Post-Test Calculation of the ROSA/LSTF Test 3-1 using RELAP5/mod3.3; NUREG/IA-409, 2012
- 9. Martínez, V; Reventós, F.; Pretel, C.; Post-Test Calculation of the ROSA/LSTF Test 3-2 using RELAP5/mod3.3; NUREG/IA-410, 2012
- 10. Martínez, V.; Reventós, F.; Pretel, C.; Sol, I.; Code Validation and Scaling of the ROSA/LSTF Test 3-1 Experiment; International Topical Meeting on Safety of Nuclear Installations TopSafe 2008; Octubre 2008; Croacia 17
NUREG/IA-0499 International Agreement Report Post-Test Calculation of the PKL-2 Test G7.1 Using RELAP5/MOD3.3 Prepared by:
V Martinez-Quiroga, F. Reventos, and C. Pretel Institute of Energy Technologies Technical University of Catalonia ETSEIB, Av. Diagonal 647, Pav. C 08028 Barcelona Spain Kirk Tien, NRC Project Manager Division of Systems Analysis Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, DC 20555-0001 Manuscript Completed: November 2017 Date Published: February 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 MATERIALS IN NRC PUBLICATIONS NRC Reference Material Non-NRC Reference Material As of November 1999, you may electronically access Documents available from public and special technical NUREG-series publications and other NRC records at libraries include all open literature items, such as books, NRCs Library at www.nrc.gov/reading-rm.html. Publicly journal articles, transactions, Federal Register notices, released records include, to name a few, NUREG-series Federal and State legislation, and congressional reports.
publications; Federal Register notices; applicant, Such documents as theses, dissertations, foreign reports licensee, and vendor documents and correspondence; and translations, and non-NRC conference proceedings NRC correspondence and internal memoranda; bulletins may be purchased from their sponsoring organization.
and information notices; inspection and investigative reports; licensee event reports; and Commission papers Copies of industry codes and standards used in a and their attachments. substantive manner in the NRC regulatory process are maintained at NRC publications in the NUREG series, NRC regulations, The NRC Technical Library and Title 10, Energy, in the Code of Federal Regulations Two White Flint North may also be purchased from one of these two sources. 11545 Rockville Pike Rockville, MD 20852-2738
- 1. The Superintendent of Documents U.S. Government Publishing Office These standards are available in the library for reference Mail Stop IDCC use by the public. Codes and standards are usually Washington, DC 20402-0001 copyrighted and may be purchased from the originating Internet: bookstore.gpo.gov organization or, if they are American National Standards, Telephone: (202) 512-1800 from Fax: (202) 512-2104 American National Standards Institute 11 West 42nd Street
- 2. The National Technical Information Service New York, NY 10036-8002 5301 Shawnee Road www.ansi.org Alexandria, VA 22312-0002 (212) 642-4900 www.ntis.gov 1-800-553-6847 or, locally, (703) 605-6000 Legally binding regulatory requirements are stated only in A single copy of each NRC draft report for comment is laws; NRC regulations; licenses, including technical available free, to the extent of supply, upon written specifications; or orders, not in NUREG-series publications.
request as follows: The views expressed in contractor prepared publications in this series are not necessarily those of the NRC.
Address: U.S. Nuclear Regulatory Commission The NUREG series comprises (1) technical and adminis-Office of Administration trative reports and books prepared by the staff (NUREG-XXXX) or agency contractors (NUREG/CR-XXXX), (2)
Multimedia, Graphics, and Storage &
proceedings of conferences (NUREG/CP-XXXX), (3) reports Distribution Branch resulting from international agreements (NUREG/IA-XXXX),
Washington, DC 20555-0001 (4) brochures (NUREG/BR-XXXX), and (5) compilations of E-mail: distribution.resource@nrc.gov legal decisions and orders of the Commission and Atomic Facsimile: (301) 415-2289 and Safety Licensing Boards and of Directors decisions under Section 2.206 of NRCs regulations (NUREG-0750).
Some publications in the NUREG series that are posted at NRCs Web site address www.nrc.gov/reading-rm/ DISCLAIMER: This report was prepared under an doc-collections/nuregs are updated periodically and may international cooperative agreement for the exchange of technical information. Neither the U.S. Government nor any differ from the last printed version. Although references to agency thereof, nor any employee, makes any warranty, material found on a Web site bear the date the material expressed or implied, or assumes any legal liability or was accessed, the material available on the date cited responsibility for any third partys use, or the results of such may subsequently be removed from the site. 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-0499 International Agreement Report Post-Test Calculation of the PKL-2 Test G7.1 Using RELAP5/MOD3.3 Prepared by:
V Martinez-Quiroga, F. Reventos, and C. Pretel Institute of Energy Technologies Technical University of Catalonia ETSEIB, Av. Diagonal 647, Pav. C 08028 Barcelona Spain Kirk Tien, NRC Project Manager Division of Systems Analysis Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, DC 20555-0001 Manuscript Completed: November 2017 Date Published: February 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
ABSTRACT The Thermalhydraulic Studies Group of Technical University of Catalonia (UPC) holds a large background in nuclear safety studies in the field of Nuclear Power Plant (NPP) code simulators.
This report analyzes with RELAP5mod3.3 the G7.1 PKL Test. This experiment is part of a Counterpart Test performed in LSTF and PKL Test Facilities within the framework of the OECD/NEA ROSA-2 and PKL-2 projects. Detailed core nodalizations and Pseudo 3D modeling have been object of study as well as the capabilities of the code for reproducing the correlation between the Core Exit Temperature (CET) and the Peak Cladding Temperature (PCT).
iii
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 80s 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 todays 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/FAP1 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 enough investigated safety issues and phenomena (e.g., boron dilution, low power and 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).
v
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 and investigation of specific safety aspects of PKL/OECD ROSA/OECD and ATLAS/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 vi
TABLE OF CONTENTS ABSTRACT ............................................................................................................................... iii FOREWORD .............................................................................................................................. v LIST OF FIGURESix LIST OF TABLES..ix EXECUTIVE
SUMMARY
.........xi ACKNOWLEDGMENTS .......................................................................................................... xiii ABBREVIATIONS AND ACRONYMS ..................................................................................... xv 1 INTRODUCTION ...............................................................................................................1 1.1 PKL-2 and ROSA-2 Counterpart Test .......................................................................1 1.2 The OECD PKL 2 Project ..........................................................................................2 2 FACILITY AND TEST DESCRIPTION ...............................................................................3 2.1 PKL III Test Facility ...................................................................................................3 2.2 Experimental Conditions ...........................................................................................4 2.3 Initial Conditions .......................................................................................................4 2.4 Test Phase ...............................................................................................................4 3 CODE INPUT MODEL DESCRIPTION ..............................................................................7 4 RESULTS ..........................................................................................................................9 5 RUN STATISTICS ........................................................................................................... 13 6 CONCLUSIONS .............................................................................................................. 15 7 REFERENCES ................................................................................................................ 17 vii
LIST OF FIGURES Figure 1 Diagram with the Different Counterpart Conditions (Courtesy of the OECD/NEA ROSA-2 Group) ....................................................................................................... 1 Figure 2 OECD PKL 2 Project Experiments (Courtesy of the OECD/NEA PKL-2 Group)....... 2 Figure 3 PKL Test Facility (Courtesy of the OECD/NEA PKL Group)..................................... 3 Figure 4 Nodalization of the Main Steam System................................................................... 7 Figure 5 Nodalization of the Bypass Heaters ......................................................................... 7 Figure 6 UPC PKL Pseudo 3D Nodalization Core Channels .................................................. 8 Figure 7 Primary Pressure ..................................................................................................... 9 Figure 8 Secondary Pressure ...............................................................................................10 Figure 9 Mass Flow Comparison ..........................................................................................10 Figure 10 Core Exit Temperature............................................................................................11 Figure 11 CET vs PCT Correlation .........................................................................................11 Figure 12 CET and PCT Versus Time.....................................................................................12 LIST OF TABLES Table 1 Chronology of the Main Events of Test G7.1 ............................................................ 5 Table 2 Initial Conditions of Test G7.1 .................................................................................. 9 Table 3 Run Statistics ..........................................................................................................13 ix
EXECUTIVE
SUMMARY
Experimental research activities are being performed in Germany by the OECD PKL 2 project with the aim of addressing thermal-hydraulic safety issues for current PWR and new PWR design concepts. These experiments are carried out at the PKL III test facility.
This report analyses the experiment G7.1 of the LSTF and PKL Counterpart Test, which was carried out as a part of the OECD/NEA ROSA-2 and PKL-2 projects. The aim of this international synergy was to analyze the effectiveness of Core Exit Temperature measurement in Accident Management strategies as well as the scaling effects that appear between counterpart transients performed at different sizes and designs.
UPC PKL Relap5mod33 nodalization, created for the simulation of the PKL boron dilution experiments, has been used and improved for this test. Two different nodalizations were prepared in order to check the capabilities of the code for simulating properly the correlation between the PCT and CET: one, the UPC PKL 1D nodalization, with just one core channel; and a second, the UPC Pseudo 3D nodalization, with concentric core channels and activating the transversal momentum equations.
Many other aspects related to the nodalization were adjusted and verified in order to improve results.
In general, the results of the simulation demonstrated that 1D nodalizations are good enough for describing the general behavior of the transient as well as the main events and phenomena. On the other hand, for simulating accurately the correlation between the CET and PCT, Pseudo 3D modeling is necessary.
xi
ACKNOWLEDGMENTS This paper contains findings that were produced within the OECD-NEA PKL-2 Project. The authors are grateful to the Management Board of the PKL-2 Project for their consent to this publication.
xiii
ABBREVIATIONS AND ACRONYMS ACC accumulator AM accident management CET core exit temperature ECCS emergency core cooling system EOP emergency operational procedure HPIS high pressure injection system HS heat structure LPIS low pressure injection system LSTF large scale test facility MS main steam MSIV main steam isolation valve MSLB main steam line break NEA Nuclear Energy Agency NPP nuclear power plant OECD Organization for Economic Cooperation and Development PCT peak cladding temperature PKL Primrkreislauf PZR pressurizer PWR pressurized water reactor RELAP reactor excursion and leak analysis program ROSA rig of safety assessment SBLOCA small break loss of coolant accident SG steam generator UP upper plenum UPC Universitat Politcnica de Catalunya (Technical University of Catalonia)
UT u-tubes xv
1 INTRODUCTION Several safety activities have been performed during the last decades under the auspices of the OECD to develop and improve computer codes. They include several experiments at integral test facilities like the Test G7.1, which forms part of the OECD PKL-2 project.
1.1 PKL-2 and ROSA-2 Counterpart Test In 2011 a Counterpart Test was performed in LSTF and PKL Test Facilities as a part of the OECD/NEA ROSA-2 and PKL-2 projects (Test 3 and test G7.1 respectively). The objective of both tests was devoted to analyze two aspects:
Core Exit Temperature measurement effectiveness in Accident Management of NPPs Scaling effects between PKL and LSTF Test Facilities In Accident Management strategies, core exit temperature measurement becomes critical for detecting core dryout and for avoiding that PCT rises until safety limits. In the Counterpart Test, the relationship between CET and PCT was object of study in order to analyze EOP set points.
The selected scenario was a hot leg SBLOCA. System failures as no high pressure safety injection and no automatic secondary-side safety cooldown were imposed. The particular test conditions for the PKL-2 Test G7.1 are described in section 2.2. The main phenomena of interest were:
Core boil-off with steam generation Steam flow toward hot-leg break Realistic pressure during core boil-off Relationship between PCT and CET Figure 1 Diagram with the Different Counterpart Conditions (Courtesy of the OECD/NEA ROSA-2 Group) 1
About the scaling issue, core dryout and AM actions were simulated under counterpart conditions at different sizes (PKL -1:145- and LSTF -1:48- ) and pressures (LSTF high and low pressure transient phases) in order to check how the scaling affects the CET and its related phenomena, and in order to sound out possible strategies for extrapolating the results to a commercial NPP (see figure one). The direct comparison of both transients was established between the whole transient of the PKL Test G7.1, and the low pressure transient phase of the LSTF Test 3.
1.2 The OECD PKL 2 Project The OECD PKL 2 is an international project carried out between 2008 and 2011 that includes several types of integral experiments (see figure two) on PKL III, ROCOM and PMK test facilities. The aim of the project is to address thermal-hydraulic safety issues for current PWR and new PWR design concepts. It also provides a wide database for the validation of computer codes and models for system integral analyses coupled with detailed analyses of local phenomena.
Begin of Project 1. PRG/MB 2. PRG/MB 3. PRG/MB 4. PRG/MB 5. PRG/MB 6. PRG/MB 7. PRG/MB 8. PRG/MB 2008 2009 2010 2011 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12 G1 (2 Tests)
PMK (2 Tests)
G2 (1 Test)
G3 (1 Test)
ROCOM Tests G4 (1 Test)
G5 (1 Test)
G6 (1 Test)
G7 (1 Test)
Figure 2 OECD PKL 2 Project Experiments (Courtesy of the OECD/NEA PKL-2 Group)
The topics of the experiments are:
G1 Systematic investigation of the heat transfer mechanisms in the SGs in presence of nitrogen, steam and water G2 Cooldown procedures with SGs isolated and emptied on the secondary side G3 Fast cooldown transients (main steam line break)
G4 Systematic study on heat transfer in the SG under reflux condenser conditions G5 Boron precipitation processes after LB - LOCA G6 RCS cooldown with void formation in RPV upper head G7 Counterpart Test with the OECD/NEA ROSA-2 Test 3 (hot leg SBLOCA) 2
2 FACILITY AND TEST DESCRIPTION 2.1 PKL III Test Facility PKL is an experimental power facility plant designed to simulate pressurized water reactors (PWR) under accidental conditions. The plant, situated in Erlangen (Germany), has been the scenario of several experiments in the last 25 years. The PKL facility replicates the entire primary system and most of the secondary system (except for the turbine and condenser) of a 1300-MW PWR plant, with elevations scaled 1:1 and diameters reduced by a factor of 12.
Therefore volumes and power have been reduced by a factor of 145. The number of rods in the core and the U-tubes in the steam generator has been divided by 145 too. The core has been modeled by 314 electrical heater rods. Unlike many experimental facilities with only two available loops (one for the break loop and one to simulate the other three loops) PKL simulates all four loops separately. This is very important in order to analyze asymmetrical transients, e.g.
with injection in two out of four loops. A diagram of the PKL test facility is shown in Figure 3.
Figure 3 PKL Test Facility (Courtesy of the OECD/NEA PKL Group)
The operating pressure of the PKL facility is limited to 45 bars on the primary side and to 56 bars on the secondary side. This allows simulation over a wide temperature range (322K to 522K) that is particularly applicable to the cooldown procedures investigated.
3
All emergency systems are represented and have a wide versatility referred to their functions and positions. There are 8 accumulators (one for each cold and hot leg). The pump injection system is available in all the hot and cold legs. Many break locations are available too.
PKL test facility has about 1500 measurement points that permit an exhaustive analysis of the tests. There are measurement devices for cladding, wall and fluid temperature, absolute and differential pressure, one and two phase mass flow, density and boron concentration.
Sixty of the measurement devices are identical to those that are used in a commercial plant to simulate what an operator would control in case of accident.
2.2 Experimental Conditions The hardware configuration of PKL is described in references [1] and [2]. Some important points are the following:
- Break: upward oriented SBLOCA (1.5 %) at Hot Leg without PZR.
- ECCs: HPIS full failure.
ACCs and LPIS set points fixed at 26 bars and 8 bars respectively.
ACC level and LPIS mass flow adjusted according to LSTF final conditions.
- Core power curve: constant at 1.8 % of the nominal power. Additional heat was added in order to compensate the differences between the environment losses of LSTF and PKL.
- SG depressurization: fully opening of 2 main steam relief valves when CET achieves 350
ºC. The depressurization of the 4 SGs was simulated connecting them via main steam header.
- Main steam relief valves: modified for having same ratio of LSTF valves.
2.3 Initial Conditions The PKL initial conditions were adjusted at its maximum pressures in order to reproduce as realistic as it can PWR SBLOCA reflux and condensation, core dry-out and Accident Management phases.
In relation to scaling, the mass inventory in the secondary side was adjusted using Kv factor in order to have the same ratio between liquid and energy storage in the SGs. In the primary system, mass inventory was tuned for having the same reflux condenser and break discharging initial conditions than the low pressure transient phase of the LSTF Test 3.
2.4 Test Phase The transient starts at t=0 seconds with the opening of the break valve. During the first 940 seconds, core is cooled under saturated conditions and reflux and condensation occurs in the UTs of the SGs. In this phase of the transient, primary pressure keeps constant over the 44 bars of the isolated secondary side. After that, core uncovery starts, and few seconds later (at 1020 seconds), primary pressure drops below secondary pressure. In this second phase, there is vapour superheating in the core, and the CET and the PCT rise above the temperature of saturation (with a delay of 270 seconds between both).
When the CET temperature achieves the AM temperature set point (1360 seconds), SGs are depressurized by MS relief valves opening. System pressures drop rapidly without a complete 4
core quenching, so CET and PCT temperatures do not drop to saturated conditions until accumulator injection system is activated (1500 seconds).
Finally, when primary pressure drops below 8 bars (2060 seconds), LPIS starts to inject water, compensating break losses and keeping constant plant parameters. At 5685 seconds, break valve is closed and the end of the transient is declared.
The main events are described in Table 1:
Table 1 Chronology of the Main Events of Test G7.1 Experimental UPC 1D UPC Pseudo-3D (s) nodalization (s) nodalization (s)
Start of the transient 0 0 0 Begin of core uncovery 940 800 940 Primary pressure below 1020 920 1010 secondary pressure Secondary side 1360 1190 1295 depressurization Start of accum. Injection 1500 1304 1450 ACC injection finished 1860 1712 1752 LPIS started 2060 1966 1993 End of the test 5685 5685 5685 5
3 CODE INPUT MODEL DESCRIPTION Two different nodalizations were prepared by making use of the UPC PKL RELAP5mod3.3 nodalization. This nodalization had been previously qualified with the simulation of the OECD/NEA PKL I tests F1.1 and F1.2, and the OECD/SETH tests E2.1, E2.2 and E2.3 (see references [3], [4] and [5] and [6]). Particular Counterpart nodalizations were:
UPC PKL 1D nodalization UPC PKL Pseudo-3D nodalization Several improvements were done for both nodalizations taking advantage of the particular specifications of the transient. In this way, main steam system was modified following the sketch shown in Figure 4. Moreover, realistic bypass heaters were added to the PZR and the SGs (as in Figure 5) in order to have a proper power balance between the environment heat losses and the power supplied by the exchangers.
Figure 4 Nodalization of the Main Steam System Figure 5 Nodalization of the Bypass Heaters 7
The differences between both nodalizations were exclusively related to the core and upper plenum (UP) modeling. UPC PKL 1D nodalization simulated them with one channel (in addition to the core bypass), having fuel and all passive heat structures (core barrel and unheated rods) linked to the same volumes. The fuel was modeled with three HSs, with the same power ratio and divided in 7 axial levels.
UPC PKL Pseudo-3D nodalization had the core and UP (until the CET thermocouple level) divided in three radial channels (see Figure 6), with one fuel HS for each channel. The HSs for the passive internal metal structures were split for each channel proportionally to the flow path of each one, and the core barrel was linked to the outer zone. The radial flow paths between cells were modeled and transversal momentum equations were activated following the recommendations of reference [7].The total number of core axial meshes was increased to 14 and the UP cell heights were adjusted so that the center of the node coincided with the elevation of the thermocouples in the test facility.
Figure 6 UPC PKL Pseudo 3D Nodalization Core Channels 8
4 RESULTS The results obtained for both nodalizations showed a quite close agreement reproducing the initial conditions (see Table 2) as well as the general behavior of the transient (see Figures 7 and 8). The Pseudo 3D nodalization provided closer results for the main events because it reduced the delay in the core uncovery (see Table 1). This was seen to be a consequence because of a higher vapor generation in the 1D nodalization during the phases of reflux and condensation and vapor superheating. It implied that, for similar break mass losses, liquid mass inventory decreased faster and core uncovery started before. In Figure 9, the differences between vapor generation and break mass flows are compared (the differences are calculated by substracting the values of the Pseudo 3D nodalization to the 1D nodalization).
Table 2 Initial Conditions of Test G7.1 Experimental UPC 1D UPC Pseudo 3D data nodalization nodalization Core power (Norm.) 1 0,996 0,996 Pressurizer pressure (Norm.) 1 1 1 Pressurizer liquid level (Norm.) 1 0,7 0,7 Secondary-side pressure (Norm.) 1 1 0.998 Secondary-side liquid level (Norm.) 1 0.998 1 Main feedwater temperature (Norm.) 1 1 1 Accumulators pressure (Norm.) 1 1 1 Accumulators temperature (Norm.) 1 1 1 LPIS pressure (initiation of system)
(Norm.) 1 1 1 LPIS temperature (Norm.) 1 1 1 Figure 7 Primary Pressure 9
Figure 8 Secondary Pressure Figure 9 Mass Flow Comparison Finally, the UPC PKL Pseudo 3D nodalization was qualified for reproducing the relationship between the CET and the PCT. This nodalization solved instabilities in the simulation of the overheated CET (Fig. 10), obtaining close results in the CETvsPCT curve (Fig. 11). The Pseudo 3D nodalization reproduced the same slope of the experimental data as well as the initial increase of the PCT. This agreement in the initial increase of the PCT is as consequence of reproducing qualitatively the reported delay between both temperatures (Fig. 12).
10
Figure 10 Core Exit Temperature Figure 11 CET vs PCT Correlation On the other hand, the Pseudo 3D nodalization was not qualified for reproducing closely the core quenching after SG depressurization action. Despite core refilling was simulated, in the calculation quench front achieved the top level of the active core before accumulators injection, showing a discrepancy with experimental results. In Fig. 10, each CET curve is associated with a vertical line that indicates the time in which the accumulators injection starts. The comparison shows that for both simulations, the temperatures dropped before accumulators injection, unlike experimental data, in which it did not occur.
11
Figure 12 CET and PCT versus time 12
5 RUN STATISTICS The calculations were performed on a Personal Computer with 3.0 GHz Intel Core Duo processor, 1.97 GB of RAM and Windows XP Service Pack 3 OS.
Table 3 shows main run statistics for all calculations performed in this report:
Table 3 Run Statistics Transient time CPU time Mass error ratio (s) (s) (emass/tmass)
UPC PKL 1D nodalization 15450.0 2325.66 7.5097*10-4 UPC PKL Pseudo-3D nodalization 15450.0 2920.2 5.6616*10-4 13
6 CONCLUSIONS The UPC PKL Relap5mod33, that was qualified for the PKL I boron dilution experiments, has been adjusted to the PKL-2 Test G7.1 proving its suitability to simulate the behavior of this transient. Results showed code and nodalization capabilities to reproduce main phenomena of the transient. Otherwise, some limitations were detected for reproducing condensation and core quenching during the fast depressurization induced by the steam generators valve opening.
Regarding core modeling, Pseudo 3D vessel modelling has shown to be a good tool for simulating core dryout and CET vs PCT correlation. On the other hand, some discrepancies have been detected with 1D modelling, as higher vapour generation and instabilities in the core exit temperatures during vapour superheating.
Closer results on CET vs PCT curve are a good starting point for later scaled plant applications.
15
7 REFERENCES
- 1. H. Kremin, H. limprecht, R. Güneysu, and K. Umminger. Description of the PKL III Test Facility. Technical report, Framatome ANP report, July 2001
- 2. S. Bernhard, and K. Umminger. Pkl II G7.1: SB-LOCA with Total Failure of HPSI (Counterpart Testing with ROSA/LSTF)-Quick Look Report-. Restricted Test Report, Areva, March 2012
- 3. Freixa, J. [et al.]. SBLOCA with Boron Dilution in Pressurized Water Reactors. Impact on Operation and Safety. "Nuclear Engineering and Design", April 2009, vol. 293, Num. 4, p.
749-760.
- 4. Freixa, J. SBLOCA with Boron Dilution in Pressurized Water Reactors. Impact to the Operation and Safety. Reventos, F. (director).Phd Thesis, Universitat Politcnica de Catalunya, DFEN, 2007.
- 5. Freixa, J. [et al.]. Boron Transport Model with Physical Diffusion for RELAP5. "Nuclear Technology", November 2007, vol. 160, Num. 2, p. 205-215
- 6. F. Reventós, J. Freixa, L. Batet, C. Pretel, D. Luebbesmeyer, D. Spaziani, J. Macek, F.
Lahovsky, F. Kasahara, K. Umminger, R. Wegner, An Analytical Comparative Exercise on the OECD-SETH PKL E2.2 Experiment, Nuclear Engineering and Design, Volume 238, Issue 4, April 2008, Pages 1146-1154,
- 7. RELAP5/MOD3 Code Manual. Volume IV: Models and Correlations. June, 1999.
SCIENTECH, Inc.
- 8. Martínez, V; Reventós, F.; Pretel, C.; Post-Test Calculation of the ROSA/LSTF Test 3-1 using RELAP5/mod3.3; NUREG/IA-409, 2012
- 9. Martínez, V; Reventós, F.; Pretel, C.; Post-Test Calculation of the ROSA/LSTF Test 3-2 using RELAP5/mod3.3; NUREG/IA-410, 2012
- 10. Martínez, V.; Reventós, F.; Pretel, C.; Sol, I.; Code Validation and Scaling of the ROSA/LSTF Test 3-1 Experiment; International Topical Meeting on Safety of Nuclear Installations TopSafe 2008; Octubre 2008; Croacia 17
NUREG/IA-0499 Post-Test Calculation of the PKL2 Test G7.1 Using RELAP/MOD3.3 February 2019 V. Martinez-Quiroga, F. Reventos, C. Pretel Technical Institute of Energy Technologies Universitat Politecnica de Catalunya ETSEIB, Av, Diagonal 647, Pav C 08028 Barcelona, Spain Division of Systems Analysis Office of Nuclear Regulatory Research U.S. Nuclear Regulatory Commission Washington, D.C. 20555-0001 K. Tien, NRC Project Manager The Thermalhydraulic Studies Group of Technical University of Catalonia (UPC) hold a large background in nuclear safety studies in the field of Nuclear Power Plan (NPP) code simulators. This report analyses with RELAP5 mod3.3 the G7.1 PKL Test. This experiment is part of a Counterpart Test performed in LSTF and PKL Test Facilities within the framework of the OECD/NEA ROSA-2 and PKL-2 projects. Detailed core nodalizations and Pseudo 3D modeling have been object of study as well as the capabilities of the code for reproducing the correlation between the Core Exit Temperature (CET) and the Peak Cladding Temperature (PCT). This report analyses the experiment F7.1 of the LSTF and PKL Counterpart Test, which was carried out as part of the OECD/NEA ROSA-2 and PKL-2 projects. The aim of this international synergy was to analyze the effectiveness of Core Exit Temperature measurement in Accident Management strategies as well as the scaling effects that appear between counterpart transients performed at different sizes and designs.
Primakrieslauf (PKL)
Committee of the Safety of Nuclear Installations (CSNI)
Consejo de Seguridad Nuclear (Spanish Nuclear Regulatory Commission, CSN)
Core Exit Temperature (CET)
Severe Accident Management Guideline/Guidance (SMAG)
Large Scale Test Facility (LSTF)
Universitat Politecnica de Cataluny (Technical University of Catalonia) UPC
Post-Test Calculation of the PKL2 Test G7.1 Using RELAP/MOD3.3 February 2019 NUREG/IA-0499