Table of Contents Author Guidelines Submit a Manuscript
Science and Technology of Nuclear Installations
Volume 2017 (2017), Article ID 5151890, 8 pages
Research Article

Analysis of NEA-NSC PWR Uncontrolled Control Rod Withdrawal at Zero Power Benchmark Cases with NODAL3 Code

1Center for Nuclear Reactor Technology and Safety, National Nuclear Energy Agency of Indonesia (BATAN), Kawasan Puspiptek, Gd. No. 80, Serpong, Tangerang Selatan 15310, Indonesia
2Nippon Advanced Information Service (NAIS Co., Inc.), 416 Muramatsu, Tokaimura, Ibaraki 319-1112, Japan
3Cooperative Major in Nuclear Engineering, Tokyo City University, 1-28-1 Tamazutsumi, Setagaya-ku, Tokyo 158-8557, Japan

Correspondence should be addressed to Tagor Malem Sembiring

Received 24 August 2017; Revised 26 October 2017; Accepted 6 November 2017; Published 23 November 2017

Academic Editor: Arkady Serikov

Copyright © 2017 Tagor Malem Sembiring et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


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.