A multi-region collision probability method for determining neutron spectra and reaction rates

The collision probability approach to neutron transport can be used to obtain the energy-dependent neutron spectrum in nuclear reactor systems as well as other quantities of interest. This method makes the approximation that the neutron distribution is constant within homogeneous regions, or cells, in the system. This assumption restricts geometries that can be modeled by the collision probability approach. The geometry modeled is typically an infinite lattice of two homogeneous cells: a fuel pin cylinder and the coolant that surrounds it. The transport of neutrons between the homogeneous cells is done using probabilities describing the chance that a neutron having a collision in one cell has its next collision in another cell. These collision probabilities can be cast in terms of escape and transmission probabilities for each cell. Some methods exist that extend the collision probability approach to systems composed of more than two homogeneous cells. In this work, we present a novel collision probability method, based on previous work by Schneider et al. (2006a), for an arbitrary number of cells. The method operates by averaging the transmission probabilities across cells of the same shape, and thus assumes a certain level of homogeneity across all cells. When using multigroup cross sections, which the collision probability approach requires, it is necessary to consider the effect that a system's geometry and composition has on those multigroup cross sections. The cross sections must be computed in a way that accounts for the resonance self-shielding that may reduce the reaction rates in the resonance region. The process of developing self-shielded cross sections in a heterogeneous system utilizes an escape cross section. We compute this escape cross section using the same collision probabilities used to obtain the energy spectrum. Results are presented for simple two-cell systems, and preliminary results for four-cell simulations are also given. An extension to the method is provided that accounts for the fact that in thermal systems the assumption of homogeneity is not always valid.