Benchmarking a diffusion code for a cylindrical-core reactor

The distorted-buckling method, proposed by us previously, allows the benchmarking of a diffusion code by comparing it with an analytic model in either 2 or 3 dimensions. Here, the method is applied to the case of a cylindrical TRIGA-type reactor to compare the fluxes predicted by an analytic model of the core and reflector, to those predicted by the code CITATION. The match is everywhere excellent. ©     

Development and application of a multi-physics and multi-scale coupling program for lead-cooled fast reactor

In this study, a multi-physics and multi-scale coupling program, Fluent/KMC-sub/NDK, was developed based on the user-defined functions(UDF) of Fluent, in which the KMC-sub-code is a sub-channel thermal–hydraulic code and the NDK code is a neutron diffusion code.The coupling program framework adopts the ‘‘master–slave’’ mode, in which Fluent is the master program while NDK and KMC-sub are coupled internally and compiled into the dynamic link library(DLL) as slave codes. The domain decomposition m...     

Verification of NETFLOW code using plant data of sodium cooled reactor and facility

The plant simulation code NETFLOW on PC applicable to the liquid-metal cooled reactors has been developed on the basis of the models developed for single-phase and two-phase light water flow systems. The functions of this code have been verified by individual tests for light water flow systems and a sodium flow system. In order to apply this code to a sodium cooled fast reactor, several extra functions were verified using the plant data obtained using 50 MW steam generators and the Monju fast breeder reactor. Finally, the turbine trip transient of the Monju was simulated and the result was compared with the measured plant data. Good agreements were obtained in these verifications. As a result of the present study, the code can be applied as an education tool for students.     

一种铅冷快堆主要构件的辐照损伤计算

为满足公众对更安全、更经济和环境更友好的核能系统的需求,提出一种铅铋合金冷却的铅冷快堆(Breeding Lead-based Economical Safe System–Demonstration,BLESS-D)。BLESS-D反应堆采用池式结构,热功率300 MW。金属材料受中子辐照时将造成材料的晶格缺陷,导致材料的宏观性能变化,改变其物理和机械性能。BLESS-D反应堆中有许多在反应堆寿期内不可更换的关键部件和设备,这些构件在反应堆运行期间如受到中子辐照损伤,将影响构件材料的性能,进而导致设备的使用寿命,限制了反应堆的寿命。本文通过计算BLESS-D反应堆主要部件和设备的原子离位数(Displacement Per Atom,DPA),评估结构材料的辐照损伤程度。利用SPECTER程序和MCNP程序进行燃料包壳、内部容器、主泵泵壳、蒸汽发生器壳和反应堆容器的DPA模拟计算,计算结果与发生材料辐照效应的DPA限值进行比较,发现内部容器的累积DPA在20年寿期内超过了材料辐照效应限值,需要进一步分析并优化设计,确保其寿期内的安全性。     

Development of a steady state analysis code for a molten salt reactor

The molten salt reactor (MSR), which is one of the ‘Generation IV’ concepts, can be used for transmutation, and production of electricity, hydrogen and fissile fuels. In this study, a single-liquid-fueled MSR is designed for conceptual research, in which no solid material is present in the core as moderator, except for the external reflector. The fuel salt flow makes the MSR neutronics different from that of conventional reactors using solid fuels, and couples the flow and heat transfer strongly. Therefore, it is necessary to study the core characteristics with due attention to the coupling among flow, heat transfer and neutronics. The standard turbulent model is adopted to establish the flow and heat transfer model, while the diffusion theory is used for the neutronics model, which consists of two-group neutron diffusion equations for fast and thermal neutron fluxes, and balance equations for six groups of delayed neutron precursors. These two models which are coupled through the temperature and heat source are coded in a microcomputer program. The distributions of the velocity, temperature, neutron fluxes, and delayed neutron precursors under the rated condition are obtained. In addition, the effects of the inflow temperature, inflow velocity, and the fuel salt residence time out of the core are discussed in detail. The results provide some valuable information for the research and design of the new generation molten salt reactors.     

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.     

Methods for calculating coupled reactor systems

Theoretical approaches to analyzing neutron kinetics in coupled reactor systems are reviewed. A systematic exposition of methods based on the use of integral parameters and an integral equation decribing neutron transport with the energy and angular variables eliminated is given. Special attention is given to a reactor-subcritical assembly as a coupled system.Translated from Atomnaya Énergiya, Vol. 97, No. 6, pp. 403–414, December, 2004.This revised version was published online in April 2005 with a corrected cover date.     

Vectorization and parallelization of a production reactor assembly code

In order to efficiently use new features of supercomputers, production codes, usually written 10 – 20 years ago, must be tailored for modern computer architectures. We have chosen to optimize the CPM-2 code, a production reactor assembly code based on the collision probability transport method. Substantional speedups in the execution times were obtained with the parallel/vector version of the CPM-2 code. In addition, we have developed a new transfer probability method, which removes some of the modelling limitations of the collision probability method encoded in the CPM-2 code, and can fully utilize parallel/vector architecture of a multiprocessor IBM 3090.     

Criticality qualification of a new Monte Carlo code for reactor core analysis

In order to accurately simulate Accelerator Driven Systems (ADS), the utilization of at least two computational tools is necessary (the thermal–hydraulic problem is not considered in the frame of this work), namely: (a) A High Energy Physics (HEP) code system dealing with the “Accelerator part” of the installation, i.e. the computation of the spectrum, intensity and spatial distribution of the neutrons source created by (p, n) reactions of a proton beam on a target and (b) a neutronics code system, handling the “Reactor part” of the installation, i.e. criticality calculations, neutron transport, fuel burn-up and fission products evolution. In the present work, a single computational tool, aiming to analyze an ADS in its integrity and also able to perform core analysis for a conventional fission reactor, is proposed. The code is based on the well qualified HEP code GEANT (version 3), transformed to perform criticality calculations. The performance of the code is tested against two qualified neutronics code systems, the diffusion/transport SCALE-CITATION code system and the Monte Carlo TRIPOLI code, in the case of a research reactor core analysis. A satisfactory agreement was exhibited by the three codes.     

Vectorization and parallelization of a production reactor assembly code

In order to efficiently use new features of supercomputers, production codes, usually written 10 – 20 years ago, must be tailored for modern computer architectures. We have chosen to optimize the CPM-2 code, a production reactor assembly code based on the collision probability transport method. Substantional speedups in the execution times were obtained with the parallel/vector version of the CPM-2 code. In addition, we have developed a new transfer probability method, which removes some of the modelling limitations of the collision probability method encoded in the CPM-2 code, and can fully utilize parallel/vector architecture of a multiprocessor IBM 3090.     

一种铅-铋合金冷却快堆的高放废物积累量研究

目前商用压水堆积累了大量的长寿命高放废物,放射毒性强,衰变时间漫长,对环境和人类构成了长期威胁,作为6种第四代核能系统堆型中的一种,铅基冷却快堆在减少长寿命高放废物产生方面具有优势。基于此本文提出了一种热功率为300 MW的铅-铋合金冷却快堆设计。利用MCNP程序对反应堆堆芯进行建模并计算了堆芯在寿期初的主要物理参数,详细分析了燃耗过程中长寿命高放核素的积累量,并与一般压水堆长寿命高放核素的积累量进行了比较。结果表明,对主要关心的次锕系核素,铅-铋合金冷却快堆的产生量远小于压水堆的,而长寿命裂变产物的产生量与压水堆的相当。总体来说,铅-铋合金冷却快堆产生的长寿命高放废物总量小于压水堆的,可看出铅-铋合金冷却快堆在减少长寿命高放废物产生方面更具有竞争性。     

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.     

Verification of FEMAXI-7 code by using irradiation test in Halden reactor for He-pressurization effect on FGR of BWR fuels under power transient

Power ramp test for He-pressurization effect on fission gas release (FGR) of about 42GWd/tUO₂ boiling water reactor (BWR) fuel rods was analyzed by the fuel performance code FEMAXI-7. The experimental data were obtained with the two rods, which were base irradiated in the Halden reactor for 12 years (IFA-409), then subjected to the power ramp tests (IFA-535) to investigate the He-pressurization effect. The FEMAXI-7 calculations were performed by inputting rod specifications and experimental conditions in both the baseand test irradiations. The results showed that the calculations reasonably followed the trends of measured cladding elongation and FGR during the power ramp test, depending on the pellet temperature and fission gas atoms diffusion rate. Based on the calculated results, the reason that no apparent He-pressurization effect was observed in the experiment was considered to be caused by insufficient gas communication during strong pellet–clad mechanical interaction (PCMI) and enhanced gap thermal conductance by the solid–solid contact due to gap closure.     

Features and performance of a coupled three-dimensional thermal-hydraulic/kinetics TRAC-PF1/NEM pressurized water reactor (PWR) analysis code

This paper summarizes the current status of the Pennsylvania State University (PSU) version of the coupled three-dimensional (3-D) thermal-hydraulic/kinetics TRAC-PF1/NEM code for pressurized water reactor (PWR) transient and accident analysis and describes applications to reactivity insertion accident (RIA) simulations as well as recent developments. The TRAC-PF1/NEM methodology utilizes closely coupled 3-D thermal-hydraulics and 3-D core neutronics transient models to simulate the vessel and a 1-D simulation of the primary system. An efficient and flexible cross-section generation procedure was developed and implemented into TRAC-PF1/NEM. These features make the coupled code capable of modeling PWR reactivity transients, including boron dilution transients, in a reasonable amount of computer time. Three-dimensional studies on hot zero power (HZP) rod ejection and main steam line break (MSLB) transients in a PWR, as well as a large break loss-of-coolant-accident (LBLOCA) and boron dilution transients, were accomplished using TRAC-PF1/NEM. The results obtained demonstrate that this code is appropriate for analysis of the space-dependent neutronics and thermal-hydraulic coupled phenomena related to most current safety issues.     

BACO (Barra Combustible), a computer code for simulating a reactor fuel rod performance

The basic philosophy and mathematical structure of the fuel performance simulation code BACO is described. This code is based on a central finite-difference quasi-bidimensional approximation. Within that approximation, the thermoelastic-plastic behaviour of a in-service fuel rod is calculated by a set of equations which are linearized and solved for each time step by a sparse matrix inversion subroutine. The numerical method is shown to be stable and to converge rapidly to physically sound results for the stresses and strains. Changes in the fuel shape due to cracking and restructuring are included in the calculation within a self-consistent mathematical frame. Code convergence and accuracy are discussed by comparing some predictions against thermoelastic and plastic analytic solutions. An example of the code predictions for the rod state during a reactor shutdown is presented and discussed.     

Comparison of two models for a high temperature reactor coupled to a gas turbine

This paper presents two independent dynamic models of a nuclear gas turbine power plant. Both the high temperature nuclear reactor (HTR) and its energy conversion system (ECS) based on a direct Brayton cycle have been modelled. One model utilises RELAP5 for the ECS, the other Aspen Custom Modeler (ACM). The reactor model used in both models is a point kinetic model derived from a detailed reactor model. The ECS model is described and compared componentwise, with an emphasis on the turbomachinery. The total plant models are compared with each other by calculating two representative transients: one load rejection transient and one transient with the system at part load.     

Qualification of the CFD code Trio_U for full scale reactor applications

The article presents a procedure to qualify the Trio_U code for the prediction of the boron concentration at the core inlet of a French 900 MWe pressurized water reactor under accidental conditions (inherent dilution problem).¹ The objective of this procedure is to ensure that the validation calculations are performed with the same modelling hypotheses as the full scale reactor analysis, for which usually no experimental data are available. A density driven ROCOM experiment as well as an UPTF Tram-C3 experiment have been used for the qualification of the Trio_U code. Both experiments present similar thermal hydraulic conditions as the reactor case. The predicted boron concentration at the core inlet of the reactor shows that the potential return to criticality might not be excluded in the case of a small break LOCA. Further neutronic calculations are necessary to confirm this result.