Application of computational codes in simulation, modeling and development of the power monitoring tools

Some developed codes and simulators for improving the power monitoring will be reviewed in this section. For example, MCNP (monte-carlo n-particle transport code) is developed for neutron detector design, or modeling a fission chamber to optimize its performance

7.1 Computational tools to conduct experimental optimization

Research reactors need a handy computational tool to predict spatial flux changes and following power distribution due to experimental requirements. Therefore it is important to get accurate and precise information ahead of any modifications. To meet this demand, flux measurements were conducted in case that a typical flux trap inside the core to be allocated. In TRR, one of standard fuel boxes, in position D6 in core configuration of the year 1999, was taken out of the core and a water trap was formed in its place. With the aid of miniature neutron detector (MND) using standard procedure, thermal neutron flux is measured inside the water trap. To calculate the flux and power theoretically, two different computational approaches such as diffusion and Monte Carlo methods were chosen. Combination of cell calculation transport code, WIMS-D5, and three-dimensional core calculation diffusion code, such as CITATION, were used to calculate neutron flux inside the whole core either in two or five energy groups. However, MCNP-4, as a Monte Carlo code, was used to calculate neutron flux again inside the whole core as well as inside the trap (Khalafi and Gharib, 1999). Figure 14 shows axial thermal flux distribution along the D6 position by measurement and computation.

It is obvious from the figure; the both calculation codes are satisfactory and a good agreement exists between detector measurements and code computations. However, diffusion method is a rational choice especially for survey calculation where the Monte Carlo approach is more time demanding. For some consideration, in order to measure spectrum, a fixed point on the midplane along D6 axis was chosen. A variety of foils of different material was selected as measuring windows to determine differential fluxes at specified energy bins. Metal foils such as Ti, Se, Mg, Ni, Al, Co, Au, In, and Fe were selected as energy windows. These foils are sensitive to a part of neutron energy spectrum starting form high energies and ending to thermal energies. Induced activity of each foil is measured based on gamma spectroscopy using high purity Germanium (HPGe) detector. By providing raw counts to SAND-II computer code, neutron energy spectrum was calculated. The measured and calculated spectrum using neutron detector, MCNP and WIMS codes is shown in Figure 15.

image507

Fig. 14. Axial thermal neutron flux distribution in trap at D6 position (Khalafi and Gharib, 1999)

image508

Fig. 15. Detector measured and MCNP and WIMS code calculated neutron spectrum(Khalafi and Gharib, 1999)

Spectrum calculations were also checked against measurements. Monte Carlo shows a better prediction while WIMS provides a fair result. It is notable that combination of WIMS/CITATION would be sufficient for neutron flux calculations while Monte Carlo technique should be reserved for the final stages of simulation. A good choice of computational tools would save time a lot in this respect and one is encouraged to perform a comprehensive simulation ahead of design and construction of irradiation facility.