Implementation of a least squares temperature fit of the thermal scattering law in mcnp 6

A. T. Pavlou,W. Ji

2015-01-01

Abstract:The accuracy of a computer simulation of neutron transport is limited by the available nuclear cross sections, which are dependent on many factors including the velocity of the neutron and the system temperature. The thermal energy range (below about 4 eV) requires many more scattering cross section datasets than at higher energies because of additional complications provided by the comparable energies of the neutron and target material. As a result, the memory burden for low-energy neutron scattering can become large. On-the-fly methods have been used to reduce the storage of cross section data by finding a temperature dependence. However, the low energy quantum, chemical and crystalline effects that can be ignored at higher energies make the temperature dependence of the double differential scattering cross section in the thermal energy range difficult to determine. A new approach has been developed by fitting the temperature dependence of energy and squared momentum transfer cumulative distribution functions using a least squares approach. The coefficients of these fits are used to directly sample the outgoing energy and angle for any temperature. This method has been implemented into the continuous-energy Monte Carlo code MCNP6 [1] and tested for the moderator material bound carbon in graphite. The secondary energy spectra and angular distributions were observed to be in good agreement by testing with Red Cullen’s ‘broomstick’ benchmark [2]. This new method only requires about 21.7 MB of total data storage which is orders of magnitude less than the current method which requires about 25 MB of data per temperature.