Search results

This paper describes the effect of cross-section data generated by several codes on calculated neutronic parameters. The Pressurized Water Reactor Mixed Oxide and Uranium Oxide (PWR MOX/UO2) Core Transient Benchmark case was chosen because it has been used widely to validate neutronic codes. The cross-section data in this study will be generated by SRAC, Serpent, and HELIOS codes. The NODAL3 code will be used to calculate neutronic parameters from each cross-section. The neutronic parameters calculated by NODAL3 are the effective multiplication factor (keff), control rod worth, critical boron concentration, and power distribution under Hot Zero Power (HZP) conditions. The Power-Weighted Error (PWE) and Error-Weighted Error (EWE), as a measure of the relative error in fuel assembly power, are less than 5 %, indicating that the calculation is consistent with DeCART as a reference. The difference in calculated radial power peaking factor for all three cross-sections to reference data reaches 6.284 % (G-3), 8.438 % (G-3), and 10.998 % (C-7), respectively, for SRAC, Serpent, and HELIOS. The axial power distribution calculated by NODAL3 at the top and bottom of the reactor core has a relative error that peaked at 16.60 %, 13.86 %, and 10.20 %, respectively, for cross-sections provided by SRAC, Serpent, and HELIOS. Further improvements are needed for NODAL3 by applying various discontinuity factors to improve its performance.

The in-house coupled neutronic and thermal-hydraulic (N/T-H) code of BATAN (National Nuclear Energy Agency of Indonesia), NODAL3, based on the few-group neutron diffusion equation in 3-dimensional geometry using the polynomial nodal method, has been verified with static and transient PWR benchmark cases. This paper reports the verification of NODAL3 code in the NEA-NSC PWR uncontrolled control rods withdrawal at zero power benchmark. The objective of this paper is to determine the accuracy of NODAL3 code in solving the continuously slow and fast reactivity insertions due to single and group of control rod bank withdrawn while the power and temperature increment are limited by the Doppler coefficient. The benchmark is chosen since many organizations participated using various methods and approximations, so the calculation results of NODAL3 can be compared to other codes' results. The calculated parameters are performed for the steady-state, transient core averaged, and transient hot pellet results. The influence of radial and axial nodes number was investigated for all cases. The results of NODAL3 code are in very good agreement with the reference solutions if the radial and axial nodes number is 2 × 2 and 2 × 18 (total axial layers), respectively.

This paper reports the results of sensitivity analysis of the multidimension, multigroup neutron diffusion NODAL3 code for the NEACRP 3D LWR core transient benchmarks (PWR). The code input parameters covered in the sensitivity analysis are the radial and axial node sizes (the number of radial node per fuel assembly and the number of axial layers), heat conduction node size in the fuel pellet and cladding, and the maximum time step. The output parameters considered in this analysis followed the above-mentioned core transient benchmarks, that is, power peak, time of power peak, power, averaged Doppler temperature, maximum fuel centerline temperature, and coolant outlet temperature at the end of simulation (5 s). The sensitivity analysis results showed that the radial node size and maximum time step give a significant effect on the transient parameters, especially the time of power peak, for the HZP and HFP conditions. The number of ring divisions for fuel pellet and cladding gives negligible effect on the transient solutions. For productive work of the PWR transient analysis, based on the present sensitivity analysis results, we recommend NODAL3 users to use2×2radial nodes per assembly,1×18axial layers per assembly, the maximum time step of 10 ms, and 9 and 1 ring divisions for fuel pellet and cladding, respectively.

A coupled neutronics thermal-hydraulics code NODAL3 has been developed based on the few-group neutron diffusion equation in 3-dimensional geometry for typical PWR static and transient analyses. The spatial variables are treated by using a polynomial nodal method while for the neutron dynamic solver the adiabatic and improved quasistatic methods are adopted. In this paper we report the benchmark calculation results of the code against the OECD/NEA CRP PWR rod ejection cases. The objective of this work is to determine the accuracy of NODAL3 code in analysing the reactivity initiated accident due to the control rod ejection. The NEACRP PWR rod ejection cases are chosen since many organizations participated in the NEA project using various methods as well as approximations, so that, in addition to the reference solutions, the calculation results of NODAL3 code can also be compared to other codes' results. The transient parameters to be verified are time of power peak, power peak, final power, final average Doppler temperature, maximum fuel temperature, and final coolant temperature. The results of NODAL3 code agree well with the PHANTHER reference solutions in 1993 and 1997 (revised). Comparison with other validated codes, DYN3D/R and ANCK, shows also a satisfactory agreement.

A nuclear reactor cooling system that has been operating for a long time can carry some debris into a fuel coolant channel, which can result in a blockage. An in-depth two-dimensional simulation of partial channel blockage can be carried out using FLUENT Code. In this study, a channel blockage simulation is employed to perform a safety analysis for the TRIGA-2000 reactor, which is converted using plate-type fuel. Heat generation on the fuel plate takes place along its axial axis. The modelling of the fuel-plate is in the form of a rectangular sub-channel with an inlet coolant temperature of 308 K with a low coolant velocity of 0.69 m/s. It is assumed that blockage is in a form of a thin plate, with the blockage area being assumed to be 60 %, 70 %, and 80 % at the sub-channel inlet flow. An unblocking condition is also compared with a steady-state calculation that has been done by COOLOD-N2 Code. The results show that a partial blockage has a significant impact on the coolant velocity. When the blockage of 80 % occurs, a maximum coolant temperature locally reaches 413 K. While the saturation temperature is 386 K. From the point of view of the safety aspect, the blockage simulation result for the TRIGA-2000 thermal-hydraulic core design using plate-type fuel shows that a nucleate boiling occurs, which from the safety aspect, could cause damage to the fuel plate.

Steady-state and transient analysis of reactor core under Reactivity-Initiated Accident (RIA) conditions are important for reactor operation safety. The reactor dynamics are influenced by neutronic and thermal-hydraulic aspects of the core. In this study, steady-state and transient analysis under RIA conditions of the RSG-GAS multipurpose reactor was carried out using MTR-DYN and EUREKA-2/RR programs. Neutronic calculations were performed using a few group cross-sections generated by Serpent 2 with the latest cross-section data ENDF/B-VIII.0. Steady-state conditions were carried out with a nominal power of 30 MW, while transient under RIA conditions occurred because the control rod was pulled too quickly while the reactor operated. These transient RIA conditions were performed for two cases, during start-up with an initial power of 1 W, and within power range with an initial power of 1 MW. Thermal-hydraulic parameters considered in this study are reactor power, the temperature of the fuel, cladding, and coolant. The calculated maximum fuel temperature at a steady state is 126.02°C. Meanwhile, the calculated maximum fuel temperature during RIA conditions at the initial power of 1 W and 1 MW are 64.38°C and 137.14°C, respectively. There are no significant differences in thermal-hydraulic parameters between each used program. The thermal-hydraulic parameters such as the maximum temperature of the coolant, cladding, and fuel under this postulated RIA condition are within the acceptable reactor operation safety limits.