Assessment of RELAP5 point kinetic model against reactivity insertion transient in the IAEA 10 MW MTR research reactor

The current study emphasizes an aspect related to the assessment of a model embedded in a computer code. The study concerns more particularly the point neutron kinetics model of the RELAP5/Mod3 code which is worldwide used. The model is assessed against positive reactivity insertion transient taking into account calculations involving thermal-hydraulic feedback as well as transients with no feedback effects. It was concluded that the RELAP5 point kinetics model provides unphysical power evolution trends due most probably to a bug during the programming process.     

Safety analysis of the IAEA reference research reactor during loss of flow accident using the code MERSAT

Using the thermal hydraulic code MERSAT detailed model including primary and secondary loop was developed for the IAEA's reference research reactor MTR 10 MW. The developed model enables the simulation of expected neutronic and thermal hydraulic phenomena during normal operation, reactivity and loss of flow accidents.Two different loss of flow accident (LOFA) have been simulated using slow and fast decrease time of core mass flow. In both cases the expected flow reversal from downward forced to upward natural circulation has been successfully simulated. The results indicate that in both accidents the limit of onset of subcooled boiling was not arrived and consequently no exceed of design limits in term of thermal hydraulic instability or DNB is observed. Finally, the simulation results show good agreement with previous international benchmark analyses accomplished with other qualified channel and thermal hydraulic system codes.     

Passive cooldown performance of a 65 MW integral reactor

Advanced integral-type pressurized water reactor with a maximum thermal power of 65 MW is under development at the Korea Atomic Energy Research Institute (KAERI). This 65 MW integral reactor incorporates a number of innovative design features. In the case of a transient, the passive residual heat removal system (PRHRS) is designed to cool the reactor coolant system (RCS) from a normal operation condition to a hot shutdown condition by a natural circulation, and the shutdown cooling system (SCS) is designed to cool the primary system from a hot shutdown condition to a refueling condition by a forced circulation. A realistic calculation has been carried out by using the TASS/SMR code and a sensitivity analysis has been performed to evaluate a passive cooldown capability for various system conditions such as natural and forced circulation conditions for the reactor coolant system or the passive residual heat removal system, and number of active PRHRS trains. The reactor coolant system and the passive residual heat removal system adequately remove the core decay heat by a natural circulation and the 65 MW integral reactor can cool the coolant to the SCS entry condition in the primary system for all the possible operational conditions.     

Analysis of loss of coolant accident in MTR pool type research reactor

In MTR research reactors, heat removal is, safely performed by forced convection during normal operation and by natural convection after reactor shutdown for residual decay heat removal. However, according to the duration time of operation at full power, it may be required to maintain the forced convection, for a certain period of time after the reactor shutdown. This is among the general requirements for the overall safety engineering features of MTR research reactors to ensure a safe residual heat removal. For instance, in safety analysis of research reactors, initiating events that may challenge the safe removal of residual heat must be identified and analyzed.In the present work, it was assumed a total loss of coolant accident in a typical MTR nuclear research reactor with the objective of examining the core behavior and the occurrence of any fuel damage.For this purpose, the IAEA 10 MW benchmark core, which is a representative of medium power pool type MTR research reactors, was chosen herein in order to investigate the evolution of cladding temperature through the use of a best estimate thermalhydraulic system code RELAP5/mod3.2.     

Experimental simulation study on start-up of the 5 MW nuclear heating reactor

An experimental simulation study on the start-up of a low temperature, natural circulation nuclear heating reactor (5 MW developed by the Institute of Nuclear Energy of Tsinghua University, Beijing) is presented. The experiment was performed on the test loop (HRTL-5), which simulates the geometry and system design of the 5 MW reactor. The manifestation of different kinds of two-phase flow instability, namely geysering, flashing instability and low steam quality density wave instability on the start-up are described. The mechanism of flashing instability, which has never been well studied in this field, is especially interpreted. Based on the study of these instabilities, it is suggested that the start-up process, from initial condition to boiling operation condition, should consist of three steps: increasing of initial pressure by means of a noncondensable gas (N₂), start-up of the reactor at this pressurized condition (single-phase regime operation), and transition to a lower pressure, boiling operation. Three transition methods are discussed. As a result of these studies, the method of transition with low heat flux and low inlet subcooling is proposed. A stable start-up process of the 5MW reactor is achieved by careful selection of the thermohydraulic parameters.     

10MW高温气冷堆的燃耗测量研究

10MW高温气冷堆的燃耗测量系统是采用非破坏性高纯锗γ谱仪在线监测来确定燃耗值,利用MCNP4A程序对测量系统的衰减因子进行计算,基于核燃料裂变核素的γ射线能谱分析,以137Cs和134Cs核素活度作为测量对象,并对燃耗测量结果进行讨论.     

Xenon transient simulation of the VVER-1000 nuclear reactor using adiabatic approximation

The neutronic behavior of very slow transients like fuel burn up and xenon studies can be performed with the sequence of instantaneous criticality calculations. Such a scheme is known as the adiabatic approximation. ¹³⁵Xe as a fission product has an enormous thermal absorption cross section, on the order of a million barns, therefore the study of xenon poisoning and its effect on flux and power distribution is very important in thermal reactors. In this work xenon transient analysis of the VVER-1000 nuclear reactor and its effect on the flux and power distribution from reactor start up to xenon saturation and the change of power from nominal to 25% of nominal is carried out using WIMS and CITATION codes. We used the WIMS cell calculation code and found some relations between xenon concentration and group constants of different FA (Fuel Assemblies); in this way we bypassed the WIMS running at each time step. Also, the CITATION code is used as a core calculation code to find the effective multiplication factor as well as flux and power distributions. In order to link WIMS and CITATION codes and facilitate numerous executions, a VISUAL BASIC program has been developed. The results have a good agreement with the safety analysis report of the reference plant such that the relative differences in most cases are less than 10%.     

Numerical simulation using CFD software of countercurrent gas-liquid flow in a PWR hot leg under reflux condensation

In order to improve the countercurrent flow model of a transient analysis code, countercurrent air-water tests were previously conducted using a 1/15 scale model of the PWR hot leg and numerical simulations of the tests were carried out using the two-fluid model implemented in the CFD software FLUENT 6.3.26. The predicted flow patterns and CCFL characteristics agreed well with the experimental data. However, the validation of the interfacial drag correlation used in the two-fluid model was still insufficient, especially regarding the applicability to actual PWR conditions. In this study, we measured water levels and wave heights in the 1/15 scale setup to understand the characteristics of the interfacial drag, and we considered a relationship between the wave height and the interfacial drag coefficient. Numerical simulations to examine the effects of cell size and interfacial drag correlations on numerical predictions were conducted under PWR plant conditions. Wave heights strongly related with the water level and interfacial drag coefficient, which indicates that the interfacial drag force mainly consists of form drag. The cell size affected the gas velocity at the onset of flooding in the process of increasing gas flow rate. The gas volumetric fluxes at CCFL predicted using fine cells were higher than those using normal cells. On the other hand, the cell size did not have a significant influence on the process of decreasing gas flow rate. The predictions for the PWR condition using a reference set of interfacial drag correlations agreed well with the Upper Plenum Test Facility data of the PWR scale experiment in the region of medium gas volumetric fluxes. The reference interfacial drag correlations employed in this study can be applied to the PWR conditions.     

Analysis of partial and total flow blockage of a single fuel assembly of an MTR research reactor core

The main aim of the following study is to perform a safety analysis of the IAEA 10 MW MTR Pool type Research Reactor [IAEA-TECDOC-233, 1980. IAEA Research Reactor Core Conversion from the use of high-enriched uranium to the use of low enriched uranium fuels Guidebook] under flow blockage of a single Fuel Assembly (FA) conditions. Such event was rarely investigated in the open literature notwithstanding the fat that it constitutes a severe accident that may lead to local dryout and eventually to loss of the FA integrity. The transients herein considered are related to partial and total obstruction of the cooling channel of a single Fuel Assembly of the reactor core. This study constitutes the first step of a larger work, which consists in performing a 3D simulation using the Best Estimate coupled code technique. However, as a first approach the instantaneous reactor power is derived through the point kinetic approach of the used thermal-hydraulic system code.     

Achievement of 1.9 MW for 300 s during the commissioning of the KSTAR ICRF transmitter

KSTAR (Korea Superconducting Tokamak Advanced Research) is a national tokamak aiming at the high beta operation based on AT (Advanced Tokamak) scenarios in Korea and ICRF (Ion Cyclotron Ranges of Frequency) is one of the essential heating and current drive tools to achieve this goal. The ICRF heating and current drive scenario requires 4 units of 2 MW transmitters with a frequency range from 25 to 60 MHz. The first KSTAR transmitter is a modified FMIT (Fusion Material Irradiation Test) transmitter consisting of four amplifier stages. An amplitude-modulated 1 mW frequency source drives a 500 W solid state wideband amplifier, which in turn drives three tuned triode/tetrode amplifier stages. The tube employed in the final power amplifier is a 4CM2500KG tetrode fabricated by CPI (Communications & Power Industries). After the fabrication of the cavity and power supply was completed in 2004, several failures of the tube during a factory and a site acceptance test occurred before eventually achieving 1.9 MW for 300 s at 33 MHz in 2007. The electrical efficiency of the FPA (Final Power Amplifier) is about 70%. Although this is a very encouraging result for the development of an ICRF transmitter for ITER (International Thermonuclear Experimental Reactor), continued efforts for a reliable operation are required to achieve the final goals of the KSTAR and ITER ICRF system.     

Experimental investigation and CFD simulation of horizontal stratified two-phase flow phenomena

For the investigation of stratified two-phase flow, two horizontal channels with rectangular cross-section were built at Forschungszentrum Dresden-Rossendorf (FZD). The channels allow the investigation of air/water co-current flows, especially the slug behaviour, at atmospheric pressure and room temperature. The test-sections are made of acrylic glass, so that optical techniques, like high-speed video observation or particle image velocimetry (PIV), can be applied for measurements. The rectangular cross-section was chosen to provide better observation possibilities. Moreover, dynamic pressure measurements were performed and synchronised with the high-speed camera system.CFD post-test simulations of stratified flows were performed using the code ANSYS CFX. The Euler–Euler two fluid model with the free surface option was applied on grids of minimum 4 × 10⁵ control volumes. The turbulence was modelled separately for each phase using the k–ω-based shear stress transport (SST) turbulence model. The results compare very well in terms of slug formation, velocity, and breaking. The qualitative agreement between calculation and experiment is encouraging and shows that CFD can be a useful tool in studying horizontal two-phase flow.     

Fully-coupled transient simulation for RCS and SPCS of an integral reactor plant

Design features of SMART such as a built-in once-through steam generator (OTSG) and a close interaction between the feedwater flow rate and steam pressure controls leads to the necessity of fully-coupled transient analysis tools of the reactor coolant system (RCS) and the steam and power conversion system (SPCS) for the purpose of a plant control system development. A fully-coupled transient simulation tool, MMS/SMART, was developed to test the capability of the plant control system for the normal load-following event and the anticipated abnormal events. The MMS/SMART was composed of several interacting MMS modules with numerical data, each of which represented a component of the SMART plant and a control logic. The RCS and the SPCS with their control logics were modeled using default modules such as a pipe, pump and tank. The developed MMS/SMART was validated by using the scaled-down experimental data and the analysis result from the TASS/SMR code. A simulation result for the 100–50–100% load-following operation with a 25%/min rate shows that the feedwater flow rate and the steam pressure are controlled well as expected, except for small-amplitudes of steam pressure fluctuation at the lower power operating region. The loss of turbine load event was also simulated and the result shows that the plant can be operated stably with the steam bypass control system.     

Design and manufacture of the fuel element for the 10 MW high temperature gas-cooled reactor

The Chinese 10 MW high temperature gas-cooled reactor (HTR-10) attained its first criticality on December 21, 2000. The fabrication of the first fuel for the HTR-10 started in February 2000 at the Institute of Nuclear Energy Technology (INET), Tsinghua University. Up to September 2000, a total of 11 721 spherical fuel elements were successfully produced. The average free uranium fraction of the first fuel-determined by the burn-leach method-was 5.0×10⁻⁵. So far, the release rate R/B of the fission gas, measured in the irradiation test, shows that not a single particle in three irradiated spherical fuel elements failed as the results of the irradiation test carried out in Russia. This paper describes the design parameter, the fabrication technology and the performance data of the HTR-10 first fuel, and the production and quality control experiences obtained from the manufacture of the first fuel for the HTR-10.     

CFD modelling of a negatively buoyant purge flow in the body of a reactor coolant circulator

The paper describes a finite element flow calculation using a two-equation turbulence model and variable properties, which analyses the effect of injecting a negatively buoyant flow stream into the primary circuit of a gas cooled reactor. This flow is intended to purge moisture-laden gas from the compartment housing the electrical motors driving the primary circuit pumps. The study highlights the sensitivity to the use or otherwise of the Boussinesq approximation and the type of boundary conditions specified at the outlet of the calculational domain. The method of calculation was benchmarked against a challenging experiment on the mixing of stably stratified saline streams, as well as the higher order CFD modelling of the same experiment.     

Effect of 10 MW/m2 plasma exposure on W-SiC/SiC dual layer tiles

W-SiC/SiC dual layer tile has many advantages as a high heat flux component (HHFC) material for fusion, in theory. However, due to insufficient data known, its high potentiality and near term availability has not been well recognized. This work provides the recent materials R&D status and the first plasma exposure test result from the world largest helical device, large helical device of National Institute for Fusion Science in Japan. Tungsten armor with SiC/SiC substrate layer survived during the LHD plasma exposure with 10 MW/m² maximum heat load for the 5.3-s operation cycle. The macro and microstructure evolution, including crack and pore formation, was analyzed and an excellent high heat load resistance was demonstrated.     

Evaluation of a numeric procedure for flow simulation of a 5 × 5 PWR rod bundle with a mixing vane spacer

The fuel assemblies of the Pressurized Water Reactors (PWR) are constituted of rod bundles arranged in a regular square configuration by spacer grids placed along its length. The presence of the spacer grids promote two antagonist effects on the core: a desirable increase of the local heat transfer downstream the grids and an adverse increase of the pressure drop due to the constriction on the coolant flow area. Most spacer grids are designed with mixing vanes which cause a cross and swirl flow between and within the subchannels, enhancing even more the heat transfer performance in the grid vicinity. The improvement of the heat transfer increases the departure from the nucleate boiling ratio, allowing higher operating power in the reactor. Due to these important thermal and fluid dynamic features, experimental and theoretical investigations have been carried out in the past years for the development of spacer grid design. More recently, the Computational Fluid Dynamics (CFD) using three dimensional Reynolds Averaged Navier Stokes (RANS) analysis has been used efficiently for this purpose. Many computational works have been performed, but the appropriate numerical procedure for the flow in rod bundle simulations is not yet a consensus. This work presents results of flow simulations performed with the commercial code CFX 11.0 in a PWR 5 × 5 rod bundle segment with a split vane spacer grid. The geometrical configuration and flow conditions used in the experimental studies performed by Karoutas et al. were assumed in the simulations. To make the simulation possible with a limited computational capacity and acceptable mesh refinement, the computational domain was divided in 7 sub-domains. The sub-domains were simulated sequentially applying the outlet results of a previous sub-domain as inlet condition for the next. In this study the k-ε turbulence model was used. The simulations were also compared with those performed by Karoutas et al. in half a subchannel and In et al. in one subchannel computational domains. Comparison between numerical and experimental results of lateral and axial velocities along of the rod bundle show good agreement for all evaluated heights downstream the spacer grid. The present numerical procedure shows better predictions than Karoutas et al. model especially further from the spacer grid where the peripheral subchannels have more influence in the average flow.     

Analysis of the IAEA research reactor benchmark problem by the RETRAC-PC code

A basic approach to perform safety analysis of a nuclear research reactor consists in using deterministic methods to verify that the established acceptance criteria related to fuel integrity are fulfilled during all the stages of the facility lifetime. These methods should be validated against a large set of experimental and postulated transients. Since measured data are not easily available in the literature, the IAEA defined typical transients in a generic 10-MW MTR nuclear reactor core as a benchmark test for computational tools verification. In this framework, an assessment study of the coupled kinetic–thermal–hydraulic RETRAC-PC code is presented herein. The considered cases include the analysis of core dynamic under ramp positive reactivity insertion, and loss of flow transients. In general, the obtained results are satisfactory and agree with results obtained by other similar codes.     

Cooling mechanism during transient reflooding of a reactor fuel bundle after loss of coolant

Cooling efficiency during transient reflooding under loss of normal coolant conditions has been examined with a 7 × 7 simulated fuel rod bundle and jet pump bypass. The bundle contains 49 electrically heated rods with 3600 mm heated length and a pseudo cosine axial power distribution. Water is injected into the lower plenum and the superheated bundle is reflooded from the bottom with some flow diverted to the simulated jet pump bypass. The results show that effective cooling can be maintained.     

10 MW高温气冷实验堆吸收球停堆系统设备热态试验

碳化硼吸收球停堆系统是10 MW高温气冷实验堆的第二停堆系统,其功能是,在控制棒失效时,吸收球落入反射层的吸收球孔道,以达到紧急停堆的目的.介绍了在堆外、空气介质、150 ℃工作温度条件下,对吸收球传动机构设备进行的热态考验、传动试验和落球试验.结果表明,吸收球停堆系统7套设备均达到了传动机构工作正常、落球时间在3 min之内、球位指示正确等设计要求.7套设备安装到10 MW高温气冷实验堆上之后,在堆上进行了氦气介质下的热态落球试验,其结果达到设计规范的要求.