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.     

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.     

Implementation of control rod movement and boron injection options by using control variables in RELAP5/PARCS V2.7 coupled code

In present neutron kinetics codes, control rods banks do not have the possibility of dynamic movement during the simulation of a transient; besides it is necessary to send the boron concentration from the thermal-hydraulic code to the neutronic code to account for changes in cross-sections due to boron dilution. For instance, control rod movements are pre-programmed with simple instructions introduced before the beginning of the calculation. Hence, control rod positions are not related to the core characteristics and the control systems at any time of the simulation. This work presents the changes introduced in RELAP5/PARCS v2.7 codes to achieve that control rods and the boron injection become more dynamic and realistic components in such kind of simulators. Furthermore, in order to test the modifications introduced in both codes, it has been analyzed a boron injection transient in a typical PWR Nuclear Power Plant. The thermal-hydraulic model includes all the primary loop components of a PWR, the core fuel assemblies modeled with PIPE components, pumps, steam generators, pressurizer, etc. The neutronic representation of the reactor has been made in a one-to-one basis fuel channel model for the whole core.     

基于修改后RELAP5/MOD3.2的氟盐高温堆实验回路分析

大型热工流体整体效应系统实验(CIET)台架是为模拟氟盐冷却高温堆(FHR)热工水力响应而设计的实验回路,采用DOWTHERM A模拟氟盐作为冷却剂。通过在RELAP5/MOD3.2程序中加入DOWTHERM A物性参数以及传热关系式,计算FHR实验回路CIET在两种工况下的热工水力行为,并与实验结果进行对比,计算工况包括强迫循环条件与自然循环条件。计算结果表明:在强迫循环条件下,堆芯热量主要靠盘管式空气换热器(CTAH)排出,堆芯进出口冷却剂温度及CTAH出口冷却剂温度与实验值符合良好,CTAH进口冷却剂温度与实验值有些微偏差;在自然循环工况中,堆芯热量主要通过DHX与堆芯辅助冷却系统(DRACS)回路的换热带走,DHX及DRACS的流量与实验值接近,相对误差在10%左右,验证了修正后RELAP5/MOD3.2的正确性。     

Integral loop test and assessment of modified RELAP5/MOD3.3 for RCS coolant inventory during a LBLOCA

This study investigates experimentally and analytically the thermal hydraulic phenomena during the transition from design basis accident (DBA) to beyond-DBA, particularly, the depletion of core coolant inventory. To investigate the overall thermal hydraulic behavior after safety injection (SI) failure during a large-break loss-of-coolant accident (LBLOCA) in a cold leg, an integral loop test was performed at the Seoul National University Integral Test Facility (SNUF), which was scaled down to 1/6.4 in length and 1/178 in area from the advanced power reactor 1400 MWe (APR1400) according to the three-level scaling method. The plant condition at 200.0 s as the base case and those at 625.0 and 1950.0 s as test cases after the initiation of LBLOCA were applied as initial conditions in the experiments, respectively. The experimental results showed that the sweepout increased the coolant flow discharged to the break depending on the steam flow rate in intact cold legs and the initial downcomer coolant level and expedited the depletion of the core coolant inventory.In the meantime, since RELAP5/MOD3.3 uses the average properties of donor volume as those of its connected junction, this scheme causes the mass and the energy flux in a junction to be calculated incorrectly if significant phase separation occurs in the donor volume such as in the downcomer during the LBLOCA. The sweepout model was developed and implemented in RELAP5/MOD3.3 to improve its calculation of coolant inventory during the LBLOCA. To assess the applicability of the modified RELAP5/MOD3.3 to the actual system, the experiments in SNUF were simulated by both the original and the modified RELAP5/MOD3.3. The original one predicted the discharge flow rate at the break larger than that of the experiment. On the other hand, the modified one calculated the discharge flow rate more similar to that of the experiment than the original one did. As a result, the modified RELAP5/MOD3.3 reduced the coolant flow discharged to the break to delay the initiation time of heater heat-up after SI failure during LBLOCA in a cold leg. This improved RELAP5/MOD3.3 will support a more realistic thermal hydraulic analysis in an integrated analysis system.     

A separate effect test on the cooling behavior in a pressurized water reactor core under a low reflooding rate condition by using ATLAS and RELAP5/MOD3.3 code assessment

A separate effect test was performed on the cooling behavior in a PWR core under a low reflooding rate condition by using the ATLAS (Advanced Thermal–Hydraulic Test Loop for Accident Simulation) which is a thermal–hydraulic integral effect test facility for the pressurized water reactors APR1400 and OPR1000. Although several integral tests for the reflood phase of a large break loss of coolant accident (LBLOCA) of APR1400 have been performed with the ATLAS, the previous integral effect tests for the reflood phase of a LBLOCA are not easily simulated by existing codes, such as the RELAP5/MOD3, due to a unique phenomena in ATLAS, that resulted from an injection of large amount of subcooled water onto the heated wall of which temperature was higher than the target value.     

Analysis of different containment models for IRIS small break LOCA, using GOTHIC and RELAP5 codes

Advanced nuclear water reactors rely on containment behaviour in realization of some of their passive safety functions. Steam condensation on containment walls, where non-condensable gas effects are significant, is an important feature of the new passive containment concepts, like the AP600/1000 ones.In this work the international reactor innovative and secure (IRIS) was taken as reference, and the relevant condensation phenomena involved within its containment were investigated with different computational tools. In particular, IRIS containment response to a small break LOCA (SBLOCA) was calculated with GOTHIC and RELAP5 codes. A simplified model of IRIS containment drywell was implemented with RELAP5 according to a sliced approach, based on the two-pipe-with-junction concept, while it was addressed with GOTHIC using several modelling options, regarding both heat transfer correlations and volume and thermal structure nodalization. The influence on containment behaviour prediction was investigated in terms of drywell temperature and pressure response, heat transfer coefficient (HTC) and steam volume fraction distribution, and internal recirculating mass flow rate. The objective of the paper is to preliminarily compare the capability of the two codes in modelling of the same postulated accident, thus to check the results obtained with RELAP5, when applied in a situation not covered by its validation matrix (comprising SBLOCA and to some extent LBLOCA transients, but not explicitly the modelling of large dry containment volumes).The option to include or not droplets in fluid mass flow discharged to the containment was the most influencing parameter for GOTHIC simulations. Despite some drawbacks, due, e.g. to a marked overestimation of internal natural recirculation, RELAP5 confirmed its capability to satisfactorily model the basic processes in IRIS containment following SBLOCA.