Calculation of the VVER-1000 coolant transient benchmark using the coupled code systems DYN3D/RELAP5 and DYN3D/ATHLET

Plant-measured data provided by the OECD/NEA VVER-1000 coolant transient benchmark programme were used to validate the DYN3D/RELAP5 and DYN3D/ATHLET coupled code systems. Phase 1 of the benchmark (V1000CT-1) refers to an experiment that was conducted during the commissioning of the Kozloduy NPP Unit 6 in Bulgaria. In this experiment, the fourth main coolant pump was switched on whilst the remaining three were running normal operating conditions. The experiment was conducted at 27.5% of the nominal level of the reactor power. The transient is characterized by a rapid increase in the primary coolant flow through the core, and as a consequence, a decrease of the space-dependent core inlet temperature. The control rods were kept in their original positions during the entire transient. The coupled simulations performed on both DYN3D/RELAP5 and DYN3D/ATHLET were based on the same reactor model, including identical main coolant pump characteristics, boundary conditions, benchmark-specified nuclear data library and nearly identical nodalization schemes. In addition to validation of the coupled code systems against measured data, a code-to-code comparison between simulation results has also been performed to evaluate the respective thermal hydraulic models of the system codes RELAP5 and ATHLET.     

Analyses of the V1000CT-1 benchmark with the DYN3D/ATHLET and DYN3D/RELAP coupled code systems including a coolant mixing model validated against CFD calculations

Plant-measured data provided within the specification of the OECD/NEA VVER-1000 coolant transient benchmark (V1000CT) were used to validate the DYN3D/RELAP5 and DYN3D/ATHLET coupled code systems. Phase 1 of the benchmark (V1000CT-1) refers to the MCP (main coolant pump) switching on experiment conducted in the frame of the plant-commissioning activities at the Kozloduy NPP Unit 6 in Bulgaria. The experiment was started at the beginning of cycle (BOC) with average core expose of 30.7 effective full power days (EFPD), when the reactor power was at 27.5% of the nominal level and three out of four MCPs were operating. The transient is characterized by a rapid increase in the primary coolant flow through the core and, as a consequence, a decrease of the space-dependent core inlet temperature. Both DYN3D/RELAP5 and DYN3D/ATHLET analyses were based on the same reactor model, including identical MCP characteristics, boundary conditions, benchmark-specified nuclear data library and nearly identical nodalization schemes. For an adequate modelling of the redistribution of the coolant flow in the reactor pressure vessel during the transient a simplified mixing model for the DYN3D/ATHLET code was developed and validated against a computational fluid dynamics calculation.

The results of both coupled code calculations are in good agreement with the available experimental data. The discrepancies between experimental data and the results of both coupled code calculations do not exceed the accuracy of the measurement data. This concerns the initial steady-state data as well as the time histories during the transient. In addition to the validation of the coupled code systems against measured data, a code-to-code comparison between simulation results has been performed to evaluate relevant thermal hydraulic models of the system codes RELAP5 and ATHLET and to explain differences between the calculation results.     

Validation of coupled codes using VVER plant measurements

A data set of five transients at different VVER type nuclear power plants was collected in order to validate neutron kinetics/thermal hydraulics codes. Two of these transients ‘drop of control rod at nominal power at Bohunice-3’ of VVER-440 type and ‘coast-down of 1 from 3 working MCPs at Kozloduy-6’ of VVER-1000 type, were then utilised for code validation. Eight institutes contributed to the validation with 10 calculations using 5 different combinations of coupled codes. The thermal hydraulic codes were ATHLET, SMABRE and RELAP5 and the neutron kinetic codes DYN3D, HEXTRAN, KIKO3D and BIPR8. The general behaviour of both the transients was quite well calculated with all the codes. Even an elementary modelling of coolant mixing in reactor pressure vessel under asymmetric transients improved correspondence to the measurements. Some differences between the calculations seem to indicate that fuel modelling and treatment of VVER-440 control rods need further consideration. The simultaneous validation interacted with the data collection effort and thus improved its quality. The complexity of data collection systems and sometimes conflicting data, however, called for compromises and interpretation guides that also taught the analysts balanced plant modelling.     

Simulation of the VVER-1000 pump start-up experiment of the OECD V1000CT benchmark with CATHARE and TRAC-PF1

In the framework of joint effort between the Nuclear Energy Agency (NEA) of OECD, the United States Department of Energy (US DOE), and the Commissariat a l'Energie Atomique (CEA), France a coupled three-dimensional (3D) thermal-hydraulics/neutron kinetics benchmark for VVER-1000 was defined. The benchmark consists of calculation of a pump start-up experiment labelled V1000CT-1 (Phase 1), as well as a vessel mixing experiment and main steam line break (MSLB) transient labelled V1000CT-2 (Phase 2), respectively. The reference nuclear plant is Kozloduy-6 in Bulgaria. The overall objective is to assess computer codes used in the analysis of VVER-1000 reactivity transients. A specific objective is to assess the vessel mixing models used in system codes. Plant data are available for code validation consisting of one experiment of pump start-up (V1000CT-1) and one experiment of steam generator isolation (V1000CT-2). The validated codes can be used to calculate asymmetric MSLB transients involving similar mixing patterns. This paper summarizes a comparison of CATHARE and TRAC-PF1 system code results for V1000CT-1, Exercise 1, which is a full plant point kinetics simulation of a reactor coolant system (RCS) pump start-up experiment. The reference plant data include integral and sector average parameters. The comparison is made from the point of view of vessel mixing and full system simulation. CATHARE used a six-sector multiple 1D vessel thermal-hydraulic model with cross flows and TRAC used a six-sector, 18-channel coarse-mesh 3D vessel model. Good agreement in terms of integral parameters and inter-loop mixing is observed.     

Validation of coupled neutron kinetic/thermal–hydraulic codes. Part 1: Analysis of a VVER-1000 transient (Balakovo-4)

Three-dimensional hexagonal reactor dynamic codes have been developed for VVER type reactors and coupled with different thermal–hydraulic system codes. In the EU Phare project SRR1/95 these codes have been validated against real plant transients by the participants from several countries. Data measured during a test in the Balakovo-4 VVER-1000 have been analysed by coupled codes. In the test, one of two working feed water pumps of the steam generators was switched off at nominal power. The steady-state assembly powers measured before and after this transient are reproduced by the codes with a maximum deviation of about 5%. The time behaviour of the most safety-relevant parameters, such as total fission power, coolant temperatures and pressures is well modelled. Thermal–hydraulic feedback effects observed in the measurement are described by the codes in a consistent manner. The analyses have shown, that an accurate treatment of the heat transfer from the fuel rods to the coolant is important. In all, the results have increased the confidence in the coupled code analyses of VVER-1000 transients.     

Simulation of mixing effects in a VVER-1000 reactor

This work has been performed in the framework of the OECD/NEA thermalhydraulic benchmark V1000CT-2. This benchmark is related to fluid mixing in the reactor vessel during a MSLB accident scenario in a VVER-1000 reactor. Coolant mixing in a VVER-1000 V320 reactor was investigated in plant experiments during the commissioning of the Unit 6 of the Kozloduy nuclear power plant. Non-uniform and asymmetric loop flow mixing in the reactor vessel has been observed in the event of symmetric main coolant pump operation. For certain flow conditions, the experimental evidence of an azimuthal shift of the main loop flows with respect to the cold leg axes (swirl) was found.Such asymmetric flow distribution was analyzed with the Trio_U code. Trio_U is a CFD code developed by the CEA Grenoble, aimed to supply an efficient computational tool to simulate transient thermalhydraulic turbulent flows encountered in nuclear systems. For the presented study, a LES approach was used to simulate turbulent mixing. Therefore, a very precise tetrahedral mesh with more than 10 million control volumes has been created.The Trio_U calculation has correctly reproduced the measured rotation of the flow when the CAD data of the constructed reactor pressure vessel where used. This is also true for the comparison of cold leg to assembly mixing coefficients. Using the design data, the calculated swirl was significantly underestimated. Due to this result, it might be possible to improve with CFD calculations the lower plenum flow mixing matrices which are usually used in system codes.     

Impact of cross-section generation procedures on the simulation of the VVER-1000 pump startup experiment in the OECD/DOE/CEA V1000CT benchmark by coupled 3D thermal-hydraulics/neutron kinetics models

In the framework of joint effort between the Nuclear Energy Agency (NEA) of OECD, the United States Department of Energy (US DOE), and the Commissariat a l'Enerige Atomique (CEA), France a coupled three-dimensional (3D) thermal-hydraulics/neutron kinetics benchmark was defined. The overall objective of OECD/NEA V1000CT benchmark is to assess computer codes used in analysis of VVER-1000 reactivity transients where mixing phenomena (mass flow and temperature) in the reactor pressure vessel are complex. Original data from the Kozloduy-6 Nuclear Power Plant are available for the validation of computer codes: one experiment of pump start-up (V1000CT-1) and one experiment of steam generator isolation (V1000CT-2). Additional scenarios are defined for code-to-code comparison. As a 3D core model is necessary for a best-estimate computation of all the scenarios of the V1000CT benchmark, all participants were asked to develop their own core coupled 3D thermal-hydraulics/neutron kinetics models using the data available in the benchmark specifications and a common cross-section library. The first code-to-code comparisons based on the V1000CT-1 Exercise 2 specifications exhibited unacceptable discrepancies between two sets of results. The present paper focuses on the analysis of the observed discrepancies. The VVER-1000 3D neutron kinetics models are based on cross-section data homogenized on the assembly level. The cross-section library, provided as part of the benchmark specifications, thus consists in a set of parameterized two group cross sections representing the different assemblies and the reflectors. The origin of the observed large discrepancies was found mainly to lie in the methods used to solve the diffusion equation. The VVER reflector properties were also found to enhance discrepancies by increasing flux gradients at the core/reflector interface thus highlighting more the difficulties in some codes to handle high exponential flux gradients. This paper summarizes the different steps applied to analyze the neutronic codes and their predictions as well as the impact of cross-section generation procedures.     

NRI Řež solution of Exercise 3 of VVER-1000 Coolant Transient Benchmark – Phase 1 with coupled codes DYN3D–ATHLET

The paper presents a solution of VVER-1000 Coolant Transient Benchmark – Phase 1 (V1000CT-1) of Exercise 3 performed with the coupled reactor dynamic code DYN3D and system code ATHLET at NRI Řež. The first part of the paper contains brief characteristics of VVER-1000 NPP input deck and describes also the applied reactor core model. The second part introduces the steady-state results and important time dependencies, compared with experimental values. The calculation results show that such type of transient can be realistically described by the coupled codes DYN3D–ATHLET.     

铅铋自然循环热工水力分析程序开发

铅铋堆内冷却剂的自然循环对于反应堆的正常运行以及事故工况下的堆芯热量导出均至关重要,相关热工水力分析工作对于支持设计及安审均有重要意义。通过对铅铋堆内一回路系统内主要部件,包括堆芯、热交换器、管道等建立热工水力物理模型,开发了适用于铅铋自然循环瞬态过程模拟的热工水力分析程序,并利用铅铋自然循环回路内开展的自然循环启动实验、功率台阶影响实验等的结果进行了程序的初步验证。结果表明,程序计算得到的结果与实验结果符合较好,能够较好模拟铅铋自然循环的瞬态过程。该程序可以为铅铋堆研发过程中自然循环热工水力分析工作提供支持。     

Comparative study of a boron dilution scenario in VVER reactors

Subsequent studies have identified many scenarios, which can lead to reactivity excursions due to boron dilution. The comparative study, presented in this paper, deals with the so-called “restart of the first reactor coolant pump’’ scenario and its reactor-dynamic consequences for both Russian designed VVER reactor types, VVER-440 and VVER-1000. The transient simulations were performed using the three-dimensional core dynamics code DYN3D. The DYN3D modeling features, including recent developments, as well as the cross-section methodology involved in these calculations, are described. The analyzed accident scenario is outlined together with the assumptions made. The results of core response in this boron dilution accident for both VVER reactors are compared within the ranges, determined by the two reactivity values of interest: the criticality limit and the reactivity initiated accident (RIA) limit.     

Modelling effects and their impact on the results of ATHLET solutions of Exercise 1 of the VVER-1000 Coolant Transient Benchmark

The VVER-1000 Coolant Transient Benchmark consists of two phases and refers to experimental data from the Kozloduy Unit 6 Nuclear Power Plant in Bulgaria. The paper describes the modelling features and their impact on the results of the Exercise 1, Phase 1 of the Benchmark obtained by two ATHLET user groups, namely GRS and NRI. The simulated transient is a main coolant pump (MCP) switching on in one loop at reduced power while three other MCPs are in operation. Particular attention is paid to the influence of the reactor vessel modelling and especially of the nodalization in the upper plenum. The comparison and discussion of the two simulation results confirm that the two solutions with the ATHLET system code achieve quite good system response of the plant transient.     

Experimental investigation of in-vessel mixing phenomena in a VVER-1000 scaled test facility during unsteady asymmetric transients

In-vessel turbulent mixing phenomena affect the time and space distribution of coolant properties (e.g., boron concentration and temperature) at the core inlet which impacts consequently the neutron kinetics response. For reactor safety evaluation purposes and to characterize these phenomena it is necessary to set and validate appropriate numerical modelling tools to improve the current conservative predictions. With such purpose, an experimental campaign was carried out by OKB Gidropress, in the framework of the European Commission Project “TACIS R2.02/02 - Development of safety analysis capabilities for VVER-1000 transients involving spatial variations of coolant properties (temperature or boron concentration) at core inlet”. The experiments were conducted on a scaled facility representing the primary system of a VVER-1000 including a detailed model of the Reactor Pressure Vessel with its internals. The simulated transients involved perturbations of coolant properties distribution providing a wide validation matrix. The main achievements of the set of experiments featuring transient asymmetric pump behaviour are presented in this paper. The potential of the obtained experimental database for the validation of thermal fluid dynamics numerical simulation tools is also discussed and the role of computational fluid dynamics in supporting the experimental data analysis is highlighted.     

The 2nd International Topical Meeting on Neutron Radiography System Design and Characterization November 12—18, 1995, the Shonan Village Center,Hayama and the Rikkyo University,Yokosuka, Japan

Framatome-ANP has developed S-RELAP5, a RELAP5/Mod2 based thermal hydraulic code, and PANBOX, a 3D core kinetics code. By coupling both codes, a powerful neutronic and thermal hydraulic plant model was developed, which is capable of calculating extremely complex transients, particularly events bearing strongly asymmetric phenomena. The capability of the code system has been tested by recalculation of several transients that occurred in Siemens built PWRs. The most complex transient was a loss of load combined with a temporary coastdown of one main coolant pump, which is presented here. Since the measured values from the data recording system of the plant were available, the calculation could be compared to measured parameters.

The key phenomenon of the transient is a highly asymmetrical neutronic condition, which was caused by:

• —the drop of 5 control rod pairs in an asymmetric pattern upon detection of “loss of load.”

• —the coastdown of one main coolant pump, due to failure to connect to the auxiliary bus, which allowed coolant in one loop to stagnate and cool. Subsequent reactivation of that pump forced a plug of cold water into one side of the core.

The physical progress of the transient is strongly dependent on this double asymmetry; therefore, a 3-D calculation is indispensable for an accurate simulation. The calculated results are in good agreement with measurements and represent an important contribution to code validation.     

A VVER-1000 LEU and MOX assembly computational benchmark analysis using the lattice burnup code EXCEL

Utilization of Mixed Uranium–Plutonium Oxide (MOX) fuel in VVER-1000 reactors envisages the core physics analysis using computational methods and validation of the related computer codes. Towards this objective, an international experts group has been established at OECD/NEA. The experts group facilitates sharing of existing information on physics parameters and fuel behaviour. Several benchmark exercises have been proposed by them with intent to investigate the core physics behaviour of a VVER-1000 reactor loaded with 2/3rd of low enriched uranium (LEU) fuel assemblies (FA) and 1/3rd of weapons grade mixed oxide (MOX) FA. In the present study an attempt is made to analyse ‘AVVER-1000LEUandMOXAssemblyComputationalBenchmark’ and predict the neutronics behaviour at the lattice level. The lattice burnup code EXCEL, developed at Light Water Reactor Physics Section, BARC is employed for this task. The EXCEL code uses the 172 energy group ‘JEFF31GX’ cross-section library in WIMS-D format. Assembly level fuel depletion calculations are performed up to a burnup of 40 MWD/kg of heavy metal (HM). Studies are made for the parametric variations of fuel and moderator temperatures, coolant density and boron content in the coolant. Both operational and off-normal states are analysed to determine the corresponding infinite neutron multiplication factor (k_∞). Pin wise isotopic compositions are computed as a function of burnup. Isotopic compositions in different annular regions of Uranium–Gadolinium (UGD) pin, fission rate distributions in UGD, UO₂ and MOX pin cells are also computed. The predicted results are compared with the benchmark mean results.     

A subchannel analysis code MATRA-LMR for wire wrapped LMR subassembly

In sodium cooled liquid metal reactors design limits are imposed on the maximum temperatures of the cladding and fuel pins. Thus an accurate prediction of the core coolant/fuel temperature distribution is essential to LMR core thermal hydraulic design. The detailed subchannel thermal hydraulic analysis code MATRA-LMR is being developed for LMFBR core design and analysis based on COBRA-IV-I and MATRA. The major modifications and improvements implemented in MATRA-LMR are as follows: sodium property calculation subprogram, sodium coolant heat transfer correlations, and most recent pressure drop correlations. To assess the development status of this code, benchmark calculations were performed with the ORNL 19 pin tests and EBR-II seven-assembly SLTHEN calculation results. The calculation results of MATRA-LMR were compared to the measurements and to the SABRE4 and SLTHEN code calculation results, respectively. Finally, the major technical results of the conceptual design for the KALIMER U-10%Zr binary alloy fueled core have been compared with the calculations of the MATRA-LMR, SABRE4 and SLTHEN codes.     

Comprehensive uncertainty and sensitivity analysis for coupled code calculations of VVER plant transients

The development of coupled codes, combining thermal-hydraulic system codes and 3D neutron kinetics codes, is an important step to perform best-estimate calculations for plant transients of nuclear power plants. For applications in safety analysis, these coupled codes should be validated by benchmark calculations and, preferably, by comparison with plant transient data from operating plants. In addition, the results should be supplemented by applying uncertainty and sensitivity analysis methods, which allow to identify relevant parameters of models and solution procedures affecting the results and to quantify their relative importance. Both objectives were part of the VALCO project. The aspect of validation is presented in [S. Mittag, et al., 2004. Neutron-Kinetic Code Validation against Measurements in the Moscow V-1000 Zero-Power Facility, in press; T. Vanttola et al., 2004. Validation of coupled codes using VVER plant measurements, in press], the aspect of a comprehensive uncertainty and sensitivity analysis for coupled code calculations is the topic of this contribution. The results and experiences obtained by the analysis for two plant transients in a VVER-440 and a VVER-1000, respectively, are presented and discussed.     

Fluid mixing and flow distribution in a primary circuit of a nuclear pressurized water reactor—Validation of CFD codes

The EU project FLOMIX-R was aimed at describing the mixing phenomena relevant for both safety analysis, particularly in steam line break and boron dilution scenarios, and mixing phenomena of interest for economical operation and the structural integrity.This report will focus on the computational fluid dynamics (CFD) code validation. Best practice guidelines (BPG) were applied in all CFD work when choosing computational grid, time step, turbulence models, modelling of internal geometry, boundary conditions, numerical schemes and convergence criteria. The strategy of code validation based on the BPG and a matrix of CFD code validation calculations have been elaborated. CFD calculations have been accomplished for selected experiments with two different CFD codes (CFX, FLUENT). The matrix of benchmark cases contains slug mixing tests simulating the start-up of the first main circulation pump which have been performed with three 1:5 scaled facilities: the Rossendorf coolant mixing model ROCOM, the Vattenfall test facility and a metal mock-up of a VVER-1000 type reactor. Before studying mixing in transients, ROCOM test cases with steady-state flow conditions were considered. Considering buoyancy driven mixing, experimental results on mixing of fluids with density differences obtained at ROCOM and the FORTUM PTS test facility were compared with calculations. Methods for a quantitative comparison between the calculated and measured mixing scalar distributions have been elaborated and applied. Based on the “best practice CFD solutions”, conclusions on the applicability of CFD for turbulent mixing problems in PWR were drawn and recommendations on CFD modelling were given. The results of the CFD calculations are mostly in-between the uncertainty bands of the experiments. Although no fully grid-independent numerical solutions could be obtained, it can be concluded about the suitability of applying CFD methods in engineering applications for turbulent mixing in nuclear reactors.     

A spatial kinetic model for simulating VVER-1000 start-up transient

An accurate prediction of reactor core behavior in transients depends on how much it could be possible to exactly determine the thermal feedbacks of the core elements such as fuel, clad and coolant. In short time transients, results of these feedbacks directly affect the reactor power and determine the reactor response. Such transients are commonly happened during the start-up process which makes it necessary to carefully evaluate the detail of process. Hence this research evaluates a short time transient occurring during the start up of VVER-1000 reactor. The reactor power was tracked using the point kinetic equations from HZP state (100 W) to 612 kW. Final power (612 kW) was achieved by withdrawing control rods and resultant excess reactivity was set into dynamic equations to calculate the reactor power. Since reactivity is the most important part in the point kinetic equations, using a Lumped Parameter (LP) approximation, energy balance equations were solved in different zones of the core. After determining temperature and total reactivity related to feedbacks in each time step, the exact value of reactivity is obtained and is inserted into point kinetic equations. In reactor core each zone has a specific temperature and its corresponding thermal feedback. To decrease the effects of point kinetic approximations, these partial feedbacks in different zones are superposed to show an accurate model of reactor core dynamics. In this manner the reactor point kinetic can be extended to the whole reactor core which means “Reactor spatial kinetic”. All required group constants in calculations are prepared using the WIMS code. In addition CITATION code was used to calculate the flux, power distribution and core reactivity inside the core. To update the last change in group constants and resultant reactivity in point kinetic equations, these neutronic codes were coupled with a developed dynamic program. This study is applied on a typical VVER-1000 reactor core to show the reactor response in short time transients caused during start-up procedure.