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.     

Super-homogenisation factors in pinwise calculations by the reactor dynamics code DYN3D

Super-homogenisation (SPH) factors were generated by a modified version of the code DYN3D for PWR fuel assemblies in hot-zero-power states defined in the OECD MOX/UO2 Benchmark. SPH factors averaged for each pin-material type and factors for each individual pin position were produced. The application of the SPH factors improves the accuracy of DYN3D calculations, especially for configurations with control rods inserted.     

DYN3D-MSR spatial dynamics code for molten salt reactors

The development of spatial dynamics code for molten salt reactors (MSRs) is reported in this paper. The graphite-moderated channel type MSR – one of the ‘Generation IV’ concepts – was selected for the numerical simulation. It has several peculiarities (e.g. the drift of delayed neutrons precursors), which disable the use of standard dynamics codes. Therefore, the own DYN3D-MSR code was developed. It is based on the light water reactor code DYN3D and it allows transients simulation by 3D neutronics and parallel channel thermal-hydraulics. The neutronics and thermal-hydraulics were modified for the MSR peculiarities, where the experience from DYN1D-MSR development was exploited. The code was validated on experimental results from the MSRE experiment done in Oak Ridge National Laboratory and by the comparison with other codes especially with the 1D version. However, by the 3D code transients can be simulated, where space-dependant efforts are relevant, like local blockage of fuel channels or local temperature perturbations.     

Safety analysis of reactivity insertion accidents in a heavy water nuclear research reactor core using coupled 3D neutron kinetics thermal-hydraulic system code technique

Nuclear power plant Safety analysis using coupled 3D neutron kinetics/thermal-hydraulic codes technique is increasingly used nowadays. Actually, the use of this technique allows getting less conservatism and more realistic simulations of the physical phenomena. The challenge today is oriented toward the application of this technique to the operating conditions of nuclear research reactors. In the current study, a three-Dimensional Neutron Kinetics and best estimate Thermal-Hydraulic model based upon the coupled PARCS/RELAP5 codes has been developed and applied for a heavy water research reactor. The objective is to perform safety analysis related to design accidents of this reactor types. In the current study two positive reactivity insertion transients are considered, SCRAM protected and self-limiting power excursion cases. The results of the steady state calculations were compared with results obtained from conventional diffusion codes, while transient calculations were assessed using the point kinetic model of the RELAP5 code. Through this study, the applicability and the suitability of using the coupled code technique with respect to the classical models are emphasized and discussed.     

基于三维MC-SN耦合方法的PWR压力容器快中子注量计算基准分析

蒙特卡罗(MC)-离散纵标(SN)耦合方法是解决同时具有复杂几何和深穿透特点的核装置屏蔽问题的有效方法。本文首次将三维MC-SN耦合方法应用于压水堆屏蔽计算。针对NUREG/CR-6115压水堆基准模型,选取热屏蔽内表面为公共交界面,将其分为几何复杂的MC模拟区和具有深穿透特点的SN模拟区。三维MC程序用于精确描述堆芯到热屏蔽精细模型,并记录穿过热屏蔽内表面的中子径迹信息。接口程序将中子径迹转换为SN计算所需的边界源,提供给三维SN程序进行热屏蔽到压力容器的计算。计算结果包括压力容器内表面、1/4壁厚处及焊缝处快中子注量(E>1.0 MeV)圆周方向分布。三维耦合方法计算结果与基准报告提供的MCNP、DORT结果符合良好,验证了该方法处理圆柱坐标系屏蔽问题的有效性和程序使用的正确性。     

Analysis of the hot gas flow in the outlet plenum of the very high temperature reactor using coupled RELAP5-3D system code and a CFD code

The very high temperature reactor (VHTR) system behavior should be predicted during normal operating conditions and postulated accident conditions. The plant accident scenario and the passive safety behavior should be accurately predicted. Uncertainties in passive safety behavior could have large effects on the resulting system characteristics. Due to these performance issues in the VHTR, there is a need for development, testing and validation of design tools to demonstrate the feasibility of the design concepts and guide the improvement of the plant components. One of the identified design issues for the gas-cooled reactor is the thermal mixing of the coolant exiting the core into the outlet plenum. Incomplete thermal mixing may give rise to thermal stresses in the downstream components. To provide flow details, the analysis presented in this paper was performed by coupling a VHTR model generated in a thermal hydraulic systems code to a computational fluid dynamics (CFD) outlet plenum model. The outlet conditions obtained from the systems code VHTR model provide the inlet boundary conditions to the CFD outlet plenum model. By coupling the two codes in this manner, the important three-dimensional flow effects in the outlet plenum are well modeled while avoiding modeling the entire reactor with a computationally expensive CFD code. The values of pressure, mass flow rate and temperature across the coupled boundary showed differences of less than 5% in every location except for one channel. The coupling auxiliary program used in this analysis can be applied to many different cases requiring detailed three-dimensional modeling in a small portion of the domain.     

Evaluation of ASME code flaw analysis procedure using the influence function method for application to PWR primary piping

This paper discusses stress intensity factor calculations and fatigue analysis for a PWR primary coolant piping system. The influence function method is applied to evaluate ASME Code Section XI Appendix A “analysis of flaw indication” for the application to a PWR primary piping. Results of the analysis are discussed in detail.     

压水堆装置三维图形仿真

利用Ｉ－ＤＥＡＳＭａｓｔｅｒＳｅｒｉｅｓ５软件的“设计”、“图面绘制”、“仿真”、“测试”、“几何资源转译器”功能,在“压水堆装置三维图形仿真”任务中完成建模的工作。校核了设计数据的可靠性。     

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.     

Analysis of the Peach Bottom flow stability test number 3 using the coupled RELAP5/PARCS code

Nowadays, the coupled codes technique, which consists in incorporating three-dimensional (3D) neutron modeling of the reactor core into system codes, is extensively used for carrying out best estimate (BE) simulation of complex transient in nuclear power plants (NPP). This technique is particularly suitable for transients that involve core spatial asymmetric phenomena and strong feedback effects between core neutronics and reactor loop thermal-hydraulics. Such complex interactions are encountered under normal and abnormal operating conditions of a boiling water reactors (BWR). In such reactors Oscillations may take place owing to the dynamic behavior of the liquid-steam mixture used for removing the thermal power. Therefore, it is necessary to be able to detect in a reliable way these oscillations. The purpose of this work is to characterize one aspect of these unstable behaviors using the coupled codes technique. The evaluation is performed against Peach Bottom-2 low-flow stability tests number 3 using the coupled RELAP5/PARCS code. In this transient dynamically complex neutron kinetics coupling with thermal-hydraulics events take place in response to a core pressure perturbation. The calculated coupled code results are herein assessed and compared against the available experimental data.     

Uncertainty and sensitivity analyses of the Kozloduy pump trip test using coupled thermal–hydraulic 3D kinetics code

The modeling of complex transients in nuclear power plants (NPP) remains a challenging topic for best estimate three-dimensional coupled code computational tools. This technique is, nowadays, extensively used since it allows decreasing conservatism in the calculation models and performs more realistic simulation and more precise consideration of multidimensional effects under complex transients in NPPs. Therefore, large international activities are in progress aiming to assess the capabilities of coupled codes and the new frontiers for the nuclear technology that could be opened by this technique. In the current paper, a contribution to the assessment and validation of coupled code technique through the Kozloduy VVER100 pump trip test is performed. For this purpose, the coupled RELAP5/3.3-PARCS/2.6 code is used. The code results were assessed against experimental data. Deviations between code predictions and measurements are mainly due to the used models for evaluating and modeling of the Doppler feedback effect. Further investigations through the use of two “antagonist” uncertainty GRS and the CIAU methods, were considered in order to evaluate and quantify the origin of the observed discrepancies. It was revealed on one hand that relative error quantification discrepancies exist between the two approaches, and further enhancements for both methods are needed.     

Pre- and post-test calculations to natural circulation experiments at the integral test facility ISB-VVER using the thermalhydraulic code ATHLET

In 1995 at the integral test facility ISB-VVER in Elektrogorsk near Moscow natural circulation experiments were performed, which were scientifically accompanied by the Forschungszentrum Rossendorf. These experiments were the first of this kind at a test facility, which models VVER-1000 thermalhydraulics. Using the code ATHLET which is being developed by ‘Gesellschaft für Anlagen und Reaktorsicherheit’, pre- and post-test calculations were done to determine the thermalhydraulic events to be expected and to define and tune the boundary conditions of the test. The conditions found for natural circulation instabilities and cold leg loop seal clearing could be confirmed by the tests. Besides the thermalhydraulic standard measuring system, the facility was equipped with needle shaped conductivity probes for measuring the local void fractions.     

Analysis of the MT-3 clad ballooning reflood test using the multi-rod coupled MATARE code

If PWR fuel rods balloon in a co-located fashion following a large break LOCA, it is possible that an un-coolable region could be formed. The MT-3 test was one of the few large-scale experiments performed to investigate this. In this paper we report the analysis of this experiment using a novel combined thermal-hydraulic and structural mechanics model, able to model multiple distinct rods in a coupled fashion. The sensitivity of the behaviour to (in particular) the channel droplet distribution makes a full ‘first principles’ validation impractical, but good agreement is achieved between observed and predicted blockage fractions, and the observed incoherence of ballooning is reproduced well.     

Assessment of the interchannel mixing model with a subchannel analysis code for BWR and PWR conditions

The influence of the interchannel mixing model employed in a traditional subchannel analysis code was investigated in this study, specifically on the analysis of the enthalpy distribution and critical heat flux (CHF) in rod bundles in BWR and PWR conditions. The equal-volume-exchange turbulent mixing and void drift model (EVVD) was embodied to the COBRA-IV-I code. An optimized model of the void drift coefficient has been devised in this study as the result of the assessment with the two-phase flow distribution data for the general electric (GE) 9-rod and Ispra 16-rod test bundles. The influence of the subchannel analysis model on the analysis of CHF was examined by evaluating the CHF test data in rod bundles representing PWR and BWR conditions. The CHFR margins of typical light water nuclear reactor (LWR) cores were evaluated by considering the influence on the local parameter CHF correlation and the hot channel analysis result. It appeared that the interchannel mixing model has an important effect upon the analysis of CHFR margin for BWR conditions.     

Validation of 3D neutronic-thermalhydraulic coupled codes RELAP5/PARCSv2.7 and TRACEv5.0P3/PARCSv3.0 against a PWR control rod drop transient

In nuclear safety field, neutronic and thermalhydraulic codes performance is an important issue. New capabilities implementation, as well as models and tools improvements are a significant part of the community effort in looking for better nuclear power plants (NPP) designs. A procedure to analyze the PWR response to local deviations on neutronic or thermalhydraulic parameters is being developed. This procedure includes the simulation of Incore and Excore neutron flux detectors signals. A control rod drop real plant transient is used to validate the used codes and their new capabilities. Cross-section data are obtained by means of the SIMTAB methodology. Detailed thermalhydraulic models were developed: RELAP5 and TRACE models simulate three different azimuthal zones. Besides, TRACE model is performed with a fully three-dimensional core, thus, the cross-flow can be obtained. A Cartesian vessel represents the fuel assemblies and a cylindrical vessel the bypass and downcomer. Simulated detectors signals are obtained and compared with the real data collected during a control rod drop trial at a PWR NPP and also with data obtained with SIMULATE-3K code.     

Survey of experiments and code development for the passive residual heat removal system of PWR in China

Three different kinds of experiments and their typical results are surveyed for the passive residual heat removal system (PRHRS) of PWR performed in Nuclear Power Institute of China (NPIC) recent ten years. The typical results of MISAP. a special code for PWR passive residual heat removal system developed and assessed by NPIC,are also described briefly in this paper.     

Analysis of the VVER1000 coolant trip benchmark using the coupled RELAP5/PARCS code

Incorporation of full three-dimensional models of the reactor core into system thermal–hydraulic transient codes allows better estimation of interactions between the core behavior and plant dynamics. Considerable efforts have been made in various countries and organizations to verify and validate the capability of the so-called coupled codes technique. For these purposes appropriate Light Water Reactor (LWR) transient benchmarks based upon programmed transients performed in Nuclear Power Plants (NPP) were recently developed on a higher ‘best-estimate’ level. The reference problem considered in the current framework is a Main Coolant Pump (MCP) switching-on transient in a VVER1000 NPP. This event is characterized by a positive reactivity addition as consequence of the increase of the core flow. In the current study the coupled RELAP5/PARCS code is used to reproduce the considered test. Results of calculation were assessed against experimental data and also through the code-to-code comparison.