Validation of coupled neutronic/thermal-hydraulic code RELAP5-3D for RBMK-1500 reactor analysis application

This paper summarizes RELAP5-3D code validation activities carried out at the Lithuanian Energy Institute, which was performed through the modeling of RBMK-1500 specific transients taking place at Ignalina NPP. A best estimate RELAP5-3D model of the INPP RBMK-1500 reactor has been developed and validated against real plant data, as well as with the calculation results obtained using the Russian STEPAN/KOBRA code. The obtained calculation results demonstrate reasonable agreement with Ignalina NPP measured data. Behaviors of the separate MCC thermal-hydraulic parameters, as well as physical processes are predicted reasonably well to the real processes, occurring in the primary circuit of RBMK-1500 reactor. The calculated reactivity and the total reactor core power behavior in time are also in reasonable agreement with the measured plant data, which demonstrates the correct modeling of the neutronic processes taking place in RBMK-1500 reactor core. The performed validation of RELAP5-3D model of Ignalina NPP RBMK-1500 reactor allowed to improve the model, which in the future would be used for the safety substantiation calculations of RBMK-1500 reactors. Future activities are discussed.     

Modeling of void, fast power and graphite temperature reactivity coefficients measurements for the validation of RELAP5-3D RBMK-1500 reactor model

This paper deals with the modeling of RBMK-1500 specific transients taking place at Ignalina NPP: measurements of void and fast power reactivity coefficients, as well as change of graphite cooling conditions transient. The simulation of these transients was performed using RELAP5-3D code model of RBMK-1500 reactor. At the Ignalina NPP void and fast power reactivity coefficients are measured on a regular basis and based on the obtained experimental results the actual values of these reactivity coefficients are determined. Graphite temperature reactivity coefficient at the plant is determined by changing graphite cooling conditions in the reactor cavity. This type of transient is unique and important from the point of view of model validation for the gap between fuel channel and the graphite bricks. The measurement results, obtained during this transient, enabled to determine the thermal conductivity coefficient for this gap and to validate the graphite temperature reactivity feedback model. The performed validation of RELAP5-3D model of Ignalina NPP RBMK-1500 reactor allowed to improve the model, which in the future would be used for the safety substantiation calculations of RBMK-1500 reactors.     

Analysis of the processes in spent fuel pools of Ignalina NPP in case of loss of heat removal

Ignalina NPP is the only nuclear power plant in Lithuania consisting of two units, commissioned in 1983 and 1987. Unit 1 of Ignalina NPP was shutdown for decommissioning at the end of 2004 and Unit 2 is to be operated until the end of 2009. Both units are equipped with channel-type graphite-moderated boiling water reactors RBMK-1500. According to the design, the spent fuel should be returned for reprocessing to Russia. However actually any fuel assembly has not been taken out from territory of the Ignalina NPP and all assemblies of spent fuel are stored in the spent fuel pools and dry on-site storage facility. Thus, the safety justification of facilities for intermediate spent fuel assemblies’ storage in Ignalina NPP is very important. This paper presents the results of loss of heat removal accidents (the most probable beyond design basis accident) in spent fuel pools of Ignalina NPP. The analysis was performed by employing best-estimate system thermal hydraulic code RELAP5 and codes for severe accidents ATHLET-CD and ASTEC. The best-estimate analysis, performed using RELAP5, allowed to investigate in the details the water evaporation, uncovering and fuel assemblies heat-up processes, when heat removal from the structures of buildings and pools are evaluated. The processes of spent fuel assemblies’ degradation due to loss of long-term heat removal were analyzed using ATHLET-CD and ASTEC codes. The results of calculations showed that the increase in water temperature in the pools from 50 °C up to 100 °C takes approximately 80-110 h, the evaporation of water volume down to uncovering of fuel assemblies takes approximately 220-260 additional hours. Later, after 200-300 h, the temperature of fuel claddings exceeds 800-1000 °C and the failures of fuel claddings occur due to cladding ballooning. The total amount of hydrogen generated up to time of complete water evaporation from spent fuel pools is about 7500-16,000 kg. These results of performed analysis were used for development of accident management guidelines for spent fuel pools of RBMK-1500.     

Benchmark analysis of main circulation pump trip events at the Ignalina NPP using RELAP5 code

The state-of-the-art code RELAP5/MOD3 was originally designed for PWRs. Because of unique RBMK designs the application of this code to RBMK-1500 encountered several problems. A successful best estimate RELAP5 model of the Ignalina NPP has been developed. This model includes the reactor main circulation circuit (MCC) and reactor control and protection system required for this kind of transient analysis. Benchmark analysis of all operating main circulation pump (MCP) trip events was performed. During the analysis the characteristics of isolation control valves and MCP throttling regulating valves were established. Comparison of calculated and measured parameters was also used to establish realistic resistances of different MCC components and realistic behaviour of the controllers of the reactor systems. Calculations performed with the RELAP5 model, which includes these modifications, compare favourably with plant data.     

Analysis of potential waterhammer at the Ignalina NPP using thermal–hydraulic and structural analysis codes

The RBMK (Russian acronym for ‘channeled large power reactor’)-1500 reactors at the Ignalina nuclear power plant (NPP) have a series of check valves in the main circulation circuit (MCC) that serve the coolant distribution in the fuel channels. In the case of a hypothetical guillotine break of pipelines upstream of the group distribution headers (GDH), the check valves and adjusted piping integrity is a key issue for the reactor safety during the rapid closure of check valve. An analysis of the waterhammer effect (i.e. the pressure pulse generated by the valves slamming closed) is needed. The thermal–hydraulic and structural analysis of waterhammer effects following the guillotine break of pipelines at the Ignalina NPP with RBMK-1500 reactors was conducted by employing the RELAP5 and PipePlus codes. Results of the analysis demonstrated that the maximum values of the pressure pulses generated by the check valve closure following the hypothetical accidents remain far below the value of pressure of the hydraulic tests, which are performed at the NPP and the risk of failure of the check valves or associated pipelines is low. Sensitivity analysis of pressure pulse dependencies on calculation time step and check valve closure time was performed. Results of RELAP5 calculations are benchmarked against waterhammer transient data obtained by employing structural mechanics code BOS fluids.     

Uncertainty and sensitivity analysis of MCPs’ trip events at Ignalina NPP

Eight main circulation pumps (MCPs) are employed for the cooling water forced circulation through the RBMK-1500 reactor at the Ignalina nuclear power plant (NPP). There have been a few events when one or more MCPs were inadvertently tripped.This paper presents investigation of a one MCP trip event and all MCPs’ trip events at Ignalina NPP. Thermal-hydraulic analysis was conducted using the best estimate system code RELAP5/MOD3.3. Uncertainty and sensitivity analysis of flow energy loss in different parts of the main circulation circuit (MCC), initial conditions and code-selected models was performed. Such analysis allows to estimate the influence of separate parameters on the calculation results and find those modelling parameters that have the largest impact on the investigated events. Uncertainty analysis indicates that natural circulation provides adequate cooling in the case of all MCPs tripped, and that the reactor is reliably cooled by forced circulation in the case of a single tripped MCP. On the basis of this analysis, recommendations for the further improvement of model are developed.     

RELAP5-3D code validation in the neutron-dynamic analysis of transient processes taking place in RBMK-1500 reactors

The paper presents an evaluation of RELAP5-3D code suitability to model-specific transients that take place during RBMK-1500 reactor operation, where the neutronic response of the core is important. Certain RELAP5-3D transient calculation results were benchmarked against calculation results obtained using the Russian complex neutronic-thermal-hydraulic code STEPAN/KOBRA, specially designed for RBMK reactor analysis. Comparison of the results obtained, using the RELAP5-3D and STEPAN/KOBRA codes, showed reasonable mutual agreement of the calculation results of both codes and their reasonable agreement with the real plant data.     

Comparison of RELAP5 calculations of VVER-1000 coolant transient benchmark phase 1 at different power

This paper provides comparisons between experimental data of “MCP switching on when the other three MCPs are in operation” and RELAP5 calculations with different initial levels of the reactor power 29.45% and 27.47% from the nominal.

The reference power plant for this analysis is Unit 6 at the Kozloduy nuclear power plant (NPP) site. RELAP5/MOD3.2 computer code has been used to simulate the investigated transient. Operational data from Kozloduy NPP have been used for the purpose of assessing how the RELAP5 model compares against plant data. During the plant-commissioning phase at Kozloduy NPP Unit 6 a number of experiments have been performed. One of them is switching on MCP when the other three MCPs are in operation.

This investigation was possible through the participation of leading specialists from Kozloduy NPP and with the support of Argonne National Laboratory, under the International Nuclear Safety Program (INSP) of the United States Department of Energy.     

Development of accident management measures for RBMK-1500 in the case of loss of long-term core cooling

Results of the Level 1 Probabilistic Safety Assessment of the Ignalina Nuclear Power Plant have shown that in the risk topography transients are dominating. Analysis has shown that failure of the long-term core cooling is the main contributor to the core damage frequency. However, the reactor core damage in the long-term indicates the potential opportunities for the accident management. The main goal of accident management is to avoid a multiple fuel channel rupture because considering the design of RBMK reactors the consequences of rupture of more than 11–16 FC at full pressure would be close to the consequences of Chernobyl accident. This paper presents a detailed thermal–hydraulic analysis of the accidents with the loss of long-term core cooling, performed using the RELAP5 model of Ignalina NPP reactor cooling circuit and safety systems. Different ways of potential accident management are discussed. On the basis of this analysis the accident management strategy was developed.     

RELAP5/MOD3.2 sensitivity calculations of loss-of-feed water (LOFW) transient at Unit 6 of Kozloduy NPP

This paper provides a comparison between the real plant data obtained by Unit 6 of Kozloduy nuclear power plant (NPP) during the loss-of-feed water (LOFW) transient and the calculation results received by RELAP5/MOD3.2 computer model of the same NPP unit.RELAP5/MOD3.2 computer model of the VVER-1000 has been developed at the Institute for Nuclear Research and Nuclear Energy-Bulgarian Academy of Sciences (INRNE-BAS) based on Unit 6 of Kozloduy NPP. This model has been used for simulation the behavior of the real VVER-1000 NPP during the LOFW transient. Several calculations have been provided to describe how the different boundary conditions reflect on the prediction of real plant parameters.This paper discusses the results of the thermal–hydraulic sensitivity calculations of loss-of-feed water transient for VVER-1000 reactor design. The report also contains a brief summary of the main NPP systems included in the RELAP5 VVER model and the LOFW transient sequences.This report was possible through the participation of leading specialists from Kozloduy NPP and with the assistance of Argonne National Laboratory (ANL) for the United States Department of Energy (US DOE), International Nuclear Safety Program (INSP).     

Benchmarking analysis of water hammer effects using RELAP5 code and development of RBMK-1500 reactor main circulation circuit model

Correct prediction of water hammer transients is of paramount importance for the safe operation of the plant. Therefore, verification of computer codes capability to simulate water hammer type transients is a very important issue at performing safety analyses for nuclear power plants. Verification of RELAP5/MOD3.3 code capability to simulate water hammer type transients employing the experimental investigations is presented. Experience gained from benchmarking analyses has been used at development of the detail RELAP5 code RBMK-1500 model for simulation of water hammer effects in reactor main circulation circuit. In RBMK-type reactors the water hammers can occur in cases of rapid check valve operation. The performed analysis using RELAP5 code RBMK-1500 model has shown that in general the maximum values of the pressure pulses due to water hammer do not exceed the permissible loads on the pipelines.     

Modelling of a supercritical CO2 power cycle for nuclear fusion reactors using RELAP5–3D

A supercritical recompression CO₂ power cycle has been simulated using the system code RELAP5–3D. This code is being developed by INL and has traditionally been used in the simulation of operational and accidental transients in fission nuclear plants. The aim of the work presented here, developed within the framework of the Spanish Fusion Technology Program Consolider TECNO_FUS, is to take advantage of the simulation capabilities of RELAP5–3D in a field where little if any experience exists in the use of the code; i.e., the simulation of the heat fluxes and the thermodynamic cycle that, in a fusion power plant, will convert thermal power from plasma into mechanical power as a previous step to electricity generation. Code capabilities that make it suitable for this purpose are, for instance, the compressor model and the libraries of fluid properties (among them CO₂ and LiPb).The reference plant for the simulation is the one being designed under TECNO_FUS, which is the Spanish proposal for DEMO. The model of the plant includes the primary coolant systems, i.e. helium and LiPb in the Spanish dual coolant modular design (doble refrigerante modular, DRM), compressors, turbine and heat exchangers (Printed Circuit type).After the model has been set-up, several steady-state calculations have been run to test the performance of the model. After designing some minimal control features and adjusting their parameters, a few transient calculations have been run in order to demonstrate the capabilities of code and model. Finally, strengths and weaknesses of code and model are highlighted, along with some conclusions on their suitability for fusion technology applications.     

Specifics of RBMK core cooling in beyond design basis accidents

The most dangerous beyond design basis accidents for RBMK reactors, leading to the worst consequences, are related to the loss of long-term heat removal from the core. Due to a specific design of RBMK, there are a few possibilities for heat removal from reactor core by non-regular means: removal of heat from graphite stack by reactor gas circuit, removal of heat from reactor core using control rods cooling circuit, depressurisation of reactor cooling system, supply of water into cooling system from low pressure water sources, etc. This paper presents the analysis of such heat removal by employing RELAP5, RELAP5-3D and RELAP/SCDAPSIM codes. The analysis was performed for Ignalina nuclear power plant with RBMK-1500 reactor. The analysis of result shows that the restoration of water supply into control rod channels enables to remove 10-30 MW of the generated heat from the reactor core. This amount of removed heat is comparable with reactor decay heat in long-term period and allows to slowdown the core heat-up process. However, the injection of water to reactor cooling system is considered as main strategy, which should be considered in RBMK-1500 accident management procedure.     

RELAP5/MOD3.2 investigation of a VVER-1000 MCP switching on Benchmark problem

This paper provides comparisons between experimental data of Kozloduy NPP “MCP switching on when the other three MCP are in operation”, with Relap5 calculations. The investigated thermal-hydraulic driven transient is characterized by spatially dependant non-symmetric processes. RELAP5/MOD3.2 computer code has been used to simulate the investigated transient. Operational data from Kozloduy NPP have been used for the purpose of assessing how the RELAP5 model compares against plant data. During the plant-commissioning phase at Kozloduy NPP Unit 6 a number of experiments have been performed. One of them is switching on MCP when the other three MCPs are in operation. The event is characterized by rapid increase in the flow through the core resulting in a coolant temperature decrease, which leads to insertion of positive reactivity due to the modeled feedback mechanisms. The main purpose of this investigation was to improve the discrepancy between the calculations and the plant data. The sensitivity calculation investigates the mixing in reactor vessel and influence of heat structure on the hot legs temperature. The areas of improvements to the Relap5 model are:

• •The non-symmetrical mixing in downcomer and reactor vessel annular exit.
• •The influence of heat structure temperature on the time delay for equipments measurements.
• •Investigation of pressurizer water level – using the hot legs temperature correction.

The RELAP5/MOD3.2 model of Kozloduy NPP VVER-1000 for investigation of operational occurrences, abnormal events, and design basis scenarios have been developed and validated in the Institute for Nuclear Research and Nuclear Energy – Bulgarian Academy of Sciences (INRNE-BAS) Sofia, and Kozloduy NPP. The model provides a significant analytical capability for the specialists working in the field of NPP safety.This investigation is a process that compares the analytical results obtained by the RELAP5 computer model of the VVER-1000 against the experimental transient data received from the Kozloduy NPP Unit 6. The comparisons indicate good agreement between the RELAP5 results and the experimental data. The sensitivity investigation improves the discrepancy between the calculation and the plant data.This investigation was possible through the participation of leading specialists from Kozloduy NPP and with the support of Argonne National Laboratory, under the International Nuclear Safety Program (INSP) of the United States Department of Energy.     

压水堆核电站堆芯物理/热工水力耦合特性研究

采用RELAP5-HD作为堆芯耦合计算程序,以秦山核电二期工程反应堆堆芯为研究对象,建立堆芯活性区的物理/热工水力耦合模型,在此基础上进行了稳态计算和掉棒事故仿真研究。结果表明,使用RELAP5-HD计算得到的结果与电厂实测值符合较好,获得的掉棒事故参数曲线能准确反映事故工况下的参数变化趋势。稳态和事故工况的计算结果均符合堆芯物理/热工水力反馈效应的理论分析,证实了所建立的堆芯耦合模型的准确性,为下一步进行核电站系统的仿真分析提供基础。     

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.     

用Simulink扩展RELAP5的控制与保护系统仿真功能

RELAP5程序本身具有模拟核电站控制与保护系统的功能，但是，由于该程序采用文本输入方式进行建模，编写复杂，可读性不强，小适合于对大型复杂控制系统进行仿真。而Simulink程序采用图形化建模方式，能够高效、便捷地对核电厂复杂控制与保护系统进行建模。因此，本文将RELAP5程序与Simulink耦合，并利用Simulink扩展RELAP5的控制与保护系统的模拟功能。为了验证两程序耦合方法的准确性，将用Simulink实现的控制与保护系统的仿真结果，与已通过验证的RELAP5实现的具有相同功能的控制与保护系统的仿真结果进行对比，结果表明二者符合较好。     

核动力装置自然循环及其过渡过程计算模型的建立

为准确分析含反应性反馈的核动力装置自然循环及其过渡过程中重要参数的响应特性,以核动力装置瞬态最佳估算程序RELAP5/MOD3为基础,采用两群三维时空中子动力学模型替代RELAP5/MOD3的点堆模型,并建立三维空间内中子物理与热工水力的耦合模型,编制相应的计算程序。利用所研制的程序对实际核动力装置的自然循环及其过渡过程进行分析计算,并与试验结果进行比较。结果表明:本文建立的时空中子动力学计算模型克服了点堆方程不能准确计算反应性反馈的缺点,计算精度高,研制的程序可作为核动力装置强迫循环与自然循环及其过渡过程的计算分析工具。     

Thermal-hydraulic analysis of accidents leading to local coolant flow decrease in the main circulation circuit of RBMK-1500

There are a few transient and loss-of-coolant accident conditions in RBMK-1500 reactors that lead to a local flow decrease in fuel channels. Because the coolant flow decreases in fuel channels (FC) leads to overheating of fuel claddings and pressure tube walls, mitigation measures are necessary. The accident analysis enabled the suggestion of the new early reactor scram actuation and emergency core cooling system (ECCS) initiation signal, which ensures the safe shutdown of the reactor and compensates the stagnation flow. Analysis of such conditions is presented in this paper. Thermal-hydraulic analysis was conducted using the state-of-the-art RELAP5 code. Results of the analysis demonstrated that, after implementation of the developed management strategy for destruction of local flow stagnation, the Ignalina nuclear power plant (NPP) would be adequately protected following accidents, leading to local coolant flow decrease in the primary circuit.     

Reactor cavity and ALS thermal-hydraulic evaluation in the case of fuel channels ruptures at Ignalina NPP

This paper presents the work analysis of the thermal-hydraulic parameters behavior in the RBMK-1500 reactor cavity (RC) and other connected volumes in the case of fuel channels ruptures. The analysis is performed with CONTAIN code using the models of accident localization system (ALS) and reactor cavity venting system (RCVS). The RCVS capacity is assessed and expressed as a number of ruptured fuel channels at which the integrity of RC is maintained. The uncertainty analysis of pressure behavior in RC during multiple fuel channel rupture was performed. The initial and boundary conditions and the code models were selected and their influence on the results is estimated.Calculation of coolant mass and energy release to the reactor cavity in case of fuel channels rupture performed using the main circulation circuit model of Ignalina NPP, which was developed by employing state-of-the-art code RELAP5/MOD3.2 [Fletcher et al., RELAP5/MOD3 code manual user’s guidelines, Idaho National Engineering Lab., NUREG/CR-5535 (1992)]. These results were applied further as the initial data for the analysis of the thermal-hydraulic parameters behavior in the affected compartments employing CONTAIN code.