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.     

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.     

Homogeneous zones definition in deterministic codes and effect on computed neutronic parameters

The design or modification and in general the analysis and control of nuclear reactors require complex calculations, which are carried out by numerical codes including neutronic and thermal-hydraulic components. Among the neutronic codes, the deterministic ones which solve the neutron transport/diffusion equation simulate the reactor core by dividing it into homogenized zones, i.e. volumes within which the macroscopic nuclear properties are considered uniform. These codes have been extensively used and tested for several decades and are shown to perform well when they analyze reactor cores containing regions with relatively homogeneous distributions of fuel, moderator and absorbing materials. In this work, the sensitivity of computed key neutronic parameters to the partitioning of the reactor core in homogenized zones is examined. Application is made for a configuration of the Greek Research Reactor (GRR-1) core, which is pool type, fueled by slab-type fuel elements. For the calculations, the neutronic code system consisting of XSDRNPM (cell-calculations) and CITATION (core analysis) is used with two different definitions of homogeneous zones for the special/control fuel assemblies. The effect on computations of neutron flux distribution, void-induced reactivity and total control rod worth is examined based on corresponding measurements. It is shown that with a more appropriate partition in homogeneous zones, the agreement of computed results with measurements can be remarkably improved concerning mainly the neutron flux, while the control rods worth is the less affected quantity.     

Introducing and validating a new method for coupling neutronic and thermal-hydraulic calculations

In this study, a new and innovative method is introduced for analyzing neutronic and thermal-hydraulic calculation. For this aim, VVR-S research reactor was selected, and the calculation procedure was performed for it. WIMS, CITATION and COBRA-EN codes were used for investigation. Calculation model consists of two sub-models: neutronic and thermo-hydraulic. The neutronic model uses WIMS and CITATION codes for neutronic simulation of the reactor core and calculating flux and power distribution over it. WIMS code simulates the fuel assemblies and CITATION models the core. The thermal-hydraulic model uses COBRA-EN code for performing the relative calculation. In this study, FORTRAN 90 program is used for linking two sub-models and performing the calculation. The proposed procedure is performed for VVR-S analysis and finally, the obtained results are compared with the experimental results that show good agreement with it.     

Thermal-hydraulic transient analysis of light water reactors using moving boundaries

This paper examines the applicability of a mathematical dynamic model developed here for the simulation of the thermal-hydraulic transient analysis for light water reactors (LWRs). The thermal-hydraulic dynamic modeling of the fuel pin and adjacent coolant channel in LWRs is based on the moving boundary concept. The fuel pin model (FUELPIN) with moving boundaries is developed to accommodate the core thermal-hydraulic model, with detailed thermal conduction in fuel elements. Some results from transient calculations are examined for the first application of the thermal-hydraulic model and the fuel pin model with moving boundaries in a boiling water reactor (BWR). An accurate minimum departure from nucleate boiling ratio (MDNBR) and its axial MDNBR boundary versus time within the fuel channel are predicted during transients. Transient analysis using a known thermal-hydraulic code, COBRA and FUELPIN linked with a PWR systems analysis code show that the thermal margin gains more by a transient MDNBR approach than the traditional quasi-steady methodology for a pressurized water reactor (PWR). The studies of the overall nuclear reactor system show that moving boundary formulation provides an efficient and suitable tool for thermal transient analysis of LWRs.     

Reactivity control technique for a pressurized water reactor with an inventive TRISO fuel particle composition

In this research paper a reactivity control technique has been suggested for the conceptual design of a compact sized pressurized water reactor (PWR) core with an inventive tristructural-isotropic (TRISO) fuel particle composition. This conceptual design is a light water cooled and moderated reactor which utilizes TRISO fuel particles in PWR technology. The use of TRISO fuel in PWR technology improves integrity of the design due to its fission fragments retention ability. The fuel provides first retention barrier within fuel itself against the release of fission fragments that makes this design concept safer and environment friendly. The suggested TRISO fuel particle composition has a small amount of Pu-240 with 2.0 w/o in the place of U-238 which acts as reactivity suppressor. Reactor codes WIMS-D/4 and CITATION have been used for simulation and core design modeling. Results reveals that the amount of excess reactivity can be reduced significantly by using a small amount of Pu-240 in TRISO fuel which in turns reduces the number of gadolinia rods in the core required for excess reactivity control and completely eliminates the requirement of soluble boron system. Therefore the effective and optimal use of reactivity suppressor and burnable poison suppresses and flattens the core excess reactivity throughout the core life and hence number of control rods can be reduced without compromising on the shutdown margin.     

FRED fuel behaviour code: Main models and analysis of Halden IFA-503.2 tests

The FRED fuel rod code is being developed for thermal and mechanical simulation of fast breeder reactor (FBR) and light-water reactor (LWR) fuel behaviour under base-irradiation and accident conditions. The current version of the code calculates temperature distribution in fuel rods, stress-strain condition of cladding, fuel deformation, fuel-cladding gap conductance, and fuel rod inner pressure. The code was previously evaluated in the frame of two OECD mixed plutonium-uranium oxide (MOX) fuel performance benchmarks and then integrated into PSI's FAST code system to provide the fuel rod temperatures necessary for the neutron kinetics and thermal-hydraulic modules in transient calculations. This paper briefly overviews basic models and material property database of the FRED code used to assess the fuel behaviour under steady-state conditions. In addition, the code was used to simulate the IFA-503.2 tests, performed at the Halden reactor for two PWR and twelve VVER fuel samples under base-irradiation conditions. This paper presents the results of this simulation for two cases using a code-to-data comparison of fuel centreline temperatures, internal gas pressures, and fuel elongations. This comparison has demonstrated that the code adequately describes the important physical mechanisms of the uranium oxide (UOX) fuel rod thermal performance under steady-state conditions. Future activity should be concentrated on improving the model and extending the validation range, especially to the MOX fuel steady-state and transient behaviour.     

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.     

Design simplification of a small nuclear reactor for large-diameter neutron transmutation doping silicon using control rods

A design concept for a small nuclear reactor for neutron transmutation doping silicon (NTD-Si) using a Pressurized Water Reactor (PWR) full-length fuel assembly was proposed in our previous work. The excess reactivity was suppressed by a combination of Gd₂O₃ and soluble boron, which results in a flatter flux profile over the core than with control rod insertion; however, the soluble boron system for reactivity control is quite complex and expensive. The removal of this system would make the design much simpler. In the current work, the removal of soluble boron is considered. Criticality, neutron transportation and core burn-up calculations were performed using the MVP/GMVP II code and MVP-BURN code. The calculation results show that the insertion of control rods in five of the nine assemblies is enough to suppress reactivity. The thermal hydraulic analysis showed that heat removal from the core was possible under 1 atm operating pressure. Silicon ingots up to 30 cm in diameter could be irradiated with sufficient uniformity in the irradiation channels.     

A methodology for the coupling of RAMONA-3B neutron kinetics and TRAC-BF1 thermal-hydraulics

The initial objective of this project was to directly couple the RAMONA and TRAC codes running on different PCs. The idea was to use the best part of each one and eliminate some of their limitations and widen the applicability of these codes to simulate different BWR and system components. It was required to try to minimize the amount of changes to present code subroutines and calculation procedures. If possible, just substitute values obtained in the parallel code. Preliminary results indicated that using a CHAN component of TRAC to model thermal-hydraulic phenomena for each neutronic channel modeled in RAMONA is rather difficult. Large amounts of CPU time consumption are obtained and lots of PCs would make this solution difficult, besides considerable large transients are introduced by the differences in thermal-hydraulic results of these codes. The substitution of the thermal-hydraulics of RAMONA, by the TRAC channel calculations, is possible but simulation of a null transient on both codes must be planed and a gradual change must be controlled by an additional supervisory subroutine. An indirect coupling of these codes, it is therefore proposed, in order to eliminate most of these limitations. In this indirect coupling, a thermal-hydraulic model of the average tube in a bundle and the global channel cooling fluid dynamics is programmed for each neutronic channel while the global reactor vessel and core is modeled by TRAC with just four channels and four rings. Results are more reliable, coupling is simpler and faster simulations are possible.     

数值反应堆堆芯通道级三维热工水力程序CorTAF开发及初步验证

堆芯是核动力系统的核心部件,其完整性是反应堆安全运行的重要前提。传统核反应堆堆芯热工水力分析方法无法满足未来先进核动力系统的高精度模拟需求。本文依托开源CFD平台OpenFOAM,针对压水堆堆芯棒束结构特点建立了冷却剂流动换热模型、燃料棒导热模型和耦合换热模型,开发了一套基于有限体积法的压水堆全堆芯通道级热工水力特性分析程序CorTAF。选取GE3×3、Weiss和PNL2×6燃料组件流动换热实验开展模型验证,计算结果与实验数据基本符合,表明该程序适用于棒束燃料组件内冷却剂流动换热特性预测。本工作对压水堆堆芯安全分析工具开发具有参考和借鉴意义。     

三维六角形组件压水堆堆芯燃料管理计算及程序系统研究

介绍所研制的WWER型压水堆堆芯燃料管理计算程序系统TPFAP-H／CSIM-H，六角形组件均匀化计算程序TPFAP-H是在压水堆正方形组件程序TPFAP的基础上，采用穿透概率法与响应矩阵方法相结合计算六角形组件内中子能谱分布，并考虑六角形栅元特点改造开发而成的CSIM-H是以先进六角形节块扩散程序为基础．参照SIMULATE程序功能而研制的物理-热工水力耦合的三维六角形节块PWR堆芯燃料管理程序两者通过接口程序LINK连接起来，可以考虑燃耗，功率、慢化剂密度变化．控制棒、氙等参数的多种反馈效应对IAEA的WWER-1000型Kalinin核电厂基准问题的校算的结果表明，临界硼浓度、功率和燃耗分布等结果与国际各研究机构的结果吻合良好，偏差均在工程要求之内。     

RIA基准题验证堆芯三维物理-热工耦合程序

以COBRA-Ⅳ和NLSANMT程序为基础,开发了堆芯三维物理-热工耦合程序C4/NK。针对两个典型的反应性引入事故（RIA）,即NEACRP弹棒基准题和提棒基准题,分别进行了验证计算。与参考值和其他程序的计算结果相比,C4/NK耦合程序具有较好的精度,能正确模拟瞬态过程中的物理-热工反馈现象。     

物理-热工耦合对超临界水堆系统特性的影响分析

物理-热工耦合是超临界水堆系统分析的关键问题之一。以日本超临界水冷热堆Super LWR的堆芯设计为例,借助Dragon编制中子截面数据库,建立双群中子扩散方程计算模块,联系同时建立的热工计算模块,得到超临界水堆的物理-热工耦合计算模型。通过对比稳态与瞬态工况下耦合前、后的热工工况,分析物理-热工耦合条件下的超临界水堆系统热工特性。结果表明：在稳态工况下,物理-热工耦合将导致内、外组件堆芯功率峰值沿轴向发生明显偏移,使得部分节点的包壳温度升高,但包壳最高温度降低;在瞬态工况下,物理-热工耦合将导致堆芯包壳最高温度的发生位置有所改变。发生给水加热丧失瞬态后,在某一时刻,外部组件的包壳最高温度将转而超过内部组件的包壳最高温度。可见,物理-热工耦合对包壳最高温度的大小和发生位置均可能产生明显影响。计算分析可为超临界水堆瞬态及安全分析提供相应理论参考。     

环形燃料热工水力性能分析程序开发及验证

本工作开发了环形燃料子通道分析程序SAAF。采用SAAF计算了西屋公司四环路压水堆所用环形燃料组件的热工水力性能,并与VIPRE-01的计算结果进行比较。结果表明,SAAF与VIPRE-01的计算结果符合较好,SAAF可用于环形燃料热工水力设计分析。     

Development of a thermal-hydraulic subchannel analysis code for motion conditions

Since the conventional subchannel analysis codes are designed for the land-based reactor core, a thermal-hydraulic subchannel analysis code was developed to evaluate thermal-hydraulic characteristics of the reactor core under motion conditions. The verification of the code was performed with experimental data and commercial codes. The ISPRA 16-rod tests were used to evaluate the steady-state prediction performance of the code, and the simulation results agree well with the test data. COBRA-EN code was applied to check the transient prediction performance of the code, and there is a good agreement between the predictions with both codes. An additional forces model for motion conditions was proposed in the code, and CFX-14.0 code was applied to verify the model. The results show that the code can be used in the thermal-hydraulic analysis of the reactor core under motion conditions. To illustrate the capabilities of the code, a fuel bundle under a complex motion condition was simulated, and the results are reasonable.     

REA analysis of the Iranian VVER-1000 core by modification of PRORIA code for annular pellets

The PRORIA code and its recent modifications are described here. The PRORIA code analyzes the transient response of the core against the reactivity increase caused by the control rod rapid withdrawal. The code solves and analyzes neutronic and thermal–hydraulic equations simultaneously. The code is designed for western PWR-type reactor performance. The equations representing thermal–hydraulic and neutronic should be modified to use the code to analyze VVER-1000 reactor core transients, because The VVER-1000 reactor fuel has a central hole in the fuel pellet. In a cylindrical solid fuel pellet, operation of an oxide fuel material at high temperature alters its morphology and the inner region is restructured to form a void at the center surrounded by a dense fuel region. Most of the restructuring occurs within the first few days of operation with slow changes afterward. Hence, the effects of a central hole in mathematical equations and in the transient are investigated. After the code modification, three accident scenarios with control rod ejection are simulated. The results are in good agreement with those reported in the plant’s FSAR. The results show that the peak fuel temperature in the hot fuel pin is lower than what the original code predicts by 150–500 °C. Furthermore, the Doppler reactivity effect, when the fuel pellet has a central hole, is higher than the solid fuel pellet.     

Development of a control rod depletion methodology for the Westinghouse NEXUS system

The NEXUS project is an effort to merge and modernize the methods employed in Westinghouse PWR and BWR steady-state reactor physics codes. The NEXUS system relies on a once-through nodal cross-section generation methodology with an innovative and efficient technique for pin power recovery. The pin power methodology overcomes a well-known limitation of existing methodologies, namely the incapacity to properly account for heterogeneity changes due to the depletion environment. The so-called control rod history problem where control rods are repeatedly inserted and withdrawn during core depletion is a good example of such a case. In addition to the control rod history impact on pin power distributions, the insertion of control rods during extended periods leads to significant control rod depletion that affects the reactivity worth of the control rods which in turn can have a significant impact on pin powers. The importance of accurately predicting pin powers, combined with the need to adequately estimate the reactivity worth and nuclear end of life of control rods in BWRs and in generation III+ PWRs, has motivated the development of a novel control rod depletion model. This methodology and its numerical qualification, initially for PWR application only, is the topic of this paper. The focus is on describing the salient features of the model and on illustrating its performance by means of numerical experiments. It is shown that together with the NEXUS pin power recovery model, the control rod depletion methodology accurately predicts the reactivity feedback from repeated control rod insertions in a PWR core.     

An assessment of a VVER-1000 core during Turbo-Generator load reduction test using RELAP5/MOD3.2 and WIMSD-5B/PARCSv2.7

Reactor dynamic tests, which are categorized as one of the power start-up test groups, are the most complex tests during commissioning of the new nuclear power plants. This paper presents the results of Turbo-Generator load reduction test as one of the reactor dynamic tests for VVER-1000/V446 unit at Bushehr Nuclear Power Plant (BNPP). In this test modeling because of the need for control rod bank worth and core reactivity coefficients, the core geometry has been modeled first by using WIMSD-5B/PARCSv2.7 codes for neutronic calculations. For performing the thermal-hydraulic analysis, the RELAP5/MOD3.2 computer code has been used. The control rod bank worth and core reactivity coefficients obtained from WIMSD-5B/PARCSv2.7 are compared with BNPP FSAR that confirm the ability and reliability of the method. Also comparison of the thermal-hydraulic core parameters obtained from RELAP5/MOD3.2 against actual plant data, indicate that this code can properly predict behavior of VVER-1000 reactor for this dynamic start-up test.     

Neutronic design of hydride fueled BWRs

This paper summarizes the neutronic part of a study of the feasibility of designing BWR cores to have enhanced power density and simplified fuel bundle by using hydride instead of oxide fuel. A 3D fuel bundle neutronic analysis is performed for a limited number of geometries to determine attainable discharge burnup, pin-by-pin power distribution, axial power distribution, reactivity coefficients, reactivity worth of control elements and burnable absorber effects. It is found that hydride fuel bundle design can be simplified by eliminating water rods and partial length fuel rods and by reducing the volume of water in-between the fuel bundles. Both an ideal and more practical bundle designs are examined. A companion study of the thermal-hydraulic and vibration characteristics of BWR cores predicts that the increase in the number of fuel rods per given core volume enables increasing the BWR power density by up to ∼30% relative to oxide fuelled core design. The net outcome is expected to be improved BWR economics even though hydride fuel requires higher uranium enrichment to compensate for its reduced uranium loading.