Comparison of Numerical Solution Techniques for Calculating Low Energy Neutrons

In this paper we investigate three numerical techniques for solving the one-dimensional straight-ahead Boltzmann equation for calculating the flux of low energy neutrons produced within a shielding material. The one-dimensional Boltzmann equation is split into a forward and backward coupled system of equations representing the production of ions of various types within a shielding material. The three numerical methods are then compared with neutron data from the Mir and ISS space station as well as Monte Carlo simulations for the production of low energy neutrons.