Research Article  Open Access
Modification of Neutron Kinetic Code for Plate Type Fuel Nuclear Reactor
Abstract
The research is conducted on the modification of neutron kinetic code for the plate type fuel nuclear reactor. REMARK is a neutron kinetic code that works only for the cylindrical type fuel nuclear reactor. In this research, our main emphasis is on the modification of this code in order to be applicable for the plate type fuel nuclear reactor. For this purpose, detailed mathematical studies have been performed and are subjected to write the program in Fortran language. Since REMARK code is written in Fortran language, so we have developed the program in Fortran and then inserted it into the source library of the code. The main emphasis is on the modification of subroutine in the source library of the code for hexagonal fuel assemblies with plate type fuel elements in it. The number of steps involved in the modification of the code has been included in the paper. The verification studies were performed by considering the small modular reactor with hexagonal assemblies and plate type fuel in it to find out the power distribution of the reactor core. The purpose of the research is to make the code work for the hexagonal fuel assemblies with plate type fuel element.
1. Introduction
Realtime multigroup advanced reactor kinetics (REMARK) code is developed by the GSE power system used for the realtime simulation of nuclear reactor core. This code utilizes twogroup, 3D timedependent diffusion theory in the form of finite difference equations to simulate nuclear reactor core. This code has two neutron energy groups to accurately simulate the characteristics of fast and thermal neutron under normal and abnormal operating conditions of the nuclear reactor. The core geometry can be modeled as 3D mesh cell structure under the limited capacity of computer resources for realtime applications. Mesh cell sizes are chosen to ensure that the material properties within cell are homogenous or heterogeneous depending on the fuel, reflector, shielding, and so forth. In the REMARK code, there is improved quasistatic solution for obtaining the flux distribution and six delayed neutrons to be calculated in each mesh cell. Reactivity feedback is based on the core thermal hydraulic conditions like fuel temperature, coolant temperature, moderator density, and void fraction. This code can calculate the reactivity effect of control rod movement and soluble boron. It can also give the exact neutron flux readings at incore and excore detector points and utilizes neutron crosssection for neutronic calculations [1, 2]. This is 3D neutron kinetic code, but it cannot examine the special kind of geometry such as hexagonal fuel assemblies. Therefore, the current research is focused on the modification of the code to make it applicable for the hexagonal fuel assemblies. The flow diagram of REMARK code is given in Figure 2.
The main concern of our research is to modify the nuclear reactor core for hexagonal fuel assemblies. The research on the new developed hexagonal fuel assembly design has already been done on the supercooled water reactor (SCWR) to show the moderation capability for the hexagonal fuel assemblies [3, 4].
The benchmark studies of different calculation schemes related to reactor design, collision probability, power reconstruction, and so forth were investigated by ten different institutes. The focus was to find out the difficulties caused by different neutron data reduction methods and for the homogenization of the fuel assemblies [5, 6].
2. Reactor Core Simulation Model
For the hexagonal fuel assemblies we consider the hexagon with sides “” and center “O” as shown in Figure 1. The distance between opposite sides is denoted by “” and the height is denoted by “” [3]. The area and volume of the hexagon can be written as The onesided area of hexagon can be justified as We have considered eight sections of adjacent block of hexagon and calculated the neutron diffusion equations. Therefore the leakage of the neutron diffusion equation integral with volume of the node is given as where can be calculated as where is neutrons flux density at the interface.
Since the net flux of neutrons at the adjacent interface is continuous, so Substituting (5) in (4) we get Similarly , , , , and can be calculated as and have the same value as . is an integral on the top, whereas is an integral on the bottom and can be written as where is neutrons flux density at the interface.
Here, the net neutrons flux at the adjacent interface is continuous; therefore Substituting (8) in (7) we get where will have the same value as .
Now we will calculate the extrapolated length, that is, , to be considered as core boundary conditions: In threedimensional hexagonal grid, is called diffusion coefficient and is the leakage term in four directions. The leakage term is calculated in four directions as For the middle section block grid, Similarly, , , , and will have the same values.
For the grid block boundary section: direction left border direction right border direction border at the bottom direction border at the top
Similarly, , , , and will have the same values: In order to solve the diffusion and convection diffusion problems, we used alternating direction implicit method. This method was first discovered by Peaceman and Rechford [4] especially to solve parabolic and 2D elliptical problems.
The reactor core is a typical case of nonuniform problem and the parameters involved are not continuous. Therefore, ADI method was considered for calculating the neutron flux of each section in the block of reactor core. Neutron flux in whole reactor core is not continuous and there is a need to divide the coarse grid of core assembly. Hence, ADI method was employed to calculate the data and to give the accuracy of results, and, consistently, an algorithm should be adjusted hereafter [5].
3. Modification of REMARK Code
In realtime simulation, the computational time of each step is less than the simulation time, so it is difficult that the nodal expansion method, the nodal green function method, and other methods meet the requirements for the realtime simulation of the reactor core.
In REMARK code, neutron diffusion timespace dynamics equations with twogroup and sixgroup delayed neutrons are utilized for the simulation [5]. These equations are given as below.
Fast Group with . We have
Thermal Group with E < 1. We have The sixgroup delayed neutron precursor concentration can be calculated as In our research, we have experienced modelling of the realtime simulation technique for hexagonal fuel assemblies especially account to solve neutron diffusion equation. In hexagonal assembly, there are eight leakage directions which change the diffusion terms in diffusion equations. For the calculation of the reactor core, there are several steps that should be taken.(i)Grid homogenization (fuel/control grid in assemblies so as to make 1D model which consists of fuel, clad, and moderator) should be performed.(ii)Results from grid homogenization are used to calculate fewgroup or twogroup diffusion calculations.(iii)The finally originating group calculations are used in diffusion equation to solve the diffusion calculation of the whole reactor core.(iv)Effective multiplication factor and power or neutron flux distribution will be calculated.(v)Neutron transport equations will be employed after a series of approximations to get the diffusion calculations.
3.1. Assumptions
There are some assumptions that need to be considered for modification in the REMARK code in terms of hexagonal assemblies.(i)Ignore transfer of neutrons from thermal to fast group; that is, transfer crosssection is zero.(ii)Delayed or prompt neutron should be considered as fast neutron.(iii)Poison concentrations like Xe and Sm from I135 and Pm149 should be considered in calculations.(iv)Taking into consideration other phenomena, fission power and decay power can be calculated.(v)Now the effective multiplication factor will be calculated by neutrons originating from fission divided by neutrons from absorption and leakage.In this paper, we have compiled all the data in separate program of hexagonal fuel assemblies into many subroutines like:(i)neutron flux distribution,(ii)group 1 and group 2 leakage,(iii)delayed neutrons,(iv)discontinuity factor,(v)method to solve diffusion equation.Since REMARK code was written in Fortran language, so we developed the program in Fortran language against each factor stated above (Figure 3). These programs are then added in the source code of REMARK code and they compile the whole subroutines. The code was then tested by considering simulation of hexagonal fuel assemblies and it gave good results, thus pointing to the accuracy of modification of REMARK code for hexagonal fuel assemblies.
4. Verification Studies of REMARK Code
For verification of REMARK code, we have considered small modular nuclear reactor of power 220 MW designed by the College of Nuclear Science and Technology, Harbin Engineering University, China [6, 7]. The main design parameters of the reactor core are illustrated in Table 1.

For the verification studies, we have considered some of the parameters of the reactor core design and utilized in the lattice physic code, that is, HELIOS [8]. This code gives neutron crosssectional data which is then used in the REMARK code for obtaining the power distribution of the reactor core. In the current research, the reactor core has hexagonal fuel assemblies with plate type fuel in it. The REMARK code in its original form cannot work for the plate type fuel nuclear reactor; therefore, the modification studies have been done in the REMARK code in order to be applicable for plate type fuel reactor. The crosssectional data from HELIOS code exchanges with modified version of REMARK code which gives the power distribution of the reactor core and this verifies the modification of the code, although REMARK code needs some thermal hydraulic data which come from the plant design data and are used to get power distribution of the reactor core. The flow diagram of the process is shown in Figure 4.
In order to generate crosssectional data, HELIOS code has been used which is twodimensional (2D) lattice physics highorder transport code and is able to solve any kind of geometry. This code utilizes neutron and gamma groups for the burnup calculations, flux distribution, and micro/macro crosssectional data. It works via two separate codes AURORA and ZENITH which are input and output processing codes. The data exchanged between the two codes are accessed by the subroutine package HERMIS [8].
The obtained power distribution of the reactor core from combination of HELIOS and modified version of the REMARK code is shown in Figure 5.
5. Conclusion
In this paper, the research is focused on the modification and verification of the neutron kinetics code, that is, REMARK code. Originally this code only works for the cylindrical type fuel elements, so we have modified this code to be applicable to plate type fuel nuclear reactor. The modification has been done by detailed mathematics and the modified files in the source library of the code have been computed. The modified code has been verified by coupling with HELIOS code to get the power distribution of the plate type fuel nuclear reactor. The result obtained from coupling analysis gives the accuracy of the modified version of the code.
Acknowledgment
The authors would like to extend their sincere appreciation to the Deanship of Scientific Research at King Saud University for its funding of this research through the Research Group Project no. RGPVPP255.
References
 REMARK Modeling Techniques Handbook, GSE Power systems, Baltimore, Md, USA.
 REMARK Code Structure Manual Users Guide, GSE Power systems, Baltimore, Md, USA.
 L. Wang, D. Qin, Q. Li, and B. Xia, “A new SCWR fuel assembly with tworow fuel rods between the hexagonal moderator channels,” Annals of Nuclear Energy, vol. 58, pp. 60–64, 2013. View at: Publisher Site  Google Scholar
 S. Kaichao, P. Seltborg, and A. Libal, MCNP modeling of hexagon VVER fuel [M.S. thesis], Reactor Physics Department, Royal Institute of Technology, Stockholm, Sweden, 2008.
 Benchmark Calculations of Power Distribution within Fuel Assemblies. Phase II: Comparison of Data Reduction and Power Reconstruction Methods in Production Codes, NEA/NSC/DOC(2003)3, 2000.
 P. Meloni and M. Polidori, “A neutronicsthermalhydraulics model for preliminary studies on TRADE dynamics,” Nuclear Engineering and Design, vol. 237, no. 15–17, pp. 1704–1717, 2007. View at: Publisher Site  Google Scholar
 Y. Zhifei and Z. Qiang, “Real time simulation research on hexagonal fuel assembly core physics,” Nuclear Safety and Simulation, vol. 1, 2010. View at: Google Scholar
 Y. Hongjun, C. Zhenhua, and C. Yipping, “Thermal conduction model and its dynamic simulation of plate type fuel element,” Chinese Journal of Nuclear Science and Engineering, vol. 22, pp. 59–62, 2002. View at: Google Scholar
Copyright
Copyright © 2013 Salah UdDin Khan 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.