Numerical Solution to Space-Angle Energy-Dependent Neutron Integral Transport Equation

Abstract

A numerical approach to the steady-state, space-, angle- and energy-dependent neutron transport equation is presented for neutron shielding calculations. The scattering integral, with anisotropic treatment of elastic scattering and isotropic treatment of inelastic scattering, is evaluated by the use of Gaussian and straightforward quadratures. A system of coupled one-group integral equations for all the energy meshes of interest, converted from the energy-dependent integral transport equation, is calculated by performing a line integration along the neutron path in the direction of motion. For this purpose the direction of neutron motion is represented by discrete-ordinate directions Ω_pq on the unit sphere.
The final presentation of the integral transport equation is derived in a difference form convenient for machine computation. A computation program PALLAS has been written in Fortran IV for IBM 360-75 computer to perform neutron transport calculations based on this approach.
Comparisons are given of the numerical solutions with analytical solutions for unscattered fluxes in various geometries such as plane, spherical and two-dimensional cylindrical, for volume sources with self absorption, and with experimental spectra for angular neutron fluxes in graphite-, polyethylene- and water-shield. Excellent agreement is obtained between the present calculations and analytical or experimental results.