Use of the Chebyshev-Legendre quadrature set in discrete-ordinate codes

The Chebyshev-Legendre (C-L) quadrature set has been implemented in the neutron transport section of the TWODANT code. The C-L quadrature set has two advantages as compared to other quadrature sets. First, it is possible to easily generate S/sub N/ orders up to S/sub 100/ with weights that are positive. TWODANT has been used to carry out calculations with as many as 10,000 discrete directions with the C-L quadrature set. Second, a product C-L set or order N will integrate exactly all the terms in the scattering source arising from P/sub N-1/ scattering. For example, the S/sub 4/ set will correctly integrate all scattering source terms if a P/sub 3/ scattering approximation is used. In TWODANT, users normally use the current S/sub 8/ quadrature set in the belief that P/sub 3/ scattering will be treated correctly. The S/sub 8/ set consists of some 40 directions and many times results in problem models that require excessive computer memory. The C-L S/sub 4/ set integrates exactly all the source terms associated with a P/sub 3/ expansion and consists of only 16 discrete directions. The results of test calculations, using both the current set and the C-L set in TWODANT, indicate that when high-ordermore » scattering is not significant, the two sets give answers that are in good agreement. 3 refs.« less