In the conventional discrete ordinates approach, the scattering source is approximated as a truncated Legendre series. In case of highly anisotropic scattering problems (e.g., incident beam problems), the truncated Legendre scattering cross sections give unphysical negative cross sections in some values of $\mu$ (cosine of scattering angle). This unphysical artifact causes negative scattering sources and negative angular fluxes. In addition to that, negative angular fluxes may cause wrong scalar flux as well.
A new method to generate non-negative scattering cross sections, which is deterministic, is proposed. The main idea of this method is to make non-negative scattering cross sections that produce equivalent scattering sources. This method does not have a practical limitation to generate non-negative scattering cross sections because the calculations of the scattering sources are performed with the conventional truncated cross section data provided from the standard processing codes.
In both the neutron and the photon/electron coupled cases, an inadequate truncated Legendre order causes problems, e.g., inaccurate angular distribution of scattering for low order, oscillations of differential scattering cross sections where cross sections are small for high order expansions, and negative differential scattering cross sections in some values of $\mu$.
In the neutron case, if the anisotropy is very high compared to the truncated Legendre order, inaccurate angular distributions of scattering occur especially in within-group scattering. Even if the truncated Legendre order is quite high to represent angular distribution of scattering in the photon/electron coupled case, it still causes oscillations where the cross sections are small.
The new method achieves both an accurate angular distribution of differential scattering cross sections and non-negative scattering cross sections for neutron and photon/electron coupled cases. The generated non-nega...