THE CASUALIDAD METHOD FOR UNCERTAINTY EVALUATION OF BEST-ESTIMATE SYSTEM THERMAL-HYDRAULICS CALCULATIONS
A. Petruzzi
Nuclear and INdustrial Engineering (NINE), Borgo Giannotti, Lucca, Italy
17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-17)
Xi'an, China, 03 - 08 September 2017
Abstract — The present paper deals with the description of the salient features of three independent approaches for estimating uncertainties associated with predictions of complex system codes. The 1st approach is the “standard” one and the most used at the industrial level: it is based upon the selection of input uncertain parameters, on assigning related ranges of variations and, possibly, PDF (Probability Density Functions) and on performing a suitable number of code runs to get the combined effect of variation on the results. In the 2nd approach the uncertainty derives from the comparison between relevant measured data and results of corresponding code calculations. The 3rd approach is based on the Bayesian inference technique and on the availability of experimental data by which computer model predictions can be improved and the ranges of variation of (in theory) „all‟ input parameters can be characterized. More details are provided in respect with the third approach that has been named CASUALIDAD (Code with the capability of Adjoint Sensitivity and Uncertainty AnaLysis by Internal Data ADjustment and assimilation).
ATUCHA-1 NPP CONTAINMENT VENTING ANALYSIS FOLLOWING SBO AND LBLOCA EVENTS BY GOTHIC CODE
A. Popa b, W. Giannottia, F. Terzuolia, A. Petruzzia, N. Forgioneb, O. Mazzantinic
a Nuclear and INdustrial Engineering (NINE), Borgo Giannotti, Lucca, Italy
b University of Pisa, Largo Lucio Lazzarino 2, 56126 Pisa, Italy
c Nucleoeléctrica Argentina Arribeños 3619, C1429BKQ, Buenos Aires, Argentina
17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-17)
Xi'an, China, 03 - 08 September 2017
Abstract — Containment behaviour plays a key role in the safety framework of a Nuclear Power Plant. The GOTHIC thermal hydraulic code has been adopted to evaluate the Atucha-1 NPP containment responses during two postulated accident scenarios, Station Black Out and Large Break Loss of Coolant Accident, while assuming the external cooling of the Reactor Pressure Vessel is carried out during the transients.
The Atucha-1 NPP has a containment designed to work at full pressure, constituted by a steel sphere enveloped by a concrete shell, and having an annular gap of air in between. The target of the analysis is the evaluation of the effects caused by the additional production of steam in the reactor cavity as a consequence of the ex-vessel cooling, which could cause an excessive pressurization of the containment, and lead to pressure values above the safety limit. The containment pressure and temperature, the distribution of hydrogen in the containment atmosphere and the water hold-up in the most relevant rooms have been monitored as target variables.
Each accident scenario was simulated using two different nodalizations, characterized by a different level of refinement. The "detailed" nodalization is meant to be the most refined nodalization according to the available computational resources; having high fidelity three dimensional details, with a high number cells. While the "coarse" nodalization was developed in order to lower the demand for computational resources without significantly compromising the global scenario response. Both nodalizations are characterized by high complexity in the representation of rooms and their connections, e.g. all doors and blow off panels have been simulated to open with the designed differential pressure logic.
Both accident transients, for each type of nodalization, were simulated for 200,000 seconds. Results showed that for the Large Break Loss of Coolant Accident pressure is predicted to reach around 5.25 bar, while the Station Black Out Scenario reaches 4.4 bar.
RELAP5-3D ANALYSIS AND QUALIFICATION OF EBR-II SHRT-17
D. De Luca, A. Petruzzi, M. Cherubini
Nuclear and INdustrial Engineering (NINE), Borgo Giannotti, Lucca, Italy
17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-17)
Xi'an, China, 03 - 08 September 2017
Abstract — Coordinated Research Project (CRP) on EBR-II Shutdown Heat Removal Tests (SHRT) was established by International Atomic Energy Agency (IAEA) in June 2012. The objective of the project is to support and to improve validation of simulation tools and projects for Sodium-cooled Fast Reactors (SFR). In the frame of this project, benchmark analysis of one of the EBR-II shutdown heat removal tests, protected loss-of-flow transient (SHRT-17), has been performed.
The aim of this paper is to present modeling of EBR-II reactor design using RELAP-3D, to show the results of the transient analysis of SHRT-17, and to discuss the results of application of the qualification process, in particular of the Fast Fourier Transform Based Method (FFTBM) to perform a quantitative accuracy evaluation of the model developed.
Complete nodalization of the reactor was made from the beginning. After achievement of acceptable steady-state results, transient analysis of the protected (with scram) loss-of-flow test was performed. Thermal-hydraulics characteristics and plant behavior was focused on prediction of natural convection cooling by evaluating the reactor core flow and temperatures and their comparison with experimental data that were provided by ANL.
Finally, the process of qualification of a system thermal-hydraulic code calculation was applied. It consists of three steps: 1) the geometrical fidelity of the nodalization, related with the evaluation and comparison of the geometrical data of the hardware respect to the estimated numerical values implemented in the nodalization; 2) the steady state level qualification, dealing with the capability of the nodalization to reproduce the steady state qualified conditions of the system; 3) the “on-transient” qualification, necessary to demonstrate the capability of the code and of the developed nodalization to reproduce the relevant thermal-hydraulic phenomena expected during the transient. The latter is a very complex step which foreseen different phases following our methodology of qualification (SCCRED, Standardized and Consolidated Calculated & Reference Experimental Database methodology). In the framework of the benchmark, the focus was only on the so called “Quantitative Accuracy Evaluation” that is performed by the FFTBM.
IAEA's Coordinated Research Project on EBR-II Shutdown Heat Removal Tests
Chirayu Batra (IAEA), Vladimir Kriventsev (IAEA), Laural L. Briggs (ANL), Stefano Monti (IAEA), Wenjun HU (CIAE), Dan-ting Sui (North ChinaElectric Power Uniersity), Guanghui Su (Xi'an Jiaotong University), Ludovic Maas (IRSN (French Institute for Nuclear Safety)), Barbara Vezzoni (KIT), Partha Sarathy UPPALA (IGCAR), Alessandro Del Nevo (ENEA), Alessandro Petruzzi (Nuclear and INdustrial Engineering (NINE)), Roberto Zanino (Politecnico di Torino-Italy), Hiroaki Ohira (JAEA), Willem F. van Rooijen (University of Fukui), K. Morita (Kyushu Univ), Chiwoong CHOI (KAERI), Andong Shin (KINS), M. Stempniewicz (Nuclear Research & Consultancy Group), Nikita Rtischev (IBRAE), Konstantin Mikityuk (Scherrer Inst), Ethan A. Bates (TerraPower)
17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-17)
Xi'an, China, 03 - 08 September 2017
Abstract — A Coordinated Research Project (CRP) on “Benchmark Analysis of EBR-II Shutdown Heat Removal Tests (SHRT)” was launched by the International Atomic Energy Agency (IAEA) in 2012. A series of transient tests were conducted on the EBR-II reactor at Argonne National Laboratory (ANL) to improve the understanding of thermal hydraulics and neutronics of fast reactors. Shutdown heat removal tests conducted in 1984 and 1986 demonstrated mechanisms by which fast reactors can survive severe accident initiators with no core damage. Two SHRT tests, SHRT-17 representing Protected Loss of Flow (PLOF) transients and SHRT-45R representing Unprotected Loss of Flow (ULOF) transients, were studied in the IAEA CRP.
The objectives of the CRP were to improve design and simulation capabilities in fast reactor thermal hydraulics, neutronics and safety analyses through benchmark analysis of these two important tests. At the first stage of the benchmark, ANL provided the input data on EBR-II geometry, as well as initial and boundary conditions for the SHRT-17 and SHRT-45R tests to perform “blind” calculations. At the second stage, ANL released the experimental observations and participants had the chance to analyze the difference and refine the models. At the third stage, a methodology to systematically analyze and compare the models and the results of each participant was applied. Nineteen organizations from eleven countries participated in the CRP, making it one of the largest CRP coordinated by the IAEA fast reactor team.
The paper provides a general CRP overview, gives the basics of the EBR-II reactor design, describes the shutdown heat removal tests, the benchmark setup, and results of numerical simulations, followed by a detailed discussion on the EBR-II CRP.
SUMMARY OF SWINTH-2016: A SPECIALISTS WORKSHOP ON ADVANCED INSTRUMENTATION AND MEASUREMENT TECHNIQUES FOR NUCLEAR REACTOR THERMAL-HYDRAULICS EXPERIMENTATION
F. Morettia, F. D'Auriab, N. Aksanb, S. Lutsanychb, K. Ummingerc, S. Guptad, D. Bestione, U. Hampelf, K.Y. Choig, H. Purhonenh, A. Gubai, D. Paladinoi j
aNuclear and INdustrial Engineering (NINE), Borgo Giannotti, Lucca, Italy
bSan Piero a Grado Nuclear Research Group (GRNSPG) – University of Pisa, Italy
cAREVA GmbH, Germany
dBecker Technologies GmbH, Germany
eCommissariat à l’Energie Atomique et aux Energies Alternatives (CEA) – Grenoble, France
fHelmholtz Zentrum Dresden Rossendorf (HZDR), Germany
gKorean Atomic Energy Research Institute (KAERI), Republic of Korea
hLappeenranta University of Technology (LUT), Finland
iHungarian Academy of Science Centre for Energy Research (MTA EK), Hungary
jPaul Scherrer Institut (PSI), Switzerland
17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-17)
Xi'an, China, 03 - 08 September 2017
Abstract — A “Specialists Workshop on Advanced Instrumentation and Measurement Techniques for Nuclear Reactor Thermal-Hydraulics (SWINTH)” was held in Italy on June 2016, promoted and organized by SILENCE, the “Significant Light and Heavy Water Reactor Thermal Hydraulic Experiments Network for the Consistent Exploitation of the Data”. The workshop covered the technology of experimentation in nuclear thermal-hydraulics, including both separate effect and integral test facilities, with coarse as well as CFD-grade measurements (although most contributions were actually devoted to the latter type), and not limited to light water-cooled reactor technology. The last international workshop with the same scope was an OECD/NEA/CSNI meeting held in Santa Barbara, CA, in 1997.
The main objectives of the workshop were to gain an up-to-date picture of the current trends in experimental programs and in technological advancement as to reactor thermal-hydraulics, to promote the technical-scientific exchange between experimentalists and scientists operating on the code development and validation side, and to help identify the existing gaps between the code validation needs and the available technology. The response of the international community to this initiative was very positive and brought valuable contribution to meeting the above-mentioned objectives through generally high-quality papers and presentations, many of which focusing on innovative solutions and state-of-the-art applications. The present paper provides an overview of the contents of the workshop papers and tries to identify lessons learned and recommendations.
In addition, questionnaires “to identify measurement needs and main gaps for further system code and CFD code development and validation” had been distributed prior to the workshop among a number of people active in this area, and the preliminary outcomes of the survey were presented at the workshop. This paper also gives an account of such outcomes. Any feedback to this paper will be taken into account for the organization of the future SWINTH workshops being planned, possibly within an OECD/NEA framework and with a scope extended to in-plant instrumentation (e.g. for Severe Accident Management).
Application of Integrated RELAP5-Transuranus Approach to a RIA Benchmark Exercise Including Uncertainty and Sensitivity
F. Moretti, M. Cherubini, A. Pop
Nuclear and INdustrial Engineering (NINE), Borgo Giannotti, Lucca, Italy
17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-17)
Xi'an, China, 03 - 08 September 2017
Abstract — A Reactivity Initiated Accident (RIA) Fuel Codes Benchmark was organized by the OECD/NEA Working Group on Fuel Safety. While the phase-I of the exercise focused on the analysis of power pulse tests performed on research reactors, the phase-II dealt with the benchmarking and comparison of participants’ modelling approaches: (task 1) to the analysis of several idealized RIA cases, and (task 2) with sensitivity analysis and uncertainty quantification (through the analysis of one particular RIA case). This paper describes the work performed by NINE as one of the participants to phase-II, and particularly on task 2. Such work involved the integrated use of two different codes, namely RELAP5 and TRANSURANUS for the thermal hydraulic and the fuel performance parts of the analysis respectively, the integration consisting in an off-line one-way sequential coupling. According to the benchmark specifications, the uncertainty analysis was performed by adopting a Monte Carlo statistical approach, requiring 200 runs with a number of uncertain input parameters randomly sampled from a specified Gaussian distribution, and by applying a 95/95 statistics to the selected output variables. Moreover, the sensitivity analysis was performed by calculating the Spearman’s rank correlation coefficients of specified output variables versus specified input parameters. The whole analysis process was handled through a Python-based “Statistical RELAP5-TRANSURANUS” (StaRT) coupled analysis programme developed on purpose. StaRT takes care of input parameter sampling, update of RELAP5 and TRANSURANUS input decks, post-processing of both codes’ results, uncertainty and sensitivity analyses. All the results produced at key steps of the process are stored in a single HDF5 file. StaRT turned out to be an effective tool for the purposes of the benchmark as well as for future applications to similar studies.
ATLAS A5.1 POST-TEST CALCULATION
Valeria Parrinello, Marco Cherubini
Nuclear and INdustrial Engineering (NINE), Borgo Giannotti, Lucca, Italy
17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-17)
Xi'an, China, 03 - 08 September 2017
Abstract — The ATLAS (Advanced Thermal-Hydraulic Test Loop for Accident Simulation) is a full pressure, full temperature integral test facility operated by KAERI (Korean Atomic Energy Research Institution). The experimental program to be conducted at the ATLAS facility in the framework of the OECD-ATLAS Project comprises about ten tests pertaining to different research topics.
Among these tests, the A5.1 experiment was performed targeting a 1% horizontal small break loss of coolant accident (SBLOCA) at the cold leg with cold leg injection of the emergency core cooling system (ECCS) and accident management (AM) action by means of secondary side depressurization. This test is the counterpart test of the SB-CL-32 test conducted using the Large Scale Test Facility (LSTF) of the Rig-of-Safety Assessment-V (ROSA-V) Program.
The two rigs are designed in accordance to different scaling approaches and the counterpart test design implied a preliminary scaling analysis for deriving the boundary conditions of the ATLAS A5.1 test by means of proper scaling parameters. The aims of ATLAS experiment are to demonstrate the progression of the cold leg line break accident and contribute to the understanding of the behavior of nuclear reactor systems and to the assessment and improvement of the existing best estimate codes.
The present paper deals with the thermal-hydraulic analysis carried out for the ATLAS A5.1 post-test calculation. The adopted system thermal-hydraulic code to perform the analysis is RELAP5Mod3.3.
A qualification procedure was pursued for developing the nodalization to be used for the analysis and demonstrating that the code results constitute a realistic approximation of the A5.1 test behavior. The bases of the procedure are embedded into the SCCRED (Standardized and Consolidated Calculated & Reference Experimental Database) methodology, a supporting tool for V&V and uncertainty evaluation of Best-Estimate system codes for licensing applications.
ATLAS A5.1 TEST BENCHMARK ACTIVITY
V. Parrinelloa, M. Cherubinia, K.Y. Choib
a Nuclear and INdustrial Engineering (NINE), Borgo Giannotti, Lucca, Italy
b Korean Atomic Energy Research Institute (KAERI), Republic of Korea
17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-17)
Xi'an, China, 03 - 08 September 2017
Abstract — The Advanced Thermal-Hydraulic Test Loop for Accident Simulation (ATLAS) was operated by Korea Atomic Energy Research Institute (KAERI) for accident simulations of advanced Pressurized Water Reactors (PWRs). In the framework of the OECD-ATLAS Project, the ATLAS A5.1 test was performed, targeting a 1% horizontal small break loss of coolant accident (SBLOCA) at the cold leg with cold leg injection of the emergency core cooling system (ECCS) and accident management (AM) action by means of secondary side depressurization. This test is the counterpart test of the SB-CL-32 test conducted using the Large Scale Test Facility (LSTF) of the Rig-of-Safety Assessment-V (ROSA-V) Program.
Fourteen institutions from ten different countries took part in the analytical benchmark on the ATLAS A5.1 test. The participants were asked to provide a total of 128 parameters and 54 time trends. This large amount of data was used to characterize the code models developed by the benchmark participants and to determine how accurately these computational models are able to represent the ATLAS behavior.
Qualitative and quantitative assessments are the two pillars of the qualification process adopted within the benchmark activity, whose bases are embedded into the SCCRED (Standardized and Consolidated Calculated & Reference Experimental Database) methodology typically adopted by NINE.
The comparison among participants’ data included a comparison of the developed nodalization features, the nodalization geometrical fidelity, and the simulation of the steady state and transient behavior of the facility. A quantitative evaluation of the results was performed too by the adoption of the FFTBM (Fast Fourier Transform Based Method).
The complex work performed by NINE for the analytical benchmark activity and the adopted qualification process will be described and the outcomes of the comparison between the experimental data and the calculation results submitted by the participants will be presented.
Defence in Depth and practical elimination of early and large releases – concepts and practice.
S. Michael Modroa , Tomislav Bajsb, Artur Lyubarskiyc
a Nuclear and INdustrial Engineering (NINE), Borgo Giannotti, Lucca, Italy
b ENCONET d.o.o., Miramarska 20, 10000 Zagreb, Croatia
c Atomenergoprojekt, Moscow, Podolskich Kursantov 1 st , Russian Federation
TopSafe 2017
Vienna, Austria, 12 - 16 February 2017
Abstract — One of the most important requirements of the International Atomic Energy Agency Safety Standard on design (SSR-2/1) is the requirement for practical elimination of accident sequences that would lead to large or early radioactive releases. This requirement is accompanied by requirements for extension of the plant design basis to include conditions considered earlier as “beyond design basis” which may include plant conditions with multiple failures of systems, including safety systems, without and with severe core damage. In this paper we are discussing the concept of practical elimination of large or early radioactive releases in the context of defence in depth and we are addressing some challenges for demonstration of this concept. The demonstration of the “practical elimination” is a complex process for which no specific guidance and criteria are available and not sufficient experience is so far generated. The “practical elimination” shall be considered as integral part of the defence in depth.