Thermal-neutron transport processes in heavy water in the temperature range of 253 to 4 K

Calculations are reported for time-dependent, space-dependent, and steady-state neutron spectra in D2O ice assemblies of different dimensions and in the temperature range of 253 to 40K. The scattering kernel used for these studies incorporates one- and two-phonon processes and is based on the Debye distribution function for the lattice vibrations of the D2O crystal. The multigroup Boltzmann diffusion equation was diagonalized to obtain transient and asymptotic spectra in assemblies at different temperatures with bucklings ranging from 0 to 0.15 cm-2. The calculated values of the effective decay constant are found to agree reasonably well with the experimental values reported by Salaita and Robeson at 2530K for waiting times of 160 to 320 μsec. The appropriate Boltzmann operator for the space-dependent problem was diagonalized to obtain its eigenvalues and eigenfunctions. By using these eigenfunctions, neutron spectra at different distances from the source plane were calculated in D2O ice at 2530K, and the diffusion lengths of neutrons were determined. Steady-state spectra in heavy ice assemblies at 77, 21, and 40K were also investigated. The results for the effective neutron temperatures and the cold-neutron fractions agree well with the experimental results of Rush et al. (14 figures, 5 tables) (U.S.)